## GLM Poisson x <- 1:20 b0 <- -2 ; b1=0.3 (lam <- exp(b0 + b1*x)) (y <- rpois(20, lambda = lam)) plot(y ~ x) plot(sqrt(y), x) cbind(x, y) ## modelagem/anĂ¡lise glm(y ~ x, family=poisson) ## GLM Binomial x <- 1:20 b0 <- -4 ; b1=0.3 nu <- b0 + b1*x (p <- exp(nu)/(1+exp(nu))) ## Binomial com n=1 (y <- rbinom(x, size=1, prob=p)) plot(x, y) lines(x, p, lty=2) ## agora com o termo espacial require(geoR) coords <- cbind(runif(50), runif(50)) plot(coords) S <- grf(grid=coords, cov.p=c(1, 0.25)) b0 <- 0 nu <- b0 + S$data (p <- exp(nu)/(1+exp(nu))) ## Binomial com n=1 (y <- rbinom(p, size=1, prob=p)) plot(coords) points(coords[y==1,], pch=19) text(coords[,1], coords[,2], as.character(round(p, dig=2)), cex=0.7, pos=3) ## agora Poisson coords <- cbind(runif(50), runif(50)) plot(coords) S <- grf(grid=coords, cov.p=c(1, 0.25)) b0 <- 0.2 nu <- b0 + S$data (lam <- exp(nu)) ## Binomial com n=1 (y <- rpois(lam, lambda=lam)) plot(coords) text(coords[,1], coords[,2], as.character(y, cex=0.7), pos=3) text(coords[,1], coords[,2], as.character(round(lam, dig=2)), cex=0.7, pos=1) mean(y); var(y)