similar to: Compare two 3d plots

Working with the built in R data set Orange, e.g. with(Orange, table(age, circumference)). How should I go about about grouping the ages and circumferences in the following ranges and having them display as such in a table? age range: 118 - 664 1004 - 1372 1582 circumference range: 30-58 62- 115 120-142 145-177 179-214 Thanks for any feedback and insights, as I hoping for an output that

Dear ExpeRts, I am trying to plot a bunch of growth curves and would like to get some more control over groups and line colors than I seem to have. Example: # make some data dat <- Orange dat$group <- ifelse(dat$Tree%in%c('1','4','5'), 'A', 'B') # plot xyplot(circumference~age, dat, groups=group) # now use lines to make the growth curve more

Hello, I want to be able to make a plot that has several series with different color and linetype. Online documentation suggest that this is possible, but I haven't found how: "We can also create redundant mappings, mapping the same variable to multiple aesthetics. This is most useful when producing a graphic for both colour and black and white display." Here's what I have to

Hi folks, I am tearing my hair out on this one. I am using an example from Pinheiro and Bates. ### this works data(Orange) mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange ) ### This works mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange, fixed = Asymp + xmid + scal ~ 1, start =

Hi List, While using 'nlme' function, I have encountered the similar problem Dr. Stevens and Dr. Graves observed (please see the posts: https://stat.ethz.ch/pipermail/r-help/2006-May/105832.html ). I have tried Dr. Stevens's original example, the problem is still there, > mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal), + data=Orange ) >

>>>>> "VW" == Viechtbauer Wolfgang (STAT) <Wolfgang.Viechtbauer at STAT.unimaas.nl> >>>>> on Fri, 26 Sep 2008 18:00:19 +0200 writes: VW> Hi all, I'm trying to fit a marginal (longitudinal) VW> model with an exponential serial correlation function to VW> the Orange tree data set. However, R crashes frequently VW>

Dear all, How to obtain the odds ratio (OR) and 95% confidence interval (CI) with 1 standard deviation (SD) change of a continuous variable in logistic regression? for example, to investigate the risk of obesity for stroke. I choose the happening of stroke (positive) as the dependent variable, and waist circumference as an independent variable. Then I wanna to obtain the OR and 95% CI with

Hello Everyone, I have three variables: Waist circumference (WC), serum triglyceride (TG) level and gender. Waist circumference and serum triglyceride is numeric and gender (male and female) is categorical. From these three variables, I want to calculate the "Lipid Accumulation Product (LAP) Index". The equation to calculate LAP is different for male and females. I am giving both

Well, something like: LAP <- ifelse(gender =='male', (WC-65)*TG, (WC-58)*TG) The exact code depends on whether your variables are in a data frame or list or whatever, which you failed to specify. If so, ?with may be useful. Cheers, Bert On Fri, Nov 3, 2023 at 3:43?AM Md. Kamruzzaman <mkzaman.m at gmail.com> wrote: > Hello Everyone, > I have three variables: Waist

df$LAP <- with(df, ifelse(G=='male', (WC-65)*TG, (WC-58)*TG)) That will do both calculations and merge the two vectors appropriately. It will use extra memory, but it should be much faster than a 'for' loop. Regards, Jorgen Harmse. ------------------------------ Message: 8 Date: Fri, 3 Nov 2023 11:10:49 +1030 From: "Md. Kamruzzaman" <mkzaman.m at gmail.com>

I found this by accident yesterday. The program that crashes is the first two lines of the example from the help page for nlmer. That example hasn't changed in a long time, so I assumed that it is an R-devel issue. It could also be a long latent nlmer bug. The second run with valgrind is puzzling. Terry T. > library(lmer) > sessionInfo() R Under development (unstable)

Thanks, posted to https://github.com/lme4/lme4/issues/412 for further discussion ... On Wed, Feb 22, 2017 at 10:03 AM, Therneau, Terry M., Ph.D. <therneau at mayo.edu> wrote: > I found this by accident yesterday. The program that crashes is the first > two lines of the example from the help page for nlmer. That example hasn't > changed in a long time, so I assumed that it is

Sometimes one might like to obtain pointwise bootstrap bias-corrected, accelerated (BCA) confidence intervals for a large number of statistics computed from a single dataset. For instance, one might like to get (so as to plot graphically) bootstrap confidence bands for the fitted values in a regression model. (Example: Chiu S et al., Early Acceleration of Head Circumference in Children with

Just a minor point in the suggested solution: df$LAP <- with(df, ifelse(G=='male', (WC-65)*TG, (WC-58)*TG)) since WC and TG are not conditional, would this be a slight improvement? df$LAP <- with(df, TG*(WC - ifelse(G=='male', 65, 58))) -----Original Message----- From: R-help <r-help-bounces at r-project.org> On Behalf Of Jorgen Harmse via R-help Sent: Friday,

Hello, I am new to R and would like to know how to label data points in the matrices of scatterplots made by the pairs() command. To be more specific, I want to assign a number to each data point, instead of the small circumference that appears as a data point. If anyone here knows if its possible to do this with R, I would greatly appreciate your help. Thank you, Mauricio Esguerra PhD

Yes, that will halve the number of multiplications. If you?re looking for such optimisations then you can also consider ifelse(G=='male', 65L, 58L). That will definitely use less time & memory if WC is integer, but the trade-offs are more complicated if WC is floating point. Regards, Jorgen Harmse. From: avi.e.gross at gmail.com <avi.e.gross at gmail.com> Date: Friday,

I'm wondering how the cex parameter is used to scale circles (i.e. does it scale the radius, diameter, area, circumference, etc.?). In my case I'm using lattice with filled circles (pch=19). Based on example, it looks like R scales the radius of the circle: library(lattice) dta<-data.frame(x=rep(1,6),y=rep(1,6),sz=c(1,2,4,8,16,32))

Hi, I have a puzzeling problem. I want to export graphics from R to TeX via postscript(). This works fine for some graphs, but for others, the eps remain empty when viewed with GSView. When such an empty eps is imported to TeX, the figure appears upside down and very small, irrespective of TeX width and height commands. If I transform the eps to pdf, the graphic shows up, but turned aqround 90°.

Hi, I have a matrix that is indexed by groups and looks something like this: 1 1 2 1 2 1 . . . 2 1 1 2 1 2 2 NA 1 . . . 3 1 NA 3 2 NA etc. The first column is the group variable and I would like to apply categorical data imputation functions to the other two columns, doing so by groups. I have tried APPLY, BY and SPLIT, but have not had much luck getting it to work. I wonder if anyone has

Hi. The following statement works fine in R-win 1.8.0, but yields a syntax error in R-win 2.0.1 (and possibly in other versions after 1.8.0): plot(c(1,2),main=expression(a==b==c)) I note that the following workaround executes successfully in both versions of R... plot(c(1,2),main=expression(a*"="*b*"="*c)) ...but I don't really understand