Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Missing one subpopulation in the plot - subsetgrp #79

Open
mullerbsf opened this issue May 17, 2021 · 1 comment
Open

Missing one subpopulation in the plot - subsetgrp #79

mullerbsf opened this issue May 17, 2021 · 1 comment

Comments

@mullerbsf
Copy link

mullerbsf commented May 17, 2021

Dear,

I am trying to plot fastSTRUCTURE results using pophelper. However, I using a list of individuals with the subpopulation label (pop_info). The first subpopulation "BOP" has 159 individuals and is not showing in the plot. I don't know how to fix it. Could you please assist me with that?

Please find below the script:

Population label

pop_info = data.frame(read.delim("WiDiv_panel_pop_info.txt", header=F,stringsAsFactors=F))
head(pop_info)
V1
1 BOP
2 Popcorn
3 Mixed
4 Sweet_corn
5 Mixed
6 BOP

K12

K12 = function(x) attr(x,"k")
spnames = paste0("K=", sapply(slist,K12))
plotQ(qlist=slist[4], splab=spnames[4], splabface="bold", splabsize=4, clustercol=c("red","blue","purple","coral","cyan","chartreuse","plum1","darkgreen","grey","tomato2","yellow","black"), height=4,
grplab=pop_info, subsetgrp=c("BOP", "IDT", "Mixed", "NSS-Mo17", "NSS-Oh43", "Popcorn", "SS-B14", "SS-B37", "SS-B73", "SS-BSSSC0", "Sweet_corn", "Tropical"), ordergrp=T, showlegend=T, legendtextsize=4,
grplabangle=-90, grplabjust=0, grplabpos=0.8, linepos=1, legendlab=c("K1","K2","K3","K4","K5","K6","K7","K8","K9","K10","K11","K12"),
showtitle=T, showsubtitle=T, titlelab="STRUCTURE", subtitlelab=expression("Population structure in" ~ italic("Zea mays") ~ "WiDiv panel"),
outputfilename="WiDiv_panel_Maize_543samples_Pop_K12", imgtype="png")

Plot

WiDiv_panel_Maize_543samples_Pop_K12

Thank you!

@mullerbsf
Copy link
Author

@royfrancis could you please help me with this issue?

Thank you!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant