Re: [unmarked] Dimension Issues

[Richard Chandler]

Hi Sarah,

Each of your observation covariates seems to have 70 columns, but your "y" matrix only has 60 columns. The number of columns should be the same. 

Hope this helps.

Richard

On Tue, Oct 23, 2012 at 2:57 PM, Sarah <fritts...@gmail.com> wrote:

I am working with a randomized complete block experimental design where I trapped 1 to 3 blocks at a time.  I am attempting to add "NA" to cells when I did not trap during those time periods.  I added NA's to the abundance cells and to the observation covariate cells.  However, I keep getting an error message that says "Error in unmarkedFrame(y, siteCovs, obsCovs, obsToY = diag(J), mapInfo = mapInfo) :
  At least one matrix in obsCovs has incorrect number of dimensions." 
Currently, all data is in one .csv file (attached). 
Here is the R code I used:

gaca.try<-read.csv("GACA_ABUN_COVS_2012_NA.csv", header=TRUE)
y=gaca.try[,10:69]
siteCovs<-gaca.try[,1:9]
ObsCovs<-list(rain=gaca.try[,c("r.1", "r.2", "r.3", "r.4", "r.5", "r.6", "r.7", "r.8", "r.9",
    "r.10", "r.11", "r.12", "r.13", "r.14", "r.15", "r.16", "r.17", "r.18", "r.19", "r.20", "r.21", "r.22", "r.23",
    "r.24", "r.25", "r.26", "r.27", "r.28", "r.29", "r.30", "r.31", "r.32", "r.33", "r.34", "r.35", "r.36", "r.37",
    "r.38", "r.39", "r.40", "r.41", "r.42", "r.43", "r.44", "r.45", "r.46", "r.47", "r.48", "r.49", "r.50", "r.51",
    "r.52", "r.53", "r.54", "r.55", "r.56", "r.57", "r.58", "r.59", "r.60", "r.61", "r.62", "r.63", "r.64", "r.65",
    "r.66", "r.67", "r.68", "r.69", "r.70")],
    date=gaca.try[,c("m.1", "m.2", "m.3", "m.4", "m.5", "m.6", "m.7", "m.8", "m.9", "m.10",
    "m.11", "m.12", "m.13", "m.14", "m.15", "m.16", "m.17", "m.18", "m.19", "m.20", "m.21", "m.22", "m.23", "m.24",
    "m.25", "m.26", "m.27", "m.28", "m.29", "m.30", "m.31", "m.32", "m.33", "m.34", "m.35", "m.36", "m.37",
    "m.38", "m.39", "m.40", "m.41", "m.42", "m.43", "m.44", "m.45", "m.46", "m.47", "m.48", "m.49", "m.50", "m.51",
    "m.52", "m.53", "m.54", "m.55", "m.56", "m.57", "m.58", "m.59", "m.60", "m.61", "m.62", "m.63", "m.64", "m.65",
    "m.66", "m.67", "m.68", "m.69", "m.70")],
    effort=gaca.try[,c("e.1", "e.2", "e.3", "e.4", "e.5", "e.6", "e.7", "e.8", "e.9", "e.10",
    "e.11", "e.12", "e.13", "e.14", "e.15", "e.16", "e.17", "e.18", "e.19", "e.20", "e.21", "e.22", "e.23", "e.24",
    "e.25", "e.26", "e.27", "e.28", "e.29", "e.30", "e.31", "e.32", "e.33", "e.34", "e.35", "e.36", "e.37",
    "e.38", "e.39", "e.40", "e.41", "e.42", "e.43", "e.44", "e.45", "e.46", "e.47", "e.48", "e.49", "e.50", "e.51",
    "e.52", "e.53", "e.54", "e.55", "e.56", "e.57", "e.58", "e.59", "e.60", "e.61", "e.62", "e.63", "e.64", "e.65",
    "e.66", "e.67", "e.68", "e.69", "e.70")])
gaca.UMF<-unmarkedFramePCount(y=y, siteCovs=siteCovs, obsCovs=ObsCovs)
summary(gaca.UMF)

If I shrink the data by deleting all NA cells, the code works fine.  However, I have different number of trapping periods for each block so I have to delete some data to do that.  Any suggestions?

Thanks,

Sarah