全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
1867 5
2013-05-12
周五发在学习版块了,好像是在这里求助更为合适


例如

一个样带包含四个样地

2010年样带物种组成信息为

plotname speciesname
plot1 sp1
plot1 sp2
plot1 sp3
plot1 sp4
plot2 sp2
plot2 sp4
plot2 sp5
plot2 sp6
plot2 sp7
plot3 sp3
plot3 sp5
plot3 sp7
plot3 sp8
plot3 sp9
plot3 sp10
plot3 sp11

2012年物种组成信息为
plotname speciesname
plot1 sp1
plot1 sp3
plot1 sp4
plot1 sp5
plot1 sp6
plot1 sp7
plot2 sp3
plot2 sp4
plot2 sp6
plot2 sp8
plot2 sp11
plot2 sp12
plot2 sp13
plot3 sp4
plot3 sp6
plot3 sp10
plot3 sp14
plot3 sp15
plot3 sp16

如何将2010年物种组成信息转化为

plotname sp1 sp2 sp3 sp4 sp5 sp6 sp7 sp8 sp9 sp10 sp11 sp12 sp13 sp14 sp15 sp16
plot1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0
plot2 0 1 0 1 1 1 1 0 0 0 0 0 0 0 0 0
plot3 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0
将2012年物种组成信息转化为
plotname sp1 sp2 sp3 sp4 sp5 sp6 sp7 sp8 sp9 sp10 sp11 sp12 sp13 sp14 sp15 sp16
plot1         1    0     1     1     1     1    1     0    0    0       0         0      0     0      0       0
plot2 0 0 1 1 0 1 0 1 0 0 1 1 1 0 0
plot3 0 0 0 1 0 1 0 0 0 1 0 0 1 1 1
即两个时期的物种组成矩阵维度相同,包含两时期出现所有物种名称及各物种在各样地的0-1信息?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-5-12 22:45:48
复制代码
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-13 15:17:27
qoiqpwqr 发表于 2013-5-12 22:45
好像是不行。

运行代码后出现如下信息

错误于mat2010[as.character(d2010[i, 1]), as.character(d2010[i, 2])] <- 1 :
  下标出界?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-13 21:28:33
应该不会有错的
复制代码
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-14 23:38:37
复制代码
供参考
> aa = read.table("clipboard",head=T)
> aa
   plotname speciesname
1     plot1         sp1
2     plot1         sp2
3     plot1         sp3
4     plot1         sp4
5     plot2         sp2
6     plot2         sp4
7     plot2         sp5
8     plot2         sp6
9     plot2         sp7
10    plot3         sp3
11    plot3         sp5
12    plot3         sp7
13    plot3         sp8
14    plot3         sp9
15    plot3        sp10
16    plot3        sp11
> aa$num = 1
> aa
   plotname speciesname num
1     plot1         sp1   1
2     plot1         sp2   1
3     plot1         sp3   1
4     plot1         sp4   1
5     plot2         sp2   1
6     plot2         sp4   1
7     plot2         sp5   1
8     plot2         sp6   1
9     plot2         sp7   1
10    plot3         sp3   1
11    plot3         sp5   1
12    plot3         sp7   1
13    plot3         sp8   1
14    plot3         sp9   1
15    plot3        sp10   1
16    plot3        sp11   1
> xtabs(num~.,data=aa)
        speciesname
plotname sp1 sp10 sp11 sp2 sp3 sp4 sp5 sp6 sp7 sp8 sp9
   plot1   1    0    0   1   1   1   0   0   0   0   0
   plot2   0    0    0   1   0   1   1   1   1   0   0
   plot3   0    1    1   0   1   0   1   0   1   1   1

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-5-14 23:39:12
wooson 发表于 2013-5-14 23:38
供参考
> aa = read.table("clipboard",head=T)
> aa
   plotname speciesname
1     plot1         sp1
2     plot1         sp2
3     plot1         sp3
4     plot1         sp4
5     plot2         sp2
6     plot2         sp4
7     plot2         sp5
8     plot2         sp6
9     plot2         sp7
10    plot3         sp3
11    plot3         sp5
12    plot3         sp7
13    plot3         sp8
14    plot3         sp9
15    plot3        sp10
16    plot3        sp11
> aa$num = 1
> aa
   plotname speciesname num
1     plot1         sp1   1
2     plot1         sp2   1
3     plot1         sp3   1
4     plot1         sp4   1
5     plot2         sp2   1
6     plot2         sp4   1
7     plot2         sp5   1
8     plot2         sp6   1
9     plot2         sp7   1
10    plot3         sp3   1
11    plot3         sp5   1
12    plot3         sp7   1
13    plot3         sp8   1
14    plot3         sp9   1
15    plot3        sp10   1
16    plot3        sp11   1
> xtabs(num~.,data=aa)
        speciesname
plotname sp1 sp10 sp11 sp2 sp3 sp4 sp5 sp6 sp7 sp8 sp9
   plot1   1    0    0   1   1   1   0   0   0   0   0
   plot2   0    0    0   1   0   1   1   1   1   0   0
   plot3   0    1    1   0   1   0   1   0   1   1   1
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群