data<-read.table("S20c.csv", header=T, sep=",")
contcat <- lm(ln.carbon~dÃas*gen,data=data)
summary(contcat)
##
## Call:
## lm(formula = ln.carbon ~ dÃas * gen, data = data)
##
## Residuals:
## Min 1Q Median 3Q Max
## -1.18079 -0.16788 0.01012 0.17739 0.75831
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.464e+00 1.962e-01 -7.462 1.19e-13 ***
## dÃas 2.398e-03 4.480e-04 5.353 9.50e-08 ***
## genTCS02 2.331e-01 2.097e-01 1.112 0.2663
## genTCS03 2.744e-01 2.092e-01 1.312 0.1898
## genTCS04 4.404e-02 2.099e-01 0.210 0.8339
## genTCS05 3.131e-01 2.092e-01 1.496 0.1347
## genTCS08 9.577e-02 2.175e-01 0.440 0.6597
## genTCS10 2.626e-01 2.092e-01 1.255 0.2096
## genTCS11 2.582e-01 2.092e-01 1.234 0.2173
## genTCS12 4.529e-01 2.092e-01 2.164 0.0305 *
## genTCS20 5.363e-01 2.092e-01 2.563 0.0104 *
## genTCS40 -4.172e-01 2.407e-01 -1.733 0.0831 .
## genTCS41 -3.718e-01 2.904e-01 -1.280 0.2006
## genTCS42 -2.812e-01 2.680e-01 -1.049 0.2942
## genTCS43 3.368e-01 2.108e-01 1.598 0.1102
## genTCS44 4.087e-01 2.095e-01 1.951 0.0512 .
## genTCS45 4.824e-01 2.097e-01 2.300 0.0215 *
## genTCS46 1.881e-01 2.092e-01 0.899 0.3688
## genTCS47 3.308e-01 2.092e-01 1.581 0.1140
## genTCS48 4.181e-01 2.095e-01 1.995 0.0461 *
## genTCS49 2.867e-01 2.092e-01 1.370 0.1707
## genTSC45 4.736e-01 4.267e-01 1.110 0.2671
## dÃas:genTCS02 -4.185e-04 4.913e-04 -0.852 0.3944
## dÃas:genTCS03 -3.958e-05 4.887e-04 -0.081 0.9354
## dÃas:genTCS04 8.549e-05 4.901e-04 0.174 0.8615
## dÃas:genTCS05 -9.895e-05 4.887e-04 -0.202 0.8395
## dÃas:genTCS08 -1.660e-04 5.129e-04 -0.324 0.7462
## dÃas:genTCS10 -1.521e-04 4.887e-04 -0.311 0.7556
## dÃas:genTCS11 1.263e-04 4.887e-04 0.258 0.7961
## dÃas:genTCS12 -4.334e-04 4.887e-04 -0.887 0.3753
## dÃas:genTCS20 -3.179e-04 4.887e-04 -0.651 0.5154
## dÃas:genTCS40 9.929e-04 5.570e-04 1.783 0.0748 .
## dÃas:genTCS41 7.254e-04 6.599e-04 1.099 0.2718
## dÃas:genTCS42 3.941e-04 6.132e-04 0.643 0.5205
## dÃas:genTCS43 -3.685e-04 4.940e-04 -0.746 0.4557
## dÃas:genTCS44 -4.214e-04 4.893e-04 -0.861 0.3892
## dÃas:genTCS45 -9.121e-04 4.901e-04 -1.861 0.0629 .
## dÃas:genTCS46 9.485e-05 4.887e-04 0.194 0.8461
## dÃas:genTCS47 -4.283e-04 4.887e-04 -0.876 0.3809
## dÃas:genTCS48 -1.899e-04 4.893e-04 -0.388 0.6980
## dÃas:genTCS49 -1.748e-04 4.887e-04 -0.358 0.7207
## dÃas:genTSC45 -2.163e-04 1.141e-03 -0.190 0.8497
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.2674 on 2371 degrees of freedom
## Multiple R-squared: 0.5068, Adjusted R-squared: 0.4983
## F-statistic: 59.42 on 41 and 2371 DF, p-value: < 2.2e-16
anova(contcat)
## Analysis of Variance Table
##
## Response: ln.carbon
## Df Sum Sq Mean Sq F value Pr(>F)
## dÃas 1 132.873 132.873 1858.3071 < 2.2e-16 ***
## gen 20 38.029 1.901 26.5928 < 2.2e-16 ***
## dÃas:gen 20 3.293 0.165 2.3031 0.0008607 ***
## Residuals 2371 169.532 0.072
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
library(emmeans)
emtrends(contcat, ~ gen, var="dÃas")
## gen dÃas.trend SE df lower.CL upper.CL
## TCS01 0.00240 0.000448 2371 0.001519 0.00328
## TCS02 0.00198 0.000202 2371 0.001584 0.00237
## TCS03 0.00236 0.000195 2371 0.001975 0.00274
## TCS04 0.00248 0.000199 2371 0.002093 0.00287
## TCS05 0.00230 0.000195 2371 0.001916 0.00268
## TCS08 0.00223 0.000250 2371 0.001742 0.00272
## TCS10 0.00225 0.000195 2371 0.001863 0.00263
## TCS11 0.00252 0.000195 2371 0.002141 0.00291
## TCS12 0.00196 0.000195 2371 0.001581 0.00235
## TCS20 0.00208 0.000195 2371 0.001697 0.00246
## TCS40 0.00339 0.000331 2371 0.002742 0.00404
## TCS41 0.00312 0.000485 2371 0.002173 0.00407
## TCS42 0.00279 0.000419 2371 0.001971 0.00361
## TCS43 0.00203 0.000208 2371 0.001621 0.00244
## TCS44 0.00198 0.000197 2371 0.001590 0.00236
## TCS45 0.00149 0.000199 2371 0.001096 0.00188
## TCS46 0.00249 0.000195 2371 0.002110 0.00288
## TCS47 0.00197 0.000195 2371 0.001586 0.00235
## TCS48 0.00221 0.000197 2371 0.001822 0.00259
## TCS49 0.00222 0.000195 2371 0.001840 0.00261
## TSC45 0.00218 0.001050 2371 0.000123 0.00424
##
## Confidence level used: 0.95
emtrends(contcat, pairwise ~ gen, var="dÃas")
## $emtrends
## gen dÃas.trend SE df lower.CL upper.CL
## TCS01 0.00240 0.000448 2371 0.001519 0.00328
## TCS02 0.00198 0.000202 2371 0.001584 0.00237
## TCS03 0.00236 0.000195 2371 0.001975 0.00274
## TCS04 0.00248 0.000199 2371 0.002093 0.00287
## TCS05 0.00230 0.000195 2371 0.001916 0.00268
## TCS08 0.00223 0.000250 2371 0.001742 0.00272
## TCS10 0.00225 0.000195 2371 0.001863 0.00263
## TCS11 0.00252 0.000195 2371 0.002141 0.00291
## TCS12 0.00196 0.000195 2371 0.001581 0.00235
## TCS20 0.00208 0.000195 2371 0.001697 0.00246
## TCS40 0.00339 0.000331 2371 0.002742 0.00404
## TCS41 0.00312 0.000485 2371 0.002173 0.00407
## TCS42 0.00279 0.000419 2371 0.001971 0.00361
## TCS43 0.00203 0.000208 2371 0.001621 0.00244
## TCS44 0.00198 0.000197 2371 0.001590 0.00236
## TCS45 0.00149 0.000199 2371 0.001096 0.00188
## TCS46 0.00249 0.000195 2371 0.002110 0.00288
## TCS47 0.00197 0.000195 2371 0.001586 0.00235
## TCS48 0.00221 0.000197 2371 0.001822 0.00259
## TCS49 0.00222 0.000195 2371 0.001840 0.00261
## TSC45 0.00218 0.001050 2371 0.000123 0.00424
##
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## TCS01 - TCS02 4.19e-04 0.000491 2371 0.852 1.0000
## TCS01 - TCS03 3.96e-05 0.000489 2371 0.081 1.0000
## TCS01 - TCS04 -8.55e-05 0.000490 2371 -0.174 1.0000
## TCS01 - TCS05 9.90e-05 0.000489 2371 0.202 1.0000
## TCS01 - TCS08 1.66e-04 0.000513 2371 0.324 1.0000
## TCS01 - TCS10 1.52e-04 0.000489 2371 0.311 1.0000
## TCS01 - TCS11 -1.26e-04 0.000489 2371 -0.258 1.0000
## TCS01 - TCS12 4.33e-04 0.000489 2371 0.887 1.0000
## TCS01 - TCS20 3.18e-04 0.000489 2371 0.651 1.0000
## TCS01 - TCS40 -9.93e-04 0.000557 2371 -1.783 0.9725
## TCS01 - TCS41 -7.25e-04 0.000660 2371 -1.099 1.0000
## TCS01 - TCS42 -3.94e-04 0.000613 2371 -0.643 1.0000
## TCS01 - TCS43 3.69e-04 0.000494 2371 0.746 1.0000
## TCS01 - TCS44 4.21e-04 0.000489 2371 0.861 1.0000
## TCS01 - TCS45 9.12e-04 0.000490 2371 1.861 0.9573
## TCS01 - TCS46 -9.49e-05 0.000489 2371 -0.194 1.0000
## TCS01 - TCS47 4.28e-04 0.000489 2371 0.876 1.0000
## TCS01 - TCS48 1.90e-04 0.000489 2371 0.388 1.0000
## TCS01 - TCS49 1.75e-04 0.000489 2371 0.358 1.0000
## TCS01 - TSC45 2.16e-04 0.001141 2371 0.190 1.0000
## TCS02 - TCS03 -3.79e-04 0.000281 2371 -1.349 0.9991
## TCS02 - TCS04 -5.04e-04 0.000283 2371 -1.779 0.9731
## TCS02 - TCS05 -3.20e-04 0.000281 2371 -1.138 0.9999
## TCS02 - TCS08 -2.53e-04 0.000321 2371 -0.787 1.0000
## TCS02 - TCS10 -2.66e-04 0.000281 2371 -0.949 1.0000
## TCS02 - TCS11 -5.45e-04 0.000281 2371 -1.940 0.9363
## TCS02 - TCS12 1.49e-05 0.000281 2371 0.053 1.0000
## TCS02 - TCS20 -1.01e-04 0.000281 2371 -0.358 1.0000
## TCS02 - TCS40 -1.41e-03 0.000388 2371 -3.641 0.0399
## TCS02 - TCS41 -1.14e-03 0.000525 2371 -2.180 0.8325
## TCS02 - TCS42 -8.13e-04 0.000465 2371 -1.748 0.9777
## TCS02 - TCS43 -5.00e-05 0.000290 2371 -0.172 1.0000
## TCS02 - TCS44 2.92e-06 0.000282 2371 0.010 1.0000
## TCS02 - TCS45 4.94e-04 0.000283 2371 1.742 0.9785
## TCS02 - TCS46 -5.13e-04 0.000281 2371 -1.828 0.9643
## TCS02 - TCS47 9.78e-06 0.000281 2371 0.035 1.0000
## TCS02 - TCS48 -2.29e-04 0.000282 2371 -0.811 1.0000
## TCS02 - TCS49 -2.44e-04 0.000281 2371 -0.868 1.0000
## TCS02 - TSC45 -2.02e-04 0.001069 2371 -0.189 1.0000
## TCS03 - TCS04 -1.25e-04 0.000279 2371 -0.449 1.0000
## TCS03 - TCS05 5.94e-05 0.000276 2371 0.215 1.0000
## TCS03 - TCS08 1.26e-04 0.000317 2371 0.399 1.0000
## TCS03 - TCS10 1.13e-04 0.000276 2371 0.407 1.0000
## TCS03 - TCS11 -1.66e-04 0.000276 2371 -0.601 1.0000
## TCS03 - TCS12 3.94e-04 0.000276 2371 1.426 0.9982
## TCS03 - TCS20 2.78e-04 0.000276 2371 1.008 1.0000
## TCS03 - TCS40 -1.03e-03 0.000384 2371 -2.686 0.4619
## TCS03 - TCS41 -7.65e-04 0.000522 2371 -1.464 0.9974
## TCS03 - TCS42 -4.34e-04 0.000462 2371 -0.939 1.0000
## TCS03 - TCS43 3.29e-04 0.000286 2371 1.152 0.9999
## TCS03 - TCS44 3.82e-04 0.000277 2371 1.377 0.9989
## TCS03 - TCS45 8.72e-04 0.000279 2371 3.130 0.1811
## TCS03 - TCS46 -1.34e-04 0.000276 2371 -0.487 1.0000
## TCS03 - TCS47 3.89e-04 0.000276 2371 1.407 0.9984
## TCS03 - TCS48 1.50e-04 0.000277 2371 0.542 1.0000
## TCS03 - TCS49 1.35e-04 0.000276 2371 0.489 1.0000
## TCS03 - TSC45 1.77e-04 0.001068 2371 0.166 1.0000
## TCS04 - TCS05 1.84e-04 0.000279 2371 0.662 1.0000
## TCS04 - TCS08 2.51e-04 0.000319 2371 0.788 1.0000
## TCS04 - TCS10 2.38e-04 0.000279 2371 0.852 1.0000
## TCS04 - TCS11 -4.08e-05 0.000279 2371 -0.146 1.0000
## TCS04 - TCS12 5.19e-04 0.000279 2371 1.861 0.9572
## TCS04 - TCS20 4.03e-04 0.000279 2371 1.447 0.9977
## TCS04 - TCS40 -9.07e-04 0.000386 2371 -2.350 0.7226
## TCS04 - TCS41 -6.40e-04 0.000524 2371 -1.222 0.9998
## TCS04 - TCS42 -3.09e-04 0.000464 2371 -0.666 1.0000
## TCS04 - TCS43 4.54e-04 0.000288 2371 1.577 0.9932
## TCS04 - TCS44 5.07e-04 0.000280 2371 1.811 0.9675
## TCS04 - TCS45 9.98e-04 0.000281 2371 3.547 0.0544
## TCS04 - TCS46 -9.36e-06 0.000279 2371 -0.034 1.0000
## TCS04 - TCS47 5.14e-04 0.000279 2371 1.843 0.9612
## TCS04 - TCS48 2.75e-04 0.000280 2371 0.984 1.0000
## TCS04 - TCS49 2.60e-04 0.000279 2371 0.934 1.0000
## TCS04 - TSC45 3.02e-04 0.001068 2371 0.282 1.0000
## TCS05 - TCS08 6.70e-05 0.000317 2371 0.211 1.0000
## TCS05 - TCS10 5.32e-05 0.000276 2371 0.193 1.0000
## TCS05 - TCS11 -2.25e-04 0.000276 2371 -0.816 1.0000
## TCS05 - TCS12 3.34e-04 0.000276 2371 1.211 0.9998
## TCS05 - TCS20 2.19e-04 0.000276 2371 0.793 1.0000
## TCS05 - TCS40 -1.09e-03 0.000384 2371 -2.841 0.3483
## TCS05 - TCS41 -8.24e-04 0.000522 2371 -1.578 0.9932
## TCS05 - TCS42 -4.93e-04 0.000462 2371 -1.067 1.0000
## TCS05 - TCS43 2.70e-04 0.000285 2371 0.944 1.0000
## TCS05 - TCS44 3.22e-04 0.000277 2371 1.163 0.9999
## TCS05 - TCS45 8.13e-04 0.000279 2371 2.918 0.2975
## TCS05 - TCS46 -1.94e-04 0.000276 2371 -0.702 1.0000
## TCS05 - TCS47 3.29e-04 0.000276 2371 1.193 0.9999
## TCS05 - TCS48 9.09e-05 0.000277 2371 0.328 1.0000
## TCS05 - TCS49 7.58e-05 0.000276 2371 0.274 1.0000
## TCS05 - TSC45 1.17e-04 0.001068 2371 0.110 1.0000
## TCS08 - TCS10 -1.39e-05 0.000317 2371 -0.044 1.0000
## TCS08 - TCS11 -2.92e-04 0.000317 2371 -0.922 1.0000
## TCS08 - TCS12 2.67e-04 0.000317 2371 0.843 1.0000
## TCS08 - TCS20 1.52e-04 0.000317 2371 0.479 1.0000
## TCS08 - TCS40 -1.16e-03 0.000415 2371 -2.795 0.3809
## TCS08 - TCS41 -8.91e-04 0.000545 2371 -1.635 0.9895
## TCS08 - TCS42 -5.60e-04 0.000488 2371 -1.149 0.9999
## TCS08 - TCS43 2.03e-04 0.000325 2371 0.623 1.0000
## TCS08 - TCS44 2.55e-04 0.000318 2371 0.803 1.0000
## TCS08 - TCS45 7.46e-04 0.000319 2371 2.337 0.7315
## TCS08 - TCS46 -2.61e-04 0.000317 2371 -0.823 1.0000
## TCS08 - TCS47 2.62e-04 0.000317 2371 0.827 1.0000
## TCS08 - TCS48 2.39e-05 0.000318 2371 0.075 1.0000
## TCS08 - TCS49 8.78e-06 0.000317 2371 0.028 1.0000
## TCS08 - TSC45 5.03e-05 0.001079 2371 0.047 1.0000
## TCS10 - TCS11 -2.78e-04 0.000276 2371 -1.008 1.0000
## TCS10 - TCS12 2.81e-04 0.000276 2371 1.018 1.0000
## TCS10 - TCS20 1.66e-04 0.000276 2371 0.600 1.0000
## TCS10 - TCS40 -1.15e-03 0.000384 2371 -2.979 0.2599
## TCS10 - TCS41 -8.78e-04 0.000522 2371 -1.680 0.9857
## TCS10 - TCS42 -5.46e-04 0.000462 2371 -1.182 0.9999
## TCS10 - TCS43 2.16e-04 0.000285 2371 0.758 1.0000
## TCS10 - TCS44 2.69e-04 0.000277 2371 0.971 1.0000
## TCS10 - TCS45 7.60e-04 0.000279 2371 2.727 0.4310
## TCS10 - TCS46 -2.47e-04 0.000276 2371 -0.894 1.0000
## TCS10 - TCS47 2.76e-04 0.000276 2371 1.000 1.0000
## TCS10 - TCS48 3.78e-05 0.000277 2371 0.136 1.0000
## TCS10 - TCS49 2.26e-05 0.000276 2371 0.082 1.0000
## TCS10 - TSC45 6.42e-05 0.001068 2371 0.060 1.0000
## TCS11 - TCS12 5.60e-04 0.000276 2371 2.026 0.9062
## TCS11 - TCS20 4.44e-04 0.000276 2371 1.608 0.9914
## TCS11 - TCS40 -8.67e-04 0.000384 2371 -2.255 0.7870
## TCS11 - TCS41 -5.99e-04 0.000522 2371 -1.147 0.9999
## TCS11 - TCS42 -2.68e-04 0.000462 2371 -0.580 1.0000
## TCS11 - TCS43 4.95e-04 0.000285 2371 1.733 0.9797
## TCS11 - TCS44 5.48e-04 0.000277 2371 1.975 0.9251
## TCS11 - TCS45 1.04e-03 0.000279 2371 3.726 0.0298
## TCS11 - TCS46 3.14e-05 0.000276 2371 0.114 1.0000
## TCS11 - TCS47 5.55e-04 0.000276 2371 2.008 0.9133
## TCS11 - TCS48 3.16e-04 0.000277 2371 1.140 0.9999
## TCS11 - TCS49 3.01e-04 0.000276 2371 1.090 1.0000
## TCS11 - TSC45 3.43e-04 0.001068 2371 0.321 1.0000
## TCS12 - TCS20 -1.15e-04 0.000276 2371 -0.418 1.0000
## TCS12 - TCS40 -1.43e-03 0.000384 2371 -3.711 0.0314
## TCS12 - TCS41 -1.16e-03 0.000522 2371 -2.218 0.8099
## TCS12 - TCS42 -8.27e-04 0.000462 2371 -1.791 0.9712
## TCS12 - TCS43 -6.48e-05 0.000285 2371 -0.227 1.0000
## TCS12 - TCS44 -1.19e-05 0.000277 2371 -0.043 1.0000
## TCS12 - TCS45 4.79e-04 0.000279 2371 1.718 0.9816
## TCS12 - TCS46 -5.28e-04 0.000276 2371 -1.913 0.9443
## TCS12 - TCS47 -5.08e-06 0.000276 2371 -0.018 1.0000
## TCS12 - TCS48 -2.43e-04 0.000277 2371 -0.878 1.0000
## TCS12 - TCS49 -2.59e-04 0.000276 2371 -0.936 1.0000
## TCS12 - TSC45 -2.17e-04 0.001068 2371 -0.203 1.0000
## TCS20 - TCS40 -1.31e-03 0.000384 2371 -3.411 0.0832
## TCS20 - TCS41 -1.04e-03 0.000522 2371 -1.997 0.9173
## TCS20 - TCS42 -7.12e-04 0.000462 2371 -1.541 0.9949
## TCS20 - TCS43 5.06e-05 0.000285 2371 0.177 1.0000
## TCS20 - TCS44 1.03e-04 0.000277 2371 0.373 1.0000
## TCS20 - TCS45 5.94e-04 0.000279 2371 2.132 0.8582
## TCS20 - TCS46 -4.13e-04 0.000276 2371 -1.495 0.9966
## TCS20 - TCS47 1.10e-04 0.000276 2371 0.400 1.0000
## TCS20 - TCS48 -1.28e-04 0.000277 2371 -0.462 1.0000
## TCS20 - TCS49 -1.43e-04 0.000276 2371 -0.518 1.0000
## TCS20 - TSC45 -1.02e-04 0.001068 2371 -0.095 1.0000
## TCS40 - TCS41 2.68e-04 0.000587 2371 0.456 1.0000
## TCS40 - TCS42 5.99e-04 0.000534 2371 1.122 0.9999
## TCS40 - TCS43 1.36e-03 0.000391 2371 3.481 0.0670
## TCS40 - TCS44 1.41e-03 0.000385 2371 3.672 0.0359
## TCS40 - TCS45 1.90e-03 0.000386 2371 4.934 0.0002
## TCS40 - TCS46 8.98e-04 0.000384 2371 2.337 0.7318
## TCS40 - TCS47 1.42e-03 0.000384 2371 3.698 0.0329
## TCS40 - TCS48 1.18e-03 0.000385 2371 3.071 0.2097
## TCS40 - TCS49 1.17e-03 0.000384 2371 3.038 0.2269
## TCS40 - TSC45 1.21e-03 0.001101 2371 1.099 1.0000
## TCS41 - TCS42 3.31e-04 0.000640 2371 0.517 1.0000
## TCS41 - TCS43 1.09e-03 0.000527 2371 2.074 0.8860
## TCS41 - TCS44 1.15e-03 0.000523 2371 2.193 0.8249
## TCS41 - TCS45 1.64e-03 0.000524 2371 3.127 0.1829
## TCS41 - TCS46 6.31e-04 0.000522 2371 1.207 0.9998
## TCS41 - TCS47 1.15e-03 0.000522 2371 2.208 0.8157
## TCS41 - TCS48 9.15e-04 0.000523 2371 1.750 0.9774
## TCS41 - TCS49 9.00e-04 0.000522 2371 1.723 0.9810
## TCS41 - TSC45 9.42e-04 0.001156 2371 0.815 1.0000
## TCS42 - TCS43 7.63e-04 0.000468 2371 1.631 0.9899
## TCS42 - TCS44 8.16e-04 0.000463 2371 1.762 0.9757
## TCS42 - TCS45 1.31e-03 0.000464 2371 2.818 0.3646
## TCS42 - TCS46 2.99e-04 0.000462 2371 0.648 1.0000
## TCS42 - TCS47 8.22e-04 0.000462 2371 1.780 0.9730
## TCS42 - TCS48 5.84e-04 0.000463 2371 1.262 0.9997
## TCS42 - TCS49 5.69e-04 0.000462 2371 1.231 0.9998
## TCS42 - TSC45 6.10e-04 0.001130 2371 0.540 1.0000
## TCS43 - TCS44 5.29e-05 0.000287 2371 0.185 1.0000
## TCS43 - TCS45 5.44e-04 0.000288 2371 1.888 0.9509
## TCS43 - TCS46 -4.63e-04 0.000285 2371 -1.623 0.9904
## TCS43 - TCS47 5.98e-05 0.000285 2371 0.209 1.0000
## TCS43 - TCS48 -1.79e-04 0.000287 2371 -0.623 1.0000
## TCS43 - TCS49 -1.94e-04 0.000285 2371 -0.679 1.0000
## TCS43 - TSC45 -1.52e-04 0.001070 2371 -0.142 1.0000
## TCS44 - TCS45 4.91e-04 0.000280 2371 1.754 0.9769
## TCS44 - TCS46 -5.16e-04 0.000277 2371 -1.862 0.9571
## TCS44 - TCS47 6.86e-06 0.000277 2371 0.025 1.0000
## TCS44 - TCS48 -2.32e-04 0.000278 2371 -0.832 1.0000
## TCS44 - TCS49 -2.47e-04 0.000277 2371 -0.890 1.0000
## TCS44 - TSC45 -2.05e-04 0.001068 2371 -0.192 1.0000
## TCS45 - TCS46 -1.01e-03 0.000279 2371 -3.613 0.0438
## TCS45 - TCS47 -4.84e-04 0.000279 2371 -1.736 0.9793
## TCS45 - TCS48 -7.22e-04 0.000280 2371 -2.581 0.5447
## TCS45 - TCS49 -7.37e-04 0.000279 2371 -2.646 0.4936
## TCS45 - TSC45 -6.96e-04 0.001068 2371 -0.651 1.0000
## TCS46 - TCS47 5.23e-04 0.000276 2371 1.894 0.9492
## TCS46 - TCS48 2.85e-04 0.000277 2371 1.027 1.0000
## TCS46 - TCS49 2.70e-04 0.000276 2371 0.976 1.0000
## TCS46 - TSC45 3.11e-04 0.001068 2371 0.291 1.0000
## TCS47 - TCS48 -2.38e-04 0.000277 2371 -0.860 1.0000
## TCS47 - TCS49 -2.54e-04 0.000276 2371 -0.918 1.0000
## TCS47 - TSC45 -2.12e-04 0.001068 2371 -0.199 1.0000
## TCS48 - TCS49 -1.51e-05 0.000277 2371 -0.055 1.0000
## TCS48 - TSC45 2.64e-05 0.001068 2371 0.025 1.0000
## TCS49 - TSC45 4.15e-05 0.001068 2371 0.039 1.0000
##
## P value adjustment: tukey method for comparing a family of 21 estimates
emmip(contcat, gen ~ dÃas, cov.reduce = range)
library(corrplot)
## corrplot 0.90 loaded
data2<-read.table("ratecor.csv", header=T, sep=",")
rownames(data2) <- data2[,1]
data2[,1] <- NULL
M=cor(data2)
corrplot(M, method="number")
library(psych)
data2<-read.table("ratecor.csv", header=T, sep=",")
pairs.panels(data2,
gap = 0,
pch=21)
pc <- prcomp(data2[,-1],
center = TRUE,
scale. = TRUE)
attributes(pc)
## $names
## [1] "sdev" "rotation" "center" "scale" "x"
##
## $class
## [1] "prcomp"
print(pc)
## Standard deviations (1, .., p=7):
## [1] 1.4816702 1.3447597 1.1415458 0.9853984 0.7175438 0.4063280 0.2053445
##
## Rotation (n x k) = (7 x 7):
## PC1 PC2 PC3 PC4 PC5 PC6
## carbon -0.24223059 -0.35146814 -0.31564321 0.57126625 -0.62483922 -0.007013014
## alfa -0.08063585 -0.06399609 -0.80190148 0.02836859 0.48909862 0.311360234
## A 0.41855885 -0.49242283 -0.20728590 -0.30933545 -0.02463908 -0.270516150
## E 0.60801197 -0.18026785 0.13824628 0.03738160 -0.18908171 0.711051890
## WUE -0.19230802 -0.64819320 0.08906996 -0.36363070 0.02735377 -0.154802451
## gsW 0.56904171 0.05169061 -0.08584105 0.45721699 0.19570337 -0.528427193
## ci.ca -0.17076677 -0.41781761 0.42419543 0.48424130 0.54309435 0.145026708
## PC7
## carbon -0.03668097
## alfa 0.09689658
## A -0.60818280
## E 0.18965453
## WUE 0.61480926
## gsW 0.37358860
## ci.ca -0.25661554
summary(pc)
## Importance of components:
## PC1 PC2 PC3 PC4 PC5 PC6 PC7
## Standard deviation 1.4817 1.3448 1.1415 0.9854 0.71754 0.40633 0.20534
## Proportion of Variance 0.3136 0.2583 0.1862 0.1387 0.07355 0.02359 0.00602
## Cumulative Proportion 0.3136 0.5720 0.7581 0.8968 0.97039 0.99398 1.00000
library(devtools)
## Loading required package: usethis
##
## Attaching package: 'devtools'
## The following object is masked from 'package:emmeans':
##
## test
library(ggbiplot)
## Loading required package: ggplot2
##
## Attaching package: 'ggplot2'
## The following objects are masked from 'package:psych':
##
## %+%, alpha
## Loading required package: plyr
## Loading required package: scales
##
## Attaching package: 'scales'
## The following objects are masked from 'package:psych':
##
## alpha, rescale
## Loading required package: grid
g <- ggbiplot(pc,
choices = c(1,2),
obs.scale = 1,
var.scale = 1,
labels = data2$gen,
ellipse = TRUE,
circle = TRUE,
ellipse.prob = 0.68)
g <- g + scale_color_discrete(name = '')
g <- g + theme(legend.direction = 'horizontal',
legend.position = 'top')
print(g)
g2 <- ggbiplot(pc,
choices = c(1,4),
obs.scale = 1,
var.scale = 1,
labels = data2$gen,
ellipse = TRUE,
circle = TRUE,
ellipse.prob = 0.68)
g2 <- g2 + scale_color_discrete(name = '')
g2 <- g2 + theme(legend.direction = 'horizontal',
legend.position = 'top')
print(g2)
