search for: andydolman

...ations 1,2,7, and 10 to be left out of the first xval sample, 3,4, and 9 out of the second, etc. Terry Therneau I can't see how this string of values, c(1,1,2,2,3,3,1,3,2,1), codes for observations 1,2,7,10 being left out of the 1st and so on. Can anyone fill me in please? Thanks, andydolman at gmail.com

Hi, I have an experiment with 2 independant factors which I have been trying to analyse in R. The problem is that there are several data points recorded on the same animal. However, no combination of treatments is repeated on the same animal. All possible combinations of treatments are done in a random order with as many points as possible being done on 1 animal before moving onto the next. The

I got a distribution function and a empirical distribution function. How do I make to Kolmogorov-Smirnov test in R. Lets call the empirical distribution function >Fn on [0,1] and the distribution function >F on [0,1] ks.test( ) thanks for the help -- View this message in context: http://www.nabble.com/Kolmogorov-Smirnov-test-tp23296096p23296096.html Sent

Dear all, I have an issue trying to install new packages (have tried with RODBC_1.3-1, gplots_2.6.1, gtools_2.7.4 packages) and get the same error message : "Error in library(gplots) : there is no package called 'gplots'" Only clue I have found so far on the Web is related to Perl (Perl modules are installed on my computer, but which one is related to gplots if any ?); no

...inference beyond this sample of regions. My problem is getting a test for the fixed effects, see http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests for background on why I want to use mcmcsamp. Thank you for your attention, Andrew Dolman - Australian Institute of Marine Science. -- andydolman@gmail.com a.dolman@aims.gov.au # required packages library('lme4') # some dummy data counts <- c(rpois(75, 2), rpois(25,4)) # number of fish in a sample [dependent variable - poisson distributed] visit <- c(rep(0, 50), rep(1,50)) # sample visit number, 0 or 1 to indicate pre and...

Hello I am trying to plot two lines of best fit on a scatter plot. I'm plotting number of visits (y) against personality (x), half the plots are blue for patchchoice one and the other plots are red for patchchoice two (these options come from one factor, see below) . So i need a line for the patch ones and another line for the patch twos. How would i go about doing this? I know i

Hi, I am trying to create a function for a goodness-of-fit test for the Pareto Distribution for some loss data that I have. So far I have the following: function(X=OTOL) { n <- length(X)-1 #calculated the number of values (extra as 0 included) i <- 2:640 #values of i j <- 1:639 #values of i-1 Y <- (n-j+1)*((X[i])-(X[j])) #First part of GoF model Y } Where OTOL is the ordered loss

hy guys i have one question :) i have two vectors markets and price market <- c(1, 5, 7, 9, 9, 6, 5, 4, 4, 3, 1, 2, 1) price <- c(100, 20, 30, 10, 50, 23, 23, 33, 96, 6, 4, 38, 96) i would like sum prices: market 1: (100+4+96), market 2: (38),..., market 9: (10+50) ao i would like get this result: (200, 38, ..., 60) and i don't wanna use while, for loops... is there any other function

Hey guys, I have a problem: I created a silly for loop without saving the results on each step. After a while I realised that it will take days to finish the loop until I get the results. Is there a way to get the data R saves in working memory or in a temporary file while runing the loop? So that stoping the loop will not result in complete data loss? Thank you very much! Greetings,

I want to make a log-log plot with a regression line, but I can't figure out what I'm doing wrong. What I'm trying is: plot(mass,area, log="xy", pch=as.numeric(food)) abline(lm(mass~area)) or plot(mass,area, log="xy", pch=as.numeric(food)) islands$logmass <- log(mass) islands$logarea <- log(area) attach(islands) abline(lm(logmass~logarea)) But that does

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

...s ceart duit dul i dteagmh?il leis an seolt?ir l?ithreach agus le mailminder[ag]justice.ie chomh maith. *********************************************************************************** ------------------------------ Message: 42 Date: Wed, 27 May 2009 15:46:31 +0200 From: Andrew Dolman <andydolman@gmail.com> Subject: Re: [R] How to write a loop? Cc: r-help@r-project.org Message-ID:     <951234ac0905270646y297cf12ai748469bf15eed1dc@mail.gmail.com> Content-Type: text/plain Try lapply(ONS, fft) and take a look here http://cran.r-project.org/doc/manuals/R-intro.html for the basics of...

Dear list members, Can anyone please point to an example of how to use glht(multcomp) with lmer objects? I am trying: summary(glht(lmerObject, linfct = mcp(x = "Tukey"))) as I would for a glm object, but with no luck. Thank you, Andy. [[alternative HTML version deleted]]