[R-es] ayuda con dos topicos
PabloEmilio.Verde en uni-duesseldorf.de
PabloEmilio.Verde en uni-duesseldorf.de
Vie Jul 17 23:18:11 CEST 2009
Hola,
El paquete "nlme" tiene implementado algunas estructuras de
covariancia espaciales.
Por ejemplo si utilizas modelos lineales con effectos aleatorios y
estructura de correlacion espacial podes usar "lme". Si usas modelos
lineales sin effectos aleatorios podes usar "glm".
Esto es un ejemplo editado desde el help:
# comienzo del ejemplo ...
fm1 <- lme(weight ~ Time * Diet, BodyWeight,
random = ~ Time)
plot(Variogram(fm1, form = ~ Time))
fm2 <- update(fm1, weights = varPower())
fm3 <- update(fm2,
correlation = corExp(form = ~ Time))
fm4 <- update(fm3,
correlation = corSpher(form = ~ Time))
anova(fm1, fm2, fm3, fm4)
# fin ...
Si miras los "Task views" de R vas a encontrar mucho mas
informacion sobre estadistica espacial.
Los "bins" son el numero de clases en un histograma. No
tengo "The R Book" asi que no se en que contexto usan la
palabra "bins".
Saludos,
Pablo
> Saludos a todos!
>
> Estoy haciendo un semivariograma y me gustaria hacer un ajuste al
> semivariograma, yo quisiera probar varios, uno gausiano, el lineal y uno
> exponencial para ver cual es mejor para mis datos. No se cuales son los
> comandos en R para hacer eso, alguien tiene un scrip que me pueda servir de
> referencia?
>
> La otra pregunta es sobre lo que significa bins?. Estoy siguiendo el libro
> de R (The R Book) para hacer los semivariogramas, pero en la pagina 776,
> aparece esto: the trend surface ..., followed by the numbers of bins....y en
> el libro aparece un 300 que no se de donde sale. Quizas alguno de ustedes
> tenga el libro y pueda darme una ayuda con esto, para entender mejor lo que
> se hizo en el ejemplo y ver como lo hago con mis datos
>
> muchas gracias!!
>
> saludos
>
> MCM
>
>
> --
> Marcia C. Muñoz
> Universidad de Puerto Rico
> Departamento de Biologia - FB246
> POBOX 23360
> San Juan, PR 00931-3360
> tel: 787-7640000, 2900
>
> [[alternative HTML version deleted]]
>
>
Más información sobre la lista de distribución R-help-es