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