[R-es] Duda interpolación (package ' gstat ')
Marcos Bermejo
markbermejo90 en hotmail.com
Jue Ago 6 17:34:03 CEST 2015
Sale plano sí.
Ya se que sin tener los datos y el código es un poco difícil, pero es que mis datos ocupan mucho, es imposible.
Seguiré mirando por internet.
Muchas gracias Rubén.
Un saludo,
> To: r-help-es en r-project.org
> From: rubenfcasal en gmail.com
> Date: Thu, 6 Aug 2015 14:21:47 +0200
> Subject: Re: [R-es] Duda interpolación (package ' gstat ')
>
> Hola Marcos,
>
> Parece que el problema es con el ajuste del variograma (sale
> plano?), sin más información no se exactamente que puede estar pasando...
>
> Si me envías el código completo y los datos lo miro con más detalle
> (e incluso te doy una alternativa no paramétrica con el paquete npsp).
>
> Un saludo, Rubén.
>
>
>
> El 04/08/2015 a las 11:24, Marcos Bermejo escribió:
> > Hola,
> >
> > # Hacemos el KED. Ver funci�n "krige()":
> > KED.rad <- krige(
> > formula=pluvPcp~layer, # covariable -> radar
> > locations=lluvia.rad.pluv.spdf,
> > newdata=radarGrid, # podr�a ser cualquier objeto Spatial
> > model=v.fit, # modelo de semivariograma.
> > maxdist=Inf
> > )
> >
> > Esta es la funci�n que me interpola los datos de lluvia. El error que me da es:
> >
> > "solve.c", line 88: singular matrix in function Usolve()
> >
> > "lufactor.c", line 208: singular matrix in function m_inverse()
> > Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, :
> > m_inverse
> > In addition: Warning message:
> > In fit.variogram(vg.aux, model = vgm(psill = 0.1, model = "Gau", :
> > Warning: singular model in variogram fit
> >
> >
> > Mi funci�n del variograma es :
> > v.fit <- fit.variogram(vg.aux, model=vgm(psill=0.15, model='Gau', range=5000,
> > nugget=0.05))
> >
> > �Alguien me podr�a ayudar?
> >
> > Gracias de antemano.
> >
> > Un saludo,
> >
> > [[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
>
>
> [[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
[[alternative HTML version deleted]]
Más información sobre la lista de distribución R-help-es