[Genabel-commits] r638 - in pkg/VariABEL: R src

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Tue Feb 8 10:53:26 CET 2011


Author: maksim
Date: 2011-02-08 10:53:26 +0100 (Tue, 08 Feb 2011)
New Revision: 638

Modified:
   pkg/VariABEL/R/var.test.gwaa.R
   pkg/VariABEL/src/tags
   pkg/VariABEL/src/var_homogeneity_test_C.cpp
   pkg/VariABEL/src/var_homogeneity_tests.cpp
   pkg/VariABEL/src/var_homogeneity_tests.h
Log:
almost ready

Modified: pkg/VariABEL/R/var.test.gwaa.R
===================================================================
--- pkg/VariABEL/R/var.test.gwaa.R	2011-02-07 19:34:13 UTC (rev 637)
+++ pkg/VariABEL/R/var.test.gwaa.R	2011-02-08 09:53:26 UTC (rev 638)
@@ -24,10 +24,9 @@
 #_____________________________________________________________________________________________________________________________
 #Available for user function for genome-wide testing variance homogeneity of trait's distribution.
 #
-"var.test.gwaa" <- function(formula, genodata, phenodata, testname="sqlm", analysis_type="AAvsABvsBB")
+"var.test.gwaa" <- function(formula, genodata, phenodata, genodata_info="", testname="sqlm", analysis_type="AAvsABvsBB")
 {
 
-
 if(testname == "bartlett") {testname <- 0}
 else if(testname == "levene") {testname <- 1}
 else if(testname == "likelihood") {testname <- 2}
@@ -47,56 +46,108 @@
 OUT <- "R" 
 FUN <- "variance_homogeneity_test_C_wrapper"
 
-design_matrix <- model.matrix(formula, data = phenodata)
-trat_name <- as.character(formula[[2]])
-trait <- phenodata[,trat_name]
-idnum <- length(trait)
-design_matrix_df <- data.frame(design_matrix)
-design_matrix_df$snp <- 0 #will be filled into iterator
 
 
+
+
+attach(phenodata, pos = 2, warn.conflicts = FALSE)
+
+idnum <- 0
+
+if(is(formula, "formula"))
+	{
+	design_matrix <- model.matrix(formula, data = phenodata)
+	trat_name <- as.character(formula[[2]])
+	trait <- phenodata[,trat_name]
+	idnum <- length(trait)
+	design_matrix_df <- data.frame(design_matrix)
+	design_matrix_df$snp <- 0 #will be filled into iterator
+	}
+	else if(is(formula, "numeric") || is(formula, "integer") || is(formula, "double"))
+	{
+	trait <- formula
+	idnum <- length(trait)
+	design_matrix_df <- data.frame(X.Intercept=rep(1,idnum), snp=rep(0,idnum))
+	}
+
+design_matrix_geno_means_df <- data.frame(X.Intercept=rep(1,idnum), snp=rep(0,idnum))
+
+
 p <- dim(design_matrix_df)[2] #covariates number + 1
 
-print(design_matrix_df)
 gtNrow <- NA
 gtNcol <- NA
 
+genodata_info_df <- data.frame()
+
+
 #Stolen from MixABEL:
 if(class(genodata) == "snp.data")
 	{
-	gtNrow <- dim(genodata)[1]
-	gtNcol <- dim(genodata)[2]
+	print("reading genotype info...")
+	genodata_info_df <- data.frame(name=genodata at snpnames,
+				 												 coding=as.character(genodata at coding),
+																 strand=as.character(genodata at strand),
+																 chromosome=as.character(genodata at chromosome),
+				 												 map=genodata at map
+																 )
+#	gtNrow <- genodata at nsnps
+#	gtNcol <- genodata at nids
+	gtNrow <- genodata at nids
+	gtNcol <- genodata at nsnps
 	genodata <- as.raw(genodata at gtps)
-	}
-	else if (class(genodata) == "matrix") 
-	{
-	gtNrow <- dim(genodata)[1]
-	gtNcol <- dim(genodata)[2]
-	storage.mode(genodata) <- "double"
-	}
+  }	
 	else if (class(genodata)=="databel") 
 	{
 	gtNrow <- dim(genodata)[1]
 	gtNcol <- dim(genodata)[2]
 	genodata <- genodata at data
+	
+	if(genodata_info != "" & file.exists(genodata_info))
+		{
+		print("reading genotype's info...")
+		genodata_info_df <- read.table(genodata_info, strings=F, header=T)
+		print(c(dim(genodata_info_df)[1], gtNcol))
+		if(dim(genodata_info_df)[1] != gtNcol) 
+			{
+			print(paste("warning: File ", genodata_info, " contains information about ", dim(genodata_info_df)[1],
+																											 " SNPs but genotype data contains ", gtNcol, ". Ignoring of the info file.", sep=""))
+			}		
+		genodata_info_df <- data.frame(name=sub("^", "snp", as.character(1:gtNcol)))
+		}
+	else
+		{
+		print(paste("warning: File ", genodata_info, "does not exist. Trying to run analysis withouth SNP's information..."))
+		genodata_info_df <- data.frame(name=sub("^", "snp", as.character(1:gtNcol)))
+		}
+
 	}
 	else
  	{
-	stop(paste("genodata class not recognised ('",class(genodata),"')",sep=""))
+	stop(paste("genodata class not recognized ('",class(genodata),"')",sep=""))
 	}
 
-														 
 
+if(testname == 1) #Levene's
+	{
+	if(is(formula, "formula"))
+		{
+		trait_regression <- lm(formula, data=phenodata)
+		trait <- trait_regression$residuals
+		}
+	}	
 
+print("Start variance analysis...")
 results_C <- .Call("iterator", genodata,
 														 as.integer(gtNrow), as.integer(gtNcol),
 													 	 as.character(FUN),
 														 as.character(OUT), 
 														 as.integer(MAR),
 														 as.integer(1),
-														 as.integer(17), #iterator additional inputa parameters number
+														 as.integer(18), #iterator additional inputa parameters number
 														 as.double(trait),
 														 as.double(data.matrix(design_matrix_df)),
+														 as.double(data.matrix(design_matrix_geno_means_df)),
 														 as.integer(p),
 														 as.integer(analysis_type),
 														 as.integer(testname),
@@ -114,19 +165,23 @@
 														 integer(idnum)
 														 )
 
+print("Variance analysis done.")
 
-if(class(genodata) == "snp.data")
-	{
-	
-	}
+print(genodata_info_df)
+print(results_C)
 
-results_df <- data.frame(results_C)
+results_df <- data.frame(genodata_info_df, results_C)
+print(results_df)
 
 cov_names <- colnames(design_matrix_df)
-print("1")
-output_column_names <- c("chisq", "df", "Intercept_effect", "Intercept_sd")
-print("2")
 
+print(cov_names)
+
+output_column_names <- c(colnames(genodata_info_df), "chisq", "df", "Intercept_effect", "Intercept_sd")
+
+
+print(output_column_names)
+
 for(i in 2:p)
 	{
 	output_column_names <- c(output_column_names, paste(cov_names[i], "_effect", sep=""))
@@ -135,7 +190,7 @@
 
 colnames(results_df) <- output_column_names
 
-print(results_df)
+results_df
 }
 
 

Modified: pkg/VariABEL/src/tags
===================================================================
--- pkg/VariABEL/src/tags	2011-02-07 19:34:13 UTC (rev 637)
+++ pkg/VariABEL/src/tags	2011-02-08 09:53:26 UTC (rev 638)
@@ -4,899 +4,21 @@
 !_TAG_PROGRAM_NAME	Exuberant Ctags	//
 !_TAG_PROGRAM_URL	http://ctags.sourceforge.net	/official site/
 !_TAG_PROGRAM_VERSION	5.6	//
-AAvsABandBB	VARlib/supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	e	enum:analysis_type_enum
-AAvsABandBB	supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	e	enum:analysis_type_enum
-AAvsABvsBB	VARlib/supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	e	enum:analysis_type_enum
-AAvsABvsBB	supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	e	enum:analysis_type_enum
-ABvsAAandBB	VARlib/supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	e	enum:analysis_type_enum
-ABvsAAandBB	supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	e	enum:analysis_type_enum
-AbstractMatrix	AbstractMatrix.h	/^class AbstractMatrix {$/;"	c
-AbstractMatrix	fvlib/AbstractMatrix.h	/^class AbstractMatrix {$/;"	c
-AbstractMatrixRFinalizer	AbstractMatrix_R.cpp	/^	static void AbstractMatrixRFinalizer(SEXP x) {$/;"	f	file:
-AbstractMatrixRFinalizer	DAlib/AbstractMatrix_R.cpp	/^	static void AbstractMatrixRFinalizer(SEXP x) {$/;"	f	file:
-BBvsAAandAB	VARlib/supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	e	enum:analysis_type_enum
-BBvsAAandAB	supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	e	enum:analysis_type_enum
-CHAR_NAN	CastUtils.cpp	/^char CHAR_NAN;$/;"	v
-CHAR_NAN	fvlib/CastUtils.cpp	/^char CHAR_NAN;$/;"	v
-COPY_AND_COMPARE	CastUtils.h	38;"	d
-COPY_AND_COMPARE	fvlib/CastUtils.h	38;"	d
-COUNTS	VARlib/supplementary_functions.h	/^	int COUNTS[GENO_TYPES_NUM];$/;"	m	struct:snp_var_data
-COUNTS	supplementary_functions.h	/^	int COUNTS[GENO_TYPES_NUM];$/;"	m	struct:snp_var_data
-CPP_dqrls	AbstractMatrix_R.cpp	/^	void CPP_dqrls($/;"	f
-CPP_dqrls	DAlib/AbstractMatrix_R.cpp	/^	void CPP_dqrls($/;"	f
-DB_CREATE	FileVector.h	26;"	d
-DB_CREATE	fvlib/FileVector.h	26;"	d
-DB_EXCL	FileVector.h	27;"	d
-DB_EXCL	fvlib/FileVector.h	27;"	d
-DB_RDONLY	FileVector.h	28;"	d
-DB_RDONLY	fvlib/FileVector.h	28;"	d
-DEBUG_LEVEL	Logger.h	17;"	d
-DEBUG_LEVEL	fvlib/Logger.h	17;"	d
-DOUBLE	const.h	13;"	d
-DOUBLE	fvlib/const.h	13;"	d
-ERROR_LEVEL	Logger.h	18;"	d
-ERROR_LEVEL	fvlib/Logger.h	18;"	d
-ErrorExit	Logger.h	/^class ErrorExit {};$/;"	c
-ErrorExit	fvlib/Logger.h	/^class ErrorExit {};$/;"	c
-FILEVECTOR_DATA_FILE_SUFFIX	FileVector.h	/^const string FILEVECTOR_DATA_FILE_SUFFIX = ".fvd";$/;"	v
-FILEVECTOR_DATA_FILE_SUFFIX	fvlib/FileVector.h	/^const string FILEVECTOR_DATA_FILE_SUFFIX = ".fvd";$/;"	v
-FILEVECTOR_INDEX_FILE_SUFFIX	FileVector.h	/^const string FILEVECTOR_INDEX_FILE_SUFFIX = ".fvi";$/;"	v
-FILEVECTOR_INDEX_FILE_SUFFIX	fvlib/FileVector.h	/^const string FILEVECTOR_INDEX_FILE_SUFFIX = ".fvi";$/;"	v
-FLOAT	const.h	12;"	d
-FLOAT	fvlib/const.h	12;"	d
-FV_AUTHORS	frversion.h	6;"	d
-FV_AUTHORS	fvlib/frversion.h	6;"	d
-FV_RELEASEDATE	frversion.h	5;"	d
-FV_RELEASEDATE	fvlib/frversion.h	5;"	d
-FV_VERSION	frversion.h	4;"	d
-FV_VERSION	fvlib/frversion.h	4;"	d
-FileHeader	frutil.h	/^    FileHeader () {$/;"	f	class:FileHeader
-FileHeader	frutil.h	/^class FileHeader$/;"	c
-FileHeader	fvlib/frutil.h	/^    FileHeader () {$/;"	f	class:FileHeader
-FileHeader	fvlib/frutil.h	/^class FileHeader$/;"	c
-FileVector	FileVector.h	/^	FileVector(char *iFilename, unsigned long cachesizeMb) : filename(string(iFilename)){$/;"	f	class:FileVector
-FileVector	FileVector.h	/^	FileVector(char *iFilename, unsigned long cachesizeMb, bool iReadOnly) : filename(string(iFilename)), readOnly(iReadOnly) {$/;"	f	class:FileVector
-FileVector	FileVector.h	/^	FileVector(string iFilename, unsigned long cachesizeMb) : filename (iFilename) {$/;"	f	class:FileVector
-FileVector	FileVector.h	/^	FileVector(string iFilename, unsigned long cachesizeMb, bool iReadOnly) : filename(iFilename), readOnly(iReadOnly) {$/;"	f	class:FileVector
-FileVector	FileVector.h	/^class FileVector: public AbstractMatrix {$/;"	c
-FileVector	fvlib/FileVector.h	/^	FileVector(char *iFilename, unsigned long cachesizeMb) : filename(string(iFilename)){$/;"	f	class:FileVector
-FileVector	fvlib/FileVector.h	/^	FileVector(char *iFilename, unsigned long cachesizeMb, bool iReadOnly) : filename(string(iFilename)), readOnly(iReadOnly) {$/;"	f	class:FileVector
-FileVector	fvlib/FileVector.h	/^	FileVector(string iFilename, unsigned long cachesizeMb) : filename (iFilename) {$/;"	f	class:FileVector
-FileVector	fvlib/FileVector.h	/^	FileVector(string iFilename, unsigned long cachesizeMb, bool iReadOnly) : filename(iFilename), readOnly(iReadOnly) {$/;"	f	class:FileVector
-FileVector	fvlib/FileVector.h	/^class FileVector: public AbstractMatrix {$/;"	c
-FilteredAndAbstractRFinalizer	DAlib/FilteredMatrix_R.cpp	/^    static void FilteredAndAbstractRFinalizer(SEXP x) {$/;"	f	file:
-FilteredAndAbstractRFinalizer	FilteredMatrix_R.cpp	/^    static void FilteredAndAbstractRFinalizer(SEXP x) {$/;"	f	file:
-FilteredMatrix	FilteredMatrix.h	/^    FilteredMatrix(AbstractMatrix &matrix) : nestedMatrix(&matrix) {$/;"	f	class:FilteredMatrix
-FilteredMatrix	FilteredMatrix.h	/^class FilteredMatrix : public AbstractMatrix {$/;"	c
-FilteredMatrix	fvlib/FilteredMatrix.h	/^    FilteredMatrix(AbstractMatrix &matrix) : nestedMatrix(&matrix) {$/;"	f	class:FilteredMatrix
-FilteredMatrix	fvlib/FilteredMatrix.h	/^class FilteredMatrix : public AbstractMatrix {$/;"	c
-FilteredMatrixRFinalizer	DAlib/FilteredMatrix_R.cpp	/^    static void FilteredMatrixRFinalizer(SEXP x) {$/;"	f	file:
-FilteredMatrixRFinalizer	FilteredMatrix_R.cpp	/^    static void FilteredMatrixRFinalizer(SEXP x) {$/;"	f	file:
-FixedChar	frutil.h	/^    FixedChar(){$/;"	f	class:FixedChar
-FixedChar	frutil.h	/^    FixedChar(string s){$/;"	f	class:FixedChar
-FixedChar	frutil.h	/^class FixedChar$/;"	c
-FixedChar	fvlib/frutil.h	/^    FixedChar(){$/;"	f	class:FixedChar
-FixedChar	fvlib/frutil.h	/^    FixedChar(string s){$/;"	f	class:FixedChar
-FixedChar	fvlib/frutil.h	/^class FixedChar$/;"	c
-GENO	VARlib/supplementary_functions.h	/^	std::string GENO[GENO_TYPES_NUM];$/;"	m	struct:snp_var_data
-GENO	supplementary_functions.h	/^	std::string GENO[GENO_TYPES_NUM];$/;"	m	struct:snp_var_data
-GENO_TYPES_NUM	VARlib/supplementary_functions.h	39;"	d
-GENO_TYPES_NUM	VARlib/var_meta_gwaa_C.cpp	79;"	d	file:
-GENO_TYPES_NUM	supplementary_functions.h	39;"	d
-GENO_TYPES_NUM	var_meta_gwaa_C.cpp	79;"	d	file:
-INT	const.h	11;"	d
-INT	fvlib/const.h	11;"	d
-INT_NAN	CastUtils.cpp	/^int INT_NAN;$/;"	v
-INT_NAN	fvlib/CastUtils.cpp	/^int INT_NAN;$/;"	v
-Logger	Logger.h	/^    Logger(int iErrorLevel) : errorLevel(iErrorLevel), enabled(true) {}$/;"	f	class:Logger
-Logger	Logger.h	/^    Logger(int iErrorLevel, bool iEnabled) : errorLevel(iErrorLevel), enabled(iEnabled) {}$/;"	f	class:Logger
-Logger	Logger.h	/^class Logger {$/;"	c
-Logger	fvlib/Logger.h	/^    Logger(int iErrorLevel) : errorLevel(iErrorLevel), enabled(true) {}$/;"	f	class:Logger
-Logger	fvlib/Logger.h	/^    Logger(int iErrorLevel, bool iEnabled) : errorLevel(iErrorLevel), enabled(iEnabled) {}$/;"	f	class:Logger
-Logger	fvlib/Logger.h	/^class Logger {$/;"	c
-MEAN	VARlib/supplementary_functions.h	/^	VARIABLE_TYPE MEAN[GENO_TYPES_NUM];$/;"	m	struct:snp_var_data
-MEAN	supplementary_functions.h	/^	VARIABLE_TYPE MEAN[GENO_TYPES_NUM];$/;"	m	struct:snp_var_data
-MESSAGE_LEVEL	Logger.h	16;"	d
-MESSAGE_LEVEL	fvlib/Logger.h	16;"	d
-MethodConvStruct	ITERlib/iterator_functions.h	/^	struct MethodConvStruct {$/;"	s
 MethodConvStruct	iterator_functions.h	/^	struct MethodConvStruct {$/;"	s
-NAMELENGTH	const.h	18;"	d
-NAMELENGTH	fvlib/const.h	18;"	d
-NA_value	VARlib/constants.h	/^const VARIABLE_TYPE NA_value = -999.999; \/\/ which numeric value is ised as NA$/;"	v
-NA_value	constants.h	/^const VARIABLE_TYPE NA_value = -999.999; \/\/ which numeric value is ised as NA$/;"	v
-NA_value_int	VARlib/constants.h	/^const int NA_value_int = -999; \/\/ which numeric value is ised as NA$/;"	v
-NA_value_int	constants.h	/^const int NA_value_int = -999; \/\/ which numeric value is ised as NA$/;"	v
-PART_SIZE	frutil.cpp	162;"	d	file:
-PART_SIZE	fvlib/frutil.cpp	162;"	d	file:
-REPORT_EVERY	convert_util.cpp	15;"	d	file:
-REPORT_EVERY	fvlib/convert_util.cpp	15;"	d	file:
-RESERVEDSPACE	const.h	5;"	d
-RESERVEDSPACE	fvlib/const.h	5;"	d
-RealHandlerWrapper	RealHandlerWrapper.h	/^    RealHandlerWrapper(): useCount(0) {}$/;"	f	class:RealHandlerWrapper
-RealHandlerWrapper	RealHandlerWrapper.h	/^class RealHandlerWrapper {$/;"	c
-RealHandlerWrapper	fvlib/RealHandlerWrapper.h	/^    RealHandlerWrapper(): useCount(0) {}$/;"	f	class:RealHandlerWrapper
-RealHandlerWrapper	fvlib/RealHandlerWrapper.h	/^class RealHandlerWrapper {$/;"	c
-ReusableFileHandle	ReusableFileHandle.h	/^    ReusableFileHandle() : isOk(false),curPos(0),$/;"	f	class:ReusableFileHandle
-ReusableFileHandle	ReusableFileHandle.h	/^    ReusableFileHandle(RealHandlerWrapper *iRealHandlerWrapper, bool iIsOk,$/;"	f	class:ReusableFileHandle
-ReusableFileHandle	ReusableFileHandle.h	/^    ReusableFileHandle(const ReusableFileHandle &rfh) :$/;"	f	class:ReusableFileHandle
-ReusableFileHandle	ReusableFileHandle.h	/^class ReusableFileHandle {$/;"	c
-ReusableFileHandle	fvlib/ReusableFileHandle.h	/^    ReusableFileHandle() : isOk(false),curPos(0),$/;"	f	class:ReusableFileHandle
-ReusableFileHandle	fvlib/ReusableFileHandle.h	/^    ReusableFileHandle(RealHandlerWrapper *iRealHandlerWrapper, bool iIsOk,$/;"	f	class:ReusableFileHandle
-ReusableFileHandle	fvlib/ReusableFileHandle.h	/^    ReusableFileHandle(const ReusableFileHandle &rfh) :$/;"	f	class:ReusableFileHandle
-ReusableFileHandle	fvlib/ReusableFileHandle.h	/^class ReusableFileHandle {$/;"	c
-SD	VARlib/supplementary_functions.h	/^	VARIABLE_TYPE SD[GENO_TYPES_NUM];$/;"	m	struct:snp_var_data
-SD	supplementary_functions.h	/^	VARIABLE_TYPE SD[GENO_TYPES_NUM];$/;"	m	struct:snp_var_data
-SHORT_INT	const.h	9;"	d
-SHORT_INT	fvlib/const.h	9;"	d
-SHORT_INT_NAN	CastUtils.cpp	/^short int SHORT_INT_NAN;$/;"	v
-SHORT_INT_NAN	fvlib/CastUtils.cpp	/^short int SHORT_INT_NAN;$/;"	v
-SIGNED_CHAR	const.h	14;"	d
-SIGNED_CHAR	fvlib/const.h	14;"	d
-SMV_CONSTANTS_H	VARlib/constants.h	3;"	d
-SMV_CONSTANTS_H	constants.h	3;"	d
-SMV_DOMETA_C_H	VARlib/inverse_variance_metaanalysis.h	19;"	d
-SMV_DOMETA_C_H	inverse_variance_metaanalysis.h	19;"	d
-SMV_GTPS_CONTAINER_H	VARlib/gtps_container.h	19;"	d
-SMV_GTPS_CONTAINER_H	gtps_container.h	19;"	d
-SMV_SUPPLEMENTARY_FUNCTIONS_H	VARlib/supplementary_functions.h	19;"	d
-SMV_SUPPLEMENTARY_FUNCTIONS_H	supplementary_functions.h	19;"	d
-SMV_VAR_HOMOGENEITY_TESTS_H	VARlib/var_homogeneity_tests.h	2;"	d
-SMV_VAR_HOMOGENEITY_TESTS_H	var_homogeneity_tests.h	2;"	d
-Snp_store_type	VARlib/supplementary_functions.h	/^typedef std::map<std::string, snp_var_data*> Snp_store_type; \/\/use std::string as key because for some reason it doesn't want to work with const char*$/;"	t
-Snp_store_type	supplementary_functions.h	/^typedef std::map<std::string, snp_var_data*> Snp_store_type; \/\/use std::string as key because for some reason it doesn't want to work with const char*$/;"	t
-Transposer	Transposer.h	/^Transposer(){ square_size =10000;};$/;"	f	class:Transposer
-Transposer	Transposer.h	/^Transposer(int opt_square_size){ square_size=opt_square_size;};$/;"	f	class:Transposer
-Transposer	Transposer.h	/^class Transposer$/;"	c
-Transposer	fvlib/Transposer.h	/^Transposer(){ square_size =10000;};$/;"	f	class:Transposer
-Transposer	fvlib/Transposer.h	/^Transposer(int opt_square_size){ square_size=opt_square_size;};$/;"	f	class:Transposer
-Transposer	fvlib/Transposer.h	/^class Transposer$/;"	c
-UNSIGNED_CHAR	const.h	15;"	d
-UNSIGNED_CHAR	fvlib/const.h	15;"	d
-UNSIGNED_CHAR_NAN	CastUtils.cpp	/^unsigned char UNSIGNED_CHAR_NAN;$/;"	v
-UNSIGNED_CHAR_NAN	fvlib/CastUtils.cpp	/^unsigned char UNSIGNED_CHAR_NAN;$/;"	v
-UNSIGNED_INT	const.h	10;"	d
-UNSIGNED_INT	fvlib/const.h	10;"	d
-UNSIGNED_INT_NAN	CastUtils.cpp	/^unsigned int UNSIGNED_INT_NAN;$/;"	v
-UNSIGNED_INT_NAN	fvlib/CastUtils.cpp	/^unsigned int UNSIGNED_INT_NAN;$/;"	v
-UNSIGNED_SHORT_INT	const.h	8;"	d
-UNSIGNED_SHORT_INT	fvlib/const.h	8;"	d
-UNSIGNED_SHORT_INT_NAN	CastUtils.cpp	/^unsigned short int UNSIGNED_SHORT_INT_NAN;$/;"	v
-UNSIGNED_SHORT_INT_NAN	fvlib/CastUtils.cpp	/^unsigned short int UNSIGNED_SHORT_INT_NAN;$/;"	v
-VARIABLE_TYPE	VARlib/supplementary_functions.h	38;"	d
-VARIABLE_TYPE	VARlib/var_meta_gwaa_C.cpp	78;"	d	file:
-VARIABLE_TYPE	supplementary_functions.h	38;"	d
-VARIABLE_TYPE	var_meta_gwaa_C.cpp	78;"	d	file:
-WRITE_SPEED_PROPORTION	AbstractMatrix.h	12;"	d
-WRITE_SPEED_PROPORTION	fvlib/AbstractMatrix.h	12;"	d
-Z	VARlib/supplementary_functions.h	/^	double Z; \/\/homogeneity test$/;"	m	struct:snp_var_data
-Z	supplementary_functions.h	/^	double Z; \/\/homogeneity test$/;"	m	struct:snp_var_data
-Z_2df	VARlib/supplementary_functions.h	/^	double Z_2df; \/\/homogeneity test 2df only$/;"	m	struct:snp_var_data
-Z_2df	supplementary_functions.h	/^	double Z_2df; \/\/homogeneity test 2df only$/;"	m	struct:snp_var_data
-_	VARlib/inverse_variance_metaanalysis.cpp	30;"	d	file:
-_	VARlib/inverse_variance_metaanalysis.cpp	33;"	d	file:
-_	VARlib/var_homogeneity_test_C.cpp	34;"	d	file:
-_	VARlib/var_homogeneity_test_C.cpp	37;"	d	file:
-_	VARlib/var_meta_gwaa_C.cpp	71;"	d	file:
-_	VARlib/var_meta_gwaa_C.cpp	74;"	d	file:
-_	inverse_variance_metaanalysis.cpp	30;"	d	file:
-_	inverse_variance_metaanalysis.cpp	33;"	d	file:
-_	var_homogeneity_test_C.cpp	34;"	d	file:
-_	var_homogeneity_test_C.cpp	37;"	d	file:
-_	var_meta_gwaa_C.cpp	71;"	d	file:
-_	var_meta_gwaa_C.cpp	74;"	d	file:
-__AbstractMatrix__	AbstractMatrix.h	2;"	d
-__AbstractMatrix__	fvlib/AbstractMatrix.h	2;"	d
-__CONST_H__	const.h	2;"	d
-__CONST_H__	fvlib/const.h	2;"	d
-__CONVERT_UTIL__	convert_util.h	2;"	d
-__CONVERT_UTIL__	fvlib/convert_util.h	2;"	d
-__CastUtils__	CastUtils.h	2;"	d
-__CastUtils__	fvlib/CastUtils.h	2;"	d
-__DAUTIL_H__	DAlib/dautil.h	2;"	d
-__DAUTIL_H__	dautil.h	2;"	d
-__FRUTIL__	frutil.h	2;"	d
-__FRUTIL__	fvlib/frutil.h	2;"	d
-__FRVECTOR__	FileVector.h	2;"	d
-__FRVECTOR__	fvlib/FileVector.h	2;"	d
-__FRVERSION__	frversion.h	2;"	d
-__FRVERSION__	fvlib/frversion.h	2;"	d
-__FilteredMatrix__	FilteredMatrix.h	2;"	d
-__FilteredMatrix__	fvlib/FilteredMatrix.h	2;"	d
-__ITERATOR_FUNCTIONS_H__	ITERlib/iterator_functions.h	2;"	d
 __ITERATOR_FUNCTIONS_H__	iterator_functions.h	2;"	d
-__LOGGER__	Logger.h	2;"	d
-__LOGGER__	fvlib/Logger.h	2;"	d
-__RSTAFF_H__	DAlib/Rstaff.h	2;"	d
-__RSTAFF_H__	Rstaff.h	2;"	d
-__RealHandlerWrapper__	RealHandlerWrapper.h	2;"	d
-__RealHandlerWrapper__	fvlib/RealHandlerWrapper.h	2;"	d
-__ReusableFileHandle__	ReusableFileHandle.h	2;"	d
-__ReusableFileHandle__	fvlib/ReusableFileHandle.h	2;"	d
-___TRANSPOSE___	Transposer.h	2;"	d
-___TRANSPOSE___	fvlib/Transposer.h	2;"	d
-addVariable	FileVector.cpp	/^void FileVector::addVariable(void *invec, string varName) {$/;"	f	class:FileVector
-addVariable	FilteredMatrix.cpp	/^void FilteredMatrix::addVariable(void * invec, string varname) {$/;"	f	class:FilteredMatrix
-addVariable	fvlib/FileVector.cpp	/^void FileVector::addVariable(void *invec, string varName) {$/;"	f	class:FileVector
-addVariable	fvlib/FilteredMatrix.cpp	/^void FilteredMatrix::addVariable(void * invec, string varname) {$/;"	f	class:FilteredMatrix
-addVariableAs	AbstractMatrix.h	/^    void addVariableAs(DT * outvec, string varname)$/;"	f	class:AbstractMatrix
-addVariableAs	fvlib/AbstractMatrix.h	/^    void addVariableAs(DT * outvec, string varname)$/;"	f	class:AbstractMatrix
-all	Logger.cpp	6;"	d	file:
-all	fvlib/Logger.cpp	6;"	d	file:
-analysis_type_enum	VARlib/supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	g
-analysis_type_enum	supplementary_functions.h	/^enum analysis_type_enum {AAvsABvsBB, AAvsABandBB, ABvsAAandBB, BBvsAAandAB};$/;"	g
-assignDoubleMatrix	AbstractMatrix_R.cpp	/^	SEXP assignDoubleMatrix(SEXP ptr, SEXP obsIndexes, SEXP varIndexes, SEXP values, SEXP direction){$/;"	f
-assignDoubleMatrix	DAlib/AbstractMatrix_R.cpp	/^	SEXP assignDoubleMatrix(SEXP ptr, SEXP obsIndexes, SEXP varIndexes, SEXP values, SEXP direction){$/;"	f
-bartlett	VARlib/supplementary_functions.h	/^enum testname_enum {bartlett, levene, likelihood, kolmogorov_smirnov};$/;"	e	enum:testname_enum
-bartlett	supplementary_functions.h	/^enum testname_enum {bartlett, levene, likelihood, kolmogorov_smirnov};$/;"	e	enum:testname_enum
-bartlett_test	VARlib/var_homogeneity_tests.cpp	/^chisq_df bartlett_test(snp_var_data * snp)$/;"	f
-bartlett_test	VARlib/var_homogeneity_tests.cpp	/^chisq_df bartlett_test(std::list<my_small_vector> * samples)$/;"	f
-bartlett_test	var_homogeneity_tests.cpp	/^chisq_df bartlett_test(snp_var_data * snp)$/;"	f
-bartlett_test	var_homogeneity_tests.cpp	/^chisq_df bartlett_test(std::list<my_small_vector> * samples)$/;"	f
-bitsPerRecord	frutil.h	/^	unsigned int bitsPerRecord;$/;"	m	class:FileHeader
-bitsPerRecord	fvlib/frutil.h	/^	unsigned int bitsPerRecord;$/;"	m	class:FileHeader
-blockWriteOrRead	RealHandlerWrapper.cpp	/^void RealHandlerWrapper::blockWriteOrRead(unsigned long iLength, char* data, bool writeAction) {$/;"	f	class:RealHandlerWrapper
-blockWriteOrRead	ReusableFileHandle.cpp	/^void ReusableFileHandle::blockWriteOrRead(unsigned long length, char* data, bool writeAction){$/;"	f	class:ReusableFileHandle
-blockWriteOrRead	frutil.cpp	/^void blockWriteOrRead(fstream& file, unsigned long length, char* data, bool writeAction){$/;"	f
-blockWriteOrRead	fvlib/RealHandlerWrapper.cpp	/^void RealHandlerWrapper::blockWriteOrRead(unsigned long iLength, char* data, bool writeAction) {$/;"	f	class:RealHandlerWrapper
-blockWriteOrRead	fvlib/ReusableFileHandle.cpp	/^void ReusableFileHandle::blockWriteOrRead(unsigned long length, char* data, bool writeAction){$/;"	f	class:ReusableFileHandle
-blockWriteOrRead	fvlib/frutil.cpp	/^void blockWriteOrRead(fstream& file, unsigned long length, char* data, bool writeAction){$/;"	f
-break_trait_up_into_groups	VARlib/supplementary_functions.cpp	/^void break_trait_up_into_groups(std::list<my_small_vector> *trait_groups, int* snp, double *trait, unsigned* nids, analysis_type_enum analys_type, int * is_trait_na)$/;"	f
-break_trait_up_into_groups	supplementary_functions.cpp	/^void break_trait_up_into_groups(std::list<my_small_vector> *trait_groups, int* snp, double *trait, unsigned* nids, analysis_type_enum analys_type, int * is_trait_na)$/;"	f
-bufToString	CastUtils.cpp	/^string bufToString(short int dataType, char *data, string nanString){$/;"	f
-bufToString	fvlib/CastUtils.cpp	/^string bufToString(short int dataType, char *data, string nanString){$/;"	f
-bytesPerRecord	frutil.h	/^	unsigned int bytesPerRecord;$/;"	m	class:FileHeader
-bytesPerRecord	fvlib/frutil.h	/^	unsigned int bytesPerRecord;$/;"	m	class:FileHeader
-cacheAllNames	FileVector.cpp	/^void FileVector::cacheAllNames(bool doCache) {$/;"	f	class:FileVector
-cacheAllNames	FilteredMatrix.cpp	/^void FilteredMatrix::cacheAllNames(bool doCache) {$/;"	f	class:FilteredMatrix
-cacheAllNames	fvlib/FileVector.cpp	/^void FileVector::cacheAllNames(bool doCache) {$/;"	f	class:FileVector
-cacheAllNames	fvlib/FilteredMatrix.cpp	/^void FilteredMatrix::cacheAllNames(bool doCache) {$/;"	f	class:FilteredMatrix
-cacheBegin	FileVector.h	/^	unsigned long cacheBegin;$/;"	m	class:FileVector
-cacheBegin	fvlib/FileVector.h	/^	unsigned long cacheBegin;$/;"	m	class:FileVector
-cacheBuffer	FileVector.h	/^	char * cacheBuffer;$/;"	m	class:FileVector
-cacheBuffer	fvlib/FileVector.h	/^	char * cacheBuffer;$/;"	m	class:FileVector
-cacheEnd	FileVector.h	/^	unsigned long cacheEnd;$/;"	m	class:FileVector
-cacheEnd	fvlib/FileVector.h	/^	unsigned long cacheEnd;$/;"	m	class:FileVector
-cache_size_Mb	FileVector.h	/^	unsigned long cache_size_Mb;$/;"	m	class:FileVector
-cache_size_Mb	fvlib/FileVector.h	/^	unsigned long cache_size_Mb;$/;"	m	class:FileVector
-cache_size_bytes	FileVector.h	/^	unsigned long cache_size_bytes;$/;"	m	class:FileVector
-cache_size_bytes	fvlib/FileVector.h	/^	unsigned long cache_size_bytes;$/;"	m	class:FileVector
-cache_size_nvars	FileVector.h	/^	unsigned long cache_size_nvars;$/;"	m	class:FileVector
-cache_size_nvars	fvlib/FileVector.h	/^	unsigned long cache_size_nvars;$/;"	m	class:FileVector
-calcCachePos	FileVector.cpp	/^void FileVector::calcCachePos(unsigned long newCenterPos, unsigned long &cacheBeginRef, unsigned long &cacheEndRef){$/;"	f	class:FileVector
-calcCachePos	fvlib/FileVector.cpp	/^void FileVector::calcCachePos(unsigned long newCenterPos, unsigned long &cacheBeginRef, unsigned long &cacheEndRef){$/;"	f	class:FileVector
-calcDataSize	frutil.cpp	/^unsigned short calcDataSize(unsigned short int type){$/;"	f
-calcDataSize	fvlib/frutil.cpp	/^unsigned short calcDataSize(unsigned short int type){$/;"	f
-calcNumLines	convert_util.cpp	/^unsigned long calcNumLines(string fileName){$/;"	f
-calcNumLines	fvlib/convert_util.cpp	/^unsigned long calcNumLines(string fileName){$/;"	f
-calcNumWordsInFirstLine	convert_util.cpp	/^unsigned long calcNumWordsInFirstLine(string fileName){$/;"	f
-calcNumWordsInFirstLine	fvlib/convert_util.cpp	/^unsigned long calcNumWordsInFirstLine(string fileName){$/;"	f
-castToAbstractMatrix	FileVector.cpp	/^AbstractMatrix* FileVector::castToAbstractMatrix(){$/;"	f	class:FileVector
-castToAbstractMatrix	FilteredMatrix.cpp	/^AbstractMatrix* FilteredMatrix::castToAbstractMatrix(){$/;"	f	class:FilteredMatrix
-castToAbstractMatrix	fvlib/FileVector.cpp	/^AbstractMatrix* FileVector::castToAbstractMatrix(){$/;"	f	class:FileVector
-castToAbstractMatrix	fvlib/FilteredMatrix.cpp	/^AbstractMatrix* FilteredMatrix::castToAbstractMatrix(){$/;"	f	class:FilteredMatrix
-check	VARlib/supplementary_functions.cpp	/^bool check(snp_var_data* snp2, snp_var_data* snp1, std::ofstream & warnings_file)$/;"	f
-check	supplementary_functions.cpp	/^bool check(snp_var_data* snp2, snp_var_data* snp1, std::ofstream & warnings_file)$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(char i){return checkNan(&i, SIGNED_CHAR);}$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(double i){return checkNan(&i, DOUBLE);}$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(float i){return checkNan(&i, FLOAT);}$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(int i){return checkNan(&i, INT);}$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(short int i){return checkNan(&i, SHORT_INT);}$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(unsigned char i){return checkNan(&i, UNSIGNED_CHAR);}$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(unsigned int i){return checkNan(&i, UNSIGNED_INT);}$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(unsigned short int i){return checkNan(&i, UNSIGNED_SHORT_INT);}$/;"	f
-checkNan	CastUtils.cpp	/^bool checkNan(void *data, int dataType){$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(char i){return checkNan(&i, SIGNED_CHAR);}$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(double i){return checkNan(&i, DOUBLE);}$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(float i){return checkNan(&i, FLOAT);}$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(int i){return checkNan(&i, INT);}$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(short int i){return checkNan(&i, SHORT_INT);}$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(unsigned char i){return checkNan(&i, UNSIGNED_CHAR);}$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(unsigned int i){return checkNan(&i, UNSIGNED_INT);}$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(unsigned short int i){return checkNan(&i, UNSIGNED_SHORT_INT);}$/;"	f
-checkNan	fvlib/CastUtils.cpp	/^bool checkNan(void *data, int dataType){$/;"	f
-checkNumBits	AbstractMatrix_R.cpp	/^	SEXP checkNumBits(){$/;"	f
-checkNumBits	DAlib/AbstractMatrix_R.cpp	/^	SEXP checkNumBits(){$/;"	f
-checkOpenForWriting	AbstractMatrix.cpp	/^void AbstractMatrix::checkOpenForWriting(const string fileName){$/;"	f	class:AbstractMatrix
-checkOpenForWriting	fvlib/AbstractMatrix.cpp	/^void AbstractMatrix::checkOpenForWriting(const string fileName){$/;"	f	class:AbstractMatrix
-checkPointer	AbstractMatrix_R.cpp	/^	void checkPointer(SEXP s) {$/;"	f
-checkPointer	DAlib/AbstractMatrix_R.cpp	/^	void checkPointer(SEXP s) {$/;"	f
-check_files_format	VARlib/supplementary_functions.cpp	/^bool check_files_format(const char** filenames, unsigned file_amount, unsigned skip_first_lines_amount, char delim)$/;"	f
-check_files_format	supplementary_functions.cpp	/^bool check_files_format(const char** filenames, unsigned file_amount, unsigned skip_first_lines_amount, char delim)$/;"	f
-chisq	VARlib/var_homogeneity_tests.h	/^		double chisq;$/;"	m	class:chisq_df
-chisq	var_homogeneity_tests.h	/^		double chisq;$/;"	m	class:chisq_df
-chisq_df	VARlib/var_homogeneity_tests.h	/^		chisq_df() $/;"	f	class:chisq_df
-chisq_df	VARlib/var_homogeneity_tests.h	/^		chisq_df(double chisq_, int df_) $/;"	f	class:chisq_df
-chisq_df	VARlib/var_homogeneity_tests.h	/^class chisq_df$/;"	c
-chisq_df	var_homogeneity_tests.h	/^		chisq_df() $/;"	f	class:chisq_df
-chisq_df	var_homogeneity_tests.h	/^		chisq_df(double chisq_, int df_) $/;"	f	class:chisq_df
-chisq_df	var_homogeneity_tests.h	/^class chisq_df$/;"	c
-chromosome	VARlib/supplementary_functions.h	/^	int chromosome;$/;"	m	struct:snp_var_data
-chromosome	supplementary_functions.h	/^	int chromosome;$/;"	m	struct:snp_var_data
-chromosome_column_name	VARlib/supplementary_functions.h	/^const std::string chromosome_column_name = "CHR";$/;"	v
-chromosome_column_name	supplementary_functions.h	/^const std::string chromosome_column_name = "CHR";$/;"	v
-close	RealHandlerWrapper.cpp	/^void RealHandlerWrapper::close(){$/;"	f	class:RealHandlerWrapper
-close	ReusableFileHandle.cpp	/^void ReusableFileHandle::close() {$/;"	f	class:ReusableFileHandle
-close	fvlib/RealHandlerWrapper.cpp	/^void RealHandlerWrapper::close(){$/;"	f	class:RealHandlerWrapper
-close	fvlib/ReusableFileHandle.cpp	/^void ReusableFileHandle::close() {$/;"	f	class:ReusableFileHandle
-closeForWriting	AbstractMatrix.cpp	/^void AbstractMatrix::closeForWriting(const string fileName){$/;"	f	class:AbstractMatrix
-closeForWriting	fvlib/AbstractMatrix.cpp	/^void AbstractMatrix::closeForWriting(const string fileName){$/;"	f	class:AbstractMatrix
-coding_array	VARlib/gtps_container.h	/^		char * coding_array;$/;"	m	class:gtps_container
-coding_array	gtps_container.h	/^		char * coding_array;$/;"	m	class:gtps_container
-compare_doubles	VARlib/supplementary_functions.cpp	/^int compare_doubles(const void *a, const void *b)$/;"	f
-compare_doubles	supplementary_functions.cpp	/^int compare_doubles(const void *a, const void *b)$/;"	f
-copyVariable	FileVector.cpp	/^void FileVector::copyVariable(char* to, char* from, int n, unsigned long * indexes ) {$/;"	f	class:FileVector
-copyVariable	fvlib/FileVector.cpp	/^void FileVector::copyVariable(char* to, char* from, int n, unsigned long * indexes ) {$/;"	f	class:FileVector
-copy_data	Transposer.cpp	/^void Transposer::copy_data(string src_data_file_name,string dest_data_file_name, unsigned long src_nvars,$/;"	f	class:Transposer
-copy_data	fvlib/Transposer.cpp	/^void Transposer::copy_data(string src_data_file_name,string dest_data_file_name, unsigned long src_nvars,$/;"	f	class:Transposer
[TRUNCATED]

To get the complete diff run:
    svnlook diff /svnroot/genabel -r 638


More information about the Genabel-commits mailing list