similar to: Density plot with lattice?

Hi, I have to analyse the number of provisioning trips to nestlings according to a number of biological and environmental factors. I was thinking of building a mixed-effects model with species and nestid as random effects, using a Poisson distribution, but the data are overdispersed (variance/mean = 5). I then thought of using a mixed-effects model with negative binomial distribution, but I have

Dear all, I need to do the same procedure on several files. But I don't know how to refer to the file name. Here is an example of what I am trying to do. List of files: file1(A,B,C, D1...Dn), file2(A,B,C,E1,...,En), file3(A,B,C,F1,...,Fn) Procedure I want to apply on each file: dft <- melt(df,id=c('A','B','C')) dft$X <- substr(dft$variable,1,3) dft$Y <-

Hi I would like to eliminate a large number of lines of the dataframe df1 The lines to delete are given here by the values of Mat (ex : 2,4,7,10). but I have a large number (300) values of Mat dput(df1) structure(list(Mat = c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 10, 10, 10, 11,

Hi all, I'm trying to get a vertical line at a specific point in a densityplot. abline seems to be what's required, but it doesn't align itself to the scale used in the plot. example: library(lattice) x<-rnorm(100) plot.new() densityplot(x) abline(v=0) ----- The line seems to use some other coordinate system. What kind of call do I use to make abline use the graph's

I'm trying to use a lattice function within a function and have problems passing the "groups" argument properly. Let's say I have a data frame d <- data.frame(x = rnorm(100), y = c("a", "b")) and want to plot variable x in a densityplot, grouped by the variable y, then I would do something like densityplot(~ x, d, groups = y) If however I wanted to

Hi, I have created two separate overlapping density plots- see example code below. What I wish now to do is combine them into one figure where they sit side by side. Any help would be great! many thanks in advance, josh. ##################### thedataA <- data.frame(x1=rnorm(100,1,1),x2=rnorm(100,3,1)) #create data thedataA.m<-melt(thedataA) densityplot(~value, thedataA.m,

Hi, I have the following script: ---- t.R --- grafica <- function() { v <- read.csv('preprocessed/komolongma.ece.uprm.edu.active',sep=',') x <- as.ts(v$active) bitmap(file="output.png") densityplot(~x,col='blue',main='Density Plot') dev.off() } grafica() ---- t.R --- When I "sourced" it from R prompt, it quietly runs.

Hi All, I'm trying to add a 'rug' representation of my data to a plot created with densityplot(). While I can do this in the simple case, I can't do it properly when I include the "groups" argument. I have an example below. I am running a reasonably new version of R. print(sessionInfo()) R version 2.12.0 Patched (2010-11-07 r53537) Platform: i686-pc-linux-gnu

Hi Dear R users, For a pair plotting, usaully we use par( ) function. Apparently it does not work anywhere. I want to have 3 plots in a single figure, like this: par(mfrow=c(3,1)) densityplot( a) densityplot(b) densityplot(c) But it does not work. How is it possible to have such a figure with densityplot( ) in a single figure? So many thanks for any help. Amir Safari

Hi folks, This is really minor but to someone not familiar with the various tentacles of the lmer package it could be really annoying. I was trying to plot the posterior density of the fixed effect parameters of a lmer model, > hr.mcmc = mcmcsamp(hr.lmer, n=50000) > densityplot(hr.mcmc, plot.points=F) There is this error, "Error in densityplot(hr.mcmc, plot.points = F) : no

Hi all, I'm interested in doing a multiple density plot on a number of columns in a dataframe. >DF lineA.1 lineA.2 lineB.1 lineB.2 r1 5.355354 6.665575 10.288498 11.74750 r2 3.643415 5.427600 11.407112 13.97065 r3 5.813674 6.438502 9.628871 11.57456 r4 5.241340 5.125049 10.456221 12.35381 r5 4.640885 8.635518 8.344704 11.98484 r6 4.559852 6.416171 10.419599 10.89247 r7

Hi, i have a script how i launch lattice to make a densityplot. in the script: jpeg(file="XXX.jpg") densityplot(~f_diametro+m_diametro+n_diametro, plot.points="rug", auto.key=T) dev.off() does'nt work and in R i dont have any output but if i launch by R row by row, runs correctly..... any idea?

Hi I would like to superimpose group-specific densityplots on top of an overall histogram using panel.histogram and panel.densityplot. Furthermore, I would like to automatically adjust the range of the y-axis to take into account the ranges of both histogram and densityplot. This last part is where I have a problem. I believe using the prepanel argument of histogram is typically the way to go,

Dear R community, I am using densityplot (lattice package) for a large dataset and wish to print it to a jpeg (the pdf is huge). R crashes consistently. Am I doing it wrong or is densityplot incompatible with jpeg? I work on a Mac, R 2.7.0. > require(lattice) > jpeg("test.jpeg") > d[1:10] [1] 0.700218 0.700175 0.700357 0.700847 0.698286 0.701520 0.698158 0.699300 0.698819

Hello, I'm doing some graphics for a paper and a need customize such with filled region above the density curve. My attempts I get something very near what I need, but I don't solve the problem of use semitransparent filled. Below a minimal reproducible code. Someone has any idea? require(lattice) # toy data... dt <- expand.grid(A=1:2, B=1:3, y=1:50) dt$y <- rnorm(nrow(dt), dt$B,

Dear People, If I define foo <-function() { x <- rnorm(500) trellis.device(postscript, file="plot.ps") densityplot(~x) dev.off() } and call foo() then plot.ps is just a blank plot. However, if I extract the lines inside foo and run them, ie. x <- rnorm(500) trellis.device(postscript, file="plot.ps") densityplot(~x) dev.off() then plot.ps is created

I am trying to plot multiple densityplots on each panel, and using panel.groups to do some additional plotting (not included in the example) as in this example. library(lattice) thedata <- data.frame(x=rnorm(1200),class=rep(1:3,each=200,times=2), group=rep(1:2,each=100,times=6)) densityplot(~x | class,groups=group, xlab='x', panel.groups =

Is there any way of adding text to a density plot? I have had a go using the text() function but I think the error is because this function doesn't work with densityplot(). Alternatively, I understand I can achieve pretty much the same result if I plot a density kernel estimate using plot() (which allows text()), but I do prefer densityplot(). Also, is it possible to specify the dimensions

I'm trying to plot a series of densities using/comparing differing weights. I see the reference to weights and subscripts, but I don't understand how to implement that. My data are of the form: I, J, Actual, Distance, Subset, Weight1, Weight2, ... I'm trying to see the effect of the distance distribution (Actual by Distance) compared to the various weighted distributions

densityplot(~x+y) does what I expect it to do if x and y have equal length. I know how to get what I want if x and y have different lengths. But what is this actually doing if x and y have different lengths? The relevant example is x=rnorm(10) y=rnorm(50,1) densityplot(~x+y)