[R] vgam cao

[romi at aon.at]

Hi!
I would like to calculate a constrained additive ordination with the package VGAM. I get the error message: 
Error in ncol(xbig.save) : object "xbig.save" not found

My data consists of 90 samples, 10 species and 3 environmental variables. The species data is continous biomass data, with values from 0 to 3000 micrograms per litre. 

I use the function:
cao (name.of.species.matrix~Temp+LF+Ptot, data=name.of.environmental.data.frame, family=gaussianff, df1.nl=3.5, Bestof=4)

After fitting model 1, I get an error message. Here is the output of model 1:

Obtaining initial values

Using initial values
    Temp    LF  Ptot
lv 0.408 0.189 -0.77

Iteration 1 
initial  value 88934920.105286 
iter  10 value 83029461.332581
iter  20 value 83028498.211367
final  value 83028488.241532 
converged
Error in ncol(xbig.save) : object "xbig.save" not found

I have tried it with fewer or more sites, fewer and more species, df1.nl from 0.5 to 5.5. I always get the same error message.
It only works when I use family=poissonff. But as I have continuous data, I think I can't use poissonff.

Can anybody help me?
Thanks,
Romana




-------------------------------------------
Versendet durch aonWebmail (webmail.aon.at)