I am having difficulty getting the trend.spatial function in geoR to
work properly. After creating a trend.spatial object with a covariate,
I try to add the command into my likfit() function as follows:
trend1.trend.spatial <- trend.spatial("1st", trend1.geodata)
trend1.spatial.EC0.1.reml <- likfit(spatial.geodata,
trend1.trend.spatial, ini.cov.pars = spatial.EC0.1.eyefit, fix.kappa =
FALSE, fix.psiA = FALSE, fix.psiR = FALSE,
+ lik.method = "REML", components = FALSE)
I receive an error stating: trend matrix has dimension incompatible with
the data
My data files have an equal number of rows or spatial locations. Any
suggestions would be useful.
Also, is there a way to add multiple covariates into the trend.spatial
function?
Thanks,
*************************
Julia L. Angstmann
Department of Botany
University of Wyoming
1000 E. University Ave.
Laramie, WY 82071
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
