similar to: Package update announcements

Displaying 20 results from an estimated 3000 matches similar to: "Package update announcements"

2005 Oct 26
2
changing memory limits to speed up lsoda
Hi All, I am running R 2.2.0 on Mac OS 10.4.2, dual G5 processors with 8 Gig RAM. I am running a simulation with lsoda that requires ~378 s to complete one set of time intervals. I need to optimize the parameters, and so need to considerably speed up the simulation. I have tried to figure out how to change the appropriate memory allocation and have search R help and Introductory
2007 May 31
2
Different fonts on different axes
Hi Folks, How do I get red bold font on my y axis and black standard font on my x axis? plot(runif(10), ylab="Red, Bold?", xlab="Black, standard?") Any pointers or examples would be great. Thanks! Hank Dr. Hank Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243
2008 Aug 12
1
tilde on a spanish keyboard?
I was trying to help someone who used a spanish keyboard on a PC running Windows. he discovered that he had no tilde key. Does anyone know of any simple work-arounds? Thanks, Hank Dr. Hank Stevens, Associate Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.cas.muohio.edu/~stevenmh/
2006 Jun 09
1
binomial lmer and fixed effects
Hi Folks, I think I have searched exhaustively, including, of course R-help (D. Bates, S. Graves, and others) and but I remain uncertain about testing fixed effects with lmer(..., family=binomial). I gather that mcmcsamp does not work with Do we rely exclusively on z values of model parameters, or could we use anova() with likelihood ratios, AIC and BIC, with (or without)
2005 Oct 10
3
Under-dispersion - a stats question?
Hello all: I frequently have glm models in which the residual variance is much lower than the residual degrees of freedom (e.g. Res.Dev=30.5, Res.DF = 82). Is it appropriate for me to use a quasipoisson error distribution and test it with an F distribution? It seems to me that I could stand to gain a much-reduced standard error if I let the procedure estimate my dispersion factor (which
2006 May 17
1
nlme model specification
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 =
2006 Jun 14
2
lmer binomial model overestimating data?
Hi folks, Warning: I don't know if the result I am getting makes sense, so this may be a statistics question. The fitted values from my binomial lmer mixed model seem to consistently overestimate the cell means, and I don't know why. I assume I am doing something stupid. Below I include code, and a binary image of the data is available at this link:
2005 Jul 01
2
Simple indexing conundrum
My apologies in advance for my thickness but I can't seem to solve the following, seemingly simple, data manipulation problem: I have a data frame that contains multiple factors and multiple continuous response variables, but duplicates of some factor combinations. The duplicates contain bad data, so I would like to eliminate the duplicates. I would like to retain the entire rows
2005 Jul 27
4
odesolve/lsoda differences on Windows and Mac
Hi - I am getting different results when I run the numerical integrator function lsoda (odesolve package) on a Mac and a PC. I am trying to simulating a system of 10 ODE's with two exogenous pulsed inputs to the system, and have had reasonably good success with many model parameter sets. Under some parameter sets, however, the simulations fail on the Mac (see error message below). The
2007 Jan 03
1
mcmcsamp and variance ratios
Hi folks, I have assumed that ratios of variance components (Fst and Qst in population genetics) could be estimated using the output of mcmcsamp (the series on mcmc sample estimates of variance components). What I have started to do is to use the matrix output that included the log(variances), exponentiate, calculate the relevant ratio, and apply either quantile or or HPDinterval to get
2004 Sep 16
1
geoR/variog4() not returning all directions
Mac OS 10.3.5, R 2.0.0 latest version of geoR I have an incomplete 5 x 20 spatial array of samples (60 out of 100 possible locations) for which I would like to calculate directional variograms using variog4(). Unfortunately, I can't get it to return all 4 directions. It returns variograms for 45, 90, and 135 degrees, omitting 0 degrees (pi/4, pi/2, 3pi/4, omitting 0). If I specify 0
2005 Oct 14
1
lattice with predicted values
Dear lattice wizards, I am trying to figure out how to plot predicted values in xyplot, where the intercept, but not the slope, varies among conditioning factor levels. I am sure it involves the groups, but I have been unsuccessful in my search in Pinhiero and Bate, in the help files, or in the archive, or in my attempts on my own. My example follows: FACT is a factor with levels a,b,c
2006 Jun 14
3
appending
All, In the function below I have 24 individuals and 6 calculations per individual. The 6 calculations are collected each time in a 1:24 loop when calculating "delta". I'd like to collect all 144 = 24*6 calculations in one vector ("delta.patient.comb"). The function works as is via indexing, but is there an easier way to collect the measurements via appendinng the 6
2007 Jan 15
1
Conflict in .Rprofile documentation FAQ vs. Help?
Hi folks, I note that in the general FAQ's we have 7.25 Why did my .Rprofile stop working when I updated R? Did you read the NEWS file? For functions that are not in the base package you need to specify the correct package namespace, since the code will be run before the packages are loaded. E.g., ps.options(horizontal = FALSE) help.start() needs to be
2006 Feb 16
2
Help to find correlation.
Respected Sir, I am trying to import excel file into R, but I need to truncate some columns from the original file. How to delete unwanted columns when I import data from excel file. How to use cor.test for the data when I want the output rowwise. How to do grouping and use cor.test on that data I need some help regarding how to calculate the correlation. I don't know whether you
2007 Aug 13
2
Error message when using zero-inflated count regression model in package zicounts
I have data on number of vines per tree for ~550 trees. Over half of the trees did not have any vines and the data is fairly skewed (median = 0, mean = 1.158, 3rd qu. = 1.000). I am attempting to investigate whether plot location (four sites), species (I'm using only the four most common species), or tree dbh has a significant influence on the number of vines per tree. When I
2006 Dec 11
2
How to write a two-way interaction as a random effect in a lmer model?
Dear All, I am working with linear mixed-effects models using the lme4 package in R. I created a model with the lmer function including some main effects, a two-way interaction and a random effect. Now I am searching how I could incorporate an interaction between the random effect and one of the fixed effects. I tried to express the interaction in:
2007 Jul 25
1
Is Rwiki down?
Dear R-friends, is Rwiki down? I've been trying to login for the past couple of days with no success. Yours sincerely, Jin ____________________________________________________________________________________ wherever you're surfing.
2005 Oct 17
1
Dunn's post hoc test
Hi Everyone. I am rather new to R and I've been trying to implement a function to carry out the above test. For a couple of days now I've been stuck on how to generate average rank differences. Say I have a vector of average ranks: averank<- c(2,5,9,12) I would like to subtract averank[1] from averank[2], averank[1] and averank[2] from averank[3] and averank[1], averank[2] and
2006 May 25
1
is possible to use update with lmer models
Hi, the question on subject is possible? example: m1 <- lmer(y~1+(1|subject)) m2 <- update(m1,.~.+x) I try but dont work, exist any other method for this? Thanks Ronaldo -- He won a Toyota now, eh? -- pal?ndromo -- |> // | \\ [***********************************] | ( ? ? ) [Prof. Ronaldo Reis J?nior ] |> V [UNIMONTES/Depto. Biologia