[Genabel-commits] r2072 - in pkg/MultiABEL: R man

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Thu Jun 15 17:47:52 CEST 2017


Author: shenxia
Date: 2017-06-15 17:47:52 +0200 (Thu, 15 Jun 2017)
New Revision: 2072

Modified:
   pkg/MultiABEL/R/MultiLoad.R
   pkg/MultiABEL/R/Multivariate.R
   pkg/MultiABEL/man/MultiLoad.Rd
   pkg/MultiABEL/man/Multivariate.Rd
Log:
Uncomment some examples

Modified: pkg/MultiABEL/R/MultiLoad.R
===================================================================
--- pkg/MultiABEL/R/MultiLoad.R	2017-06-15 15:35:20 UTC (rev 2071)
+++ pkg/MultiABEL/R/MultiLoad.R	2017-06-15 15:47:52 UTC (rev 2072)
@@ -34,7 +34,7 @@
 #' \code{\link{Multivariate}}
 #' 
 #' @examples 
-#' \dontrun{
+#' 
 #' ## loading example gwaa.data in GenABEL
 #' require(GenABEL)
 #' data(ge03d2ex.clean)
@@ -43,6 +43,7 @@
 #' loaded <- MultiLoad(gwaa.data = ge03d2ex.clean, trait.cols = c(5, 6, 8), 
 #'                     covariate.cols = c(2, 3))
 #' 
+#' \dontrun{
 #' ## converting the same dataset into DatABEL format files
 #' require(DatABEL)
 #' write.table(phdata(ge03d2ex.clean), 'pheno.txt', col.names = TRUE, row.names = TRUE, 

Modified: pkg/MultiABEL/R/Multivariate.R
===================================================================
--- pkg/MultiABEL/R/Multivariate.R	2017-06-15 15:35:20 UTC (rev 2071)
+++ pkg/MultiABEL/R/Multivariate.R	2017-06-15 15:47:52 UTC (rev 2072)
@@ -30,17 +30,18 @@
 #' \code{\link{MultiLoad}}
 #' 
 #' @examples 
-#' \dontrun{
+#' 
 #' ## loading example gwaa.data in GenABEL
+#' require(GenABEL)
 #' data(ge03d2ex.clean)
 #' 
 #' ## running multivariate GWAS for 3 traits: height, weight, bmi
 #' loaded <- MultiLoad(gwaa.data = ge03d2ex.clean, trait.cols = c(5, 6, 8), 
 #'                     covariate.cols = c(2, 3))
 #' 
-#' ## running the multivariate GWAS again
+#' ## running the multivariate GWAS 
 #' res <- Multivariate(loaded)
-#' }
+#' 
 #' @aliases Multivariate, multivariate
 #' @keywords multivariate
 #' 
@@ -75,7 +76,7 @@
 	scan <- .Fortran('MultiSummaryLoopPrecise', k = as.integer(k), m = as.integer(m), nn = as.numeric(n), 
 			varg = x$cvg, betamat = betamat, R = x$cor.pheno, invR = solve(x$cor.pheno), D = matrix(0, m, m),
 			sdY = diag(sqrt(x$var.pheno)), invsdY = diag(1/sqrt(x$var.pheno)), sY = sqrt(x$var.pheno), 
-			b = betamat, s = betamat, pil = numeric(k), coef = betamat, PACKAGE = "MultiABEL")
+			b = betamat, s = betamat, pil = numeric(k), coef = betamat, ss = betamat, PACKAGE = "MultiABEL")
 	cat(' OK\n')
 	res$Beta.S <- scan$pil
 	scan$fstat <- res$Beta.S/m/(1 - res$Beta.S)*(n - m - 1)

Modified: pkg/MultiABEL/man/MultiLoad.Rd
===================================================================
--- pkg/MultiABEL/man/MultiLoad.Rd	2017-06-15 15:35:20 UTC (rev 2071)
+++ pkg/MultiABEL/man/MultiLoad.Rd	2017-06-15 15:47:52 UTC (rev 2072)
@@ -46,7 +46,7 @@
 DatABEL format input, individual IDs in \code{phenofile} and \code{genofile} have to match!
 }
 \examples{
-\dontrun{
+
 ## loading example gwaa.data in GenABEL
 require(GenABEL)
 data(ge03d2ex.clean)
@@ -55,6 +55,7 @@
 loaded <- MultiLoad(gwaa.data = ge03d2ex.clean, trait.cols = c(5, 6, 8), 
                     covariate.cols = c(2, 3))
 
+\dontrun{
 ## converting the same dataset into DatABEL format files
 require(DatABEL)
 write.table(phdata(ge03d2ex.clean), 'pheno.txt', col.names = TRUE, row.names = TRUE, 

Modified: pkg/MultiABEL/man/Multivariate.Rd
===================================================================
--- pkg/MultiABEL/man/Multivariate.Rd	2017-06-15 15:35:20 UTC (rev 2071)
+++ pkg/MultiABEL/man/Multivariate.Rd	2017-06-15 15:47:52 UTC (rev 2072)
@@ -34,18 +34,19 @@
 DatABEL format input, individual IDs in \code{phenofile} and \code{genofile} have to match!
 }
 \examples{
-\dontrun{
+
 ## loading example gwaa.data in GenABEL
+require(GenABEL)
 data(ge03d2ex.clean)
 
 ## running multivariate GWAS for 3 traits: height, weight, bmi
 loaded <- MultiLoad(gwaa.data = ge03d2ex.clean, trait.cols = c(5, 6, 8), 
                     covariate.cols = c(2, 3))
 
-## running the multivariate GWAS again
+## running the multivariate GWAS 
 res <- Multivariate(loaded)
+
 }
-}
 \references{
 Xia Shen, ..., Jim Wilson, Gordan Lauc, Yurii Aulchenko (2015).
 Multi-omic-variate analysis identified novel loci associated with 



More information about the Genabel-commits mailing list