require(geoR) data(ca20) summary(ca20) ## .. analises descriticvas aqui!!! par(mfrow=c(2,2)) plot(variog(ca20, max.dist=500)) plot(variog(ca20, trend=~area,max.dist=500)) plot(variog(ca20, trend=~area+altitude,max.dist=500)) plot(variog(ca20, trend=~area+poly(altitude, 2),max.dist=500)) ## M1 ml1 <- likfit(ca20, ini=c(110, 150), nug=30) ml2 <- likfit(ca20, trend=~area, ini=c(70, 100), nug=30) ml3 <- likfit(ca20, trend=~area+altitude, ini=c(70, 100), nug=30) ml4 <- likfit(ca20, trend=~area+poly(altitude,2), ini=c(60, 100), nug=30) ml1 ml2 ml3 ml4 logLik(ml1) logLik(ml2) logLik(ml3) logLik(ml4) c(logLik(m1),logLik(m3),logLik(m3),logLik(m4)) ## ... mas.... precisa transformar?? ## reajustar com fix.lambda=F ml1l <- likfit(ca20, ini=c(110, 150), nug=30, fix.lam=F) ml2l <- likfit(ca20, trend=~area, ini=c(70, 100), nug=30, fix.lam=F) ml3l <- likfit(ca20, trend=~area+altitude, ini=c(70, 100), nug=30, fix.lam=F) ml4l <- likfit(ca20, trend=~area+poly(altitude,2), ini=c(60, 100), nug=30, fix.lam=F) c(logLik(ml1),logLik(ml3),logLik(ml3),logLik(ml4)) c(logLik(ml1l),logLik(ml2l),logLik(ml3l),logLik(ml4l)) require(MASS) par(mfrow=c(2,2)) boxcox(ca20) boxcox(ca20,trend=~area) boxcox(ca20,trend=~area+altitude) boxcox(ca20,trend=~area+poly(altitude,2)) ##... ##...