Skip to content

Commit

Permalink
Updating tutorial
Browse files Browse the repository at this point in the history
  • Loading branch information
lvclark committed Jan 9, 2018
1 parent 72ae496 commit 2a098de
Show file tree
Hide file tree
Showing 4 changed files with 45 additions and 1 deletion.
2 changes: 1 addition & 1 deletion DESCRIPTION
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
Package: polyRAD
Version: 0.1
Date: 2019-01-02
Date: 2019-01-09
Title: Genotype Calling with Uncertainty from Sequencing Data in Polyploids and Diploids
Authors@R: c(person("Lindsay V.", "Clark", email = "[email protected]",
role = c("aut", "cre")),
Expand Down
4 changes: 4 additions & 0 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,10 @@ install.packages("polyRAD_0.x.zip", repos = NULL)

replacing the file name with the one you just downloaded.

## Tutorial

The tutorial document for the package is available [on Github](https://github.com/lvclark/polyRAD/raw/master/vignettes/polyRADtutorial.pdf).

## Funding

This material is based upon work supported by the National Science Foundation under Grant No.
Expand Down
40 changes: 40 additions & 0 deletions vignettes/polyRADtutorial.Rnw
Original file line number Diff line number Diff line change
Expand Up @@ -32,6 +32,46 @@ A single ``RADdata'' object contains the entire dataset of read depth and
locus information, as well as parameters that are estimated during the
course of analysis.

\section{Summary of available functions}

For any function named in this section, see its help page for more information.
(For example by typing \texttt{?VCF2RADdata} into the R console.)

Several functions are available for import of read depth data and (optionally)
alignment information into a RADdata object:

\begin{itemize}
\item \texttt{VCF2RADdata}
\item \texttt{readTagDigger}
\item \texttt{readStacks1}
\item \texttt{readHMC}
\end{itemize}

More generally, the \texttt{RADdata} function is used for constructing RADdata
objects; see the help page for that function for more information on what
data are needed.

Several pipelines are available for genotype estimation, depending on how the
population is structured (i.e. what the genotype prior probabilities should be.):

\begin{itemize}
\item \texttt{PipelineMapping2Parents}
\item \texttt{IterateHWE}
\item \texttt{IteratePopStruct}
\end{itemize}

Lastly, for exporting the estimated genotypes to other software:

\begin{itemize}
\item \texttt{ExportGAPIT}
\item \texttt{Export\_rrBLUP\_Amat}
\item \texttt{Export\_rrBLUP\_GWAS}
\end{itemize}

If you need continuous numerical genotypes exported in some other format, see
\texttt{GetWeightedMeanGenotypes}. Also, \texttt{GetLikelyGen} returns the
most likely genotypes for a single sample.

\section{Estimating genotype probabilities in a mapping population}

In this example, we'll import some data from an F1 mapping population
Expand Down
Binary file modified vignettes/polyRADtutorial.pdf
Binary file not shown.

0 comments on commit 2a098de

Please sign in to comment.