similar to: multiple plots problem

Hello. I am new to R. And, I want to perform a multiple nested anova on a large datasets (with 9448 observations). Under the helps from R-Sig-ecology mailing list, I have gained many progresses. But I still have some confusions. I want to ask for some helps here. my dataset("SeedL.txt") was not attached. Data are not sorted by factors. In this dataset, "SpecN"

I am new to the R language. I am trying to plot multiple figures on one page through a loop, but the code just produce one graph on one page. Can someone show some light on what's wrong? Here is my code: library("quantreg") tcdata<-read.table("mydata.txt",header=TRUE) postscript("myfigure.ps") basins<-

Hello, R-help, Although this is my first post to the mailing list, I have been a subscriber for a while now, and have found all of the posts to be very informative in my quest to attempt to master R. Thank you for that in advance. This question involves the use of the par() variables pin, oma, and mar to resize "sub-plots" in an output window. I am aware from ?plot that the use of

Dear R users. I'm new with layout and I can't figure how to teak my graphs. I have the following code: layout(matrix(c(1,2,3), 3, 1, byrow = TRUE), heights=c(0.3,0.3,0.6)) boxplot(rnorm(100), horizontal=TRUE, axes=FALSE) boxplot(rnorm(100), horizontal=TRUE, axes=FALSE) hist(rnorm(100)) Is it possible to have the two horizontal boxplot closer to each other? With regards, Phil

Dear R-users, I often stack plots that have the same x-axis. To save space and have the plots themselves as large as possible I like to minimize the margins between the plots to zero. I use the "mfrow" and "mar" parameters to achieve this. However, the different margin settings for the individual plots lead to the inner plots being higher than the two outer plots. To make the

Hola a todos, en efecto, lo crucial es el new, he utilizado estas líneas para mi código: par(new=T, mfrow=c(1,1), mar=c(0,4,2,2), oma=c(1,1,0,1)) plot(-1,type="n",axes=F,xlab='',ylab='')

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings around each. Each tree has a TreeID, which I use as the random effect (blocking factor). There are two

Dear R gurus, I have XY data giving the locations of tree seedlings that were surveyed during a 210 meter belt transect. This belt transect was taken by stretching a line across the field, then measuring all seedlings within 1 meter on either side of the line. The end result was XY coordinates and height for ~1,300 seedlings. I would like to use that data to calculate density of

Dear all I was trying to put a title for my plot, but i got this result, > x11(width=10,height=5,title="seedling");par(mfrow=c(1,2))Error in x11(width = 10, height = 5, title = "seedling") : unused argument(s) (title = "seedling")> The title is not taking How can i give a title for the plot (where in need to make 2 plots within a window of x11() ) thanking

Hi everybody, I am trying to fit a model with the lmer function for mixed effects. I have an experimental design consisting of 5 field plots. Each plot is divided in 12 subplots where the influence of three factors on the growing of tree seedlings is tested: (1) seed (1 = presence; 0 = absence); (2) seedling species (oak holm vs. pine); (3) treatment (three different treatments). In each of

Hi, I ran the example on pp. 799-800 from Machael Crawley's "The R Book" using package survival v. 2.36-5, R 2.13.0 and RStudio 0.94.83. The model is a Cox's Proportional Hazards model. The result was quite different compared to the R Book. I have compared my code to the code in the book but can not find any differences in the function call. My results are attached as well as a

Full_Name: Charles J. Geyer Version: 1.8.0 OS: i686-pc-linux-gnu (Suse 8.2) Submission from: (NULL) (134.84.86.22) Two bugs (perhaps related, perhaps independent) revealed by the same Poisson regression with offset mydata <- read.table(url("http://www.stat.umn.edu/geyer/5931/mle/seeds.txt")) out.fubar <- glm(seedlings ~ burn01 + vegtype * burn02 + offset(log(totalseeds)),

Hi, I'm new to R and GLMMs, and I've been unable to find the answers to my questions by trawling through the R help archives. I'm hoping someone here can help me. I'm running an analysis on Seedling survival (count data=Poisson distribution) on restoration sites, and my main interest is in determining whether the Nutrients (N) and water absorbing polymer Gel (G) additions to the

Hello, I am trying to use the projection.matrix( ) function and am following the example given. I have my data formatted very similar to the test.census example. > str(AsMi05mat) `data.frame': 1854 obs. of 6 variables: $ Tag : num 501 502 503 504 505 506 507 508 509 510 ... $ Year : int 1995 1995 1995 1995 1995 1995 1995 1995 1995 1995 ... $ Length : num 34

Hi there, I've looked through the very helpful advice about adding fitted lines to plots in the r-help archive, and can't find a post where someone has offered a solution for my specific problem. I need to plot logistic regression fits from three differently-sized data subsets on a plot of the entire dataset. A description and code are below: I have an unbalanced dataset

Hello, After putting together interaction code that worked for a single pair of interactions, when I try to evaluate two pairs of interactions( flowers*gopher, flowers*rockiness) my computer runs out of memory, and the larger desktop I use just doesn't go anywhere after about 20 minutes. Is it really that big a calculation? to start: mle2(minuslogl = Lily_sum$seedlings ~ dnbinom(mu = a,

Could someone explain the meaming of this error message from augPred: > augPred(area3.pen.nlme, primary=~day) Error in predict.nlme(object, value[1:(nrow(value)/nL), , drop = FALSE], : Levels 1,2,3 not allowed for block > predict.nlme(area3.pen.nlme) does not produce an error. area3.pen.nlme was created with: > area3.pen.nlme <- nlme(area ~ SSlogis(day, Asym, xmid, scal),

Hola a todos, Mi pregunta es si hay alguna manera (sencilla/rápida o no) de colocar una leyenda común para dos gráficos colocados en paralelo con par(mfrow=c(1,2)). He probado con par(xpd=NA) y con par("usr") pero no encuentro la solución. Está claro que podría poner la misma leyenda para cada gráfica, pero quedaría mejor con una global para ambos. Muchas gracias de antemano por la

Dear R users, I have been having problems getting believable estimates from anova on a model fit from lmer. I get the impression that F is being greatly underestimated, as can be seen by running the example I have given below. First an explanation of what I'm trying to do. I am trying to fit a glmm with binomial errors to some data. The experiment involves 10 shadehouses, divided between

Hi, I'm new to R and GLMMs, and I've been unable to find the answers to my questions by trawling through the R help archives. I'm hoping someone here can help me. I'm running an analysis on Seedling survival (count data=Poisson distribution) on restoration sites, and my main interest is in determining whether the Nutrients (N) and water absorbing polymer Gel (G) additions to the