similar to: a new-bie question about obtaining certain value from the print out

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

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 =

Hi, how can I extract the t-values after running a pairwise.t.test? The output just list the p-values. Many thanks for your help. Cheers Guido ____________________________________ Guido J. Parra School of Tropical Environment Studies and Geography James Cook University Townsville Queensland 4811 Phone: 61 7 47815824 Fax: 61 7 47814020 Mobile: 0437630843 e-mail:

It seems that MASS suggest to judge on the basis of sum(residuals(mode,type="pearson"))/df.residual(mode). My question: Is there any rule of thumb of the cutpoiont value? The paper "On the Use of Corrections for Overdispersion" suggests overdispersion exists if the deviance is at least twice the number of degrees of freedom. Are there any further hints? Thanks. -- Ronggui

hi: still experimenting. is there a function that tells me what an S object is, or how it is constructed? s <- cor.test ( x, y ); s$estimate$name = 'correlation' ; <- try to rename 'cor' to 'correlation' fails. obviously, name is not a part here. s$newfield = c("another info field", 3.2 ) ; <- this is not congruous so what.is(s) #tells

Good afternoon everybody, I'm writing a little function to visualise hypothesis testing. Therefore I need to extract the confidence level of a t-test. Here a little example: x <- str(t.test(1:10) gives List of 9 $ statistic : Named num 5.74 ..- attr(*, "names")= chr "t" $ parameter : Named num 9 ..- attr(*, "names")= chr

Dear all, I coded a function called u.glm u.glm <- function (x,ahi,age,bmiz,gender) { library(nortest) lil.rslt <- lillie.test(x) if (lil.rslt$p.value >0.05) { cat("Logtrans=0, lillie=",lil.rslt$p.value,"\n") xmodel<-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel) confint(xmodel) } else { cat("Logtrans=1,

Dear list, I have a series of data points which I want to approximate with exactly two linear functions. I would like to choose the intervals so that the total deviation from my fitted lines is minimal. How do I best do this? Thanks! Kamila The information transmitted in this electronic communication...{{dropped}}

What's the difference between t.test(x, y) and pairwise.t.test()? Is it just that the former takes two vectors, whereas the latter takes a vector and a factor?

Rusers, I am trying to append multiple .csv files of different dimensions (but # of columns are the same for all .csv files). I do have .csv files whose names are CA1.csv ~ CA100.csv. CA means california and 1 means the first file. So what I have been doing (after googling how to append by row multiple files) was: cleanup_data<-function(state,count) { out<-matrix() for (i in 1:

Hello All, I have a doubt regarding the fault tolerance of the rsync tool. What happens if the network goes down for some time and comes back after few seconds? Will rsync retry ? If yes, how many times/ how much time it will keep retrying ? How to ensure the completion of the rsync operation when there are network break ups? thanks for any help. -- Regards, Ravi.Thati -------------- next

Sorry guys, but I chose a really stupid name before (no "reference classes"). Hope it's okay to re-post. Cheers, Janko >>> ORIGINAL MESSAGE <<< Dear list, Is it possible to update or reassign '.self' with an image of '.self' (e.g. a locally stored .Rda file) from within a method? I know that this might sound akward, but here's the use

Hi list, I tried to load the package SAGx and failed because it complains it's looking for the Biobass which is not there. Then I looked up the package list and Biobass is not found. I'm trying to run the Jonckheere-Terpstra test and from what I see in the R archive, SAGx is the only place it's been implemented. > library(SAGx) Loading required package: multtest Loading required

Hi all, I am very new to this group and also the Unix platform. I am having some questions related to Throtling. I have visited lot of sites but not satisfied except www.docum.org I am from India. Working on *nix platform & love to work on RedHat Linux. Thanks in Adv. to Group BY __________________________________ Do you Yahoo!? SBC Yahoo! - Internet access at a great low price.

hi guys, I''m reading up on Rspec, Mocha and some related material to put BDD in my new rails app. I have also checked out Ryan Bates'' railscasts on rspec (that''s how I got to know about Mocha). Reading up on the Rspec''s main site, the main example in http://rspec.rubyforge.org/rspec/1.3.0/ does not show any use of assert_equals. Rather it just uses the

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

When I do t.test on two distributions (see example below), it outputs numerous data about the t.test. What I'd like to do is individually capture some of this data and assign it to other variables. However, I am unable to find anything in the help section. In the example below, the t value is -4.0441 and the p-value is 0.006771 How can I assign these values to two variables, let's

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

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20081112/618073fe/attachment.pl>

I'd like to learn the process of revising R functions & packages and then submitting proposed patches to the R Core team. Would someone be willing to mentor me through one example? For starters, consider an example. I'd like to revise the t.test function to return the stderr value to the user. We only need to change the "rval" in the third-from-the end line of