similar to: box() and hist()

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

Hi, I'm confused about how to specify random and fixed factors in an aov() term. I tried to reproduce a textbook example: One fixed factor (Game, 4 levels) and one random factor (Store, 12 levels), response is Points. The random factor Store is nested in Game. I tried > str(kh.df) `data.frame': 48 obs. of 4 variables: $ Subj : Factor w/ 48 levels

Hi everyone. I have an histogram like this one: http://imageshack.us/photo/my-images/442/dfsdfsdj.jpg/ I would like to remove the extra space under the histogram so it "touch" the x axis. Is it possible? Here's my code. hist(X, freq = F, col = 'gray', axes = FALSE) abline(v = cyano_euk_min, col = 'red', lty = 3, lw = 2) abline(v = cyano_euk_max, col =

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, 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

Hi, Does anybody encounter the same problem when we overlap histogram and density     that the density line seem to shift to the right a little bit?           If you do have the same problem, what should we do to correct that?           Thank you.           par(mar=c(4,4,2,1.2),oma=c(0,0,0,0))     hist(datobs,prob=TRUE, main ="Volume of a catchment from four    

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

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,

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:

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 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 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, I want to get a histogram with the legend for my data. I drew a normal density curve and kernel density curve in the histogram, and I also label mean and median in the X axis. From the code, I got two legend: One shows "Normal Density" and "Kernel Density" and their corresponding lines, the other shows "Mean = value" and "Median = value" and their

hi all, I'm trying for quite some time to have an x- and y-axis, but no entire box. >plot(..,axes=F) >axis(1) >axis(2) Gives this, but their axes do not go to the origin. Quite a number of people find this gap between the two axes disturbing. Has anyone an idea how to let these axes go to the origin? thank you in advance [[alternative HTML version deleted]]

Hi, I apologise for the stupid question but how to you get rid of the box around a plot in the package zoo? I can't seem to find an equivalent for bty="l" i.e. just x and y axis. Cheers James -- View this message in context: http://www.nabble.com/Zoo-Plot-Can%27t-Get-Rid-of-Box-tp21468452p21468452.html Sent from the R help mailing list archive at Nabble.com.

How about -- d is a vector as you described it, then define d1 <- d[1:(length(d)-1)] d2 <- d[2:length(d)] d3 <- d1!=0 & d2!=0 & (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

Hi all, I don't understand the following output. I've created a data subset from a data frame by > p1.sub <- subset(p1.dat, vp!="p1") this is ok. But > attach(p1.sub) > vp [1] p1ab p1ab p1ab p1ab p1ab p1br p1br p1br p1br p1br p1kf p1kf p1kf p1kf p1kf [16] p1mg p1mg p1mg p1mg p1mg p1mw p1mw p1mw p1mw p1mw Levels: p1 p1ab p1br p1kf p1mg p1mw shows me that the

Hello there, Is there a functionality or command that generates axes in the shape of a capital Greek letter gamma (upside down L) ? I use these for making sediment profiles. I have read through help-lists, tried different things and asked several people but never got a satisfactory result. I've tried the following: xx <- c(2,4,6,5.8,5,4,3,1.5,0, 0, 0, 0) yy <- (0:-11) plot(NULL,