扫描(文件...在convert.inp中读取时的错误和不同的观察数量

问题描述 投票:0回答:1

我试图在包RMark中使用convert.inp()读取数据集。

我的数据集包含个人ID号,重新获取历史(即0101..00101 ......),站点,物种,年龄,性别,体重和寄生虫负担。我将数据格式化为.inp,网站作为组 - 紧跟重新捕获历史记录 - 并将物种作为第二个.inp文件中的组。这两个文件都有列标题的/* */注释和唯一的ID号。作为组文件的站点读取了213个22个变量的观察结果。作为组文件的Species会产生以下错误:

Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec,  : 
line 2 did not have 37 elements

如果我删除/* */列标题注释,该文件将读入但只有211个观察13个变量。

为什么我首先得到这个错误?为什么我只需删除列标题就可以减少第二个文件中的观察次数?

这是我尝试过的代码:

文件1网站作为组:

SiteGroup <- convert.inp("SiteGroup.inp", group.df = data.frame(Site= c("Albany Pine Bush", "SChroon Lake")), covariates = c("Peromyscus", "Tamias striatus", "Myodes gapperi", "Blarina brevicauda", "Sorex cinereus", "Tamiasciurus hudsonicus", "Glacomys volans", "Glacomys sabrinus", "Napeaozapus insignis", "Mustela frenata", "Sorex hoyi", "Male", "Female", "Juvenile", "Adult", "Weight (g)", "dWeight (g)", "Length (cm)", "Avg Tick Burden"), use.comments = TRUE)

*通过coersion引入警告NAs - 我缺少一个人的体重

将2种文件作为具有列标题的组:

SpeciesGroup <- convert.inp("SpeciesGroup.inp", group.df = data.frame(Species = c("Peromyscus", "Tamias striatus", "Myodes gapperi", "Blarina brevicauda", "Sorex cinereus", "Tamiasciurus hudsonicus", "Glacomys volans", "Glacomys sabrinus", "Napeaozapus insignis", "Mustela frenata", "Sorex hoyi")), covariates = c("Albany Pine Bush", "Schroon Lake", "Male", "Female", "Juvenile", "Adult", "Weight (g)", "dWeight (g)", "Length (cm)", "Avg Tick Burden"), use.comments = TRUE)

文件2输出:

Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec,  : 
line 2 did not have 37 elements
r file rmark
1个回答
0
投票

在将转换后的数据框导出为.csv之后,我能够将物种组文件中的每一行与源文件进行比较,发现2个人在物种代码的所有字段中都有0。

© www.soinside.com 2019 - 2024. All rights reserved.