Não foi possível enviar o arquivo. Será algum problema com as permissões?
Essa é uma revisão anterior do documento!
## comandos de list xydata <- data.frame(x = 1:20, y = c(0.9,5.1,2.4,8.1,4.2,7.1,5.6,7.6,5.9,7.7, 11.8,6.9,9.3,10.9,8.4,11.6,13.0,13.8,13.1,9.3) ) with(xydata, plot(x,y)) attach(xydata) plot(x,y) detach(xydata) plot(xydata$x, xydata$y) reg <- lm(y ~ x, data = xydata) is.list(reg) names(reg) reg$call summary(reg) regS <- summary(reg) is.list(regS) names(regS) sqrt(sum(reg$residuals^2)/reg$df.res) regS$sigma ## dg <- c(9, 5, 2, 3, 0, 3, 2, 4, 14, 3, 4, 1, 0, 6, 1) fp <- function(p) sum(log(p * (1-p)^dg)) ps <- seq(0.01, 0.99, length=100) lps <- sapply(ps, fp) plot(ps, lps , ty="l") est <- optimize(fp, c(0, 1), maximum=T) est abline(v = est$max) fp <- function(par) sapply(par, function(p) sum(log(p) + dg*log(1-p))) plot(ps, fp(ps), ty="l") fp <- fp(par) sapply(par, function(p) dgeom(dg, prob=p, log=T))