Closed BerzhanKurmanov closed 6 years ago
Dear Berzhan, I can't say, because I can't reproduce your error without some more information from you. It would good if you can send me the code and the tree / data you used, so that I can reproduce the error (see also https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example). Kind regards, Klaus
This can send me the data also of list to: klaus dot schliep at gmail dot com
Dear Klaus, Thank you for quick response! Here is the code. Also I attached my data file in csv format.
First <- read.csv("…\MLVA_repeat_data.csv") rownames(first) <- first[,1] second <- first[,-1] third <- data.frame(t(second)) yersinia <- as.phyDat(third, type ='USER', levels = c(1:47)) dm.yersinia = dist.hamming(yersinia) UPGMA.yersinia <- upgma(dm.yersinia) NJ.yersinia <- nj(dm.yersinia) plot(NJ.yersinia) plot(UPGMA.yersinia) parsim.all <- optim.parsimony(NJ.yersinia, data = yersinia)
After this I received that error notification.
Thank you in advance!
Best regards, Berzhan
From: Klaus Schliep notifications@github.com Sent: Monday, August 27, 2018 8:38 PM To: KlausVigo/phangorn phangorn@noreply.github.com Cc: Kurmanov,Berzhan b.kurmanov@ufl.edu; Author author@noreply.github.com Subject: Re: [KlausVigo/phangorn] Issue with maximum parsimony analysis (#60)
Dear Berzhan, I can't say, because I can't reproduce your error without some more information from you. It would good if you can send me the code and the tree / data you used, so that I can reproduce the error (see also https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-examplehttps://urldefense.proofpoint.com/v2/url?u=https-3A__stackoverflow.com_questions_5963269_how-2Dto-2Dmake-2Da-2Dgreat-2Dr-2Dreproducible-2Dexample&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=Unqeu_jfAebaglEj1biE4JUYVit_v71q9rnh0Z0vfcg&s=u3keymNLX4gyOOmVKqVSgBjqGLx8CB4wlHMMIyugPBg&e=). Kind regards, Klaus
This can send me the data also of list to: klaus dot schliep at gmail dot com
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_KlausVigo_phangorn_issues_60-23issuecomment-2D416412792&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=Unqeu_jfAebaglEj1biE4JUYVit_v71q9rnh0Z0vfcg&s=-h-Nq_CyuYI_IL8NVVBq7EbbVP_s1sawvGSzfEqIe1E&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_Al-5FdC1jN0sydIeci8nlGtivlruSCit-5F5ks5uVJD9gaJpZM4WOEk4&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=Unqeu_jfAebaglEj1biE4JUYVit_v71q9rnh0Z0vfcg&s=UGPJ8n10n2OyPEmEpE9IsxM9OWJd9V4Tw_qqbtd1IqQ&e=.
Hi @BerzhanKurmanov can you send me the csv data to my email? Regards, Klaus
To which address?
From: Klaus Schliep notifications@github.com Sent: Tuesday, August 28, 2018 1:31 PM To: KlausVigo/phangorn phangorn@noreply.github.com Cc: Kurmanov,Berzhan b.kurmanov@ufl.edu; Mention mention@noreply.github.com Subject: Re: [KlausVigo/phangorn] Issue with maximum parsimony analysis (#60)
Hi @BerzhanKurmanovhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_BerzhanKurmanov&d=DwMCaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=9_k8e7nG77mJcDKLtsY1eiciuWla7Qy5h__Flf0ZPoo&s=FpRdCqSOXTvHUWj1KUlbXTCD55IwY_QxnCzeKazeo1w&e= can you send me the csv data to my email? Regards, Klaus
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_KlausVigo_phangorn_issues_60-23issuecomment-2D416673260&d=DwMCaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=9_k8e7nG77mJcDKLtsY1eiciuWla7Qy5h__Flf0ZPoo&s=UgVsY7TUPsB5e2_RltzO2mF7yU5WtKTaqExbOSv_X5c&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_Al-5FdC0oxbnG0gTxTUuzW7T4pWzHNmrp-2Dks5uVX46gaJpZM4WOEk4&d=DwMCaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=9_k8e7nG77mJcDKLtsY1eiciuWla7Qy5h__Flf0ZPoo&s=URmu9XlpUjsHKA68via-nt3jkzRXw67C4ckBlwna0v0&e=.
klaus.schliep@gmail.com
Ok! Thank you!
From: Klaus Schliep notifications@github.com Sent: Tuesday, August 28, 2018 1:35 PM To: KlausVigo/phangorn phangorn@noreply.github.com Cc: Kurmanov,Berzhan b.kurmanov@ufl.edu; Mention mention@noreply.github.com Subject: Re: [KlausVigo/phangorn] Issue with maximum parsimony analysis (#60)
klaus.schliep@gmail.commailto:klaus.schliep@gmail.com
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_KlausVigo_phangorn_issues_60-23issuecomment-2D416674763&d=DwMCaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=5ZYq7PwqYpasFlt4X6QDK174L_ZOWcBWT5xSlkuSfAA&s=zwqmi4IP_woAk1S04GjdnAlF4WF_K47OphgEF3WkUZs&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_Al-5FdC-2Dr3SXnHp8DLpPwO44HFMKBXOQYsks5uVX9LgaJpZM4WOEk4&d=DwMCaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=5ZYq7PwqYpasFlt4X6QDK174L_ZOWcBWT5xSlkuSfAA&s=-ADVZW_vp-M80_5tG-o-zBp0avY1T9rw78gWttEF1jY&e=.
Hi @BerzhanKurmanov
at the moment it only possible 31 states for the fitch algorithm (which is faster), but can use the
sankoff algorithm which at the moment only uses a NNI rearrangement.
This code works for me:
parsim.all <- optim.parsimony(NJ.yersinia, data = yersinia, method = "sankoff")
Some of the sequences are identical
unique(yersinia) 58 sequences with 23 character and 23 different site patterns. The states are 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 yersinia 74 sequences with 23 character and 23 different site patterns. The states are 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 The latest version of phangorn will return a tree with multifurcations. Regards, Klaus
Oh, thank you very much! That works for me as well!
Best regards, Berzhan
From: Klaus Schliep notifications@github.com Sent: Tuesday, August 28, 2018 2:22 PM To: KlausVigo/phangorn phangorn@noreply.github.com Cc: Kurmanov,Berzhan b.kurmanov@ufl.edu; Mention mention@noreply.github.com Subject: Re: [KlausVigo/phangorn] Issue with maximum parsimony analysis (#60)
Hi @BerzhanKurmanovhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_BerzhanKurmanov&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=fkFr7u1n2hsdDiK2cGROn84bbP6hNph1IFjUF1E3luc&s=mNhJfrtApSEdP1bKBDzo45WJ3wMkpBOsoAfTvuWX3zI&e= at the moment it only possible 31 states for the fitch algorithm (which is faster), but can use the sankoff algorithm which at the moment only uses a NNI rearrangement. This code works for me: parsim.all <- optim.parsimony(NJ.yersinia, data = yersinia, method = "sankoff")
Some of the sequences are identical
unique(yersinia) 58 sequences with 23 character and 23 different site patterns. The states are 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 yersinia 74 sequences with 23 character and 23 different site patterns. The states are 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 The latest version of phangorn will return a tree with multifurcations. Regards, Klaus
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_KlausVigo_phangorn_issues_60-23issuecomment-2D416690129&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=fkFr7u1n2hsdDiK2cGROn84bbP6hNph1IFjUF1E3luc&s=_Z9kCKpc_19dG8TnVrRtB950aJNxCnvB2tNmGpAFExg&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_Al-5FdCzDJr-5FOszsa7LzbMWmCm9bMvJ1ISks5uVYovgaJpZM4WOEk4&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=fkFr7u1n2hsdDiK2cGROn84bbP6hNph1IFjUF1E3luc&s=SoIzNp5x5eMU5N9VlqxTR5fihA6gB3DWbWqtFb38Vuc&e=.
Dear Klaus,
Could you please help me with one issue? When I use bootstrap.pml function in phangorn package, I receive this notification:
Error in random.addition(x) : NAs in foreign function call (arg 1) In addition: Warning messages: 1: In if (!is.na(tmp)) { : the condition has length > 1 and only the first element will be used 2: In if (tmp == 1) { : the condition has length > 1 and only the first element will be used 3: In if (tmp == 2) do_rearr <- extras$rearrangement %in% c("NNI", : the condition has length > 1 and only the first element will be used 4: In prepareDataFitch(data) : NAs introduced by coercion to integer range
What can be the reason for this and how can I solve this problem? If you need I can send you my dataframe and the code. looking forward to your reply!
Best regards, Berzhan.
From: Klaus Schliep notifications@github.com Sent: Tuesday, August 28, 2018 2:22 PM To: KlausVigo/phangorn phangorn@noreply.github.com Cc: Kurmanov,Berzhan b.kurmanov@ufl.edu; Mention mention@noreply.github.com Subject: Re: [KlausVigo/phangorn] Issue with maximum parsimony analysis (#60)
Hi @BerzhanKurmanovhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_BerzhanKurmanov&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=fkFr7u1n2hsdDiK2cGROn84bbP6hNph1IFjUF1E3luc&s=mNhJfrtApSEdP1bKBDzo45WJ3wMkpBOsoAfTvuWX3zI&e= at the moment it only possible 31 states for the fitch algorithm (which is faster), but can use the sankoff algorithm which at the moment only uses a NNI rearrangement. This code works for me: parsim.all <- optim.parsimony(NJ.yersinia, data = yersinia, method = "sankoff")
Some of the sequences are identical
unique(yersinia) 58 sequences with 23 character and 23 different site patterns. The states are 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 yersinia 74 sequences with 23 character and 23 different site patterns. The states are 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 The latest version of phangorn will return a tree with multifurcations. Regards, Klaus
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_KlausVigo_phangorn_issues_60-23issuecomment-2D416690129&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=fkFr7u1n2hsdDiK2cGROn84bbP6hNph1IFjUF1E3luc&s=_Z9kCKpc_19dG8TnVrRtB950aJNxCnvB2tNmGpAFExg&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_Al-5FdCzDJr-5FOszsa7LzbMWmCm9bMvJ1ISks5uVYovgaJpZM4WOEk4&d=DwMFaQ&c=pZJPUDQ3SB9JplYbifm4nt2lEVG5pWx2KikqINpWlZM&r=RbfmuMpRAj-rHk7fjLsRQjkYUzwT44SUG_4gLkxQUdw&m=fkFr7u1n2hsdDiK2cGROn84bbP6hNph1IFjUF1E3luc&s=SoIzNp5x5eMU5N9VlqxTR5fihA6gB3DWbWqtFb38Vuc&e=.
Hi @BerzhanKurmanov, what's your script looking like? Klaus
This is my script: mlva<-read.csv("\\file.epi.ufl.edu\home\b.kurmanov\MY FOLDER\ MLVA_data_repeats.csv") rownames(mlva)<-mlva[,1] vntrs<-mlva[,-1] vntrst<-data.frame(t(vntrs)) analysisA<-as.phyDat(vntrst, type='USER', levels=c(0:54)) distA<-dist.hamming(analysisA) upgma.test<-upgma(distA) nj.test<-nj(distA) parsim.all<-optim.parsimony(nj.test, data=analysisA, method="sankoff") fit.4 = pml(nj.test, data=analysisA) fitGTR.4 <- optim.pml(fit.4, optInv=TRUE, optGamma=TRUE, rearrangement = "NNI", control = pml.control(trace = 1)) fitGTR.4 <- optim.pml(fitGTR.4, rearrangement = "stochastic", ratchet.par = list(iter = 5L, maxit = 5L, prop = 1/3),trace=0) fitGTR.4 <- optim.pml(fitGTR.4, rearrangement = "stochastic", multicore=T,trace=0) bs.pall.ref <- bootstrap.pml(fitGTR.4, bs=100, optNni=TRUE, control = pml.control(trace=0))
From: Klaus Schliep notifications@github.com Sent: Friday, July 12, 2019 3:09 PM To: KlausVigo/phangorn phangorn@noreply.github.com Cc: Kurmanov,Berzhan b.kurmanov@ufl.edu; Mention mention@noreply.github.com Subject: Re: [KlausVigo/phangorn] Issue with maximum parsimony analysis (#60)
Hi @BerzhanKurmanovhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_BerzhanKurmanov&d=DwMCaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=sr-lh5wEgoMxtlbsTKjNgurfkzUYoWYj-29N4qSIw2c&m=8SgZ8xqLb5rwT9RmE389PijglljM6ugj8N3_RJIfHA4&s=13yp5VUMtLiB60etj8mMimPyNE2G6qMnAsNWjzBFLF0&e=, what's your script looking like? Klaus
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_KlausVigo_phangorn_issues_60-3Femail-5Fsource-3Dnotifications-26email-5Ftoken-3DAJP52C3DR3FZ33GSTMR36DTP7DJDTA5CNFSM4FRYJE4KYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGODZ2T4QA-23issuecomment-2D511000128&d=DwMCaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=sr-lh5wEgoMxtlbsTKjNgurfkzUYoWYj-29N4qSIw2c&m=8SgZ8xqLb5rwT9RmE389PijglljM6ugj8N3_RJIfHA4&s=wQxxLE618wE0Ha-HGKJvU1_j5FrAenot1UNhb9WaOUQ&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_AJP52C6PVC7CA3KBQGDMNNTP7DJDTANCNFSM4FRYJE4A&d=DwMCaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=sr-lh5wEgoMxtlbsTKjNgurfkzUYoWYj-29N4qSIw2c&m=8SgZ8xqLb5rwT9RmE389PijglljM6ugj8N3_RJIfHA4&s=ZH9ITBjIuG8-XyNvY64Sdk3is1nAm9gJ-_KDg27us1A&e=.
Good morning dear Klaus, Could you tell me, please, what was my mistake? Does my script wrong?
Best regards, Berzhan.
From: Kurmanov,Berzhan Sent: Friday, July 12, 2019 3:12 PM To: 'KlausVigo/phangorn' reply@reply.github.com; KlausVigo/phangorn phangorn@noreply.github.com Cc: Mention mention@noreply.github.com Subject: RE: [KlausVigo/phangorn] Issue with maximum parsimony analysis (#60)
This is my script: mlva<-read.csv("\\file.epi.ufl.edu\home\b.kurmanov\MY FOLDER\ MLVA_data_repeats.csvfile://file.epi.ufl.edu/home/b.kurmanov/MY%20FOLDER/%20MLVA_data_repeats.csv") rownames(mlva)<-mlva[,1] vntrs<-mlva[,-1] vntrst<-data.frame(t(vntrs)) analysisA<-as.phyDat(vntrst, type='USER', levels=c(0:54)) distA<-dist.hamming(analysisA) upgma.test<-upgma(distA) nj.test<-nj(distA) parsim.all<-optim.parsimony(nj.test, data=analysisA, method="sankoff") fit.4 = pml(nj.test, data=analysisA) fitGTR.4 <- optim.pml(fit.4, optInv=TRUE, optGamma=TRUE, rearrangement = "NNI", control = pml.control(trace = 1)) fitGTR.4 <- optim.pml(fitGTR.4, rearrangement = "stochastic", ratchet.par = list(iter = 5L, maxit = 5L, prop = 1/3),trace=0) fitGTR.4 <- optim.pml(fitGTR.4, rearrangement = "stochastic", multicore=T,trace=0) bs.pall.ref <- bootstrap.pml(fitGTR.4, bs=100, optNni=TRUE, control = pml.control(trace=0))
From: Klaus Schliep notifications@github.com<mailto:notifications@github.com> Sent: Friday, July 12, 2019 3:09 PM To: KlausVigo/phangorn phangorn@noreply.github.com<mailto:phangorn@noreply.github.com> Cc: Kurmanov,Berzhan b.kurmanov@ufl.edu<mailto:b.kurmanov@ufl.edu>; Mention mention@noreply.github.com<mailto:mention@noreply.github.com> Subject: Re: [KlausVigo/phangorn] Issue with maximum parsimony analysis (#60)
Hi @BerzhanKurmanovhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_BerzhanKurmanov&d=DwMCaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=sr-lh5wEgoMxtlbsTKjNgurfkzUYoWYj-29N4qSIw2c&m=8SgZ8xqLb5rwT9RmE389PijglljM6ugj8N3_RJIfHA4&s=13yp5VUMtLiB60etj8mMimPyNE2G6qMnAsNWjzBFLF0&e=, what's your script looking like? Klaus
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_KlausVigo_phangorn_issues_60-3Femail-5Fsource-3Dnotifications-26email-5Ftoken-3DAJP52C3DR3FZ33GSTMR36DTP7DJDTA5CNFSM4FRYJE4KYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGODZ2T4QA-23issuecomment-2D511000128&d=DwMCaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=sr-lh5wEgoMxtlbsTKjNgurfkzUYoWYj-29N4qSIw2c&m=8SgZ8xqLb5rwT9RmE389PijglljM6ugj8N3_RJIfHA4&s=wQxxLE618wE0Ha-HGKJvU1_j5FrAenot1UNhb9WaOUQ&e=, or mute the threadhttps://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_notifications_unsubscribe-2Dauth_AJP52C6PVC7CA3KBQGDMNNTP7DJDTANCNFSM4FRYJE4A&d=DwMCaQ&c=sJ6xIWYx-zLMB3EPkvcnVg&r=sr-lh5wEgoMxtlbsTKjNgurfkzUYoWYj-29N4qSIw2c&m=8SgZ8xqLb5rwT9RmE389PijglljM6ugj8N3_RJIfHA4&s=ZH9ITBjIuG8-XyNvY64Sdk3is1nAm9gJ-_KDg27us1A&e=.
Dear Klaus, I've tried to perform a maximum parsimony analysis using phangorn package, but have obtained such a warning: Error in optim.fitch(tree = tree, data = data, trace = trace, rearrangements = rearrangements, : NAs in foreign function call (arg 1) In addition: Warning message: In prepareDataFitch(data) : NAs introduced by coercion to integer range
What is the reason for this? How can I solve this problem?
Best regards, Berzhan