[R] Dotchart and its arguments

2019-02-20 Thread Robert Zimbardo
Hi all I was recently trying to customise a dotchart of a matrix dats <- matrix(1:6, nrow=2, dimnames=list(R=letters[1:2], C=letters[14:16])) dotchart(dats) with pch and pt.cex and noticed some irregularities, namely that R doesn't use the values in the positions it uses for plotting also for th

[R] R memory limits on table(x, y) (and bigtabulate)

2017-07-03 Thread Robert Zimbardo
I have two character vectors x and y that have the following characteristics: length(x) # same as length(y) # 872099 length(unique(x)) # 47740 length(unique(y)) # 52478 I need to crosstabulate them, which would lead to a table with 47740*52478 # 2505299720 cells, which is more than 2^31 # 2

[R] Error message from allEffects(model) / effect(model)" ‘range’ not meaningful for factors"

2015-08-19 Thread Robert Zimbardo
Hi I cannot figure out why the effects package throws me error messages with the following simple code: rm(list=ls(all=TRUE)); set.seed(1); library(effects) # set up data x <- factor(rep(letters[1:3], each=100)) y <- c(rnorm(100, 3, 3), rnorm(100, 4, 3), rnorm(100, 5, 3)) # fit linear model m

[R] Two-tailed exact binomial test with binom.test and sum(dbinom(...))

2014-12-13 Thread Robert Zimbardo
Hi R experts, I have a few related questions that are actually a combination of an R and a hopefully not too trivial (?) statistics question, namely regarding the computation of an exact two-tailed binomial test. Let's assume the following scenario: - number of trials = 10 - p of success = 0.6

[R] anova comparisons

2013-02-23 Thread Robert Zimbardo
I have several linear models on the same data: m1 <- lm(y ~ poly(x,1)) m2 <- lm(y ~ poly(x,2)) m3 <- lm(y ~ poly(x,3)) What I don't understand is why anova(m1, m2, m3, test="F") - yields the same RSS and SS values, but a different p-value from anova(m1, m2, test="F") - when it also yields the S

Re: [R] Failing to print mer object in an RData image

2009-04-24 Thread Robert Zimbardo
t on a different computer. Question though: does that mean that the mer object is now inaccessible and the only way to get the data back is to run the analyses again? The changes are not backward compatible?? Thx, RZ > > > 2009/4/23 Robert Zimbardo : >> Hi all >> >> I hav

[R] Failing to print mer object in an RData image

2009-04-23 Thread Robert Zimbardo
Hi all I have problems in accessing a mer object called model.01 from a workspace that was created with R 2.8.1 and saved with save into an .RData file (on Windows XP or Ubuntu 8.10, don't remember anymore). Now I want to open it in R 2.9.0 on Ubuntu 8.10. I use # load workspace load("name.RData"

[R] R 2.8.0 for Ubuntu Intrepid Ibex i386

2008-11-05 Thread Robert Zimbardo
Hi all Is there going to be a binary of R 2.8.0 for Ubuntu Intrepid Ibex i386 - i.e., not just amd 64bit? Thx, RZ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do r