Hi everyone,
I'm having trouble running a PGLS model with the package "AICmodavg". I
continuously get the error of false convergence with certain Lambda
values (even
when trying to run the model with different Lambda values) and for other La
mbda values I run into "error in eigen(val) : infinite or missing values in
'X' ". I've tried several optimizers and removing some outlier values but
the same errors keep on popping up.. Does anyone know how to solve this
problem?
Below a part of my script with the specific dependent and independent varia
bles and I've also attached files with the relevant data and phylogenetic
tree information.
Cand.models = list()
niter = 100
for (i in 1:niter) {
Cand.models[[i]] = gls(Maxgrowthrate ~ log.Forrang+log.Weight_av), data =
d, method= "ML", na.action=na.omit
correlation = corPagel(value=0.4, trees[[i]]))
}
Thank you in advance for any help, it's much appreciated!
Kind regards,
Hannah van Noort
"Maxgrowthrate" "Forrang" "Weight_av"
Diomedea_exulans 0.6600 1500.0000 8.6500
Diomedea_epomophora 0.9200 1250.0000 9.0000
Thalassarche_chrysostoma 1.8400 250.0000 3.3250
Thalassarche_bulleri na 1200.0000 2.9100
Thalassarche_melanophrys 2.1700 470.0000 3.5700
Phoebetria_palpebatra 1.8400 1350.0000 2.9500
Macronectes_giganteus 1.7300 200.0000 4.5650
Fulmarus_glacialis 5.0800 500.0000 0.6150
Bulweria_bulwerii na 200.0000 0.9900
Procellaria_aequinoctialis 1.8000 7000.0000 1.3350
Procellaria_parkinsoni na 200.0000 0.7100
Puffinus_tenuirostris 2.2400 50.0000 0.5850
Puffinus_puffinus 2.8600 450.0000 0.4200
Oceanites_oceanicus 8.0000 250.0000 0.0350
Oceanodroma_leucorhoa 3.5600 615.0000 0.0450
Phaethon_rubricauda 2.2400 200.0000 0.6700
Pelecanus_occidentalis 3.6000 20.0000 3.4500
Fregata_magnificens 0.9000 500.0000 1.4750
Morus_bassanus 3.0000 300.0000 3.0000
Sula_dactylatra 1.7000 50.0000 1.9900
Sula_leucogaster 2.0000 40.0000 1.3500
Sula_sula 1.7000 40.0000 1.1500
Phalacrocorax_carbo 3.6500 7.0000 2.7500
Phalacrocorax_aristotelis 3.1100 9.0000 1.7700
Larus_dominicanus 2.0000 10.0000 1.0250
Larus_occidentalis 2.5000 75.0000 1.1000
Larus_argentatus 3.2000 21.0000 1.1350
Rissa_tridactyla 5.0000 36.0000 0.4000
Sterna_bergii 3.0000 40.0000 0.3600
Sterna_sandvicensis 4.5000 20.0000 0.2200
Sterna_dougallii 6.4000 25.0000 0.1250
Sterna_hirundo 6.2000 15.2000 0.1350
Sterna_paradisaea 5.9000 20.0000 0.1180
Sterna_fuscata 2.2000 500.0000 0.2000
Anous_stolidus 3.6000 180.0000 0.1725
Anous_minutus 5.0000 10.0000 0.1210
Cepphus_columba 4.0000 7.0000 0.5000
Cepphus_grylle 2.9000 4.0000 0.4700
Uria_lomvia 1.0500 170.0000 1.1000
Alca_torda 3.1000 12.0000 0.7200
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
