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