search for: chidata

...(vdata$X08-vdata$Y08) Z09 <- (vdata$X09-vdata$Y09) I would like to use chisq.test for each "row" and output the p-value for each in a stored variable. I don't know how to do it. Can you help? so far I have done it for one row (but I want it done automatically for all my data): chidata=rbind(c(vdata$Y08[1],Z08[1]),c(vdata$Y09[1],Z09[1])) results <- chisq.test(chidata) results$p.value I tried removing the [1] and the c() but that didn't work... Any ideas? THANKS!

...ll, I am trying to use the scan function in an R script that I am calling from the command line on a Mac; at the shell prompt I type: $ Rscript get_q_values.R LRT_codeml_output in the hope that LRT_codeml_output will get passed to the get_q_values R script. The first line of that script is: chidata <- scan(file="") which, as I understand how scan works, will read the contents of the file from the command line into the object chidata. I did this a few times and it worked like a charm. And then, it stopped working. Now, every time I try to do this, I get "Read 0 items" a...