[R] Splitting dataset for Tuning Parameter with Cross Validation

2009-07-12 Thread Tim
Hi, My question might be a little general. I have a number of values to select for the complexity parameters in some classifier, e.g. the C and gamma in SVM with RBF kernel. The selection is based on which values give the smallest cross validation error. I wonder if the randomized splitting of

Re: [R] Splitting dataset for Tuning Parameter with Cross Validation

2009-07-13 Thread Tim
g smallest validation error. It appears to me the splitting is once for all choices of the complexity parameter. Thanks! --- On Sun, 7/12/09, Tim wrote: > From: Tim > Subject: [R] Splitting dataset for Tuning Parameter with Cross Validation > To: r-h...@stat.math.ethz.ch > Date: Sunda

[R] Get (feature, threshold) from Output of rpart() for Stump Tree

2009-05-08 Thread Tim
f the stump tree , and 13.5, which I believe to be the threshold on the feature.  If I am able to get these two out, then I will be able to further process them or write them into a file. Any hint? Thanks and regards! -Tim [[alternative HTML version deleted]]

Re: [R] Get (feature, threshold) from Output of rpart() for Stump Tree

2009-05-08 Thread Tim
feature ID. Is this the correct way to get the row label of fit$splits[1,]? Regards, - Tim --- On Fri, 5/8/09, Terry Therneau wrote: From: Terry Therneau Subject: Re: Get (feature, threshold) from Output of rpart() for Stump Tree To: "Tim" Cc: r-help@r-project.org Date: Friday, May

[R] ROCR: auc and logarithm plot

2009-05-12 Thread Tim
to have a more distinguishable plot of the ROCs. So is it possible to have a logarithm FP axis which might probably separate them well? Or zoom in the part close to the leftup corner of ROC plot? Or any other ways to make the ROCs more separate? Thanks and regards! --Tim [[alterna

Re: [R] ROCR: auc and logarithm plot

2009-05-12 Thread Tim
quot;, the color of bars are still black, the default one, instead of red "2". --Tim --- On Tue, 5/12/09, Tobias Sing wrote: From: Tobias Sing Subject: Re: [R] ROCR: auc and logarithm plot To: timlee...@yahoo.com, r-help@r-project.org Date: Tuesday, May 12, 2009, 5:54 AM > 1. I hav

[R] #Keeping row names when using as.data.frame.matrix

2013-05-17 Thread Tim
#question I have the following data set: Date<-c("9/7/2010","9/7/2010","9/7/2010","9/7/2010","9/7/2010","9/7/2010","9/8/2010") EstimatedQuantity<-c(3535,2772,3279,3411,3484,3274,3305) ScowNo<-c("4001","3002","4002","BR 8","4002","BR 8","4001") dataset<- data.frame(EstimatedQuantity,Date,ScowNo)

[R] need help reshaping table using aggregate

2012-06-20 Thread Tim
I am trying to learn how to reshape my data set. I am new to R, so please bear with me. Basically, I have the following data set: site<-c("A","A","B","B") bug<-c("spider","grasshopper","ladybug","stinkbug") count<-c(2,4,6,8) myf <- data.frame(site, bug, count) myf site bug count 1

[R] Installation: not creating necessary directories

2008-10-29 Thread Tim
file or directory make[5]: *** [Rdfiles] Error 127 make[5]: Leaving directory `/home/USERACCOUNT/mybin/R-2.8.0/src/library' file ../../library/base/R-ex cannot be opened at ../../share/perl/massage-Examples.pl line 136. ... Is there any advice on what might be happening or on what I

[R] CRAN I-MR / Xbar-R / Xbar-S control chart package ?

2017-07-08 Thread Tim Smith
Hi, I've had a quick look through the package list, and unless I've missed something, I can't seem to find anything that will do I-MR / Xbar-R / Xbar-S control charts ? Assuming there is something out there, can anyone point me in the right direction ?

Re: [R] CRAN I-MR / Xbar-R / Xbar-S control chart package ?

2017-07-09 Thread Tim Smith
Interesting, thanks for that. I came accross qcc but my quick scan of the docs is that it only did xbars but maybe I need to re-read the docs, I guess it does the individual plot versions (I-MR) too. Tim On 8 July 2017 at 20:53, Rui Barradas wrote: > Hello, > > I have no experience

[R] Interesting behavior of lm() with small, problematic data sets

2017-09-05 Thread Glover, Tim
perfect fit, but a p-value of 0.3112 doesn’t seem reasonable. As a side note, the various r^2 values seem odd too. Tim Glover Senior Scientist II (Geochemistry, Statistics), Americas - Environment & Infrastructure, Amec Foster Wheeler 271 Mill Road, Chelmsford, Massachusetts, USA 01824-4105 T

[R] rtmvnorm {tmvtnorm} seems broken

2014-12-09 Thread Tim Benham
General linear constraints don't seem to work. I get an error message if I have more constraint equations than variables. E.g. executing the following code print(R.version) library('tmvtnorm') cat('tmvtnorm version ') print(packageVersion('tmvtnorm')) ## Let's constrain our sample to the dwarfed h

[R] read.table with missing data and consecutive delimiters

2015-02-11 Thread Tim Victor
34 22 NA NA5 3 NA NA NA6 What am I missing? Thanks, Tim [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEAS

[R] problem with labeling plots, possibly in font defaults

2014-08-07 Thread Tim Blass
Hello, I am using R 3.1.1 on a (four year old) MacBook, running OSX 10.9.4. I just tried making and labeling a plot as follows: > x<-rnorm(10) > y<-rnorm(10) > plot(x,y) > title(main="random points") which produces a scatter plot of the random points, but without the title and without any numbe

[R] How can I overwrite a method in R?

2014-10-08 Thread Tim Hesterberg
How can I create an improved version of a method in R, and have it be used? Short version: I think plot.histogram has a bug, and I'd like to try a version with a fix. But when I call hist(), my fixed version doesn't get used. Long version: hist() calls plot() which calls plot.histogram() which fa

Re: [R] How can I overwrite a method in R?

2014-10-09 Thread Tim Hesterberg
thz.ch > > Subject: Re: [R] How can I overwrite a method in R? > > > > This is usually ill-advised, but I think it's the right solution for > > your problem: > > > > assignInNamespace("plot.histogram", function(...) plot(1:10), "graphics") > > hi

[R] big datasets for R

2014-10-28 Thread Tim Hoolihan
Another source of large datasets is the Public Data Sets on AWS http://aws.amazon.com/public-data-sets/ Tim Hoolihan @thoolihan http://linkedin.com/in/timhoolihan On Oct 28, 2014, at 7:00 AM, r-help-requ...@r-project.org wrote: > -- > > Message: 2 > D

[R] Download showing as exploit

2015-09-21 Thread Tim Kingston
Hi , I work for the NHS, and our IT service has been unable to download as its anti-virus software says it contains an exploit. Is this normal? Is there a way around this? Kind regards, Tim Kingston Sent from my HTC [[alternative HTML version deleted

[R] Problem installing/loading packages from Africa

2016-04-22 Thread Tim Werwie
insight would be very helpful. See error messages below. Thanks -- Tim - When attempting to install, the print-out I get in the Console in RStudio is any length of: > install.packages("swirl") % Total% Received % Xferd Average Speed TimeTime Time Current

[R] [R-pkgs] CRAN updates: rpg and odeintr

2017-02-20 Thread Tim Keitt
rpg is a package for working with postgresql: https://github.com/thk686/rpg odeintr is a package for integrating differential equations: https://github.com/thk686/odeintr Cheers, THK http://www.keittlab.org/ [[alternative HTML version deleted]] __

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

2008-07-27 Thread Tim Hesterberg
rap is used for inference (bias, standard error, confidence intervals), not improving on ThetaHat. Tim Hesterberg >Hi Doran, > >Maybe I am wrong, but I think bootstrap is a general resampling method which >can be used for different purposes...Usually it works well when you do not >have a

[R] Vista problem -- can't type commands at prompt

2008-08-07 Thread Tim Calkins
thanks for any help you might be able to provide. tim -- Tim Calkins 0406 753 997 __ 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 pr

[R] Question about type conversion in read.table with columns that contain "+" and "-" in R > 2.7

2008-09-01 Thread Tim Beissbarth
or "-" and would like to convert these to characters or factors, without having to specify the individual column types manually. Is there any way to still do so in a new version of R? Many thanks and best wishes, Tim __ R-help@r-project.o

[R] Counting the number of non-NA values per day

2009-08-11 Thread Tim Chatterton
rements for each day to check that the mean and max are likely to be valid (for example I need at least 18 hourly measurements to calculate a valid daily mean) Try as I might I have not found a simple way of doing this. Can anybody help please? Many thanks, Tim. -- __

[R] REMOVE ME

2009-08-12 Thread Tim Paysen
This mailing list is too intrusive.  Remove my name. [[alternative HTML version deleted]] __ 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.ht

[R] Finding minimum of time subset

2009-08-13 Thread Tim Clark
quot; "12:30:01" "12:45:01" "13:00:00" "13:15:02" This gives me the data taken at quarter hour intervals (removes 12:20:00) but I am still left with multiple values at the quarter hours. I would like to obtain: "12:00:00" "12:15:05&quo

Re: [R] Finding minimum of time subset

2009-08-14 Thread Tim Clark
ontex, or what the number one at the end does. Thanks for you help, Tim Tim Clark Department of Zoology University of Hawaii --- On Fri, 8/14/09, jim holtman wrote: > From: jim holtman > Subject: Re: [R] Finding minimum of time subset > To: "Tim Clark" > Cc: r-help@r

Re: [R] Finding minimum of time subset

2009-08-14 Thread Tim Clark
Jim, Got it! Thanks for the explanation and the example. Always nice to learn new tricks on R. Aloha, Tim Tim Clark Department of Zoology University of Hawaii --- On Fri, 8/14/09, jim holtman wrote: > From: jim holtman > Subject: Re: [R] Finding minimum of time subset >

Re: [R] Finding minimum of time subset

2009-08-14 Thread Tim Clark
Thanks for everyones help and for the alternate ways of doing this. I am always amazed at how many solutions this list comes up with for things I get stuck on! It really helps us non-programmers learn R! Aloha, Tim Tim Clark Department of Zoology University of Hawaii --- On Fri, 8/14/09

[R] Trying to rename spatial pts data frame slot that isn't a slot()

2009-08-30 Thread Tim Clark
slotNames(). > slotNames(mymcp) [1] "data""polygons""plotOrder" "bbox" "proj4string" What am I missing here? Is Slot "ID" not a slot? Can I export the ID's with the shapefiles to qGIS? Can I rena

[R] Zoomable graphs with multiple plots

2009-09-02 Thread Tim Shephard
tell, neither can handle the task. Does anyone know anything else that might work? Or generally know different? Cheers, Tim. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://

Re: [R] Zoomable graphs with multiple plots

2009-09-03 Thread Tim Shephard
This is what I've done. Just capture two identifies and then replot. If the identifies are right left, I zoom out. Works quite well. Still, can't wait for iplot xtreme. Cheers, Tim. On Thu, Sep 3, 2009 at 5:03 AM, Jim Lemon wrote: > Tim Shephard wrote: >> >&

Re: [R] Zoomable graphs with multiple plots

2009-09-03 Thread Tim Shephard
On Thu, Sep 3, 2009 at 1:25 PM, Jim Porzak wrote: > Tim, > > I've had success (& user acceptance) simply plotting to a .pdf & > passing zoom functionality to Acrobat, or whatever. > > Worked especially well with large US map with a lot of fine print annotation. >

[R] Amazon SimpleDB and R

2009-09-19 Thread Tim Shephard
As far as I know there isn't anything available for this, but I thought I'd check before working up something of my own. Is there a way to query Amazon SimpleDB and import the data results directly into R? Cheers, Tim. __ R-help@r-project.o

[R] ROCR.plot methods, cross validation averaging

2009-09-23 Thread Tim Howard
st entry with only two values p...@x.values[[1]] <- c(0,1) p...@y.values[[1]] <- c(0,1) p...@alpha.values[[1]] <- c(Inf,0) plot(perf, avg="threshold") ##output results in an error with this message # Error in if (from == to) rep.int(from, length.out) else as.vector(c(from, : #

[R] Find each time a value changes

2010-02-10 Thread Tim Clark
,y1,y1,y1,y2,y1,y2,y1,y2) xy<-cbind(x,y) I would like to know each time the numbers change. Correct answer should be: x=1,11,21,51,61,71,81,91 I would appreciate any help or suggestions. It seems like it should be simple but I’m stuck! Thanks, Tim Tim Clark Department of Zoology Unive

Re: [R] Find each time a value changes

2010-02-10 Thread Tim Clark
Thanks everyone! I would have been banging my head around for quite a while and still wouldn't have come up with either solution. The function rle() is a good one to know! Aloha, Tim Tim Clark Department of Zoology University of Hawaii --- On Wed, 2/10/10, Ben Tupper wrote: &

Re: [R] Find each time a value changes

2010-02-11 Thread Tim Clark
# SEE WHAT RLE GIVES test <- rle(xy[,2]) print(str(test) # USE JIM"S TRICK OF CUMULATIVE SUMMING # TO GET THE LOCATIONS result <- cumsum(c(1,rle(xy[,2])$lengths)) Tim Clark Department of Zoology University of Hawaii --- On Wed, 2/10/10, Ben Tupper wrote: > From:

[R] Formatting question for separate polygons

2010-02-11 Thread Tim Clark
format (x1)? #Sample data x<-c(1,2,3,4,5,6) y<-c(1,2,2,1) #I need to format the data like this x1<-c(1,1,2,2,NA,3,3,4,4,NA,5,5,6,6,NA) y1<-rep(c(1,2,2,1,NA),length(x)/2) #Final plot plot(c(1,6), 1:2, type="n") polygon(x1,y1,density=c(40)) Thanks, Tim Tim Clark Department

Re: [R] Formatting question for separate polygons

2010-02-12 Thread Tim Clark
Thanks Uwe! And Peter for the correction. I would never have come up with that! Tim Tim Clark Department of Zoology University of Hawaii --- On Fri, 2/12/10, Uwe Ligges wrote: > From: Uwe Ligges > Subject: Re: [R] Formatting question for separate polygons > To: "Peter

[R] Non-monotonic spline using splinefun(method = "monoH.FC")

2010-02-15 Thread Tim Heaton
Hi, In my version of R, the stats package splinefun code for fitting a Fritsch and Carlson monotonic spline does not appear to guarantee a monotonic result. If two adjoining sections both have over/undershoot the way the resulting adjustment of alpha and beta is performed can give modified values

Re: [R] Non-monotonic spline using splinefun(method = "monoH.FC")

2010-02-15 Thread Tim Heaton
direction and moving one knot would no longer guarantee that the curve changes in only a finite region) or possibly shrink the coefficients twice (which would create a flatter spline than necessary but would give a finite effect and the same curve if the data were entered in the opposite dir

[R] r help date format changes with c() vs. rbind()

2010-02-19 Thread Tim Clark
n't c give me the correct answer? Thanks, Tim > a<-as.POSIXct("2000-01-01 12:00:00") > a [1] "2000-01-01 12:00:00 HST" > b<-rbind(a,a) > b [,1] a 946764000 a 946764000 > c<-as.POSIXct(b,origin="1970-01-01") > c [1] "2000

Re: [R] r help date format changes with c() vs. rbind()

2010-02-19 Thread Tim Clark
Glad to know it isn't just me! I couldn't use Phil's data.frame method since my real problem went from a POSIXct object to a large matrix where I used rbind and then back to POSIXct. Jim's function worked great on converting the final product back to the proper date. Than

Re: [R] (quite possibly OT) Re: how to make R running on a Linux server display a plot on a Windows machine

2010-02-23 Thread Glover, Tim
One additional option on the Linux-on-the-local-machine option described below. Consider running Microsoft's free virtual machine with a copy of Linux in there. Now you have the advantages of dual-boot with BOTH operating systems available with a toggle key sequence. Tim Glover S

Re: [R] Plotting 15 million points

2010-02-25 Thread Glover, Tim
Have you considered taking a random subset and plotting that? I'd bet you can get a really impression of the distribution with a few hundred thousand points at most. Tim Glover Senior Environmental Scientist - Geochemistry Geoscience Department Atlanta Area MACTEC Engineering and Consu

[R] reading data from web data sources

2010-02-27 Thread Tim Coote
o that it is easier to consume by machines, while being easy to read for humans. The early records are a serious machine parsing challenge as they are tiff images of old notebooks ;-) tia Tim Tim Coote t...@coote.org vincit veritas __ R-he

Re: [R] reading data from web data sources

2010-02-27 Thread Tim Coote
documentation for. What is it that's removing the -999 and -888 values in this code -they seem to be gone, but I cannot see why. Phil's reads in the data, but interleaves rows with just a year and all other values as NA. Tim On 27 Feb 2010, at 17:33, Gabor Grothendieck wrote:

Re: [R] Best Hardware & OS For Large Data Sets

2010-02-27 Thread Tim Coote
self-contained, reproducible code. Tim Coote t...@coote.org +44 (0)7866 479 760 __ 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 com

Re: [R] turn character string into unevaluated R object

2010-03-02 Thread Tim Calkins
3.2 1.3 0.2 setosa > 4 4.6 3.1 1.5 0.2 setosa > 5 5.0 3.6 1.4 0.2 setosa > 6 5.4 3.9 1.7 0.4 setosa > > Liviu > > ____

[R] ggplot2 rose diagram

2010-03-09 Thread Tim Howard
size of each slice, which makes for a different presentation than ggplot2. Any suggestions on how to use this calculation method in the ggplot framework? Thanks in advance for your help. Tim Howard New York Natural Heritage Program __ R-help@r-projec

[R] Question on passing in parameter to Cox hazard

2010-03-09 Thread Tim Smith
Hi, I wanted to do the cox model using a matrix. The following lines illustrate what I want to do: dat <- matrix(rnorm(30), ncol=3,dimnames = list(1:10,letters[1:3])) Survival <- rexp(10) Status <- ifelse(runif(10) < .7, 1, 0) mat <- as.data.frame(cbind(dat,Surviv

Re: [R] ggplot2 rose diagram

2010-03-10 Thread Tim Howard
geom_bar(width=1) + opts(axis.text.x = theme_text(size = 18)) smp2 + coord_polar() Cheers, Tim >>> Tim Howard 3/9/2010 9:25 AM >>> Dear R gurus - consider this plot: library(ggplot2) dat <- sample(1:8,100,replace=TRUE) smp <- ggp

Re: [R] ggplot2 rose diagram

2010-03-10 Thread Tim Howard
.g. sqrt(countInBin/totalCount) . My goal is for area of the slice to correlate with number of records in each bin (rather than area expanding at a faster rate). Thanks, Tim >>> hadley wickham 3/10/2010 9:14 AM >>> For Q2 you can use opts(legend.position = c(0.9, 0.9)). For Q3, y

Re: [R] ggplot2 rose diagram

2010-03-10 Thread Tim Howard
Got it. Thanks again so much for the help. Best, Tim >>> hadley wickham 3/10/2010 2:46 PM >>> > By Q2 I was trying to refer to the Y-axis labels. For the polar plot, the > Y-axis labels reside left of the panel. I was looking for a way to get the > Y-axis lab

[R] heatmap.2 - ColSideColors question

2010-03-16 Thread Tim Smith
Hi, I wanted to make more than one side color bar. For example, I can make one side color bar (col1) with the following code: --- library(gplots) mat <- matrix(sample(1:100,40),nrow=5) class1 <- c(rep(0,4),rep(1,4)) col1 <- ifelse(class1 == 0,"blue","red") class2

[R] abline on heatmap

2010-03-20 Thread Tim Smith
Hi, Is there a way I can draw an abline on a heatmap? I try the abline function, but don't get the line. My sample code is: mat <- matrix(sample(1:100,40),nrow=5) heatmap(mat,col=greenred(75),trace="none", dendrogram = "column",labCol = NULL) abline(h=5,v=4) thanks! [[alternati

[R] Calculating distance between spatial points

2009-06-25 Thread Tim Clark
iff<-diff(SP.utm$Lat) d=(long.diff^2+lat.diff^2)^.5 speed=d/date.diff Aloha, Tim Tim Clark Department of Zoology University of Hawaii __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the po

[R] How should I denormalise a data frame list of lists column?

2009-07-01 Thread Tim Slidel
t I'd want there is something like: > f1_denorm c1c2 Split SplitDenorm 1 0 A,B,C A, B, C A 2 0 A,B,C A, B, C B 3 0 A,B,C A, B, C C 4 1 A,EA, E A 5 1 A,EA, E E 6 2 F,G,H F, G, H F 7 2 F,G,H F, G, H G 8 2 F,G,H F, G, H H I thought perhap

[R] How should I denormalise a data frame list of lists column?

2009-07-01 Thread Tim Slidel
#x27;t find a way to do this, the result I'd want there is something like: > f1_denorm c1c2 Split SplitDenorm 1 0 A,B,C A, B, C A 2 0 A,B,C A, B, C B 3 0 A,B,C A, B, C C 4 1 A,EA, E A 5 1 A,EA, E E 6 2 F,G,H F, G, H F 7 2 F,G,H F, G

[R] Getting the month out of my date as a number not characters

2009-07-02 Thread Tim Chatterton
quot;],"%H") and hf$month <- as.POSIXct(strptime(hf$date, format = "%m")) but they are still coming out as characters. Any ideas please? Thanks, Tim. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] Getting the month out of my date as a number not characters

2009-07-02 Thread Tim Chatterton
;hour"],"%H") and hf$month <- as.POSIXct(strptime(hf$date, format = "%m")) but they are still coming out as characters. Any ideas please? Thanks, Tim. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/lis

Re: [R] Collinearity in Linear Multiple Regression

2009-07-21 Thread Tim Paysen
Actually, the CI index and VIF are just a start.  It is best to look at what they call a matrix of "variance proportions" (found in SAS and a few other places...)--which hardly anyone understands (including the SAS folks).  It is a matrix of estimates of what the variences of the regression coef

[R] Duplicated date values aren't duplicates

2009-07-23 Thread Tim Clark
(07/30/02 00:00:00) (08/06/03 17:45:00) Is there a reason different dates would come up as duplicate values? I would prefer not to have to delete them if I don't have to. Any suggestions on how to get R to realize they are different? Thanks, Tim Tim Clark Department of Zoolo

Re: [R] Duplicated date values aren't duplicates

2009-07-24 Thread Tim Clark
Don and Jim, Thanks! I got it! Duplicated is only returning one of the two duplicated dates (the second date). It all makes sense now! Tim Tim Clark Department of Zoology University of Hawaii --- On Fri, 7/24/09, Don MacQueen wrote: > From: Don MacQueen > Subject: Re: [R] Dupl

[R] Maximizing values in subsetted dataframe

2009-07-28 Thread Tim Clark
choose between these and obtain the one with the fewest missing values. In my real dataset I have multiple variables collected every second and I would like to subsample it every 5, 10, and 15 minutes. I appreciate your help. Tim Tim Clark Department of Zoology Univer

[R] creating MS Access query objects using RODBC

2009-08-05 Thread Tim Calkins
write sql to make a table or that I could run the sql, extract it to R, and then use sqlSave to save the dataframe as a table in the db. thanks in advance, -- Tim Calkins 0406 753 997 [[alternative HTML version deleted]] __ R-help@r-project.o

[R] xts off by one confusion or error

2010-04-08 Thread Tim Coote
nux-gnu system i386, linux-gnu status major 2 minor 10.0 year 2009 month 10 day26 svn rev50208 language R version.string R version 2.10.0 (2009-10-26) And from installed.packages (): xtsNA NA "GPL-3"

Re: [R] xts off by one confusion or error

2010-04-08 Thread Tim Coote
x Jan 2010 6.08 > zoo (aggregate (vv, as.yearmon(index(vv)), mean)) x Feb 2010 6.08 On 8 Apr 2010, at 15:53, Tim Coote wrote: Hullo I may have missed something blindingly obvious here. I'm using xts to handle some timeseries data. I've got daily measurements for 100 y

Re: [R] xts off by one confusion or error

2010-04-08 Thread Tim Coote
x Jan 2010 6.08 > zoo (aggregate (vv, as.yearmon(index(vv)), mean)) x Feb 2010 6.08 On 8 Apr 2010, at 15:53, Tim Coote wrote: On 8 Apr 2010, at 15:53, Tim Coote wrote: Hullo I may have missed something blindingly obvious here. I'm using xts to handle some timeseries d

[R] Best subset of models for glm.nb()

2010-04-20 Thread Tim Clark
r a valid alternative to narrowing down models when you have a lot of data and a large number of variables your interested in. Any thoughts on either of these methods? Or should I be doing something else? Thanks for your help, Tim Tim Clark Department of Zoology University

[R] Periodic regression - lunar percent cover

2010-04-29 Thread Tim Clark
, 0.88, 0.62, 0.19, 0.21, 0.01, 0.70, 1.00, 0.88, 0.04, 0.70, 0.93, 0.23, 0.99, 0.19, 0.79, 1.00, 0.03, 0.01, 0.00) g1<-glm(y~cos((2*pi*lp)/29.530589)+sin((2*pi*lp)/29.530589)) Thanks, Tim Tim Clark Department of Zoology University of Hawaii __

[R] Estimating theta for negative binomial model

2010-05-03 Thread Tim Clark
theta for the glm model? The two models are: mod.nb<-glm.nb(mantas~site,data=mydata) mod.glm<-glm(mantas~site,data=mydata, family=negative.binomial(mod.nb$theta)) How else can I get theta for the family=negative.binomial(theta=???) Thanks! Tim Tim Clark Department of Z

[R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

2009-09-24 Thread Tim Howard
n't use this method in ROCR library(ROCR) data(ROCR.xval) RCdat <- plot(perf, avg="threshold") RCdat ## output: NULL Does anyone have any tricks for piping or extracting these data? Or, perhaps for steering me in another direction? Thanks, Tim From: "Tim Howard" Sub

Re: [R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

2009-09-24 Thread Tim Howard
Whoops, sorry. Here is the full set with the missing lines: library(ROCR) data(ROCR.xval) pred <- prediction(ROCR.xval$predictions, ROCR.xval$labels) perf <- performance(pred,"tpr","fpr") RCdat <- plot(perf, avg="threshold") RCdat Thanks. Tim >>> D

Re: [R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

2009-09-24 Thread Tim Howard
ll see if I can get in there and tweak it. Thanks to both of you for the help. Tim >>> David Winsemius 9/24/2009 9:43 AM >>> On Sep 24, 2009, at 9:09 AM, Tim Howard wrote: > All, > I'm trying again with a slightly more generic version of my first > questio

Re: [R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

2009-09-24 Thread Tim Howard
Yes, that's exactly what I am after. Thank you for clarifying my problem for me! I'll try to dive into the plot.performance function. Best, Tim >>> Tobias Sing 9/24/2009 9:57 AM >>> Tim, if I understand correctly, you are trying to get the numerical values of aver

[R] col headers in read.table()

2009-09-24 Thread Tim Smith
Hi, I was trying to read in a file test.txt, which has the following data: norm normnormclass class class a 1 2 3 4 5 6 b 3 4 5 6 7 8 c 5 6 7 8 9 10 in my R code, I do the f

[R] Data formatting for matplot

2009-09-27 Thread Tim Clark
ting a multiple line plot. Thanks, Tim Tim Clark Department of Zoology University of Hawaii __ 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-gu

Re: [R] Data formatting for matplot

2009-09-27 Thread Tim Clark
? I don't have to use matplot(), but would like to understand its use. Thanks, Tim Tim Clark Department of Zoology University of Hawaii --- On Sun, 9/27/09, Henrique Dallazuanna wrote: > From: Henrique Dallazuanna > Subject: Re: [R] Data formatting for matplot > To: "

Re: [R] Data formatting for matplot

2009-09-28 Thread Tim Clark
Thanks for everyones help. It is great to have a number of options that result in the same graph. Aloha, Tim Tim Clark Department of Zoology University of Hawaii --- On Mon, 9/28/09, Henrique Dallazuanna wrote: > From: Henrique Dallazuanna > Subject: Re: [R] Data formatti

[R] xyplot help - colors and break in plot

2009-09-28 Thread Tim Clark
Data Point", ylab = "MCP Area", ylim = c(0,150), scales = list(tck = c(1, 0)), #Removes tics on top and r-axis pch=19,cex=.4, auto.key = list(title = "Mantas", x = .05, y=.95, corner = c(0,1),border = TRUE)) #Legend Thanks, Tim

Re: [R] xyplot help - colors and break in plot

2009-09-29 Thread Tim Clark
n you have data that groups in clusters separated by a large range. Aloha, Tim Tim Clark Department of Zoology University of Hawaii --- On Mon, 9/28/09, Felix Andrews wrote: > From: Felix Andrews > Subject: Re: [R] xyplot help - colors and break in plot > To: "Tim Cla

[R] bwplot scales in alphabetical order

2009-09-30 Thread Tim Clark
Did not work! id<-factor(id,levels=c(7,8,9,4,5,6,1,2,3),lables=letters[1:9]) Your help would be greatly appreciated. Tim Tim Clark Department of Zoology University of Hawaii __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

Re: [R] bwplot scales in alphabetical order

2009-09-30 Thread Tim Clark
Peter, Thanks, that did it! Tim Tim Clark Department of Zoology University of Hawaii --- On Wed, 9/30/09, Peter Ehlers wrote: > From: Peter Ehlers > Subject: Re: [R] bwplot scales in alphabetical order > To: "Tim Clark" > Cc: r-help@r-project.org > Date: Wednesd

[R] Paste a character to an object

2009-10-03 Thread Tim Clark
nd to "100". So I am stuck. What I want is: homerange[[1]]$polygons$"100" What I need is something like the following, but that works: percent<-c("100","75","50") p=1 paste(homerange[[1]]$polygons$,percent[p],sep="") Thanks fo

Re: [R] Paste a character to an object

2009-10-03 Thread Tim Clark
at is what I am stuck on. The result should be a vector the values for the "100","75",and "50" levels, but not the "90" level. Aloha, Tim Tim Clark Department of Zoology University of Hawaii --- On Sat, 10/3/09, David Winsemius wrote: > From:

Re: [R] Paste a character to an object

2009-10-04 Thread Tim Clark
ode and couldn't figure it out. I knew there had to be an easier way. I greatly appreciate all your help, Tim Tim Clark Department of Zoology University of Hawaii --- On Sat, 10/3/09, David Winsemius wrote: > From: David Winsemius > Subject: Re: [R] Paste a character to an

[R] Satellite ocean color palette?

2009-10-09 Thread Tim Clark
similar to what is found on the CoastWatch web page: http://oceanwatch.pifsc.noaa.gov/imagery/GA2009281_2009282_sst_2D_eddy.jpg Thanks! Tim Tim Clark Department of Zoology University of Hawaii __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] Satellite ocean color palette?

2009-10-09 Thread Tim Clark
Thanks! The colorRampPalette() did just what I need. Tim Tim Clark Department of Zoology University of Hawaii --- On Fri, 10/9/09, Barry Rowlingson wrote: > From: Barry Rowlingson > Subject: Re: [R] Satellite ocean color palette? > To: "Tim Clark" > Cc: r-hel

[R] Bezier interpolation

2009-10-29 Thread Tim Clark
between min(x) and max(x)). Mark Hindell asked the same question in 2006 (http://tolstoy.newcastle.edu.au/R/e2/help/06/12/7034.html). I contacted him and he never found a workable function. Has one been developed since then? Thanks, Tim Tim Clark Department of Zoology University of Hawaii

[R] Time Series methods

2009-11-01 Thread Tim Bean
within the series? (So for a monthly data set, calculate median for all Januarys, all Februarys, all Marchs, etc.) Thanks, Tim [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Time Series methods

2009-11-02 Thread Tim Bean
at 4:09 PM, Kjetil Halvorsen < kjetilbrinchmannhalvor...@gmail.com> wrote: > introduce a factor variable with the months and then use tapply? > > Kjetil > > On Sun, Nov 1, 2009 at 9:07 PM, Tim Bean wrote: > > Hello, I have a quick question about time series methodology. If

[R] Discontinuous graph

2009-11-16 Thread Tim Smith
Hi, I wanted to make a graph with the following table (2 rows, 3 columns): a b c x 1 3 5 y 5 8 6 The first column represents the start cordinate, and the second column contains the end cordinate for the x-axis. The third column contains the y-axis co-ordinate. For example, the first row in the m

[R] Lattice plot

2009-11-17 Thread Tim Smith
Hi, I was trying to get a graph in lattice with the following data frame (7 rows, 5 cols): chr start1 end1 meth positive 1 1 10 20 1.5y 2 2 12 18 -0.7n 3 3 22 34 2.0y 4 1 35 70 3.0y 5 1120 140 -1.3n 6 1180 19

Re: [R] Discontinuous graph

2009-11-17 Thread Tim Smith
David Winsemius : > > On Nov 16, 2009, at 12:40 PM, Tim Smith wrote: > >> Hi, >> I wanted to make a graph with the following table (2 rows, 3 columns): >> a b c >> x 1 3 5 >> y 5 8 6 >> The first column represents the start cordinate, and the second c

[R] Adding columns to lower level of list

2009-11-21 Thread Tim Clark
ke my final result to have two more columns in cut.ltr[[1]]. One with speed, and the second with burst. Thanks in advance for your help. Tim > cut.ltr[1:5] *** List of class ltraj *** Type of the traject: Type II (time recorded) Irregular traject. Variable time lag between

[R] Removing objects from a list based on nrow

2009-11-29 Thread Tim Clark
er=c("A","B"),number=c(1,2)) df3<-data.frame(letter=c("A","B","C","D","E"),number=c(1,2,3,4,5)) df4<-data.frame(letter=c("A","B","C","D","E"),number=c(1,2,3,4,5)) lst<-li

Re: [R] Removing objects from a list based on nrow

2009-11-29 Thread Tim Clark
Linlin, Thanks! That works great! Tim Tim Clark Department of Zoology University of Hawaii --- On Sat, 11/28/09, Linlin Yan wrote: > From: Linlin Yan > Subject: Re: [R] Removing objects from a list based on nrow > To: "Tim Clark" > Cc: r-help@r-project.org > Da

  1   2   3   4   >