交互表 - 與寵物和房屋的案例我有一份房屋清單和一份動(dòng)物物種清單。houses = c(1,1,2,3,4,4,4,4,5,6,5)animals = c('cat','dog','cat','dog','rat', 'cat', 'spider', 'snake', 'cat', 'cat', 'rat')我正在嘗試創(chuàng)建一個(gè)函數(shù)，該函數(shù)返回一個(gè)上三角表，指示每只寵物，它被觀(guān)察到的生活在同一個(gè)房子中的次數(shù)比其他動(dòng)物物種。是否有意義？對(duì)于上面的例子，表應(yīng)該是這樣的（希望沒(méi)有錯(cuò)?。?nbsp; dog rat spider snakecat 1 2 1 1 dog 0 0 0rat 1 1spider 1注意：此函數(shù)應(yīng)適用于任何兩個(gè)相同長(zhǎng)度的向量，無(wú)論它們是否包含數(shù)字或字符串
2 回答
料青山看我應(yīng)如是
TA貢獻(xiàn)1772條經(jīng)驗(yàn) 獲得超8個(gè)贊
使用table和crossprod：
out <- crossprod(table(houses, animals))
out[lower.tri(out, diag=TRUE)] <- NA
out
# animals
# animals cat dog rat snake spider
# cat NA 1 2 1 1
# dog NA NA 0 0 0
# rat NA NA NA 1 1
# snake NA NA NA NA 1
# spider NA NA NA NA NA
由于輸出是a，matrix您可以NA直接在print以下位置禁止打印值：
print(out,na.print="")
# animals
# animals cat dog rat snake spider
# cat 1 2 1 1
# dog 0 0 0
# rat 1 1
# snake 1
# spider
梵蒂岡之花
TA貢獻(xiàn)1900條經(jīng)驗(yàn) 獲得超5個(gè)贊
只需添加一個(gè)稀疏的替代方案，以防萬(wàn)一在制表期間可以減少任何內(nèi)存占用：Matrix::crossprod(xtabs(~ houses + animals, sparse = TRUE))
- 2 回答
- 0 關(guān)注
- 579 瀏覽
添加回答
舉報(bào)
0/150
提交
取消