Você pode fazer a conta manualmente: list(stD = with(bur1, kst - outer(ks, kt)), stR = with(bur1, kst - outer(ks, kt))/bur1se) (Usando o exemplo do help da função: stdiagn) On Fri, May 25, 2012 at 9:14 AM, Ana Lúcia Souza Silva Mateus <anawilson.ufam@gmail.com> wrote:
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.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O