similar to: Circular statistics with (direction,size) data

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 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 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 <-

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

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:

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

I am plotting subsets of my data, using ylim. This works fine, but the outer margin line widths of the plot are thin, due to clipping. If I include > trellis.par.set(clip=list(panel = "off")) then the outer margin line widths are fine, but the outlying data is visible. Is there any way of achieving both correct margin line widths and clipping of outlying data? Thanks, Dan Bebber

I received the following error when I attempted to update my packages: updating HTML package descriptions Warning messages: 1: DESCRIPTION file of package 'file28862' missing or broken in: packageDescription(p, lib = lib, fields = pkgFlds) 2: number of columns of result not a multiple of vector length (arg 2) in: rbind(retval, c(p, lib, desc)) 3: DESCRIPTION file of package

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

Greetings, you kind and wonderful people. I've spent the past few hours or so (read: all day) attempting to get a basic PDC up and running, with no success. I've even gone to the extent of copying John Herpstra's config, which contains all sorts of weird and wonderful stuff, in the hope that he had something I missed. He didn't. I've got one box running debian testing, with

Could someone please give a brief explanation, or pointer to an explanation, of the following error: > arima(ts.growth, order = c(1,0,0),include.mean=T) Error in arima(ts.growth, order = c(1, 0, 0), include.mean = T) : non-stationary AR part from CSS and why it does not arise with > arima0(ts.growth, order = c(1,0,0)) Many thanks ____________________________ Dr. Daniel P. Bebber

I need to set the content to 0 because actually is not using malloc(), it uses calloc() that does set the memory to 0. Best Regards --------------------------------------------- Helios Cuevas L?pez Embedded Software Engineer Castle Care-tech Ltd. First Floor, 6 Bracknell Beeches, Old Bracknell Lane West, Bracknell, Berkshire, RG12 7BW Office - 01344 469470 Direct Dial - 01344 469485 FAX - 01344

Specifying xlim or ylim in the mapCountryData function of the rworldmap library do not alter the plotting region on my system. #Using the example from rworldmap library(rworldmap) mapCountryData() #uses the sample dataset mapCountryData(ylim = c(-45,45)) #makes no difference to the plot R 2.14.1 on Linux Mint 11 64bit. Thanks, Dan ________ Dr Dan Bebber Head of Climate Change Research,

Yanick & Jean-Marc, First of all thank you very much for taking your time to try to help me. Hi Yanick, This is what I am trying to do, but for example the function SpeexPreprocessState *speex_preprocess_state_init has lots of pointers that are initialize with static inline void *speex_alloc (int size). So I am trying to convert all this pointers into arrays but it is getting messy trying

Dear all, I am trying to use the Speex library pre-processor, I only need the VAD and Echo Canceller. I am using a STM32F103 processor, but I am using a tiny OS and it doesn't allow to use malloc. Is there a version of the library that will work without the need of using malloc? Obviously I will use all the time the same parameters, like 8kHz sampling rate, 160 frame... I will appreciate

Dear R-users, As noted by Paul Murrell < p.murrell at auckland.ac.nz > there is errors in the code for polar plotting I send to R-help under the title "Polar plot, circular plot (angular data)" at Thu Oct 17 2002 - 12:18:20 CEST. Thanks! I have reorganized the code into a structure ('pp'). This allows plots to be modified to a greater extent by passing arguments by ...

Dear R-users, Hereby a polar plot function for plotting angular data. I hope it will be usefull for some of you. I had a need to plot frequencies of wind-directions. The not-that-cheap SigmaPlot software did not allow me to change the orientation of the angular axis to clockwise orientation (what is used for meteorological observations). I even tried the latest version availible at the time

Dear list-members I am new to R and a statistics beginner. I really like the ease with which I can extract and manipulate data in R, and would like to use it primarily. I've been learning by checking analyses that have already been run in SAS. In an experiment with Y being a response variable, and group a 2-level between-subject factor, and time a 5-level within-subject factor. 2

I want to get the minimum arc (in degrees) needed to include a set of compass directions. I would like to use the range.circular() function of the package circular, because that package understands a compass-type of angle convention, but it gives results I don't understand. Howver, I can get the correct answer in the example below, 90 degrees, using the CircStats package. How can I make the

Dear All, Ulric Lund and me are developing a new package called 'circular' that as soon as possible will substitute 'CircStats' package. A pre-released is available at: http://www.dst.unive.it/~claudio/R/circular_0.1.tar.gz I kindly ask everybody is using CircStats to check this new package and send me comments, suggestions and bugs to claudio at unive.it and not to the list. As