Re: [R] Simple (?) subset problem

# Trim white space (leading and/or trailing). Includes tabs. trim <- function(str)gsub('^[[:space:]]+', '', gsub('[[:space:]]+$', '', str)) # levels(SurveyData$direction_) <- trim(levels(SurveyData$direction_)) > Date: Thu, 14 Aug 2008 15:55:38 -0700 > From: [EMAIL PROTECTED] > To: r-help@r-proj

[R] paste: multiple collapse arguments?

Dear R-helpers, I have a numeric vector, like: x <- c(1,2,3,4,5,6) I make this into a string for output to a text file, separated by \n: paste(x, collapse="\n") Is there a way to alternate the collapse argument? So between the first two elements of x, I want to separate by " ", then by "\n",

Re: [R] Plot(x,y)

What is the problem exactly? Remko - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109

[R] RGoogleDocs/RCurl through proxy

server and the port number. (Windows XP). thanks, Remko - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University

Re: [R] RGoogleDocs/RCurl through proxy

pt out of bounds If I get the Google Docs authorization by turning the proxy off; options( RCurlOptions = list(verbose = TRUE, proxy = "")) And then turn it back on, I can connect just fine. thanks, Remko --------

[R] How to set default plotting colors by treatment?

' function to find the 'treatment' variable (after all, it finds the other variables just fine). Can someone explain why this is, and how to fix this? thanks for your help Remko - Remko Duursma Post-Doctoral Fellow Centre for Pla

Re: [R] How to set default plotting colors by treatment?

> col=c("blue","red")mydfr$[treatment] Yes, but I would like to use the function for lots of other dataframes as well, so embedding 'mydfr' in the function is not the ideal solution... remko --------- Remko Duursma Post-D

Re: [R] How to set default plotting colors by treatment?

The example is reproducible! Did you see the first post? remko - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie

[R] Problem with as.POSIXct on dates object

ected? If so, can you explain why? thanks for your help, Remko - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie Univ

Re: [R] Problem with as.POSIXct on dates object

mko ----- Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia Mobile: +61 (0)422 096908 On Mon, Jul 20, 2009 at 1:41 PM, Gabor Grothendi

Re: [R] Assigning rank based on total count

fr$ POBR <- c(nstates:1)[table(testdfr$POB)][testdfr$POB] greetings, Remko - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie Univer

Re: [R] combine venn diagram and pie chart

Search on http://finzi.psych.upenn.edu/search.html gives lots of hits on "pie chart". See for example ?pie3D in the plotrix package. Or, be more specific about what you are looking for. Remko ----- Remko Duursma Post-Doctoral Fellow

Re: [R] Build a dataframe row by row?

rame(do.call("rbind", resultlis)) greetings, Remko ----- Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie Universi

Re: [R] Build a dataframe row by row?

If that is not what you mean, can you give a specific example? remko > -N > > On 8/4/09 4:57 PM, Remko Duursma wrote: >> >> Hi Noah, >> >> there are a few ways to do this. Easiest is to keep adding an element >> to a list, and then make it into a datafr

[R] SHLIB not working (Win Vista)

mko ----- Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia Mobile: +61 (0)422 096908 www.remkoduursma.

Re: [R] Data sets with usage in documentation object 'data' but not in the code

d luck, remko ----- Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia Mobile: +61 (0)422 096908 www.rem

Re: [R] R user defined language file for NotePad++

it. Remko - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia Mobile: +61 (0)422

Re: [R] Possible bug in plot.POSIXct regarding x axis

; > -- > Karl Ove Hufthammer > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, sel

[R] system() or shell() with python script

run the script fine from a command window just fine, from the same directory. Any pointers? thanks, Remko - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 D

Re: [R] spliting a word into letters

strsplit("hello","")[[1]] r --------- Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University

Re: [R] 19900501 into 1990-05-01

Something like: format(as.Date("20090501",format="%Y%m%d"), "%Y %m %d") r ----- Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 27

Re: [R] Print methods

Because print.foo is not defined if you only include the function "g" in your namespace. remko ----- Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 275

Re: [R] 3D plot, rotatable and with adjustable symbols

Try the rgl package. r - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia

Re: [R] Concave hull

See the function 'convhulln' in the 'geometry' package. It uses this algorithm : http://www.qhull.org/ remko --------- Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury

Re: [R] Concave hull

Oh right I think I did not catch that *because of* the caps. Sorry. r - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science

Re: [R] vector as data.frame element?

umbers: index <- c(0,1) and then map like this: > mylis[[which(index == 0)]] [1] "x" > mylis[[which(index == 1)]] [1] "y" "z" Is that what you want? Remko - Remko Duursma Post-Doctoral Fellow Centre for Plan

Re: [R] Is there a function to test if all the elements in a vector are unique

any(duplicated(c(1,2,2))) r - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia

Re: [R] Export R output to Word/RTF?

http://en.wikipedia.org/wiki/Software_release_life_cycle#Beta r - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie

[R] Box-counting dimension and package 'fdim'

nd what is what there. Does anyone know of an implementation of this algorithm in R, or elsewhere, or uses the fdim package successfully and can help me out? thanks, Remko --------- Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Scienc

Re: [R] processing all files with certain extension in a directory

Albert, try something like this: extfiles <- list.files(pattern=".ext") for(f in extfiles){ process.data(f) #etc } greetings, Remko ----- Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of W

Re: [R] How to delete matrix rows based on NA frequency?

tings, Remko ----- Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia Mobile: +61 (0)422 096908 www.remkoduursma.com 20

Re: [R] manipulating longitudinal data in r

You can use ave(), like this: ke$maxa <- ave(ke$a, as.factor(ke$patid), FUN=max) greetings, remko -- View this message in context: http://r.789695.n4.nabble.com/manipulating-longitudinal-data-in-r-tp4649855p4650138.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Arrange values on a timeline

(a$timestamp,nrow(b ) greetings, Remko --------- Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde

Re: [R] Storing CA Results to a Data Frame?

Please provide a reproducible example, like: library(ca) data(author) p <- ca(author) # now look at this: str(p) to find that this object of class 'ca' has lots of different results in it - it is up to you to decide which ones you make into a dataframe. cheers, Remko -- View this message i

Re: [R] matrix help

I think you want ?solve ... Remko -- View this message in context: http://r.789695.n4.nabble.com/matrix-help-tp2714378p2714896.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/ma

Re: [R] subtraction based on two groups in a dataframe

Try something like this: dfr <- read.table(textConnection("plate.id well.id Group HYB rlt1 1 P1 A1 Control SKOV3hyb 0.190 2 P1 A2 Control SKOV3hyb 0.210 3 P1 A3 Control SKOV3hyb 0.205 4 P1 A4 Control SKOV3hyb 0.206 5 P

Re: [R] fixing the dispersion parameter in glm

How about: y[y[,2] %in% x.samp[,2],] gives you the subset of y where values in the second column are restricted to your sample from x. You can then sample from this matrix, if you need to... greetings, Remko -- View this message in context: http://r.789695.n4.nabble.com/drawing-samples-bas

[R] DLL not found

not attached): [1] grid_2.12.0 tools_2.12.0 --------- Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Mobile: +61 (0)422 096908 www.remkoduursma.com ___

Re: [R] DLL not found

e error message. Remko --------- Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Mobile: +61 (0)422 096908 www.remkoduursma.com On Wed, Oct 20, 2010 at 10:00 AM, Duncan Murdoch wrote: > On

Re: [R] DLL not found

led in the latter I deleted it from the second library, but it still gives me the same error message. Remko --------- Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2

[R] R-Fortran question (multiple subroutines)

call g(x,y,z) end subroutine g(x,y,z) z = x*y end calling this from R shows that subroutine g is not called. The code compiled as executable works fine. thanks, Remko ----- Remko Duursma Research Lecturer Centre for Plants and the E

Re: [R] R-Fortran question (multiple subroutines)

Hi Berwin and all others who replied: that did the trick, thanks for your help! remko > Actually, it turns out that this example is simplified enough. :) > > I put this snippet into a file, compiled it via "R CMD SHLIB", loaded > it into R and then was very surprised about the result of ".Fort

Re: [R] calculate area between intersecting polygons

Dear Jonas, if you can write the difference in y-values between your polygons as a function, you can use integrate() to get the area between the polygons. It sounds like perhaps your x-values will not match between the polygons because they come from different sources, so you probably have to do

Re: [R] calculate area between intersecting polygons

Here is a different solution: library(gpclib) p1 <- as(poly1, "gpc.poly") p2 <- as(poly2, "gpc.poly") area.poly(p2) + area.poly(p1) - area.poly(union(p1,p2)) I.e., take areas of both polygons and subtract the union (check plot(union(p1,p2)) ) to get the area of the intersection. greetings, R

Re: [R] datetime objects

like this: format(strptime(datetime, "%m/%d/%Y %H:%M"), "%m/%d/%Y %H:%M %Z") remko -- View this message in context: http://r.789695.n4.nabble.com/datetime-objects-tp3013417p3013604.html Sent from the R help mailing list archive at Nabble.com. __ R-

Re: [R] calculate area between intersecting polygons

I don't know why I forgot that you can do this as well : area.poly(intersect(p1,p2)) ... a bit more straightforward. greetings, Remko -- View this message in context: http://r.789695.n4.nabble.com/calculate-area-between-intersecting-polygons-tp3012980p3014581.html Sent from the R help mailin

Re: [R] writing a table

This probably fixes it, although an example would be nice : write.csv(dfr, row.names=FALSE) you change the names with the ?names function before writing it to a csv file. remko -- View this message in context: http://r.789695.n4.nabble.com/writing-a-table-tp3014821p3014839.html Sent from the

Re: [R] clinical data analysis

?prop.test -- View this message in context: http://r.789695.n4.nabble.com/clinical-data-analysis-tp3014811p3014856.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] One silly question about "tapply output"

Or, tap <- tapply(rate, rating, mean) data.frame(Rating=as.factor(rownames(tap)), Mean=as.vector(tap)) remko -- View this message in context: http://r.789695.n4.nabble.com/One-silly-question-about-tapply-output-tp3015202p3015274.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] How Plot, Axis and Rect work together?

Hi Alex, is this what you want? par(xaxs="i", yaxs="i", cex.axis=0.3) plot(1, type='n', ann=FALSE, axes=FALSE,xlim=c(0,100),ylim=c(0,100)) axis(1, at=1:100) axis(2, at=1:100) abline(h=1:100,col="grey") abline(v=1:100,col="grey") X <- sample(1:100, 10) Y <- sample(1:100, 10) rect(X,Y,X+1,Y+1,col=

Re: [R] how to make multiple curves in one plot

Dear Karena, Please read the posting guide : your example cannot be reproduced because we don't know in which package the 'running' function resides. Generally speaking, use 'points' or 'lines' to add points or lines to the current device. greetings, remko -- View this message in context: h

[R] rgl.snapshot() : no longer works?

Hi all, > library(rgl) > plot3d(1,1,1) > snapshot3d("somefile.png") Error in rgl.snapshot(...) : pixmap save format not supported in this build Why does this no longer work? thanks, Remko > sessionInfo() R version 2.12.0 (2010-10-15) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_

Re: [R] rgl.snapshot() : no longer works?

Ok, thanks. I just found the new (?) function rgl.postscript() , which works better for me anyway. Remko Duncan Murdoch-2 wrote: > > On 02/11/2010 8:24 PM, Remko Duursma wrote: >> >> Hi all, >> >>> library(rgl) >>> plot3d(1,1,1) >>>

Re: [R] how to add frequencies to barplot

Hi Casper, try this: b <- barplot(table(x2),density=4, ylim=c(0,100)) tx2 <- table(x2) text(b,tx2+5, as.character(tx2)) greetings, Remko -- View this message in context: http://r.789695.n4.nabble.com/how-to-add-frequencies-to-barplot-tp3051923p3052067.html Sent from the R help mailing lis

Re: [R] Change the text size of the title in a legend of a R plot.

uot;b"),title=" ") legend("topright",c(" "," "),title="Legend",cex=0.6, bty='n', title.adj=0.15) A bit of a hack but it works If you want the title larger, it will probably not fit the box, which you can omit by setting bty='

Re: [R] Exporting interactive 3D plots with axes and labels

Does ?snapshot3d not do what you want? Or what do you mean by 'export' ? remko ----- Remko Duursma Research Lecturer Hawkesbury Institute for the Environment University of Western Sydney Hawkesbury Campus, Richmond Mobile: +61 (0)

Re: [R] complex search between dataframes

Here is another solution, using apply() and lapply(). # list of sequences corresponding to each row in data1: seqlist <- apply(data1, 1, function(x)seq(from=x[1], to=x[2])) # Find which rows of data2 are inside those sequences: rowsInSeq<- unlist(lapply(seqlist , function(x)which(data2$position

[R] Plot polygon in 3D with rgl

.tri(m, delaunayn(m))) plot3d(m[,1],m[,2],m[,3], type='l', col="black", size=2) rgl.triangles(m[tm, 1], m[tm, 2], m[tm, 3], col="green") --------- Remko Duursma Research Lecturer Centre for Plants and the Environment Uni

Re: [R] Plot polygon in 3D with rgl

Duncan, thanks for the tip! I actually saw this 2D function but had no idea how to use it in 3D. Works great. Remko - Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond

[R] Faster matrix operation?

ks, Remko ----- Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia Mobile: +61 (0)4

[R] Faster union of polygons?

soutline, leaves[[i]]) }) # about 1sec here. # Check it: plot(leavesoutline) thanks! Remko --------- Remko Duursma Research Lecturer Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biolog

Re: [R] Faster union of polygons?

Baptiste, thanks for the tip but this would give me an approximate union, and I really need a (nearly) exact one. I am also not sure how to set the alpha parameter in a non-arbitrary way. Remko - Remko Duursma Research Lecturer Centre for

Re: [R] Faster union of polygons?

duce('union', x) > leaveout <- Unionall(leaves) > > > On Tue, Jun 1, 2010 at 9:53 PM, Remko Duursma > wrote: >> >> Dear R-helpers, >> >> thanks for yesterday's speeding-up tip. Here is my next query: >> >> I have lots of polygons

Re: [R] Faster union of polygons?

PM, Remko Duursma wrote: >> Thanks for the tip - this cleans up the code a lot! >> Unfortunately, there is no gain in speed. > > Playing a little bit dirty, > > punion <- >    function(...) > { >    n <- nargs() >    if (0L == n) new("gpc.poly") >

[R] Sweave : allowing errors in R code?

Dear R-helpers, sorry if this is obvious, but I can't find this in the documentation. I am using Sweave, and have some code that does not actually work - but I want to include it anyway, including the error message that R produces. But on running Sweave() on my .Rnw file, it simply halts when i

Re: [R] Sweave : allowing errors in R code?

> Hi Remko, > > How about ?try > > Hope it helps, > > Tsjerk > Not quite, because then the code chunk in the final doc. will have this try() around it as well - not too pretty for a user manual. remko > On Aug 9, 2011 5:30 AM, "Remko Duursma" wrote: &g

Re: [R] Sweave : allowing errors in R code?

> You can modify the behaviour upon errors by setting the error option. Check > ?option and ?stop I looked at the 'error' settings in ?options, but I am not sure how that would be helpful, since I am trying to avoid Sweave from halting after an error occurs. I still want the error to be printed,

Re: [R] Sweave : allowing errors in R code?

oject.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. On 11-08-08 9:48 PM, Remko Duursma wrote: > Dear R-helpers, > > > sorry if this is obvious, but I can't find this in the documentation. > > I am using Sweave, and have some code that

Re: [R] Sweave : allowing errors in R code?

Works like a charm, thanks! remko -- View this message in context: http://r.789695.n4.nabble.com/Sweave-allowing-errors-in-R-code-tp3728790p3729802.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list http

[R] Sweave : some comments disappear

print=TRUE and keep.source=TRUE, but neither seem to affect this behavior. thanks, Remko --------- Remko Duursma Research Lecturer Hawkesbury Institute for the Environment University of Western Sydney Hawkesbury Campus, Richmond Mobile: +6

Re: [R] Sweave : some comments disappear

on of R. remko --------- Remko Duursma Research Lecturer Hawkesbury Institute for the Environment University of Western Sydney Hawkesbury Campus, Richmond Mobile: +61 (0)422 096908 www.remkoduursma.com On Tue, Sep 6, 2011 at 6:55 PM, Prof Brian Ripley wrote: > On Tue, 6 Sep 2

[R] Read text into a nested list

Dear R-list, I am looking for a way to efficiently read a text file (see example below) into a nested list. I can probably find an ugly way to do this, but I would appreciate if anyone has (ideas for) efficient solutions. It seems to me this is a bit similar to parsing XML, could tools from the

[R] gam confidence interval (package mgcv)

s(X), data=dfr, family=Gamma(link=log)) # I am interested in the relationship: # Y(treatment =="B") / Y(treatment=="A") as a function of X, with a confidence interval! Do I just do a bootstrap here? Or is there a more appropriate method? Thanks a lot for any help. greetin

Re: [R] gam confidence interval (package mgcv)

am without a problem, but I don't know how to get the CI for the ratio! thanks, remko --------- Remko Duursma Research Lecturer Hawkesbury Institute for the Environment University of Western Sydney Hawkesbury Campus, Richmond Mobile: +61 (0)422 096908 w

[R] How to 'mute' a function (like confint())

Dear R-helpers, I am using confint() within a function, and I want to turn off the message it prints: x <- rnorm(100) y <- x^1.1+rnorm(100) nlsfit <- nls(y ~ g0*x^g1, start=list(g0=1,g1=1)) > confint(nlsfit) Waiting for profiling to be done... 2.5%97.5% g0 0.4484198 1.143761 g1 1.038

Re: [R] Searching for the first two characters of a string

dfr <- data.frame(txt= c("abab","ghghg","ththt","dfdfdf"), yvar=1:4) ind <- substr(dfr$txt,1,2) dfr[ind == "ab",] greetings, Remko -- View this message in context: http://r.789695.n4.nabble.com/Searching-for-the-first-two-characters-of-a-string-tp4089862p4090042.html Sent from the R help ma

[R] Align expression and text in mtext()

Dear R-helpers, I have trouble aligning an expression with a subscript, and text, in margin text: par(mar=c(6,6,1,1)) b <- barplot(1:3) mtext(c("A","B","C"), at=b, side=1, line=1, cex=1.3) mtext(expression(italic(C)[a]~(more~text)), at=0, line=1, side=1,cex=1.3) As explained in the help files

Re: [R] strsplit question

Hi Erin, this is one way: Block <- c("5600-5699","6100-6199","9700-9799","9400-9499","8300-8399") splBlock <- strsplit(Block,"-") sapply(splBlock, "[", 1) greetings, Remko -- View this message in context: http://r.789695.n4.nabble.com/strsplit-question-tp3896847p3896850.html Sent from th

[R] deparse(substitute(x)) on an object with S3 class

ass = "testclass") # But this does : > testlist2 <- unclass(testlist) > print.testclass(testlist2) Your object name is testlist2 thanks, Remko Duursma -- View this message in context: http://r.789695.n4.nabble.com/deparse-substitute-x-on-an-object-with-S3-class-

Re: [R] deparse(substitute(x)) on an object with S3 class

Thanks for that - too bad there isn't a simple workaround! greetings Remko -- View this message in context: http://r.789695.n4.nabble.com/deparse-substitute-x-on-an-object-with-S3-class-tp4605592p4607815.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] plot: abline() - define line length

a <- 0 b <- temp } else { a <- temp[1] b <- temp[2] } } segments(x0=from,x1=to, y0=a+from*b,y1=a+to*b,...) } --------- Remko Duursma Post-Doctoral Fellow Cent

Re: [R] Goodness of fit for gamma distributions

Remko - Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109 Australia Mobile: +61 (0)422 096908 On Tue, Jan 27, 2009 at 3:38 AM, Dan31415

Re: [R] Goodness of fit for gamma distributions

ese estimated parameters: shapeest <- 1.101 rateest <- 2.49 h <- hist(mydata, breaks=50, plot=FALSE) plot(h$mids, log(h$density)) curve(log(dgamma(x, shape=shapeest, rate=rateest)), add=TRUE) #Remko ----- Remko Duursma Post-Doctoral Fellow Cen

Re: [R] Ignore text when reading data

r[setdiff(1:length(r),notread)] # make it into a dataframe (I think this can be done prettier, but whatever) z <- paste(mydata, collapse="\n") read.table(textConnection(z)) greetings Remko - Remko Duursma Post-Doctoral Fellow Centre for

Re: [R] Ignore text when reading data

Sorry, forgot this line after the textConnection bit: r <- readLines(myfile) ----- Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Scie

Re: [R] pause in function execution

Have you tried flush.console() ? r - Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW 2109

[R] How to tell if lattice is current device?

() or another lattice function? dev.cur() does not help me, it just tells me "windows 2". (Windows XP). thanks Remko ----- Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science University of Western Sydney Hawkesbury Campus Ri

Re: [R] How to tell if lattice is current device?

Thanks for the tip. I did some reading on setHook, etc.. but I am terribly confused. Do you (or someone else!) have a quick example of what the call will look like? thanks remko - Remko Duursma Post-Doctoral Fellow Centre for Plant and Food

[R] Running examples failed (but there are none).

'Maeswrap-Ex.R' file myself, there are no problems. Any pointers? thanks Remko ----- Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Scie

Re: [R] R Help

x <- list(integer(0),1,2) x[sapply(x, length) > 0] - Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University Nort

[R] interpSpline with dates?

estdfr$yvar <- rnorm(10) sp <- interpSpline(yvar ~ Date, testdfr) preddfr <- data.frame(Date=as.POSIXct("2008-08-02")) predict(sp, preddfr) thanks, Remko ----- Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science Univ

Re: [R] interpSpline with dates?

Thanks for that... the error message threw me off. r - Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North

Re: [R] reading scanned graphs

Check out http://www.datathief.org/ Remko - Remko Duursma Post-Doctoral Fellow Centre for Plant and Food Science University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie University North Ryde NSW

[R] Help on vectorizing

d, so that we have the vector r, and want to find x. # I have a (very) ugly solution: chars <- paste(r,collapse="")zeros <- strsplit(chars,"1")x <- nchar(zeros[[1]]) Thanks for your help! Remko Duursma [[alternative HTML version deleted]]

Re: [R] scanning text file to extract particular word

Hi Sit, # Read the data: dat <- scan(textConnection(" 0.293290E-05 0.117772E-05 -0.645205 *rs2282755* 0.307521E-05 0.000314 0.412997 *rs1336838* 0.484017E-05 0.218311 0.188669 *rs2660664 rs967785* 0.977861E-05 0.070474 0.294653

[R] multiple "next" in for loops?

Dear R list, I was wondering if there is an easy fix to this problem (there are workarounds, as always): Within a for loop, we can use "next" to skip to the next index, but how can we skip the next n indices? So, I would like something that looks like; for(i in 1:10){ if(somecondition)n