similar to: R-alpha: new class for chisquare tests?

From: Martin Maechler <maechler@stat.math.ethz.ch> To: Kurt.Hornik@ci.tuwien.ac.at CC: R-devel@stat.math.ethz.ch Subject: Re: R-alpha: class for chisquare tests; Thoughts on print & summary Kurt, I think we hopefully are coming to an agreement that 1) your ctest collection should make its way into 'R core' ((and you are the one who can make it happen now ..)) By all

The newest version of R for Unix (version 0.50 alpha-1) is now (or will soon be) available from the following sites. NORTH AMERICA: http://lib.stat.cmu.edu/R/Alpha EUROPE: ftp://ftp.stat.math.ethz.ch/R/ ftp://statlab.uni-heidelberg.de/pub/mirrors/auckland/R/ JAPAN: ftp://ftp.u-aizu.ac.jp/pub/lang/R/ NEW ZEALAND: ftp://stat.auckland.ac.nz/pub/R/

Hi all, I just used power.t.test() in Version 1.9.0 beta (2004-03-31) today under FC1. The following is the output: > power.t.test(delta = .5, power = .95, type = "paired") $n [1] 53.94062 $delta [1] 0.5 $sd [1] 1 $sig.level [1] 0.05 $power [1] 0.95 $alternative [1] "two.sided" $note [1] "n is number of *pairs*, sd is std.dev. of *differences* within

In this years biostat teaching I will include Rcommander (it indeed simplifies syntax problems that makes students frequently miss the core statistical problems). But I could not find how to make a simple chisquare comparison between observed frequencies and expected frequencies (eg in genetics where you expect phenotypic frequencies corresponding to 3:1 in standard dominant/recessif

Hello! First of all I must appologize if this has been raised previously, but search provided by Robert King at the University of Newcastle seems to be down these days. Additionally let me know if such a question should be sent to R-help. I did a contribution to function hwe.hardy in package 'gap' during the weekend. That functions performs Hardy-Weinberg equilibrium test using MCMC. The

Hello, I've posted a question about this subject yesterday, but since there was no R code to comment, no one did. I'm trying to have the print method for class 'htest' print some extra information common in some test, like the time series linearity related tests. Many of them have an 'order' parameter, representing a lag or embedding dimension, and it would be a nice

Hello, everyone, I don't quite understand the following message: > TTT <- t.test(1:10, y=c(7:20)) > class(TTT) [1] "htest" > TTT@p.value Error: Trying to get slot "p.value" from an object whose class ("htest") is not defined > TTT$p.value [1] 1.855282e-05 Why the message says the class of TTT is not defined while class(TTT) gets

I have a dataset called "raw-data" . I am trying to use the following code - col_name<-names(raw_data) for (i in 1:(length(names(raw_data))-2)) { tbl=table(raw_data$Pay.Late.Dummy, raw_data$col_name[i]) chisqtest<-chisq.test(tbl) } Say the 1st column of my raw_data is Column1. The idea is when i=1 then raw_data$col_name[i] will automatically become raw_data$Column1 ,

Hi all .. there is no doubt a simple answer to this, but it eludes me. In the first session below ( with jarque.bera.test) you will see that p.value prints with a name of X-squared . This is easily fixed by changing the source to assign a more appropriate name - no name is assigned in the source listing below (the original source code of jarque.bera.test() from tseries).. but what I

In the past months (maybe even years), I have received a variety of suggestions on improving/extending/... my ctest package. Some of the wishes have already been fulfilled. In the future, I will at least: * provide a formula-based interface * provide an option for controlling the behavior of 2-sided p-values (sum of more extreme, tail balancing a la S-PLUS, double one-sided, mid p-value,

Dear R developers, sorry about the empty bug-report that I seem to have filed earlier. (I wrote it via `bug.report()' from within R using `emacsclient'. I wasn't aware that I have to save the file after pressing C-x # and before sending it.) Hopefully this time it will work...... I noticed the following problem with the on-line help system: 1) I start an R session and type

Hi, New to both R and SEM, so this may be a very simple question. I am trying to run a very simple path analysis using the sem package. There are 2 exogenous (FARSCH, LOCUS10) and 2 endogenous (T_ATTENT, RMTEST) observed variables in the model. The idea is that T_ATTENT mediates the effect of FARSCH and LOCUS10 on RMTEST. The RAM specification I used is FARSCH -> T_ATTENT, y1x1, NA

Hi, I am trying to test for pairwise associations between genotypes ( Rows=individuals, Columns =genes, data are up to 4 genotypes per gene, some with 2,3 or 4) where each chisquare comparison is different depending on the genes tested. The test is the observed multilocus (across columns for each individual) genotypes vs the expectation, which is the product of the individual frequency for each

o <- (structure(list(sand.silt = structure(list(statistic = structure(185, .Names = "W"), parameter = NULL, p.value = 0.0478835773838087, null.value = structure(0, .Names = "location shift"), alternative = "two.sided", method = "Wilcoxon rank sum test with continuity correction", data.name = ".column by site"), .Names =

Hello, I have a question concerning the R-function chisq.test. For example, I have some count data which can be categorized as follows class1: 15 observations class2: 0 observations class3: 3 observations class4: 4 observations I would like to test the hypothesis whether the population probabilities are all equal (=> Test for discrete uniform distribution) If you have a small sample size

Hi everyone, I am quite new to R. my command is t.test(c(1:5,7:11), y=c(1:10),alternative = c("two.sided"), paired = TRUE) The output is Paired t-test data: c(1:5, 7:11) and c(1:10) t = 3, df = 9, p-value = 0.01496 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: 0.1229738 0.8770262 sample estimates: mean of the differences

Here are a handful of oddities in pre0.50-7 (at some risk of reprisals from people waiting for the real 0.50...): I was mucking about trying to figure out how to know which tags has been used in a function call, and some weird stuff happened. ---------- A call with tagged arguments is something like a list, the tags can be used to access elements, but the names attribute is absent, until the

Hello all you helpful people out there! I am stil R Beginner using R 2.5.1 on a Apple Power Book G4 with Mac OS X 10.4.10 . I have two tabels like this: MalTabChi X1 X4 X6 X8 X10 X14 X21 X24 X29 X38 X43 X50 X67 X76 X78 X80 X82 X84 Anth_cap 1 1 1 1 6 5 1 45 12 4 12 6 56 5 2 4 1 1 Anth_crin 12 1 1 2 76 5 1

# Why does expressing one function require(ctest) t.test # return only function (x, ...) UseMethod("t.test") <environment: namespace:ctest> # but expressing another function shapiro.test # returns more complete code? function (x) { DNAME <- deparse(substitute(x)) x <- sort(x[complete.cases(x)]) n <- length(x) if (n < 3 || n > 5000)

Hi all, I'm getting contradictory results from bptest and ncvTest on a model calculated by GLS as: olslm = lm(log(rr)~log(aloi)*reg*inv, data) varlm = lm(I(residuals(olslm)^2)~log(aloi)*reg*inv, data) glslm = lm(log(rr)~log(aloi)*reg*inv, data, weights=1/fitted(varlm)) Testing both olslm and glslm with both ncvTest and bptest gives: > ncvTest(olslm) Non-constant Variance Score Test