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