search for: bto

...nd especially PROC SQL, are invaluable for this. But once I get to datasets manageable for R/S-Plus, then I ship to these tools for the programming and graphics. This seems to work great. Walt Paczkowski Data Analytics Corp. -----Original Message----- >From: Rob Robinson <rob.robinson at bto.org> >Sent: Jan 17, 2008 4:31 AM >To: r-help at stat.math.ethz.ch >Subject: Re: [R] things that are difficult/impossible to do in SAS or SPSSbut simple in R > > >I wonder if those who complain about SAS as a programming environment have >discovered SAS/IML which provides a p...

....data2) Question 2: Any hints on what might be causing this? Am I fitting the wrong (or too complicated a model)? Btw if it is relevant I am using mgcv_1.3-19 and R 2.3.1. Many thanks for any hints on where I should start digging Cheers rob *** Want to know about Britain's birds? Try www.bto.org/birdfacts *** Dr Rob Robinson, Senior Population Biologist British Trust for Ornithology, The Nunnery, Thetford, Norfolk, IP24 2PU Ph: +44 (0)1842 750050 E: rob.robinson at bto.org Fx: +44 (0)1842 750030 W: http://www.bto.org eSafe scanned this email for viruses, vandals and ma...

...preciated. I'm getting quite stressed as I have 10 files with approx 100,000 records in each to analyse so manipulating data becomes a pain. Here is the code below, not sure it makes much sense without seeing the dataset:- chaff<-read.spss("/Users/Kat/Desktop/papers in progress/btopaper/ edited BTO data/fatnewchaff.sav", to.data.frame=TRUE) chafffat<-read.spss("/Users/Kat/Desktop/papers in progress/btopaper/ edited BTO data/fatmethods.sav") chaffmerge2<-merge(chaff, chafffat, by.x=c("RINGNO", "FAT", "FATMTD"), by.y=c(&quo...

Hi there I have an interesting problem: My csv file is of array dimensions [12,50], but it was saved the wrong way: there should be only 11 colums. What happens now if I read it into R is that the whole data set is shifted ( in the first row, the last column contains already the first value of the supposed second row and so on...) how can I tell R to switch after 11 read values to the next row,

...ixmap and read.jpeg() from rgl. the pnm example provided works fine x <- read.pnm(system.file("pictures/logo.ppm", package="pixmap")[1]) plot(x) #draws an R logo, nice :-) But my code returns an error x<-read.pnm("C:/Documents and Settings/simonp/My Documents/Simon BTO/RELU/GIS data/ten km areas in analysis.bmp") Error in pm.readmagicnumber(con) : Not a PNM format file I have tried changing the dpi when I export the image but still no cigar. I have dredged the help forum for related files and read the help files. My Os is windows XP and i'm running 2.8....

R-users, Is there a function for ceiling to the nearest ten? a <- 1:10*4 a [1] 4 8 12 16 20 24 28 32 36 40 The resulting vector should look like this ("ceiling to the nearest ten") [1] 10 10 20 20 20 30 30 40 40 40 Thanks in advance Lauri

...25, June 2001, especially the figures showing the partial predictions of Mackerel egg densities. Using those graphs as an example, I would like to extract the partial predictions for e.g. "s(b.depth)", given the inclusion of the other variables. -- Staffan Roos, PhD Research Ecologist BTO Scotland

I am interested to take a course in R. Someone know of some course in europe for the first semester of the next year?. Happy Year Thank you, Xavier [[alternative HTML version deleted]]

Dear users, a very simple question: Given two vectors x and y x<-as.character(c("A","B","C","D","E","F")) y<-as.factor(c("1","2","3","4","5","6")) i want to combine them into a single vector z as A1, B2, C3 and so on. z<-x*y is not working, i tried several others

I am using bagging to perform Bagged Regression Trees on count data (bird abundance in Britain and Ireland, in relation to climate and land cover variables). Predictions from the final model are visually believable but I would really like a diagnostic equivalent to classification success that can be used to decide if a model is adequate. Whereas with classification data an error rate is returned,

Hi all. I'm migrating from a small OpenLDAP setup and currently have users' password hashes in {SSHA} and {CRYPT}$5$.16s format. Can I just ldbedit or ldbmodify user's supplementalCredentials fields in /var/lib/samba/private/sam.ldb.d/DC%3DAD%2CDC%3DEXAMPLE%2CDC%3DCOM.ldb to migrate passwords? Provided that I could get the data structure right. (Documentations about

Hi all, Suppose I had the below example where a survey was carried out recording the number of each type of pet in each house count<-c(2,1,2,1,2,3,4) house<-c("house1","house1","house2","house3","house4","house4","house4")

Hi all, This is partly a statistical question as well as a question about R, but I am stumped! I have count data from various sites across years. (Not all of the sites in the study appear in all years). Each site has its own habitat score "habitat" that remains constant across all years. I want to know if counts declined faster on sites with high "habitat" scores. I can

hi all, I often have a data frame like this example data.frame(sq=c(1,1,1,2,2,3,3,3,3),area=c(1,2,3,1,2,3,1,2,3),habitat=c("garden","garden","pond","field","garden","river","garden","field","field")) for each "sq" I have multiple "habitat"s each with an associated "area". I

Hi all, This could be naivety/stupidity on my part rather than a problem with model output, but here goes.... I have fitted a fairly simple model m1<-glm(count~siteall+yrs+yrs:district,family=quasipoisson,weights=weight,data=m[x[[i]],]) I want to know if yrs (a continuous variable) has a significant unique effect in the model, so I fit a simplified model with the main effect ommitted...

Hi all I have produced a biplot for a PCA (see attached pdf) that I ran however the names of the variables which are placed at the end of the arrows overlap and are thus unreadable. Similarly some of the numbered points overlap. I was wondering if there was a way to edit the biplot to move the label names and if not what the best alternative is. Thanks Anna pca<-biodata[,3:10]

...21209<D6>^@^_<<D7>^@^B^_><D8>^@^A<B9>3:42:47 -06<A8>^Ac(from -^B^A<BA>^@Ohost<AC>^@^\<86>/Submit)<A8>^@!cd<A8>^@N7; 4^B <8C>^@^_DZ^B^D &^@<F2>^KMessage-Id: <200411220942.^Y^A^Dq^@^O<F1>^B <FF>^BTo: Renaud AllardB^A^E<B2> Then, I randomly see errors like this in maillog: truncated lz4 chunk at 555395 Given I am still able to read the mails using imap, I think there might be a way to recover them to a proper format. However, using dsync to make a mirror throws a lot of errors like...

...are. .Since I am in India , the cost differences can cause almost a digital divide in who can and who cant use sophisticated software. Thanks to the Angels here....Yes we Can R............... Regards, Ajay www.decisionstats.com On Wed, Jan 7, 2009 at 8:03 PM, Simon Pickett <simon.pickett@bto.org> wrote: > I would like to add that I would have spent many more years doing my PhD if > it wasnt for R! all data management, statistics and graphics were conducted > using it. This was the direction my university and many more research > institutes appear to be heading. > &gt...

Dear list members. I have problems with the usage of cv.glm from the boot package. Here are some parts of the script I wanted to use: data <- read.table("selected_2D.csv", header=TRUE, sep=",") ? glm.fitted <- glm("ydata$ y ~ 1 + density + vsurf_ID6 + vsurf_S ", data=data) error <- cv.glm(data=data, glm.fitted, K=6) ydata$y is a separate data set, where

Is R an appropriate tool for data manipulation and data reshaping and data organizing? I think so but someone who recently joined our group thinks not. The new recruit believes that python or another language is a far better tool for developing data manipulation scripts that can be then used by several members of our research group. Her assessment is that R is useful only when it comes to data