I would like to calculate some summary statistics such as FIS, FST, Ho, He for "each population" using a hierfstat object with 217 samples and 38,000 sites. Although the basic.stats() function estimated FIS, Ho, Hs for each population, FST and other statistics are not estimated for each population but only with the overall values. Is it possible to calculate the population-specific FST for each population using hierfstat in R?
Thanks for any responses and suggestions!
I tried separating samples based on population first and calculated FST and FIS for each population, but I am not sure if it makes sense.
How to calculate basic.stats of each population using hierfstat in R
54 Views Asked by Pei-Wei Sun At
1
There are 1 best solutions below
Related Questions in R
- How to make an R Shiny app with big data?
- How do I keep only specific rows based on whether a column has a specific value?
- Likert scale study - ordinal regression model
- Extract a table/matrix from R into Excel with same colors and stle
- How can I solve non-conformable arguments in R netmeta::discomb (Error in B.matrix %*% C.matrix)?
- Can raw means and estimated marginal means be the same ? And when?
- Understanding accumulate function when .dir is set to "backwards"
- Error in if (nrow(peaks) > 0) { : argument is of length zero Calls: CopywriteR ... tryCatch -> tryCatchList -> tryCatchOne -> <Anonymous> Execution ha
- How to increase quality of mathjax output?
- Convert the time intervals to equal hours and fill in the value column
- How to run an R function getpoints() from IPDfromKM package in an R shiny app which in R pops up a plot that utilizes clicks to capture coordinates?
- Replace NA in list of dfs in certain columns and under certain conditions
- R and text on Cyrillic
- The ts() function in R is returning the correct start and frequency but not end value which is 1 and not 179
- TROUBLING with the "DROP_NA" Function
Related Questions in FST
- How to calculate basic.stats of each population using hierfstat in R
- Deserializing byte array into Java object in a different JVM
- Pass string into map_dfr to filter using tidyfst
- Load R data objects' attributes without loading object from file?
- read.fst() crashes R : workaround needed to detect corrupted file before read.fst()
- union two FSA's with different symbol lists
- How to understand this formula in Direct Construction of Minimal Acyclic Subsequential Transducers?
- Does the R arrow package have anything like the random access capability of the fst package?
- Saving and readings lubridate intervals to/from disk
- There is a set of pairs (set_pairs) and I want to create a set (set_fsts) of first elements of those pairs. How can I do this in Isabelle?
- Fastest way to search for a long list of patterns in a text
- Detect cyclic feeding interactions without applying XFST replace rules to lexicon
- Is there any downside/risk to using fst to compress a serialized list?
- Migrate lucene FST files from 5.1.0 to 8.9.0
- Kaldi error while running script prepare_lang.sh -- GenericRegister::GetEntry: No such file or directory
Trending Questions
- UIImageView Frame Doesn't Reflect Constraints
- Is it possible to use adb commands to click on a view by finding its ID?
- How to create a new web character symbol recognizable by html/javascript?
- Why isn't my CSS3 animation smooth in Google Chrome (but very smooth on other browsers)?
- Heap Gives Page Fault
- Connect ffmpeg to Visual Studio 2008
- Both Object- and ValueAnimator jumps when Duration is set above API LvL 24
- How to avoid default initialization of objects in std::vector?
- second argument of the command line arguments in a format other than char** argv or char* argv[]
- How to improve efficiency of algorithm which generates next lexicographic permutation?
- Navigating to the another actvity app getting crash in android
- How to read the particular message format in android and store in sqlite database?
- Resetting inventory status after order is cancelled
- Efficiently compute powers of X in SSE/AVX
- Insert into an external database using ajax and php : POST 500 (Internal Server Error)
Popular # Hahtags
Popular Questions
- How do I undo the most recent local commits in Git?
- How can I remove a specific item from an array in JavaScript?
- How do I delete a Git branch locally and remotely?
- Find all files containing a specific text (string) on Linux?
- How do I revert a Git repository to a previous commit?
- How do I create an HTML button that acts like a link?
- How do I check out a remote Git branch?
- How do I force "git pull" to overwrite local files?
- How do I list all files of a directory?
- How to check whether a string contains a substring in JavaScript?
- How do I redirect to another webpage?
- How can I iterate over rows in a Pandas DataFrame?
- How do I convert a String to an int in Java?
- Does Python have a string 'contains' substring method?
- How do I check if a string contains a specific word?
If you are still looking for an answer, here's what worked for me: After creating a geno_object with the df2genind function, do the following;
This will give you a pairwise Fst.
Check this link for more details.