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

On 2024-11-28 8:36 a.m., Michael Dewey wrote: > To make the plot clearer I have removed the axes but I wish to remove > all the boz except the bottom horizontal line. Using the bty parameter > does not seem to enable me to just leave the horizontal line at the > foot. I can get the "l" version to remove everything except the left > hand side and the base. I could also

Thank you Duncan, I will try that next. Michael On 28/11/2024 13:52, Duncan Murdoch wrote: > On 2024-11-28 8:36 a.m., Michael Dewey wrote: >> To make the plot clearer I have removed the axes but I wish to remove >> all the boz except the bottom horizontal line. Using the bty parameter >> does not seem to enable me to just leave the horizontal line at the >> foot. I can

That's pretty similar to call axis function, i.e. ?plot(rnorm(100), yaxt="n", bty="n") ?usr <- par("usr") ?## lines(usr[c(1,2)], usr[c(3,3)], xpd = TRUE) ?axis(1, ...) # where 1 is bottom, 2 left, 3 top, and 4 right sides of the plot box) Best Fer On 11/28/24 14:52, Duncan Murdoch wrote: > On 2024-11-28 8:36 a.m., Michael Dewey wrote: >> To

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

To make the plot clearer I have removed the axes but I wish to remove all the boz except the bottom horizontal line. Using the bty parameter does not seem to enable me to just leave the horizontal line at the foot. I can get the "l" version to remove everything except the left hand side and the base. I could also remove the entire box but that looks odd. Am I reading the

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