Search Archives

Hi, The spatial library is already in rw1041\library: library(spatial) library(help=spatial) but do check Ribeiro's library geoR at CRAN too. For an introduction, see http://www.maths.lancs.ac.uk/~ribeiro/geoR/ Cheers, Tord At 15...

Dear all, I try to plot include expression(alpha) in a legend(): x<- seq(0,150) par(mfrow=c(1,1)) matplot(x, outer(x, c(0.001, 0.01, 0.05, 0.1, 0.5), function(x, k) exp...

Hi, I would like to plot a map of US counties using different colors. map() seems to be the function to use, e.g. library(maps); map('usa'); map('county', 'colorado', add=T,fill = T, col=c(1:5)) plots...

Dear Duncan, At 19:58 2003-10-20 -0400, you wrote: >On Mon, 20 Oct 2003 18:50:22 +0200 (MET DST), you wrote: > > >>I have read CHANGES and can't find changes about return(). > >CHANGES is the list of...

Dear all, I want to aggregate() my data based of the classes of a factors, just as in the help file, but how do I define the "Region varable" in the ?aggregate example? Now, the column is just like any...

Dear all, I get an error message when I use polr() in MASS. These are my data: skugg grupp frekv 4 1 gr3 0 5 2 gr3 3 6 3 gr3 6 10 1 gr5 1 11 2 gr5 12...

Dear all, I try to load the nlme library for Linear and nonlinear mixed effects models by Jose Pinheiro, Douglas Bates and Saikat DebRoy but I get the following error message: > library(nls) > library(nlme) Error in dyn.load(x...

Dear all, I would like to standardise my data. Here is an exemple for two columns: > onorm<- read.table("onorm.txt", header=T) > names(onorm) [1] "on1" "on2" "on3" "on4" "on5" > norm<- data.frame((onorm$on1-mean(onorm$on1))/stdev...

Dear all, I've used glm(family=binomial(link="logit")) several times, but now I think that a log link is more appropriate. I want to fit a model for probability of tree fall (TF)), with tree diameter (dbh) and...

Dear all, Thanks for, as always the quick replies! I have a suggestion: add the upper.tri() and lower.tri() under "See Also" in the diag() help page. I searched for a such function there, and I can see that...

Dear all, Sorry to bother you with one more newbie question. I have a dataobject with several hundreds of columns. I want to remove columns with a mean of the column values below a certain value: > a<- c(1,2...

Dear all, I have a tree data matrix. For some trees I lack info about tree species, but I want to set them to be spruce. For some reason the tree species names on the remaining (non-NA) rows are...

Dear all, A question related to the following has been asked on R-help before, but I could not find any answer to it. Input will be much appreciated. I got an unexpected sign of the "slope" parameter associated with...

Dear all, I have written a function for calculating the volume of a tree (=trad) or snag (=h?gst). The included volume regreesion model includes ten parameter values, which are tree species specific. bj?rk.formh?jd.pars is an...

Thanks to you all for prompt replies! As evident from the replies pasted below there are different solutions to this problem. cheers, Tord ------------------------------------------------------------------------------ "There are some tools avaiable to transform maps to \pkg{GeoBUGS} format: the function \code{sp2WB} of...

To whom it may concern, I get the following message when I run my function: Warning message: multi-argument returns are deprecated in: return(call.fn, repl, time, from, to, last.year, occup.m, ant.occ.m, > version platform i386...

Dear all, Does anyone have a suggestion on how to make this code more efficient? I have a large data frame that I want to merge with another data frame b.art.temp<- b.art[, c("BalteSFyndID", "RUBINKOD")] I first...

Dear all, I've worked with optim before but never encountered this error message: Nelder-Mead direct search function minimizer 0.23 0Error: subscript out of bounds The error seems to depend on the initial parameter values. However, strangely (I...

Dear all, I try extracting information from a list with several levels, but I would be happy for recommendation on writing more efficient code: > h0<- seq(0,100, by = 20); expo<- seq(0.1, 0.5, l = 5) > grid<- expand...

Dear all, I use the code below for my barplot and it looks fine on my screen. But when I paste into my MSWord file, after having copied it by presing the camera, the colours change. The same actually happens...