[Rd] Pl. provide and Input for var.test (PR#716)

Kurt Hornik Kurt.Hornik@ci.tuwien.ac.at
Thu, 2 Nov 2000 08:10:18 +0100 (CET)


>>>>> 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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._