################################################## # FinalPhiAltModelSelection.R # Sean Connolly # Final models for the cubic transofrmation of log(synchrony) ######################################################### ##################################### # Start with the full synthetic model (without CDI) DivDDclSEM.AltPhi <- psem( lm(cvpop ~ latitude+CrossShelfPos+mHC+sdHC+estRichness+estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ latitude+CrossShelfPos+mHC+sdHC+estRichness+estUnevenness+DD.cublog, ecol.table.new), lm(estUnevenness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(DD.cublog ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi) #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # altphi ~ cvpop + ... coef 30 -2.0964 0.0446 * # estUnevenness ~ estRichness + ... coef 33 -1.1009 0.2789 # DD.cublog ~ estRichness + ... coef 33 -0.1328 0.8952 # DD.cublog ~ estUnevenness + ... coef 33 1.3542 0.1849 # #-- #Global goodness-of-fit: # #Chi-Squared = 8.865 with P-value = 0.065 and on 4 degrees of freedom #Fisher's C = 12.372 with P-value = 0.135 and on 8 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop latitude 0.0033 0.0130 31 0.2526 0.8023 0.0605 # cvpop CrossShelfPos -0.0093 0.0982 31 -0.0947 0.9252 -0.0181 # cvpop mHC -0.0322 0.0675 31 -0.4768 0.6368 -0.0816 # cvpop sdHC 0.0574 0.0568 31 1.0104 0.3201 0.2017 # cvpop estRichness -0.0513 0.1759 31 -0.2914 0.7727 -0.0642 # cvpop estUnevenness 0.5050 0.3937 31 1.2826 0.2092 0.3133 # cvpop DD.cublog 0.0055 0.0023 31 2.3933 0.0229 0.4284 * # altphi latitude 0.5729 0.2907 31 1.9708 0.0577 0.4382 # altphi CrossShelfPos -6.8960 2.2017 31 -3.1321 0.0038 -0.5556 ** # altphi mHC -0.7941 1.5140 31 -0.5245 0.6036 -0.0832 # altphi sdHC -0.6546 1.2735 31 -0.5140 0.6109 -0.0952 # altphi estRichness -0.5873 3.9429 31 -0.1489 0.8826 -0.0305 # altphi estUnevenness 17.9212 8.8247 31 2.0308 0.0509 0.4603 # altphi DD.cublog -0.0500 0.0515 31 -0.9703 0.3394 -0.1612 # estUnevenness latitude -0.0245 0.0036 34 -6.7870 0.0000 -0.7288 *** # estUnevenness CrossShelfPos -0.0107 0.0327 34 -0.3271 0.7456 -0.0335 # estUnevenness mHC 0.0196 0.0270 34 0.7237 0.4742 0.0798 # estUnevenness sdHC 0.0281 0.0196 34 1.4336 0.1608 0.1590 # estRichness latitude 0.0255 0.0079 34 3.2437 0.0026 0.3756 ** # estRichness CrossShelfPos 0.3811 0.0712 34 5.3568 0.0000 0.5917 *** # estRichness mHC -0.1549 0.0589 34 -2.6289 0.0128 -0.3128 * # estRichness sdHC 0.0906 0.0427 34 2.1243 0.0410 0.2541 * # DD.cublog latitude 0.6698 0.6071 34 1.1034 0.2776 0.1588 # DD.cublog CrossShelfPos 4.6898 5.4975 34 0.8531 0.3996 0.1171 # DD.cublog mHC -7.2873 4.5514 34 -1.6011 0.1186 -0.2368 # DD.cublog sdHC 14.6458 3.2967 34 4.4426 0.0001 0.6606 *** AIC(DivDDclSEM.AltPhi,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 84.372 1416.372 36 39 ############## # Now start successively dropping terms: ################# DivDDclSEM.AltPhi.trim1 <- psem( lm(cvpop ~ mHC+sdHC+estRichness+estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ latitude+CrossShelfPos+mHC+sdHC+estUnevenness+DD.cublog, ecol.table.new), lm(estUnevenness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(DD.cublog ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim1) #--- #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # cvpop ~ latitude + ... coef 32 0.2905 0.7733 # cvpop ~ CrossShelfPos + ... coef 32 -0.1666 0.8687 # altphi ~ cvpop + ... coef 30 -2.0964 0.0446 * # estUnevenness ~ estRichness + ... coef 33 -1.1009 0.2789 # DD.cublog ~ estRichness + ... coef 33 -0.1328 0.8952 # altphi ~ estRichness + ... coef 31 -0.1489 0.8826 # DD.cublog ~ estUnevenness + ... coef 33 1.3542 0.1849 # #-- #Global goodness-of-fit: # #Chi-Squared = 9.007 with P-value = 0.252 and on 7 degrees of freedom #Fisher's C = 13.418 with P-value = 0.494 and on 14 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop mHC -0.0329 0.0640 33 -0.5138 0.6108 -0.0833 # cvpop sdHC 0.0560 0.0542 33 1.0339 0.3087 0.1968 # cvpop estRichness -0.0607 0.1263 33 -0.4810 0.6337 -0.0761 # cvpop estUnevenness 0.4188 0.2539 33 1.6494 0.1085 0.2599 # cvpop DD.cublog 0.0056 0.0021 33 2.6418 0.0125 0.4401 * # altphi latitude 0.5641 0.2802 32 2.0130 0.0526 0.4314 # altphi CrossShelfPos -7.1164 1.6050 32 -4.4340 0.0001 -0.5733 *** # altphi mHC -0.7092 1.3809 32 -0.5136 0.6111 -0.0743 # altphi sdHC -0.7124 1.1943 32 -0.5965 0.5550 -0.1036 # altphi estUnevenness 18.1682 8.5339 32 2.1289 0.0411 0.4666 * # altphi DD.cublog -0.0501 0.0507 32 -0.9889 0.3302 -0.1617 # estUnevenness latitude -0.0245 0.0036 34 -6.7870 0.0000 -0.7288 *** # estUnevenness CrossShelfPos -0.0107 0.0327 34 -0.3271 0.7456 -0.0335 # estUnevenness mHC 0.0196 0.0270 34 0.7237 0.4742 0.0798 # estUnevenness sdHC 0.0281 0.0196 34 1.4336 0.1608 0.1590 # estRichness latitude 0.0255 0.0079 34 3.2437 0.0026 0.3756 ** # estRichness CrossShelfPos 0.3811 0.0712 34 5.3568 0.0000 0.5917 *** # estRichness mHC -0.1549 0.0589 34 -2.6289 0.0128 -0.3128 * # estRichness sdHC 0.0906 0.0427 34 2.1243 0.0410 0.2541 * # DD.cublog latitude 0.6698 0.6071 34 1.1034 0.2776 0.1588 # DD.cublog CrossShelfPos 4.6898 5.4975 34 0.8531 0.3996 0.1171 # DD.cublog mHC -7.2873 4.5514 34 -1.6011 0.1186 -0.2368 # DD.cublog sdHC 14.6458 3.2967 34 4.4426 0.0001 0.6606 *** AIC(DivDDclSEM.AltPhi.trim1,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 79.418 528.218 33 39 DivDDclSEM.AltPhi.trim2 <- psem( lm(cvpop ~ sdHC+estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ latitude+CrossShelfPos+sdHC+estUnevenness+DD.cublog, ecol.table.new), lm(estUnevenness ~ latitude+mHC+sdHC,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(DD.cublog ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim2) #--- #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # cvpop ~ latitude + ... coef 34 0.2665 0.7915 # cvpop ~ CrossShelfPos + ... coef 34 -0.4239 0.6743 # estUnevenness ~ CrossShelfPos + ... coef 34 -0.3271 0.7456 # cvpop ~ mHC + ... coef 34 -0.3837 0.7036 # altphi ~ mHC + ... coef 32 -0.5136 0.6111 # altphi ~ cvpop + ... coef 32 -2.0909 0.0446 * # estRichness ~ cvpop + ... coef 31 -0.2914 0.7727 # DD.cublog ~ estUnevenness + ... coef 33 1.3542 0.1849 # estRichness ~ estUnevenness + ... coef 33 -1.1009 0.2789 # estRichness ~ DD.cublog + ... coef 33 -0.1328 0.8952 # estRichness ~ altphi + ... coef 31 -0.1489 0.8826 # #-- #Global goodness-of-fit: # #Chi-Squared = 9.89 with P-value = 0.54 and on 11 degrees of freedom #Fisher's C = 16.669 with P-value = 0.781 and on 22 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop sdHC 0.0410 0.0468 35 0.8755 0.3873 0.1439 # cvpop estUnevenness 0.4374 0.2235 35 1.9570 0.0584 0.2714 # cvpop DD.cublog 0.0058 0.0020 35 2.9344 0.0059 0.4484 ** # altphi latitude 0.5623 0.2771 33 2.0295 0.0505 0.4301 # altphi CrossShelfPos -7.0682 1.5842 33 -4.4616 0.0001 -0.5695 *** # altphi sdHC -0.9490 1.0895 33 -0.8710 0.3900 -0.1381 # altphi estUnevenness 17.3693 8.2968 33 2.0935 0.0441 0.4461 * # altphi DD.cublog -0.0424 0.0479 33 -0.8862 0.3819 -0.1369 # estUnevenness latitude -0.0243 0.0035 35 -6.8857 0.0000 -0.7245 *** # estUnevenness mHC 0.0206 0.0265 35 0.7787 0.4414 0.0842 # estUnevenness sdHC 0.0266 0.0188 35 1.4139 0.1662 0.1506 # estRichness latitude 0.0255 0.0079 34 3.2437 0.0026 0.3756 ** # estRichness CrossShelfPos 0.3811 0.0712 34 5.3568 0.0000 0.5917 *** # estRichness mHC -0.1549 0.0589 34 -2.6289 0.0128 -0.3128 * # estRichness sdHC 0.0906 0.0427 34 2.1243 0.0410 0.2541 * # DD.cublog latitude 0.6698 0.6071 34 1.1034 0.2776 0.1588 # DD.cublog CrossShelfPos 4.6898 5.4975 34 0.8531 0.3996 0.1171 # DD.cublog mHC -7.2873 4.5514 34 -1.6011 0.1186 -0.2368 # DD.cublog sdHC 14.6458 3.2967 34 4.4426 0.0001 0.6606 *** AIC(DivDDclSEM.AltPhi.trim2,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 74.669 268.002 29 39 DivDDclSEM.AltPhi.trim3 <- psem( lm(cvpop ~ estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ latitude+CrossShelfPos+estUnevenness, ecol.table.new), lm(estUnevenness ~ latitude+sdHC,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(DD.cublog ~ latitude+mHC+sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim3) #--- #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # cvpop ~ latitude + ... coef 35 0.1558 0.8771 # estUnevenness ~ CrossShelfPos + ... coef 35 -0.4203 0.6769 # DD.cublog ~ CrossShelfPos + ... coef 34 0.8531 0.3996 # cvpop ~ CrossShelfPos + ... coef 35 -0.3338 0.7405 # altphi ~ sdHC + ... coef 34 -1.5724 0.1251 # cvpop ~ sdHC + ... coef 35 0.8755 0.3873 # estUnevenness ~ mHC + ... coef 35 0.7787 0.4414 # altphi ~ mHC + ... coef 34 -0.6415 0.5255 # cvpop ~ mHC + ... coef 35 -0.0235 0.9814 # DD.cublog ~ estUnevenness + ... coef 34 1.3060 0.2003 # estRichness ~ estUnevenness + ... coef 33 -1.1009 0.2789 # altphi ~ DD.cublog + ... coef 32 -0.9889 0.3302 # estRichness ~ DD.cublog + ... coef 33 -0.1328 0.8952 # estRichness ~ altphi + ... coef 32 -0.1695 0.8665 # cvpop ~ altphi + ... coef 33 -2.2306 0.0326 * # cvpop ~ estRichness + ... coef 31 -0.2914 0.7727 # #-- #Global goodness-of-fit: # #Chi-Squared = 15.886 with P-value = 0.461 and on 16 degrees of freedom #Fisher's C = 28.348 with P-value = 0.652 and on 32 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop estUnevenness 0.5111 0.2064 36 2.4755 0.0181 0.3171 * # cvpop DD.cublog 0.0067 0.0016 36 4.0666 0.0002 0.5208 *** # altphi latitude 0.4639 0.2683 35 1.7291 0.0926 0.3548 # altphi CrossShelfPos -7.8572 1.5463 35 -5.0811 0.0000 -0.6330 *** # altphi estUnevenness 12.3412 7.9475 35 1.5528 0.1295 0.3170 # estUnevenness latitude -0.0250 0.0034 36 -7.3770 0.0000 -0.7457 *** # estUnevenness sdHC 0.0310 0.0178 36 1.7359 0.0911 0.1755 # estRichness latitude 0.0255 0.0079 34 3.2437 0.0026 0.3756 ** # estRichness CrossShelfPos 0.3811 0.0712 34 5.3568 0.0000 0.5917 *** # estRichness mHC -0.1549 0.0589 34 -2.6289 0.0128 -0.3128 * # estRichness sdHC 0.0906 0.0427 34 2.1243 0.0410 0.2541 * # DD.cublog latitude 0.6072 0.6003 35 1.0115 0.3187 0.1440 # DD.cublog mHC -7.7552 4.5006 35 -1.7231 0.0937 -0.2520 # DD.cublog sdHC 15.2948 3.1952 35 4.7868 0.0000 0.6899 *** AIC(DivDDclSEM.AltPhi.trim3,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 76.348 162.062 24 39 DivDDclSEM.AltPhi.trim4 <- psem( lm(cvpop ~ estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ latitude+CrossShelfPos, ecol.table.new), lm(estUnevenness ~ latitude+sdHC,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(DD.cublog ~ mHC+sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim4) #--- #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # DD.cublog ~ latitude + ... coef 35 1.0115 0.3187 # cvpop ~ latitude + ... coef 35 0.1558 0.8771 # estUnevenness ~ CrossShelfPos + ... coef 35 -0.4203 0.6769 # DD.cublog ~ CrossShelfPos + ... coef 35 0.7227 0.4747 # cvpop ~ CrossShelfPos + ... coef 35 -0.3338 0.7405 # altphi ~ sdHC + ... coef 35 -1.0126 0.3182 # cvpop ~ sdHC + ... coef 35 0.8755 0.3873 # estUnevenness ~ mHC + ... coef 35 0.7787 0.4414 # altphi ~ mHC + ... coef 35 -0.3089 0.7592 # cvpop ~ mHC + ... coef 35 -0.0235 0.9814 # DD.cublog ~ estUnevenness + ... coef 34 1.3060 0.2003 # altphi ~ estUnevenness + ... coef 34 1.9686 0.0572 # estRichness ~ estUnevenness + ... coef 33 -1.1009 0.2789 # altphi ~ DD.cublog + ... coef 33 -0.4886 0.6283 # estRichness ~ DD.cublog + ... coef 33 -0.1328 0.8952 # estRichness ~ altphi + ... coef 33 -0.5103 0.6132 # cvpop ~ altphi + ... coef 33 -2.2306 0.0326 * # cvpop ~ estRichness + ... coef 31 -0.2914 0.7727 # #-- #Global goodness-of-fit: # #Chi-Squared = 19.608 with P-value = 0.355 and on 18 degrees of freedom #Fisher's C = 32.815 with P-value = 0.621 and on 36 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop estUnevenness 0.5111 0.2064 36 2.4755 0.0181 0.3171 * # cvpop DD.cublog 0.0067 0.0016 36 4.0666 0.0002 0.5208 *** # altphi latitude 0.1329 0.1660 36 0.8001 0.4289 0.1016 # altphi CrossShelfPos -7.8807 1.5763 36 -4.9995 0.0000 -0.6349 *** # estUnevenness latitude -0.0250 0.0034 36 -7.3770 0.0000 -0.7457 *** # estUnevenness sdHC 0.0310 0.0178 36 1.7359 0.0911 0.1755 # estRichness latitude 0.0255 0.0079 34 3.2437 0.0026 0.3756 ** # estRichness CrossShelfPos 0.3811 0.0712 34 5.3568 0.0000 0.5917 *** # estRichness mHC -0.1549 0.0589 34 -2.6289 0.0128 -0.3128 * # estRichness sdHC 0.0906 0.0427 34 2.1243 0.0410 0.2541 * # DD.cublog mHC -8.9323 4.3489 36 -2.0539 0.0473 -0.2903 * # DD.cublog sdHC 14.6581 3.1336 36 4.6777 0.0000 0.6611 *** AIC(DivDDclSEM.AltPhi.trim4,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 76.815 140.065 22 39 DivDDclSEM.AltPhi.trim5 <- psem( lm(cvpop ~ estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ CrossShelfPos, ecol.table.new), lm(estUnevenness ~ latitude,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(DD.cublog ~ mHC+sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim5) #--- #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # estUnevenness ~ CrossShelfPos + ... coef 36 -0.0587 0.9535 # DD.cublog ~ CrossShelfPos + ... coef 35 0.7227 0.4747 # cvpop ~ CrossShelfPos + ... coef 35 -0.3338 0.7405 # DD.cublog ~ latitude + ... coef 35 1.0115 0.3187 # altphi ~ latitude + ... coef 36 0.8001 0.4289 # cvpop ~ latitude + ... coef 35 0.1558 0.8771 # estUnevenness ~ mHC + ... coef 36 1.2429 0.2219 # altphi ~ mHC + ... coef 36 -0.5658 0.5751 # cvpop ~ mHC + ... coef 35 -0.0235 0.9814 # estUnevenness ~ sdHC + ... coef 36 1.7359 0.0911 # altphi ~ sdHC + ... coef 36 -1.2045 0.2363 # cvpop ~ sdHC + ... coef 35 0.8755 0.3873 # DD.cublog ~ estUnevenness + ... coef 34 1.3060 0.2003 # altphi ~ estUnevenness + ... coef 35 1.5528 0.1295 # estRichness ~ estUnevenness + ... coef 33 -1.1009 0.2789 # altphi ~ DD.cublog + ... coef 34 -0.3973 0.6936 # estRichness ~ DD.cublog + ... coef 33 -0.1328 0.8952 # estRichness ~ altphi + ... coef 33 -0.5103 0.6132 # cvpop ~ altphi + ... coef 34 -2.0649 0.0466 * # cvpop ~ estRichness + ... coef 31 -0.2914 0.7727 # #-- #Global goodness-of-fit: # #Chi-Squared = 23.43 with P-value = 0.268 and on 20 degrees of freedom #Fisher's C = 38.593 with P-value = 0.534 and on 40 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop estUnevenness 0.5111 0.2064 36 2.4755 0.0181 0.3171 * # cvpop DD.cublog 0.0067 0.0016 36 4.0666 0.0002 0.5208 *** # altphi CrossShelfPos -8.0789 1.5491 37 -5.2151 0.0000 -0.6509 *** # estUnevenness latitude -0.0268 0.0033 37 -8.0529 0.0000 -0.7979 *** # estRichness latitude 0.0255 0.0079 34 3.2437 0.0026 0.3756 ** # estRichness CrossShelfPos 0.3811 0.0712 34 5.3568 0.0000 0.5917 *** # estRichness mHC -0.1549 0.0589 34 -2.6289 0.0128 -0.3128 * # estRichness sdHC 0.0906 0.0427 34 2.1243 0.0410 0.2541 * # DD.cublog mHC -8.9323 4.3489 36 -2.0539 0.0473 -0.2903 * # DD.cublog sdHC 14.6581 3.1336 36 4.6777 0.0000 0.6611 *** AIC(DivDDclSEM.AltPhi.trim5,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 78.593 125.26 20 39 DivDDclSEM.AltPhi.trim6 <- psem( lm(cvpop ~ estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ CrossShelfPos, ecol.table.new), lm(estUnevenness ~ latitude,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos+mHC+sdHC,ecol.table.new), lm(DD.cublog ~ sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim6) #--- #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # estUnevenness ~ CrossShelfPos + ... coef 36 -0.0587 0.9535 # DD.cublog ~ CrossShelfPos + ... coef 36 0.8733 0.3883 # cvpop ~ CrossShelfPos + ... coef 35 -0.3338 0.7405 # DD.cublog ~ latitude + ... coef 36 1.4688 0.1506 # altphi ~ latitude + ... coef 36 0.8001 0.4289 # cvpop ~ latitude + ... coef 35 0.1558 0.8771 # estUnevenness ~ mHC + ... coef 36 1.2429 0.2219 # DD.cublog ~ mHC + ... coef 36 -2.0539 0.0473 * # altphi ~ mHC + ... coef 36 -0.5658 0.5751 # cvpop ~ mHC + ... coef 35 -0.0235 0.9814 # estUnevenness ~ sdHC + ... coef 36 1.7359 0.0911 # altphi ~ sdHC + ... coef 36 -1.2045 0.2363 # cvpop ~ sdHC + ... coef 35 0.8755 0.3873 # DD.cublog ~ estUnevenness + ... coef 35 1.0302 0.3100 # altphi ~ estUnevenness + ... coef 35 1.5528 0.1295 # estRichness ~ estUnevenness + ... coef 33 -1.1009 0.2789 # altphi ~ DD.cublog + ... coef 35 -0.3428 0.7338 # estRichness ~ DD.cublog + ... coef 33 -0.1328 0.8952 # estRichness ~ altphi + ... coef 33 -0.5103 0.6132 # cvpop ~ altphi + ... coef 34 -2.0649 0.0466 * # cvpop ~ estRichness + ... coef 31 -0.2914 0.7727 # #-- #Global goodness-of-fit: # #Chi-Squared = 27.752 with P-value = 0.147 and on 21 degrees of freedom #Fisher's C = 45.611 with P-value = 0.324 and on 42 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop estUnevenness 0.5111 0.2064 36 2.4755 0.0181 0.3171 * # cvpop DD.cublog 0.0067 0.0016 36 4.0666 0.0002 0.5208 *** # altphi CrossShelfPos -8.0789 1.5491 37 -5.2151 0.0000 -0.6509 *** # estUnevenness latitude -0.0268 0.0033 37 -8.0529 0.0000 -0.7979 *** # estRichness latitude 0.0255 0.0079 34 3.2437 0.0026 0.3756 ** # estRichness CrossShelfPos 0.3811 0.0712 34 5.3568 0.0000 0.5917 *** # estRichness mHC -0.1549 0.0589 34 -2.6289 0.0128 -0.3128 * # estRichness sdHC 0.0906 0.0427 34 2.1243 0.0410 0.2541 * # DD.cublog sdHC 12.2766 3.0352 37 4.0448 0.0003 0.5537 *** AIC(DivDDclSEM.AltPhi.trim6,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 83.611 123.611 19 39 DivDDclSEM.AltPhi.trim7 <- psem( lm(cvpop ~ estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ CrossShelfPos, ecol.table.new), lm(estUnevenness ~ latitude,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos+mHC,ecol.table.new), lm(DD.cublog ~ sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim7) #--- #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # estUnevenness ~ CrossShelfPos + ... coef 36 -0.0587 0.9535 # DD.cublog ~ CrossShelfPos + ... coef 36 0.8733 0.3883 # cvpop ~ CrossShelfPos + ... coef 35 -0.3338 0.7405 # DD.cublog ~ latitude + ... coef 36 1.4688 0.1506 # altphi ~ latitude + ... coef 36 0.8001 0.4289 # cvpop ~ latitude + ... coef 35 0.1558 0.8771 # estUnevenness ~ mHC + ... coef 36 1.2429 0.2219 # DD.cublog ~ mHC + ... coef 36 -2.0539 0.0473 * # altphi ~ mHC + ... coef 36 -0.5658 0.5751 # cvpop ~ mHC + ... coef 35 -0.0235 0.9814 # estUnevenness ~ sdHC + ... coef 36 1.7359 0.0911 # altphi ~ sdHC + ... coef 36 -1.2045 0.2363 # estRichness ~ sdHC + ... coef 34 2.1243 0.0410 * # cvpop ~ sdHC + ... coef 35 0.8755 0.3873 # DD.cublog ~ estUnevenness + ... coef 35 1.0302 0.3100 # altphi ~ estUnevenness + ... coef 35 1.5528 0.1295 # estRichness ~ estUnevenness + ... coef 34 -0.5269 0.6017 # altphi ~ DD.cublog + ... coef 35 -0.3428 0.7338 # estRichness ~ DD.cublog + ... coef 33 -0.1328 0.8952 # estRichness ~ altphi + ... coef 34 -0.8068 0.4254 # cvpop ~ altphi + ... coef 34 -2.0649 0.0466 * # cvpop ~ estRichness + ... coef 32 0.0173 0.9863 # #-- #Global goodness-of-fit: # #Chi-Squared = 32.613 with P-value = 0.068 and on 22 degrees of freedom #Fisher's C = 50.705 with P-value = 0.226 and on 44 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop estUnevenness 0.5111 0.2064 36 2.4755 0.0181 0.3171 * # cvpop DD.cublog 0.0067 0.0016 36 4.0666 0.0002 0.5208 *** # altphi CrossShelfPos -8.0789 1.5491 37 -5.2151 0.0000 -0.6509 *** # estUnevenness latitude -0.0268 0.0033 37 -8.0529 0.0000 -0.7979 *** # estRichness latitude 0.0228 0.0081 35 2.8005 0.0083 0.3357 ** # estRichness CrossShelfPos 0.4160 0.0726 35 5.7287 0.0000 0.6459 *** # estRichness mHC -0.1152 0.0586 35 -1.9654 0.0573 -0.2327 # DD.cublog sdHC 12.2766 3.0352 37 4.0448 0.0003 0.5537 *** AIC(DivDDclSEM.AltPhi.trim7,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 86.705 120.905 18 39 DivDDclSEM.AltPhi.trim8 <- psem( lm(cvpop ~ estUnevenness+DD.cublog, ecol.table.new), lm(altphi ~ CrossShelfPos, ecol.table.new), lm(estUnevenness ~ latitude,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos,ecol.table.new), lm(DD.cublog ~ sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim8) #--- #Tests of directed separation: # # Independ.Claim Test.Type DF Crit.Value P.Value # estUnevenness ~ CrossShelfPos + ... coef 36 -0.0587 0.9535 # DD.cublog ~ CrossShelfPos + ... coef 36 0.8733 0.3883 # cvpop ~ CrossShelfPos + ... coef 35 -0.3338 0.7405 # DD.cublog ~ latitude + ... coef 36 1.4688 0.1506 # altphi ~ latitude + ... coef 36 0.8001 0.4289 # cvpop ~ latitude + ... coef 35 0.1558 0.8771 # estUnevenness ~ sdHC + ... coef 36 1.7359 0.0911 # altphi ~ sdHC + ... coef 36 -1.2045 0.2363 # estRichness ~ sdHC + ... coef 35 1.2588 0.2164 # cvpop ~ sdHC + ... coef 35 0.8755 0.3873 # DD.cublog ~ estUnevenness + ... coef 35 1.0302 0.3100 # altphi ~ estUnevenness + ... coef 35 1.5528 0.1295 # estRichness ~ estUnevenness + ... coef 35 -0.8824 0.3836 # altphi ~ DD.cublog + ... coef 35 -0.3428 0.7338 # estRichness ~ DD.cublog + ... coef 34 0.5182 0.6077 # estRichness ~ altphi + ... coef 35 -0.6755 0.5038 # cvpop ~ altphi + ... coef 34 -2.0649 0.0466 * # cvpop ~ estRichness + ... coef 33 0.0240 0.9810 # #-- #Global goodness-of-fit: # #Chi-Squared = 25.605 with P-value = 0.109 and on 18 degrees of freedom #Fisher's C = 38.468 with P-value = 0.358 and on 36 degrees of freedom # #--- #Coefficients: # # Response Predictor Estimate Std.Error DF Crit.Value P.Value Std.Estimate # cvpop estUnevenness 0.5111 0.2064 36 2.4755 0.0181 0.3171 * # cvpop DD.cublog 0.0067 0.0016 36 4.0666 0.0002 0.5208 *** # altphi CrossShelfPos -8.0789 1.5491 37 -5.2151 0.0000 -0.6509 *** # estUnevenness latitude -0.0268 0.0033 37 -8.0529 0.0000 -0.7979 *** # estRichness latitude 0.0283 0.0079 36 3.5598 0.0011 0.4165 ** # estRichness CrossShelfPos 0.4233 0.0754 36 5.6181 0.0000 0.6573 *** # DD.cublog sdHC 12.2766 3.0352 37 4.0448 0.0003 0.5537 *** AIC(DivDDclSEM.AltPhi.trim8,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 72.468 101.611 17 39 DivDDclSEM.AltPhi.trim9 <- psem( lm(cvpop ~ DD.cublog, ecol.table.new), lm(altphi ~ CrossShelfPos, ecol.table.new), lm(estUnevenness ~ latitude,ecol.table.new), lm(estRichness ~ latitude+CrossShelfPos,ecol.table.new), lm(DD.cublog ~ sdHC,ecol.table.new), data=ecol.table.new ) summary(DivDDclSEM.AltPhi.trim9) AIC(DivDDclSEM.AltPhi.trim9,aicc=TRUE,AIC.type="dsep") # AIC AICc K n #1 80.935 105.662 16 39 # Worse. trim8 is the best model