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@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
