similar to: anova question

Dear UseRs, I'm making box plots from a data set that looks like this: Chr Start End GeneDensity ReadCount_Explant ReadCount_Callus ReadCount_Regen 1 1 1 10000 107.82 1.243 1.047 1.496 2 1 10001 20000 202.50 0.835 0.869 0.456 3 1 20001 30000 158.80 1.813 1.529 1.131

Dear all, I have a quasipoisson glm for which I need confidence bands in a graphic: gm6 <- glm(num_leaves ~ b_dist_min_new, family = quasipoisson, data = beva) summary(gm6) library('VIM') b_dist_min_new <- as.numeric(prepare(beva$dist_min, scaling="classical", transformation="logarithm")). My first steps for the solution are following: range(b_dist_min_new)

Hello R experts, I haven't found anything like the S function 'varcomp' as described in W.N. Venables & B.D. Ripley's 'Modern Applied Statistics ...' for R; does something comparable exist for R, or is planned for future releases? More generally, are there libraries with post-anova test procedures, like Student-Newman-Keuls? Or do those of you who frequently use

Hello R-Experts, 2 questions: does anybody have example functions for determining the power of, say, a t-test or a single factor anova with fixed or random factors? and: is there an implementation of Cochran's test for heterogeneity of variances (i.e., is there a frequency distribution of Cochran's C)? Many thanks for your attention ... Ulf (please CC: me, I'm not on the

I'm using ssplot for drawing a map of Austria and colour the nine provinces regarding their share of employment. Now I wanted to add the figures in each province and failed miserably. Using the locator() and text() function caused the error message "invalid graphics state". I try to show you what I have done below, maybe you can find a general fault in my codes. I know that it's

Dear R users, I want to group numerical values in classes with different size and count the values for each classes. My data is in this forma: d 15 12,5 30,4 20,5 80,4 100,5 8,2 40,5 33 21 11 And I want the group them in classes with 4 (5,....etc) cm size like this: class d 16 16 32 24 84 104 12 44 36 24 12 and final to count values for each class: class d

Is it possible to change the limits on plots that are already on the screen? In particular, is there any R equivalent to the sequence of matlab commands plot(1:10,1:10) hold on plot(2:12,5:15) I know I can use points and lines to add points and lines to plots, but the limits of the plot do not change when I do this. Looking at various examples, it seems that the answer is "no", but I

I noticed that congigure reported libpng as "not found" although a recent version exists. Further looking into this issue told me, that the test failed because the test from R.m4 AC_CHECK_LIB(png, png_create_write_struct, [ BITMAP_LIBS="${BITMAP_LIBS} -lpng -lz" AC_DEFINE(HAVE_PNG) ], , ${LIBS}) ], AC_MSG_RESULT([no])) ]) did

I would like to use gridBase to place four separate heatmaps (actually, a stripped-down heatmap.2 from ght gregmisc package that contains only the "image" part) into four different viewports. I can get the placement correct, but I keep 'losing' the previous plot. Any suggestions? Here is some quick example code trying to put a heatmap into the left viewport and then

Hello out there, in their "R Complements" to Modern Applied Statistics W. Venables and B. Ripley say that there is a lack of time series related functions in R compared to S plus. Looking at CRAN I didn't find much, too. As I have to learn something about +-basic ts analysis, e.g. spectral analysis, I would like to know whether somebody has developed a "unpublished" time

Hello out there, my first try to install a package (e.g. ctest from CRAN) didn't succeed; looking at the error messages I think it must be a configuration problem of my machine, but as I'm not very experienced I would be glad if someone could give me a hint ... R's version is 0.63.0, my i386-Linux is kernel 2.0.35. I hadn't had problems to compile R, and R itself seems to work.

I''m trying to use environments and seem to be failing. Right now I have 4 defined environments: production, cat, development, beta They are defined as follows on my puppetmaster: cat /etc/puppet/puppet.conf [main] pluginsync = true vardir = /var/lib/puppet manifest = /etc/puppet/environments/production/site.pp modulepath = /etc/puppet/environments/production/modules [master] reports =

Dear community, I have a question regarding the standard error of "marginal hazard ratio" estimates from a cox proportional hazard model with an interaction term. I have to dichotomous variables ADMA and ALA (both low/high) where I have divided the continuous variables at a certain level. I am aware about the contra arguments categorizing a continous variable and I could as well use

Thanks for your help, Rui! That works and will save me a lot of trouble. --Kelly -----Original Message----- From: Rui Barradas [mailto:ruipbarradas at sapo.pt] Sent: Tuesday, July 10, 2012 2:24 AM To: Vining, Kelly Cc: r-help at r-project.org Subject: Re: [R] boxplot with "cut" Hello, Maybe this iss what you're looking for. GD is your data.frame. multi.boxplot <-

Hi, I am trying to generate plots on our unix (no X) server to be included in our web-pages -> specifically pngs. Without X, png() obviously doesn't work, so I have been trying to define the graphics print device bitmap as described in the help: > bitmap(file="plot.png", type = "png256", height = 6, width=6,res=72) > pie(c(12,5)) >

Excuse me, if the following is FAQ or nonsense, I'm new to R (and the list)! After trying to import a textfile of 184897 Bytes, R (v0.63.0 on Linux 2.0.35/DLD 5.4) gives the following error message: > read.table ("db.out") Error: heap memory (1953 Kb) exhausted [needed 150 Kb more] > Is there a possibility to increase

Dear useRs We have a data-set (comma delimited) with 12Millions of rows, and 5 columns (in fact many more, but we need only 4 of them): id, factor 'a' (5 levels), factor 'b' (15 levels), date-stamp, numeric measurement. We run R on suse-linux 9.1 with 2GB RAM, (and a 3.5GB swap file). on average we have 30 obs. per id. We want to aggregate (eg. sum of the measuresments under

I have a similar question, I know about --vsize, but I use R under Emacs using ESS. I know there is an easy way to call R from ESS with command line options but I don't remember how. I can't find help on this in ESS documentation or in R documentation. Could someone please remind me how this is done. ------------------------------------ | Robert Denham | |

How's wine support coming with translucent windows? Guild wars 2 has a new launcher that makes extensive use of transparency to mimic the in-game interface, but of course, it turns into a massive black square here. Screenshot of the launcher, all that is black is supposed to be transparent: http://i.imgur.com/BqNo7.jpg I know it's nitpicky but google implied some progress on this and I

Anyone have a user configured auto attendant setup? Something that can be used without the * admin helping to make changes. Something where the operator can record the message like 'press 1 for john, 2 for bill, 3 for jean' and then the operator can enter the extension that gets dialed when the caller presses 1 or 2 or 3? This would be useful if Bill leaves the company, the operator can