我正在尝试使用包RMark中的convert.inp()读取数据集。
我的数据集包含个人ID号,重新捕获的历史记录(即0101..00101 ...),位置,物种,年龄,性别,体重和寄生虫负担。我将数据格式化为.inp,将站点作为组-紧随重获历史之后-将物种作为第二个.inp文件中的组。这两个文件的列标题和唯一ID号都有/* */
条注释。站点作为组文件可以很好地读取22个变量的213个观察值。 “物种作为组”文件会产生以下错误:
Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, :
line 2 did not have 37 elements
如果删除/* */
列标题注释,该文件将被读入,但仅包含13个变量的211个观察值。
为什么我首先遇到此错误?为什么只通过删除列标题在第二个文件中观察到的内容更少?
这是我尝试过的代码:
将文件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)
*由于胁迫而发出警告NAs-我缺少1个人的体重
将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
答案 0 :(得分:0)
将转换后的数据帧导出为.csv后,我能够比较物种组文件和源文件中的每一行,发现在物种代码的所有字段中有2个2个人有0。