search for: ytrafo

...",ylab="breeding strategy",col="gold",lty=1) ### Kruskal-Wallis test kruskal_test(breeding ~ habitat, data = mydata, distribution = approximate(B = 9999)) ### Nemenyi-Damico-Wolfe-Dunn test (joint ranking) NDWD <- oneway_test(breeding ~ habitat, data = mydata, ytrafo = function(data) trafo(data, numeric_trafo = rank), xtrafo = function(data) trafo(data, factor_trafo = function(x) model.matrix(~x - 1) %*% t(contrMat(table(x), "Tukey"))), teststat = "max", distribution = approximate(B = 900000)) ### global p-value...

Is it possible to implement the Jonckheere-Terpstra test for ordered alternatives using the coin package: Conditional Inference Procedures in a Permutation Test Framework? I found jonckheere.test{clinfun}, but it uses a normal approximation when ties are present in the data. To make this concrete, I've include a small dataset. Thanks. --Dale Hollander and Wolfe, 1999 Table 6.6, pg 205

...vival data. I am using function "ctree" from the "party" package. I came up with this command: test <- ctree(Surv(time, event)~., data =data.test, controls=ctree_control(teststat="max", testtype="Bonferroni", mincriterion=0.95,savesplitstats = TRUE), ytrafo = function(data)trafo(data, numeric_trafo = rank), xtrafo=function(data)trafo(data, surv_trafo=logrank_trafo(data, ties.method = "logrank")) ) which works well but as I am not a statistician it is quite confusing and i might not run it properly. My technical problem is that I would lik...

...ch classes are different. i find this code in a R help: http://stackoverflow.com/questions/2478272/kruskal-wallis-test-with-details-on-pairwise-comparisons class <- m.class.length.lf var <- m.class.l dft <- data.frame(class,var) NDWD <- oneway_test(var ~ class, data = dft, ytrafo = function(data) trafo(data, numeric_trafo = rank), xtrafo = function(data) trafo(data, factor_trafo = function(x) model.matrix(~x - 1) %*% t(contrMat(table(x), "Tukey"))), teststat = "max", distribution = approximate(B=1000)) ### global p-value...