search for: wtest

Hi, i need to use pspline. In this pspline function coxph.wtest was used. When I try to make some change to this function by pulling out the pspline function, it turns out R gave me an error msg, saying coxph.wtest cannot be found. Even if i dont change anything in pspline and just rename it and run the function, it did not work out. Can any one help me with th...

Hello, I use t.test for normal distributed and wilcox.test for non-normal distributed samples. It is easy to write a function for t.test that calculates the effect size, because all parts of the formula are available from the t.test result: r = sqrt(t*t / (t*t + df)) However, for Wilcoxon tests, the formula for effect sizes is: r = Z / sqrt(N) I wonder how I can calculate the Z-score in R for

Hi, I'm trying to configure a Domain trust between Samba 3.5.6 (TEST domain) and Windows 2008 r2 (WTEST Domain). Samba is the trusting side and Windows is the trusted side. I created the "incoming trust" in the W2K8 called TEST. Then I executed the "net rpc trustdom establish WTEST" and got: Enter TEST$'s password: Could not connect to server W2K8SERVER Could not initialise...

...= TRUE) pw.test[i] <- (test2$p.value < 0.05) test3 <- ks.test(samp1, samp2,alternative = c("two.sided"), exact = TRUE) pks.test[i] <- (test3$p.value < 0.05) } pow.ttest[j] <- sum(pt.test)/sim.size # Calculate powers for jth setting pow.wtest[j] <- sum(pw.test)/sim.size pow.kstest[j] <- sum(pks.test)/sim.size } # End of j loop Error in pt.test[i] <- (test1$p.value < 0.05) : object "pt.test" not found Chat on a cool, new interface. No download required. Go to http://in.messenger.yahoo.com/webme...

I propose a patch to ifelse that leverages anyNA(test) to achieve an improvement in performance. For a test vector of length 10, the change nearly halves the time taken and for a test of length 1 million, there is a tenfold increase in speed. Even for small vectors, the distributions of timings between the old and the proposed ifelse do not intersect. The patch does not intend to change the

Hi, Would be very glad for help on this problem. Using this code: temp<-function(x, bins, tot) { return(as.numeric(lapply(split(x, bins), wtest, tot))); } wtest <- function(x, y) { return(wilcox.test(x,y)$p.value); } rs <- function(x, bins) { binCount <- length(split(x[,1], bins)); tot <- as.numeric(x); result<-matrix(apply(x, 2, temp, bins, tot), nrow=binCount, byrow=F); rownames(result)<-names...

...lcox.test on a loop, when the test fail, I would like to continue doing the next iteration and getting the p-value. I tried to do tryCatch or try but I cannot retrieve the p-value if the test is not fail. sample code: test2=list(numeric(0),c(10,20)); test1=list(c(1),c(1,2,3,4)); for (i in 1:2){ wtest=wilcox.test(test1[[i]],test2[[i]]) } i=1 will fail, I want to ignore this and get the pvalue for i=2. Thanks, Lin [[alternative HTML version deleted]]

...A, -53L)) ?str(winter) ?w_melt <- melt(winter,id="week",variable="year");str(w_melt) ?#? Reorder DOESN'T WORK ?w_melt <- reorder(w_melt$year,c("BY2005","BY2009","BY2006","BY2008","BY2007","BY2010")) pdf("wtest.pdf") ggplot(w_melt,aes(week,value/1000,colour=year,order= - as.numeric(year))) + geom_line(size=.75)+? theme_bw() + opts(title="Cumulative",axis.text.x = theme_text(angle=45,hjust=1)) + ?labs(y="Number of individuals X 1,000",x="week") dev.off() ? Felipe D. Carr...

Hi, I want to fit multiple cox models using the coxph() function. To do this, I use a for-loop and save the relevant results in a separate matrix. In the example below, only two models are fitted (my actual matrix has many more columns), one gives a warning message, while the other does not. Right now, I see all the warning message(s) after the for-loop is completed but have no idea which model

Dear R-forumites, I want to apply a Wilcoxon test on subsets of the data frame mydata, splitted using the myindice variable. When I send : wilcoxtest <- by(mydata, mydata$myindice, function(x) {wilcox_test(x$value~x$fact)}) I get : Error in eval(expr, envir, enclos) : object "x" not found whereas wilcoxtest <- by(mydata, mydata$myindice, function(x)

Dear R-users, I have noticed small discrepencies in the reported estimate of the variance of the frailty by the print method for survreg() and the 'theta' component included in the object fit: # Examples in R-2.6.2 for Windows library(survival) # version 2.34-1 (2008-03-31) # discrepancy fit1 <- survreg(Surv(time, status) ~ rx + frailty(litter), rats) fit1 fit1$history[[1]]$theta

Dear List, How do I extract the approximate Wald test for the frailty (in the following example 17.89 value)? What about the P-values, other Chisq, DF, se(coef) and se2? How can they be extracted? ######################################################> kfitm1 Call: coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id, dist = "gauss"), data = kidney)