Gracias por la respuestas!

Yo estoy usando library(geoR) y tengo el semivariograma ya, el problema es
que no me lo da ajustado

Los comandos para hacer el ajuste los reconozco (e.g. corSpher o corRatio),
no se como escribir la formula...ahi es donde me siento perdida.

Gracias por su ayuda

saludos

MCM

El 17 de julio de 2009 18:11, Kjetil Halvorsen <kjetil1001@gmail.com>escribió:

> 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@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@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
>



-- 
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]]

