Saudações,<br><br>Tenho a seguinte dúvida, se algúem poder me responder ficaria grata.<br>Estou trabalhando com a seguinte função no R (stdiagn), apresentada abaixo:<br><br><i>> stdiagn<br>function (pts, stkh, stse, stmc = 0, Dzero = FALSE) <br>
{<br> oldpar <- par(mfrow = c(2, 2))<br> pointmap(as.points(pts))<br> title("Data map")<br> oprod <- outer(stkh$ks, stkh$kt)<br> <span style="color:rgb(204,0,0)"> </span><u style="color:rgb(204,0,0)"><b>st.D</b></u> <- stkh$kst - oprod<br>
<b><u style="color:rgb(255,0,0)">st.R</u> </b><- st.D/(stse)<br> if (!Dzero) {<br> persp(stkh$s, stkh$t, st.D, xlab = "Distance", ylab = "Time", <br> zlab = "D", expand = 0.5, ticktype = "detailed", <br>
theta = -30, shade = 0.4, cex = 0.7)<br> title("D plot")<br> }<br> else {<br> persp(stkh$s, stkh$t, st.D/oprod, xlab = "Distance", <br> ylab = "Time", zlab = "Dzero", expand = 0.5, ticktype = "detailed", <br>
theta = -30, shade = 0.4, cex = 0.7)<br> title("Dzero plot")<br> }<br> plot(outer(stkh$ks, stkh$kt), st.R, xlab = "K(s)K(t)", ylab = "R")<br> abline(h = 0)<br> title("Residual Plot")<br>
if (length(stmc) > 1) {<br> hist(c(stmc$t0, stmc$t), nclass = 50, include.lowest = TRUE, <br> xlab = "Test statistic", main = "")<br> y.75 <- par()$usr[4] * 0.75<br> y.8 <- par()$usr[4] * 0.8<br>
lines(rep(stmc$t0, 2), c(0, y.75), lwd = 2)<br> text(stmc$t0, y.8, "Data Statistic")<br> title("MC results")<br> }<br> else {<br> plot(c(0, 1), c(0, 1), type = "n", axes = FALSE, xlab = "", <br>
ylab = "")<br> ad <- par()$adj<br> par(adj = 0.5)<br> text(0.5, 0.5, "No MC available")<br> par(adj = ad)<br> }<br> par(oldpar)<br>}<br><environment: namespace:splancs><br>
> </i><br><br>Gostaria de obter os valores de <i><span style="color:rgb(204,0,0)"> </span><u style="color:rgb(204,0,0)"><b>st.D</b></u> e <b><u style="color:rgb(255,0,0)">st.R</u> </b>.<br>Desde já agradeço pelo auxilio.<br>
<br><br><br>
</i><br><br>