Displaying 4 results from an estimated 4 matches for "numeric_trafo".
2013 Jan 11
0
Error with looping through a list of strings as variables
Dear R users:
I have been trying to figure out how to include string variables in a for
loop to run multiple random forests with little success. The current code
returns the following error:
Error in trafo(data = data, numeric_trafo = numeric_trafo, factor_trafo =
factor_trafo, :
data class character is not supported
In addition: Warning message:
In storage.mode(RET@predict_trafo) <- "double" : NAs introduced by coercion
The code runs fine with the data before I add the > for (h in varlist){ loop.
Loops i...
2012 Jan 09
2
Unexpected results using the oneway_test in the coin package
...quot;,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
print(pvalue(NDWD))
### sites-by-site p va...
2012 Aug 23
0
party package: ctree - survival data - extracting statistics/predictors
...uot;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 like to extract the statistics output from my...
2012 Mar 26
0
Different result with "kruskal.test" and post-hoc analysis with Nemenyi-Damico-Wolfe-Dunn test implemented in the help page for oneway_test in the coin package that uses multcomp
...his 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
print(pvalue(NDWD))
[1] 0.074
99 percent...