shalabhgandhi@usa.net
2000-Nov-02 05:41 UTC
[Rd] Pl. provide and Input for var.test (PR#716)
I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of mine for which I am mentioning 2 cases CaseA and CaseB. CaseA is working but CaseB is not. (CaseA) Following are 4 commands I gave: R library(ctest) x<-rnorm(50) y<-rnorm(30) var.test(x,y) The output generated by R was: F test to compare two variances data: x and y F = 0.9505, num df = 39, denom df = 29, p-value = 0.8707 alternative hypothesis: true ratio of variances is not equal to 1 95 percent confidence interval: 0.4675884 1.8646602 sample estimates: ratio of variances 0.950451 (CaseB) Now I want R to give me the same output as above, when I read it from a file. The commands I gave are: R library(ctest) x<-read.table("input1.txt", header = FALSE) y<-read.table("input2.txt", header = FALSE) var.test(x,y) The output generated by R is an error: "not enough x observations" I want to know why is this error coming, since I am creating the input1.txt from rnorm(50) and input2.txt from rnorm(30),which should mean the same? Pl. see if you can provide any help on this. regards shalabh gandhi ____________________________________________________________________ Get free email and a permanent address at http://www.netaddress.com/?N=1 ____________________________________________________________________ Get free email and a permanent address at http://www.netaddress.com/?N=1 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
>>>>> shalabhgandhi writes:> I want to use var.test function of ctest library. I was able to > generate the output from R on the input given by rnorm. But when I > store the same data set provided by rnorm into a file, and then read a > file into a dataframe, then on using the dataframe as a parameter to > var.test function - it gives an error "not enough x observations".> I'll explain the above mentioned problem of mine for which I am > mentioning 2 cases CaseA and CaseB. CaseA is working but CaseB is not.> (CaseA) Following are 4 commands I gave:> R > library(ctest) > x<-rnorm(50) > y<-rnorm(30) > var.test(x,y)> The output generated by R was: > F test to compare two variances> data: x and y > F = 0.9505, num df = 39, denom df = 29, p-value = 0.8707 > alternative hypothesis: true ratio of variances is not equal to 1 > 95 percent confidence interval: > 0.4675884 1.8646602 > sample estimates: > ratio of variances > 0.950451> (CaseB) Now I want R to give me the same output as above, when I read it from > a file. The commands I gave are:> R > library(ctest) > x<-read.table("input1.txt", header = FALSE) > y<-read.table("input2.txt", header = FALSE) > var.test(x,y)> The output generated by R is an error: > "not enough x observations"> I want to know why is this error coming, since I am creating the > input1.txt from rnorm(50) and input2.txt from rnorm(30),which should > mean the same?> Pl. see if you can provide any help on this.First, why send 8 *DIFFERENT* bug reports with the same content? Second, why is this a bug report? Have you read the FAQ first? You do not give enough information to help us reproduce your problem, as you do not describe what your "input1.txt" and "input2.txt" look like. I would assume that reading x back from your input file failed somehow, which should explain the problem. -k -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._