similar to: lattice: clipping data, not plot margins

Dear list, If anyone knows how to assign dimnames to matrices or arrays I would be most grateful for help. I've tried various permutations of likely-looking code but get error messages every time. I could find no example in the documentation. Many thanks, Dan Bebber Department of Plant Sciences University of Oxford South Parks Road Oxford OX1 3RB UK Tel. 01865 275000

Dear R list I have an m * n matrix P and a vector V of length n containing indices for rows in P. For each of the m columns I want to extract the value in the row specified by V, and put these values into a new vector W of length n. At present I am doing this with a for.... loop, but I imagine there is a faster way that doesn?t involve loops. If anyone knows the way I would be most grateful.

Dear R list: I modelled changes in a variable (mconc) over time (d) for individuals (replicate) given one of three treatments (treatment) using: mconc.lme <- lme(mconc~treatment*poly(d,2), random=~poly(d,2)|replicate, data=my.data) summary(mconc.lme) shows that the linear coefficient of one of the treatments is significantly different to zero, viz. Value Std.Error

Dear list, there are now >400 packages available on CRAN. Would it be useful to classify these packages according to what they do (e.g. classification, graphics, spatial statistics), to assist the user in finding the appropriate package for their problem? Or perhaps the search facility is enough. I would attempt such a classification, but my knowledge of statistical methods isn't good

Hello, I wrote a wrapper for symbols() that produces a bivariate bubble plot, for use when plot(x,y) hides multiple occurrences of the same x,y combination (e.g. if x,y are integers). Circle area ~ counts per bin, and circle size is controlled by 'scale'. Question: how can I automatically make the smallest circle the same size as a standard plot character, rather than having to

Hello, My data are numbers of trees in plots sampled in a number of forest stands. Some stands were subjected to a treatment, others not. Several plots were sampled per stand to get a better idea of what the stand means were, but replication is really at the stand level. Therefore I think this is a split-plot design. I would like to know whether the treatment affected the number of trees, so:

Some notes on installing rgeos in R 2.14.0 on a Linux Mint 11 (x86_64-pc-linux-gnu) machine 1. rgeos 0.1-15 will not run with GEOS 3.2.0-1, which is the version currently available on Synaptic package manager 2. I installed GEOS 3.3.1 (the latest version) from http://trac.osgeo.org/geos/, but rgeos will not run with this either 3. I then install GEOS 3.2.3, and now rgeos works fine! It seems

I would like to extract predicted failure times from a coxph model in library(survival). However, none of the prediction options ("lp", "risk", "expected", "terms") seem to bear any relationship to failure time. Perhaps I am asking the wrong question, but can coxph provide predicted failure times? Thanks, Dan Bebber Department of Plant Sciences University

Hello, when I try to use the filled circle (pch = 16) in a Windows metafile, it appears highly pixelated rather than as a smooth vector. The other filled circles (pch = 19 and 20) are vector circles, filled with pixels. Results are the same whether I use windows() and save or copy as a metafile, or win.metafile() Also, some of the other symbols are misshapen, e.g. the diamond pch = 23 is not

Hello, I've just located the illuminating explanation by Douglas Bates on degrees of freedom in mixed models. The take-home message appears to be: don't trust the p-values from lme. Questions: Should I give up hypothesis testing for fixed effects terms in mixed models? Has my time spent reading Pinheiro & Bates been in vain? Is there a publication on this issue? Thanks, Dan Bebber

I noticed that clipping occurs in lossy audio compression even if there is no clipping in the original file (though the original file has peaks that are just below the maximum). I know that this happens due to all the filtering involved during compression, but I'm wondering just how audible this clipping is because I don't hear anything wrong. I mean, as things are now, is it safe to say

Dear List, I have circular data comprising direction and size for a number of observations. I wish to calculate statistics such as the mean direction and concentration for these data. I note that the CircStats package only deals with angular data (i.e. "Each observation is treated as a unit vector or a point on the unit circle". There appears to be no facility for weighting

Motivated by the discovery of 'loglet analysis' (http://phe.rockefeller.edu/LogletLab/) that allows one to decompose growth curves into a series of logistic equations, I attempted to do the same thing in R. SIMULATED DATA Time <- 1:200 pop.size <- SSlogis(Time,10,20,5) + SSlogis(Time,20,100,20) + rnorm(length(Time)) MY ANALYSIS results <- nls(size ~ SSlogis(Time, Asym1, xmid1,

Dear list, Apologies, I have sent this message before but received no replies so I'm trying again just in case... Motivated by the discovery of 'loglet analysis' (http://phe.rockefeller.edu/LogletLab/) that allows decomposition of growth curves into a series of logistic equations, I attempted to do the same thing in R. #SIMULATED DATA Time <- 1:200 pop.size <-

I have obtained shapefiles for Indian states from here: http://www.maptell.com/index.php?option=com_remository&Itemid=159&func=fileinfo&filecatid=115&parent=category Problem: I want to extract centroid coordinates for each State, but there is some coding problem with the shapefiles that prevents this. #Code: #After extracting the shapefiles from the india_state.zip file, then:

Hello, I am having trouble specifying a suitable nlme model. My data structure is described by gd <- groupedData(ppath ~ lcut | exp, outer = ~ bait, inner = ~ weight, data = d) i.e. the response (ppath) of several subjects (sub) was measured at levels of a continuous variable (lcut). Subjects were given either of one level of a factor (bait), and all subjects were measured at two levels

Hello! When I draw points on grid's viewport with clip set to TRUE clipping does not have effect on some symbols. It happens only when drawing over gtkDevice (GTK ver. 1.2) (with postscript device clipping works OK - take a look at this screenshot http://www.hppi.troitsk.ru/Kondrin/clipbug.png) and only with symbols that are drawn with commands other than gdk_draw_polygon() (i.e. circles

Hi, with the Lame MP3 encoder I reduce the amplitude of the input signal prior to encoding with the --scale option to prevent clipping that can be introduced due to resampling errors. Is Vorbis prone to this sort of clipping? <p>Headless --- >8 ---- List archives: http://www.xiph.org/archives/ Ogg project homepage: http://www.xiph.org/ogg/ To unsubscribe from this list, send a

I'm baffled as to how the Lattice package achieves clipping. Would someone mind explaining this to me? Firstly, my attempt using "just" the grid package: x<-seq(0,3,by=0.3)/2.8 y<-seq(0,1,by=0.1) grid.newpage() grid.rect(gp=gpar(fill="pink")) vp<-viewport(width=0.8, height=0.8) pushViewport(vp) grid.rect(gp=gpar(fill="white")) grid.xaxis()

Hello, Below is a function (test.map) that permits drawing the same map using three different devices. The "pdf" device doesn't clip polygons to the plot region as I see it does by both the native device (in my case "Quartz") and the "png" device. test.map("pdf") # produces "test-map.pdf" with no clipping test.map("png") #