Я запустил анализ разделения дисперсии, используя r, и мне нужно найти значение значимости (p-значение) для каждого компонента в ResultTable?
Вот коды, которые я использовал до сих пор:
#first make 4 tables of explanatory variables
Nutrients=c("DIN","DRP")
Sediments=c("SIS","SOS")
Habitat=c("RHA")
thevarpart = varpart(Y = All.data[, "MCI"],
log10(All.data[, Nutrients]),
All.data[,Sediments],
All.data[,Habitat])
x11(); plot(thevarpart, cutoff = 0, digits =
1,bg=c("hotpink","skyblue","yellow"),Xnames=c('Nutrients','Fine Sediment','Habitat'))
x11(); showvarparts(3,bg=c("hotpink","skyblue","yellow"))
ResultsTable <- rbind.data.frame(thevarpart$part$fract, thevarpart$part$indfract); ResultsTable
Мои данные выглядят так:
MCI SIS SOS DIN DRP RHA
118 282.5595 23.90908052 0.0335 0.002 60
90 315.9466376 32.12509173 0.675 0.006 55
97 573.8064879 71.16698102 0.9655 0.026 53.5
119 123.8694668 28.89448716 0.03 0.002 72
120 39.27030208 9.465458598 0.0175 0.002 71
97.5 475.0791544 103.2212162 0.915 0.0145 52
120 588.576473 107.5916448 0.21 0.01 80.5
101 319.4977937 45.97594364 0.465 0.0045 48.5
89 615.1995445 85.89353893 0.6475 0.0105 53
120 349.7319968 28.82643698 1.546 0.002 70.5
102 381.8583689 44.76326487 0.078 0.004 51
87 656.9546267 87.62515067 3.62 0.032 49.5
97 1186.128223 116.1077253 0.7205 0.012 48.5
83 330.235734 32.63877857 1.1915 0.014 51
106 121.9643987 12.22626485 0.385 0.002 63.5
99 528.7900893 24.33578253 0.93 0.006 51
90 1546.765404 126.9419995 0.915 0.0085 56
115 625.869261 31.09234308 0.37 0.005 59.5
97 516.7899503 70.82898363 0.815 0.011 48.5
121 446.7218383 116.7642469 0.1585 0.011 83
74 1245.867382 151.7759208 0.302 0.071 48.5
100 641.1141972 81.31781456 1.0575 0.007 58
93 518.8630661 86.74102689 0.51 0.009 49
119.5 170.5972465 15.07056712 0.63 0.002 61
117 80.33983969 10.6866821 0.038 0.002 55
98 541.6111975 56.74776662 0.972 0.005 61
101 600.0329388 68.70973658 0.832 0.0095 58
95 784.6798499 116.2767808 0.919 0.0235 53.5
100 1071.666097 108.9452369 1.0485 0.021 47
83 1385.849876 112.2573191 2.11 0.029 47
78 796.8454586 99.79436559 1.531 0.036 41
75 1573.820802 196.5013325 1.568 0.0275 37