[R-es] ayuda con dos topicos
Kjetil Halvorsen
kjetil1001 en gmail.com
Sab Jul 18 00:11:39 CEST 2009
Hola!
Para dar buen ayuda necesitas en realidad explicar algo más lo que
quieres hacer.
la cosa es que , en R, o mejor, en CRAN, hay muchas maneras alternativas para
estimar (semi)variogramas.
Un paqute útil es gstat
install.packages("gstat", dep=TRUE)
library(gstat)
library(help=gstat)
vignette("gstat")
?fit.variogram
Una alternativa interesante, especialmente si vas a trabajar con GLM's
en el contexto
espacial, es geoR:
install.packages(c("geoR","geoRglm"), dep=TRUE)
Si vas a trabajar con conjuntos de datso muy grandes (más que un poco
de centenares
de puntos), puede ser de auyda fields:
install.packages("fields", dep=TRUE)
En adición hay otras alternativas!!
Kjetil
2009/7/17 <PabloEmilio.Verde en uni-duesseldorf.de>:
> 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]]
>>
>>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
--
"... an entire human genome would fit on a music CD."
--- www.thinkgene.com
Más información sobre la lista de distribución R-help-es