search for: myowelt

Hi there, I tried this code from homepage: http://myowelt.blogspot.de/2008/04/beautiful-correlation-tables-in-r.html <http://myowelt.blogspot.de/2008/04/beautiful-correlation-tables-in-r.html> corstarsl <- function(x){ require(Hmisc) x <- as.matrix(x) R <- rcorr(x)$r p <- rcorr(x)$P ## define notions for significance levels; spacing...

...t/interpret a main effect of stereotype threat in the confirmed presence of the interaction effect GENDER:STTHREAT, because a main effect of stereotype threat would actually be caused by the interaction (an error-bar plot illustrating this can be found here if one scrolls a little downwars: http://myowelt.blogspot.com/2008/05/obtaining-same-anova-results-in-r-as-in.html) . I thus tend to use Type-II SSs and calculate my ANOVA with > library(car) > Anova(lm(MWP ~ GENDER * STTHREAT), type="II") Anova Table (Type II tests) Response: MWP Sum Sq Df F value Pr(&gt...

Hi! I searched some time ago a way to get the Wilcoxon test results as a table more or less formatted. Nobody told me any solution and I found nothing on the Internet. Recently I came across this link ( http://myowelt.blogspot.com/2008/04/beautiful-correlation-tables-in-r.html), which helped me to find a solution. Here's the solution (I'm using R Commander): W <- as.matrix(lapply(Dataset[2:11], function(x) wilcox.test(x ~ GrFac, alternative="two.sided", data=Dataset)$statistic)) P <- as...

Dear all, I have been trying to investigate the behaviour of different weights in weighted regression for a dataset with lots of missing data. As a start I simulated some data using the following: library(MASS) N <- 200 sigma <- matrix(c(1, .5, .5, 1), nrow = 2) sim.set <- as.data.frame(mvrnorm(N, c(0, 0), sigma)) colnames(sim.set) <- c('x1', 'x2') # x1 & x2 are