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]]
PabloEmilio.Verde en uni-duesseldorf.de
2009-Jul-17 21:18 UTC
[R-es] ayuda con dos topicos
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]] > >