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

[noreply at r-forge.r-project.org]

Author: maksim
Date: 2011-02-08 11:55:30 +0100 (Tue, 08 Feb 2011)
New Revision: 639

Modified:
   pkg/VariABEL/R/var.test.gwaa.R
   pkg/VariABEL/src/var_homogeneity_test_C.cpp
Log:
ready

Modified: pkg/VariABEL/R/var.test.gwaa.R
===================================================================
--- pkg/VariABEL/R/var.test.gwaa.R	2011-02-08 09:53:26 UTC (rev 638)
+++ pkg/VariABEL/R/var.test.gwaa.R	2011-02-08 10:55:30 UTC (rev 639)
@@ -107,7 +107,7 @@
 		{
 		print("reading genotype's info...")
 		genodata_info_df <- read.table(genodata_info, strings=F, header=T)
-		print(c(dim(genodata_info_df)[1], gtNcol))
+#		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],
@@ -167,20 +167,20 @@
 
 print("Variance analysis done.")
 
-print(genodata_info_df)
-print(results_C)
+#print(genodata_info_df)
+#print(results_C)
 
 results_df <- data.frame(genodata_info_df, results_C)
-print(results_df)
+#print(results_df)
 
 cov_names <- colnames(design_matrix_df)
 
-print(cov_names)
+#print(cov_names)
 
 output_column_names <- c(colnames(genodata_info_df), "chisq", "df", "Intercept_effect", "Intercept_sd")
 
 
-print(output_column_names)
+#print(output_column_names)
 
 for(i in 2:p)
 	{

Modified: pkg/VariABEL/src/var_homogeneity_test_C.cpp
===================================================================
--- pkg/VariABEL/src/var_homogeneity_test_C.cpp	2011-02-08 09:53:26 UTC (rev 638)
+++ pkg/VariABEL/src/var_homogeneity_test_C.cpp	2011-02-08 10:55:30 UTC (rev 639)
@@ -125,8 +125,11 @@
 								/*auxiliary variables:*/ qty, jpvt, qraux, work, v, x_for_ch2inv);
 
 	//___________________________________________
-	
 
+
+//	for(int i=0 ; i<p ; i++) {std::cout<<"betas["<<i<<"]="<<betas[i]<<", se["<<i<<"]="<<se[i]<<"\n";}
+
+
 	*chi2 = (betas[p-1]/se[p-1])*(betas[p-1]/se[p-1]);
 	*df = 1;