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))