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