是否有辦法將觀察到的魚的原始數量（"num"）和每個環境資料列（"temp"、"do "等）之間的相關統計數字按物種（"group"）制成一個檔案？
*以及num與env.因素的平均值和中位數之間的相關性？
我還希望能夠選擇使用哪種相關方法（Pearson相關、Kendall等級相關、Spearman相關等）
。我的資料：
zeros < -結構(list(年 =結構(c(1L。 1L。 1L, 1L, 1L。 1L, 1L, 1L,
1L,/span> 1L。 1L, 1L。 1L, 1L。 1L, 1L。 1L, 1L。 1L, 1L。 1L, 2L。 2L, 2L, 2L,
2L,/span> 2L。 2L, 2L。 2L, 2L。 2L, 2L。 2L, 2L。 2L, 2L。 2L, 2L。 2L, 2L, 2L,
2L), 。 標簽 = c("2019"/span>。 "2020"）。 class = " factor"）。 季節=結構(c(1L。
1L, 1L。 1L, 1L。 1L, 1L。 1L, 1L。 1L, 2L。 2L, 2L。 2L, 2L。 2L, 2L, 2L,
2L,/span> 2L。 2L, 1L。 1L, 1L。 1L, 1L。 1L, 1L。 1L, 1L。 1L, 2L。 2L, 2L, 2L,
2L,/span> 2L。 2L, 2L。 2L, 2L。 2L), 。 標簽 = c("dry"。 "wet"）。 class = " factor"),
網站=結構(c(1L。 1L, 2L, 2L。 3L, 3L, 4L。 4L, 5L。 5L,
1L, 1L。 2L, 2L。 3L, 3L。 4L, 4L。 5L, 5L。 1L, 1L。 2L, 2L。 3L,
3L,/span> 4L。 4L, 5L。 5L, 1L, 1L。 2L, 2L。 3L, 3L。 4L, 4L。 5L,/span> 5L
), 。 標簽 = c("1"/span>。 "2", "3"。 "4", "5"）。 class = " factor"),
組=結構(c(1L。 2L, 1L, 2L。 1L, 2L, 1L。 2L, 1L。 2L, 2L,
1L, 2L。 1L, 2L。 1L, 2L。 1L, 2L。 1L, 2L。 1L, 2L。 1L, 2L。 1L,
2L, 1L。 2L, 1L。 2L, 1L, 2L。 1L, 2L。 1L, 2L。 1L, 2L。 1L, 2L
), 。 標簽 = c("hardhead silverside"。 "Sailfin molly"）。 class = " factor"),
num = c(0。 8, 0。 9, 0。 13, 0。 9, 0。 10, 0, 0。 0, 0。 0, 0, 0,
0, 0,, 0, 0, 0, 0, 0, 0, 0。 0, 0。 0, 0, 4。 0,/span> 7。 0, 2,
0,/span> 3。 0,/span> 13。 0）。 temp = c(23L。 36L, 35L, 34L。 30L, 28L。 18L,
19L, 33L,, 25L, 20L, 33L。 23L, 36L, 32L。 28L, 17L。 34L,
31L, 26L， 34L。 26L, 35L, 15L。 25L, 26L, 20L。 18L, 14L。 23L,
17L, 26L,。 17L, 19L, 29L。 31L, 18L。 15L）。 sal = c（12. 5,。
25.5, 8.5, 15. 5, 17.5, 27.5, 9. 5, 31.5, 1。 5, 34.5, 25.5,
21.5, 10.5, 8. 5, 32.5, 19.5, 6． 5, 5.5, 15.5, 28. 5, 6.5, 3.5,
29.5, 13.5, 7． 5, 16.5, 3.5, 28. 5, 22.5, 5.5, 9. 5, 12.5, 29.5,
24.5, 8.5, 32. 5, 37.5, 3.5, 12。 5, 19. 5),做=c（9. 66, 7.66,
1.66, 14.66, 15. 66, 1.66, 14. 66, 15.66, 0. 66, 5.66, 10.66,
11.66, 4.66, 0。 66, 13.66, 1. 66, 13.66, 6. 66, 6.66, 10.66,
9.66, 15.66, 9. 66, 15.66, 4. 66, 13.66, 1. 66, 11.66, 6.66,
8.66, 12.66, 0。 66, 6.66, 0. 66, 9.66, 16. 66, 1.66, 10.66,
15.66, 10。 66）。 深度= c（120L。 161L。 52L, 52L。 43L, 105L,
165L, 23L, 79L, 136L, 41L, 59L。 65L, 118L, 122L。 69L, 137L,
88L, 152L, 105L。 108L, 79L， 96L。 80L, 22L, 110L。 157L, 118L,
126L, 93L, 156L, 64L, 74L, 24L。 111L, 113L, 157L。 78L, 121L,
130L））。 class = "data. frame",行。 names = c(NA。 -40L））
uj5u.com熱心網友回復：
你問題的第一部分很直接：
zeros. <- split(zeros, zeros$group)
zeros.cors <- sapply(zeros. spl。 function(x) cor(x[。 "num"]。 x[, 6。 9]））
dimnames(zeros. cors)[[1]] < -colnames(zeros)[6。 9]
zeros.cors
# Hardhead silverside Sailfin molly
# temp -0.3080334 0.36174046
# sal 0.1393580 0.47095129;/span>
# does 0.2544695 -0.06646818
# depth 0.1296208 0.08777425
t（zeros.cors）
# temp sal do depth
# Hardhead silverside -0.3080334 0.1393580 0.25446948 0.12962078。
# Sailfin molly 0.3617405 0.4709513 -0.06646818 0.08777425。
使用write.csv(zeros.cors, file="results.csv")或write.csv(t(zeros.cors), file="results.csv")取決于你希望行/列是什么。
第二個問題不清楚。一個組的平均值/中間值將是一個單一的值，所以你不能把它與環境變數相關聯。你可以通過聚合來計算各組的平均值：
aggregate(zeros[。 5。 9]/span>。 by=list(zeros$group)。 "mean")
# Group.1 num temp sal do depth
# 1 Hardhead silverside 1.45 25.95 15.35 8.51 105.20
# 2 Sailfin molly 2.45 25.00 18.90 9.06 90.25
總計(zeros[, 5。 9]/span>。 by=list(zeros$group)。 "median"）
# Group.1 num temp sal do depth
# 1 Hardhead silverside 0 26 11.5 9.66 115.5
# 2 Sailfin molly 0 24 19.5 10.66 90.5
轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/312344.html
標籤：