R help - Jan 2009 - Trend.spatial function in geoR

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