search for: garbad

Hi all, how do I change the magnification of xaxis annotation? I need a smaller axis text size in some of my plots. I tried cex.axix=0.5 in plot(), but this doesn't work. Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hi all, I don't understand the behavior of is.integer(): > x <- integer() > is.integer(x) [1] TRUE > x <- 10 > is.integer(x) [1] FALSE > x <- 1:10 > is.integer(x) [1] TRUE Why is.interger() returns FALSE if x has only one element? And how can someone check if x is an integer but contains only one value? (R 1.5.0 on Linux i386) Thanks, Sven

Hi all, there is a vector v with several NAs. I want to create a new vector n of the same length as v and the same NAs as in v and tried this: n <- vector(length=length(v), mode="numeric") replace(n, which(is.na(v)), NA) but this does't work, all values in n are 0. What went wrong? Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Hi all, suppose I've got a vector y with some data (from a repeated measure design) observed given the conditions in f1 and f2. I've got a model with two unknown fix constants a and b which tries to predict y with respect to the values in f1 and f2. Here is an exsample # "data" y <- c(runif(10, -1,0), runif(10,0,1)) # f1 f1 <- rep(c(-1.4, 1.4), rep(10,2)) # f2 f2 <-

Hi there, if I copy an x11() graphics device to an eps-file (with dev.copy2eps()) the font in the legend is very large and doesn't fit to the legend box in the eps-file (same with a postscript file). I'm not sure if this is a R problem rather than a ghostscript one. But is there a way to solve this problem in R or depends this on my ghostscript installation? System: R Version 1.2.3 on

...(d1-d2)!=0 I think d3 is the final vector you wanted. I didn't test it, so I'm not sure my logic is right, but I think this is definitely a path you can take... -----Original Message----- From: owner-r-help at stat.math.ethz.ch [mailto:owner-r-help at stat.math.ethz.ch]On Behalf Of Sven Garbade Sent: Wednesday, September 18, 2002 9:35 AM To: r-help at stat.math.ethz.ch Subject: [R] comparing elements in a vector Hi all, there is a vector d which contains three different valus: 0, 1 or 2. I wrote a loop which gives me T, when the values in d changed from 1 to 2 or from 2 to 1, otherwis...

Hallo R-users, I'm a R (release 0.65) novice, but have read some about S. Is there a similar function in R like the S function "by()" ? This function groups data typically from a matrix by a variable and applies a function on the grouped data. Example in S: > by(data, year, summary) where "data" is a matrix containing the variable "year". After grouping

Dear R-list, is it possible to convert a threedimensional array to a matrix? The array has the form , , -4 -4 -3 -2 -1 0 1 2 3 4 -80 1588 NA 1171 NA 1121 NA 852 NA 1580 0 1497 NA 1311 NA 762 NA 1185 NA 1466 80 1960 NA 1257 NA 941 NA 1435 NA 1636 [...] I need a matrix like -80 -4 -4 1588 -80 -3 -4 NA -80 -2 -4 1171 etc. I?ve tried write.table() with no succes:

Morning, I've plotted multiple figures on one page, setup the page with n2mfrow. This works fine, but I can't figure out how to change to a specific plot (e.g. first row, second column) to add some details. Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hi there, I want to read all data from a binary file (with package Rstreams). I tried the following function (small piece of code): readdata <- function (file) { # open file s <- openstream(file, "read") # read int values while (s$position <= s$size) { readint(s, 1, 2) -> type [...] } return() } and got the following error: > Error in

Hi there, when plotting a time series I got this error message: Error in xy.coords(x, y, xlabel, ylabel, log) : x and y lengths differ In addition: Warning messages: 1: longer object length is not a multiple of shorter object length in: vx.1[int.lr] + (-1) * vx.2[int.rl] vx.1[int.rl] and vx.2[int.rl] have the same length. Does anybody know what the warning message means? Thanks, Sven

Thanks for the advice, folks! Shravan Vasishth <vasishth at ling.ohio-state.edu> Try pico. There's no substitute for/comparison with emacs, though, IMHO. but just to be sure: you don't have to have ESS to use (x)emacs for editing R scripts. Sven Garbade <garbade at psy.uni-muenchen.de> ESS can be installed via rpm on Mandrake. baron at cattell.psych.upenn.edu (Jonathan Baron) ESS is the real answer to your question, and anything else has got to be more difficult. Nels Tomlinson <tomlinso at purdue.edu> Installing ESS on debian was...

Hi all, I can't figure out how to compute Wilks Lambda in a one way repeated measure design. My matrix looks like: > t2.m Blank ECR ENC UEA UED 1 -0.15 0.14 0.16 0.09 0.14 2 0.30 0.08 0.14 0.14 0.14 [...] where each row is a case and the columns are levels of one factor (named trial): > t2.fit <- manova(t2.m ~ 1) > summary(t2.fit, intercept=T,

Try vcov() from the MASS library. -----Original Message----- From: Sven Garbade [mailto:garbade at psy.uni-muenchen.de] Sent: Monday, June 24, 2002 11:19 PM To: R-Help Subject: [R] Covarianz matrix from an aov model Hi all, how can I extract the covariance matrix from an object of class aov? Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

Hi all, I computed a Wilks Lambda Test with manova: > df.man <- manova(df.mul ~ 1, na.rm=TRUE) > summary(df.man, intercep=T, test="Wilks") Df Wilks approx F num Df den Df Pr(>F) (Intercept) 1 0.0002824 393.3 9 1 0.03911 * Residuals 9 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*'

Hi, does anybody know why the following is not working: > hist(rnorm(200)) > box(bty="o") gives me a box without rounded corners. System: > R.Version() $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status [1] "" $major [1] "2" $minor [1] "0.0" $year

Dear list, how can I compute the inverse of the X'WX matrix ("inverse of the weighted sum of squares and crossproducts matrix") from an object of class "lm" from a weigthed linear regression? Thanks, Sven

Hi list, can anybody point me to the trick how glm is computing the dispersion parameter in quasi-poisson regression, eg. glm(...,family="quasipoisson")? Thanks &regards, Sven

...ur request was to reposition to the beginning of the file, and that is what happened. This was already fixed in R-patched: from NEWS there: o seek() on file connections was ignoring the `origin' argument. It's always worth looking at the current patched version. On Tue, 29 Jan 2002 garbade@psy.uni-muenchen.de wrote: > Full_Name: Sven Garbade > Version: 1.4.0 > OS: Linux i386 > Submission from: (NULL) (129.187.254.49) > > > It seems that seek() can't reposition a file connection, but it can report the > actual file position correctly if the file connectio...

Dear list members, I have problems to interpret the coefficients from a lm model involving the interaction of a numeric and factor variable compared to separate lm models for each level of the factor variable. ## data: y1 <- rnorm(20) + 6.8 y2 <- rnorm(20) + (1:20*1.7 + 1) y3 <- rnorm(20) + (1:20*6.7 + 3.7) y <- c(y1,y2,y3) x <- rep(1:20,3) f <- gl(3,20,