[Genabel-commits] r1624 - pkg/ProbABEL/doc

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Wed Feb 26 22:43:42 CET 2014


Author: lckarssen
Date: 2014-02-26 22:43:42 +0100 (Wed, 26 Feb 2014)
New Revision: 1624

Modified:
   pkg/ProbABEL/doc/ProbABEL_manual.tex
Log:
Small changes to the section "Basic Analysis options" in the ProbABEL manual. The example command lines are now a bit simpler and also include the "gtdata" subdirectory we have in the examples subdir.


Modified: pkg/ProbABEL/doc/ProbABEL_manual.tex
===================================================================
--- pkg/ProbABEL/doc/ProbABEL_manual.tex	2014-02-23 22:16:53 UTC (rev 1623)
+++ pkg/ProbABEL/doc/ProbABEL_manual.tex	2014-02-26 21:43:42 UTC (rev 1624)
@@ -337,33 +337,29 @@
 However, for a simple run only three options are mandatory, which
 specify the necessary files needed to run the regression analysis.
 
-These options are
-\texttt{--dose} (or \texttt{-d}),
-specifying the genomic predictor/MLDOSE file described in section \ref{ssec:dosein};
-\texttt{--pheno} (or \texttt{-p}),
-specifying the phenotypic data file described in section \ref{ssec:phenoin}; and
-\texttt{--info} (or \texttt{-i}),
-specifying the SNP information file described in section \ref{ssec:infoin}.
+These options are \texttt{--dose} (or \texttt{-d}), specifying the
+genomic predictor/MLDOSE file described in section \ref{ssec:dosein};
+\texttt{--pheno} (or \texttt{-p}), specifying the phenotypic data file
+described in section \ref{ssec:phenoin}; and \texttt{--info} (or
+\texttt{-i}), specifying the SNP information file described in section
+\ref{ssec:infoin}.
 
 If you change to the \texttt{examples} directory you can run
 an analysis of height by running
 \begin{verbatim}
-user at server:~/ProbABEL/examples/$ ../bin/palinear -p height.txt \
-                                 -d test.mldose -i test.mlinfo
+palinear -p height.txt -d gtdata/test.mldose -i gtdata/test.mlinfo
 \end{verbatim}
-Output from the analysis will be directed to the
+Output from the analysis will be stored in the
 \texttt{regression.out.csv} file.
-
 The analysis of a binary trait (e.g.~chd) can be run with
 \begin{verbatim}
-user at server:~/ProbABEL/examples/$ ../bin/palogist -p logist_data.txt \
-                                 -d test.mldose -i test.mlinfo
+palogist -p logist_data.txt -d gtdata/test.mldose \
+    -i gtdata/test.mlinfo
 \end{verbatim}
-
 To run a Cox proportional hazards model, try
 \begin{verbatim}
-user at server:~/ProbABEL/examples/$ ../bin/pacoxph -p coxph_data.txt \
-                                 -d test.mldose -i test.mlinfo
+pacoxph -p coxph_data.txt -d gtdata/test.mldose \
+    -i gtdata/test.mlinfo
 \end{verbatim}
 
 Please have a look at the shell script files \texttt{example\_qt.sh},
@@ -372,13 +368,22 @@
 
 To run an analysis with MLPROB files, you need specify the MLPROB file
 with the \texttt{-d} option and also specify that there are two
-genetic predictors per SNP, e.g.~you can run linear model with
+genetic predictors per SNP, e.g.~you can run a linear model with
 \begin{verbatim}
-user at server:~/ProbABEL/examples/$ ../bin/palinear -p height.txt \
-                                 -d test.mlprob -i test.mlinfo \
-                                 --ngpreds=2
+palinear -p height.txt -d gtdata/test.mlprob -i gtdata/test.mlinfo \
+    --ngpreds=2
 \end{verbatim}
 
+When using genomic predictor files (dosages or probabilities) stored
+in filevector (a.k.a.~DatABEL) format (i.e.~a combination of
+\texttt{.fvi} and \texttt{.fvd} files) you can specify these like you
+would with ordinary text files. This is how the previous example would
+change:
+\begin{verbatim}
+palinear -p height.txt -d gtdata/test.mlprob.fvi -i gtdata/test.mlinfo \
+    --ngpreds=2
+\end{verbatim}
+
 \subsection{Advanced analysis options}
 The option \texttt{--interaction} allows you to include interaction
 between SNPs and any covariate. If for example your model is



More information about the Genabel-commits mailing list