similar to: problem with putting objects in list

Hi I have a table called npl containing results of simulations. It contains about 19000 entries and the structure looks like this: NoPlants sim run year DensPlants 1 6 lng_cs99_renosterbos 1 4 0.00192 . . . it has 43 different entries for sim and year goes from 1 to 100, and run from 1 to 5. I would like to calculate the mean of DensPlants for each

Hi I have the following code which I would like to simplify. Id does linear regressions and returns the r-squares, and the coefficients. It runs slow, as it is doing the regressions for each - is it possible to get the values in a dataframe which looks as follow: expert | xx | seeds | r.squared | slope | intercept Thanks in advance, Rainer library(reshape) rsqs <- as.data.frame(

Hi is it possible to specify the shape of the point to be used in ggplot (as with pch in plot)? I couldn't find anything in the help. Thanks Rainer -- Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation Biology (UCT) Department of Conservation Ecology and Entomology University of Stellenbosch Matieland 7602 South Africa Tel: +27 - (0)72 808 2975 (w) Fax: +27 - (0)86 516 2782 Fax:

Hi I heard so much about Emacs and ESS that I decided to try it out - but I am stuck at the beginning. Is there anywhere a beginners manual for Emacs & ESS to be used with R? even M-x S tells me it can't start S-Plus - obviously - but I want it to start R... Any help welcome (otherwise I will be stuck with Eclipse and R) Rainer -- Rainer M. Krug, Dipl. Phys. (Germany), MSc

Hi I asked about if there is any way of opening URLs from the help browser in the same window of the same dillo browser - here is the answer. Just to reiterate: dillo is for me the perfect browser for the help of R when you use ?... Rainer -------- Original Message -------- Subject: Re: [Dillo-dev] Opening new URL in same instance and -s option Date: Thu, 21 Sep 2006 21:31:57 -0400 From:

Hi I generate a new dataframe by doing: npl.agg <- aggregate(npl$DensPlants, list(year=npl$year, sim=npl$sim), mean, na.rm=TRUE ) Now I want to plot it by using coplot(npl.agg$x ~ npl.agg$year | npl.agg$sim, type="l") but, as npl.agg$year is seen as a factor, the order of the points on the x-axis (time axis) does not follow the numerical sorting 1...100, but rather the text

Hi I have (8 measures * 96 groups) = 768 datasets for which I did linear regressions using lm(). Now I want to compare the slopes for each of the 8 measures in each of the 96 groups. As I understand , I can not use > anova(lm1, ..., lm8) as the lm1 ... lm8 are based on different datasets. I also read in previous discussions in this list, that I can see if the slope +- stddev(slope)

Hi the function density() does normally integrate to one - I've checked it and it works and I also read the previous threads. But I realised that it does not integrate to one if I use from, to or cut. My scenario: I simulated densities of a plants originating from an sseed source at distance zero. Therefore the density of the plants will be highest close to zero. Is there anything I can do

Hi I have problems installing R 2.1.1. from rpm (R-base-2.2.0-1.i586.rpm) on SuSE 10. It installs fine, all dependances are OK (at least Yast does not complain) but when I try to run R, it gives the following error message: symbol lookup error: /usr/lib/libblas.so.3: undefined symbol: _gfortran_filename and R quits. blas version: 3.0-926 Any help appreciated, Rainer -- -- Rainer M.

Hi I installed R 2.2.0 from source and want to use install.packages but it doesn't work. http_proxy is set to http://proxy.sun.ac.za:3128 but it still can't connect to the repository. The mirror is available, I can connect to it via the internet. Any help welcome, Rainer -- NEW TELEPHONE NUMBER Tel: +27 - (0)72 808 2975 (w) Rainer M. Krug, Dipl. Phys. (Germany), MSc

Hi I have the following loop and would like to vectorize it. Any ideas if it is possible? Thanks, Rainer Tha Loop: for (i in 2:Result$NoSims) { ppp <- runifpoint(Result$NoPlants) K <- Kest(ppp) Result$LSim[i,] <- sqrt(K$iso / pi) - K$r CM <- (Result$LSim[i,] * Result$LSim[i,]) / abs(K$r[2] - K$r[1]) Result$SigCM[i] <- sum(CM, na.rm=TRUE) print(i) flush.console() }

Hi I am using R 2.2.1 under Linux (SuSE 10) and would like to know if it is possible to create graphs, i.e. jpeg(filename=fn) try( coplot( mor$I_Morisita ~ mor$Year | mor$RunStr, show.given=FALSE) ) dev.off() from a text console? It gives me an error message on the jpeg() command: Error in X11(..snip..) unable to start device jpeg In addition: warning message: unable to open connection to

Hi I try to run WineTools, It starts fine, but when I click any button (e.g. Basic Setup), it gives me the following message and returns to the main menu: ------------------------ ckrug@consecol-ckrug:~> wt detecting Wine version... done. Drive C: is /home/ckrug/.wine/drive_c Wine 0.9.3 wine is executed as wine Parameters are --noexit Browser is /usr/bin/firefox. WINEVER is

Dear R-help, I have a loop, which is set to take about 26 hours to run at the rate it's going - this is ridiculous and I really need your help to find a more efficient way of loading up my array gpcc.array: #My data is stored in a table format with all the data in one long column #running though every longitute, for every latitude, for every year. The #original data is sotred as

This code below produces empty XXX.png files - if I use plot(), it works and if I enter the commands png(filename=fn) bubble(positions, do.sqrt=FALSE, main=q) dev.off() manually, it works as well. I am lost - any help appreciated. The weird thing is that it worked before I made some changes... (I don't have a copy of the working version...). R version: platform i686-pc-linux-gnu arch

Hi I am using R 2.2.0 under SuSE 10 I want to use lm() to get the slope and intercept for several daatasets and store them in a database. So far so good - but how do I extract the slope and the intercept from the result from lm()? my code looks like this: lmNNDens <- lm(log(DensNN$MeanNN) ~ log(DensNN$MeanDensity)) anovaLM <- anova(lmNNDens) Results$slope[No] <- ???lmNNDens???

R 2.3.0 Linux, SuSE 10.0 Hi I have two problems with mle - probably I am using it the wrong way so please let me know. I want to fit different distributions to an observed count of seeds and in the next step use AIC or BIC to identify the best distribution. But when I run the script below (which is part of my original script), I get one error message for the first call of mle: Error in

Rainer M Krug writes: > Thorsten Wiegand used in his paper Wiegand T., and K. A. Moloney 2004. > Rings, circles and null-models for point pattern analysis in ecology. > Oikos 104: 209-229 a statistic he called O-Ring statistic which is > similar to Ripley's K, only that it uses rings instead of circles. > > http://www.oesa.ufz.de/towi/towi_programita.html#ring

Dear list members I am pleased to annonce that I have just finished a native Excel reader/writer. It's wrapped up in two packages: either "xlsReadWrite" (open source) or the slightly beefed-up "xlsReadWritePro" (shareware). Working with Excel data is now as easy as writing read.xls and write.xls. Some more details: - Infos and download: http://treetron.googlepages.com -

Hello All, We would like to call quantile() function from the R-package STATS in a Delphi program. If this is possible, could anyone provide us with an example? Thanks in advance. --Anna ----------------------------------------- Anna Belova Abt Associates Inc. 4800 Montgomery Ln, St 600 Bethesda, MD-20814 phone: 301-347-5304 fax: 301-652-7530 http://www.abtassociates.com/environment