Re: [R] A question about list

2011-03-17 Thread jim holtman
use 'unlist' > x <- list(1:10, 20:27, 30:50) > x [[1]] [1] 1 2 3 4 5 6 7 8 9 10 [[2]] [1] 20 21 22 23 24 25 26 27 [[3]] [1] 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 > quantile(unlist(x)) 0% 25% 50% 75% 100% 1.0 15.0 31.0 40.5 50.0 > On Thu, Mar 17, 2011

Re: [R] Using barplot() with zoo -- names.arg not permitted?

2011-03-17 Thread Gabor Grothendieck
On Thu, Mar 17, 2011 at 4:54 PM, David Wolfskill wrote: > On Thu, Mar 17, 2011 at 10:23:33AM -0400, Gabor Grothendieck wrote: >> On Thu, Mar 17, 2011 at 9:38 AM, David Wolfskill wrote: >> ... >> > But the X-axis labels show up as "large integers" -- the POSIXct values >> > are apparently treated

Re: [R] Strange R squared, possible error

2011-03-17 Thread Berwin A Turlach
G'day Gabor, On Thu, 17 Mar 2011 11:36:38 -0400 Gabor Grothendieck wrote: > The idea is that if you have a positive quantity that can be broken > down into two nonnegative quantities: X = X1 + X2 then it makes sense > to ask what proportion X1 is of X. For example: 10 = 6 + 4 and 6 is > .6 of

Re: [R] Using barplot() with zoo -- names.arg not permitted?

2011-03-17 Thread David Wolfskill
On Thu, Mar 17, 2011 at 07:52:29PM -0400, Gabor Grothendieck wrote: > ... > > I hope to learn what I'm misunderstanding. > ... > > This line in the code above produces an invalid object: > class(df$time) <- "POSIXct" > > It should be: > class(df$time) <- c("POSIXct", "POSIXt") Ah. Thank you, Si

Re: [R] Histograms with strings,

2011-03-17 Thread David Wolfskill
On Thu, Mar 17, 2011 at 10:06:21PM +, Khanvilkar, Shashank wrote: > Hello, > Thanks in advance for any help, > > I have read a CSV file in which there is a column for an IP addr as in: > > tmpInFile$V2 > [1] "74.125.224.38" "74.125.224.38" "129.46.71.19" "129.46.71.19" > [5] "129.46.71.1

Re: [R] Strange R squared, possible error

2011-03-17 Thread Gabor Grothendieck
On Thu, Mar 17, 2011 at 8:22 PM, Berwin A Turlach wrote: > G'day Gabor, > > On Thu, 17 Mar 2011 11:36:38 -0400 > Gabor Grothendieck wrote: > >> The idea is that if you have a positive quantity that can be broken >> down into two nonnegative quantities: X = X1 + X2 then it makes sense >> to ask wh

Re: [R] Help with plotting a line that is multicoloured based on levels of a factor

2011-03-17 Thread Tóth Dénes
Hi! Not an elegant solution, but seems to work: date <- c(1:300) flow <- sin(2*pi/53*c(1:300)) levels <- factor(rep(c("high","med","low"),100)) data <- cbind.data.frame(date, flow, levels) colours <- as.numeric(levels)+1 # interpolate resolution <- 0.001 appres <- approx(date,flow,seq(min(date

Re: [R] Strange R squared, possible error

2011-03-17 Thread Berwin A Turlach
G'day Gabor, On Thu, 17 Mar 2011 20:38:21 -0400 Gabor Grothendieck wrote: > > Or am I missing something? O.k., because the residuals don't add to zero, there may be a non-zero correlation between residuals and fitted values, which messes up the equation at the variance level. > Try it on an ex

Re: [R] Help with plotting a line that is multicoloured based on levels of a factor

2011-03-17 Thread baptiste auguie
Hi, because each colour is defined on non-consecutive points, you'll probably need to cut the intervals to define segments around each point. One approach might be the following, d = transform(data, start = date - c(0, diff(date)/2), end = date + c(0, diff(date)/2) ) d$start.y = approx(d$date, d

Re: [R] Help with Time Series Plot‏

2011-03-17 Thread Dennis Murphy
Hi: Try this: plot(ts(rnorm(100), start = 2004, freq = 12), xaxt = 'n', yaxt = 'n', xlab = '', ylab = '') axis(1, at = c(2004:2012), cex.axis = 0.7) axis(2, cex.axis = 0.7) title(xlab = 'My X lab', ylab = 'RQI', cex.lab = 0.1) The illegible dots in the region where the axis labels would nor

Re: [R] Help with plotting a line that is multicoloured based on levels of a factor

2011-03-17 Thread Tóth Dénes
Indeed, I forgot about the segments function. with(d,plot(date,flow,type="n")) with(d,segments(start,start.y,end,end.y,col=colour)) > Hi, > > because each colour is defined on non-consecutive points, you'll > probably need to cut the intervals to define segments around each > point. One approa

Re: [R] Help with Time Series Plot‏

2011-03-17 Thread Peter Ehlers
On 2011-03-17 16:37, Axel Urbiz wrote: Dear List, This is an embarrassing question, but I can seem to make this work…How do I change the font size on the xlab and on the numbers shown in the x-axis on the time series plot below. The arguments cex.lab and cex.axis do not seem to be 'passing' to

Re: [R] R² for non-linear model + comparing linear + non-linear models

2011-03-17 Thread David Winsemius
This mostly seemed like a condensation of what had been said before ... threads from the past which were readily available through searching. http://finzi.psych.upenn.edu/R/Rhelp02/archive/20370.html http://finzi.psych.upenn.edu/R/Rhelp02/archive/73750.html http://finzi.psych.upenn.edu/R/Rhe

Re: [R] Histograms with strings,

2011-03-17 Thread jim holtman
names(which.max(table(x))) This is a 'named' vector On Thu, Mar 17, 2011 at 8:58 PM, Khanvilkar, Shashank wrote: > Thanks > >> which.max(table(x)) > 129.46.71.19 >          10 > > How do I get only "129.46.71.19" back as a str... > > > -Original Message- > From: jim holtman [mailto:jholt

Re: [R] calculating AUCs for each of the 1000 boot strap samples

2011-03-17 Thread Frank Harrell
Taby, At the end of your note are you referring to the bootstrap confidence intervals in the "external validation" case, i.e., not corrrected for overfitting? If so you can get that without the bootstrap (e.g., Hmisc package rcorr.cens function). You can get bootstrap overfitting-corrected ROC a

Re: [R] A question about list

2011-03-17 Thread David Winsemius
On Mar 17, 2011, at 5:44 PM, ufuk beyaztas wrote: Hi dear all, It may be a simple question, i have a list output with different number of elements as following; [[1]] [1] 0.86801402 -0.82974691 0.3974 -0.98566707 -4.96576856 -1.32056754 [7] -5.54093319 -0.07600462 -1.34457280 -1.04

Re: [R] Histograms with strings,

2011-03-17 Thread David Winsemius
On Mar 17, 2011, at 6:06 PM, Khanvilkar, Shashank wrote: Hello, Thanks in advance for any help, I have read a CSV file in which there is a column for an IP addr as in: tmpInFile$V2 [1] "74.125.224.38" "74.125.224.38" "129.46.71.19" "129.46.71.19" [5] "129.46.71.19" "129.46.71.19" "129

Re: [R] Help with plotting a line that is multicoloured based on levels of a factor

2011-03-17 Thread David Winsemius
On Mar 17, 2011, at 6:33 PM, Pamela Allen wrote: Hi All, I'm trying to plot data that is a time series of flows that are associated with a specific level, and I would like each level to represent a colour in a line plot. Here is some data that approximates what I'm using: date=c(1:300)

Re: [R] Spatial cluster analysis of continous outcome variable

2011-03-17 Thread Jon Toledo
I read it but it said PDF file and Ps, didn´t specify which other files, so I attached a csv file, which I thought would work.I have uploaded the file in rapidhare (second option was putting on the web):http://rapidshare.com/files/453101614/Coordinates_and_values.csv Hope this works.Thanks for

[R] pattern of panel waves using R

2011-03-17 Thread Sebastián Daza
Hi everyone, Is there any command to identify the pattern of responses of a database with this format: year id 20081 20091 20082 20092 20083 20093 20084 20094 20104 I just need the frequency of the patterns grouped by id: 2008 2009 2010 = 80 2009

[R] quantmod Some Single Letter Tickers Not getFin

2011-03-17 Thread Sparks, John James
Hi, I have been learning the quantmod package over the last several days. I went to check some of my data pulls against other sources and was surprised to find that a few tickers that have single characters do not successfully scrape from Google Finance using getFin(). Particularly require(quan

[R] trouble in call of "texteval"

2011-03-17 Thread André Camargo
  Hi all. I'm having a little trouble with the  function "texteval" (session package). I have used "texteval" in the construction of the function "ExpandData1". "ExpandData1" does not work as expected. However, when I run only the inside code of "ExpandData1" I get the right result. Apparently "te

Re: [R] Histograms with strings,

2011-03-17 Thread Khanvilkar, Shashank
Thanks > which.max(table(x)) 129.46.71.19 10 How do I get only "129.46.71.19" back as a str... -Original Message- From: jim holtman [mailto:jholt...@gmail.com] Sent: Thursday, March 17, 2011 4:39 PM To: Khanvilkar, Shashank Cc: r-help Subject: Re: [R] Histograms with strings,

[R] Replace split with regex for speed ?

2011-03-17 Thread rivercode
Have timestamp in format HH:MM:SS.MMM.UUU and need to remove the last "." so it is in format HH:MM:SS.MMMUUU. What is the fastest way to do this, since it has to be repeated on millions of rows. Should I use regex ? Currently doing it with a string split, which is slow: >head(ts) [1] 09:30:00

[R] help please: put output into dataframe

2011-03-17 Thread Ram H. Sharma
Dear R community members I have been struggling on this simple question, but never get appropriate solution. So please help. # my data, though I have a large number of variables var1 <- rnorm(500, 10,4) var2 <- rnorm(500, 20, 8) var3 <- rnorm(500, 30, 18) var4 <- rnorm(500, 40, 20) datafr1 <- da

[R] confirmatory factor analysis program in R

2011-03-17 Thread rvohen
Does someone have confirmatory factor analysis program in R,which includes factor loading and some tests?thank you! -- View this message in context: http://r.789695.n4.nabble.com/confirmatory-factor-analysis-program-in-R-tp3386133p3386133.html Sent from the R help mailing list archive at Nabble.

[R] OT: Trends in statistical computing?

2011-03-17 Thread Joshua Wiley
Dear List, First off, this is completely off topic, but I thought others might find it interesting. The good people at IBM are apparently providing written interpretations of results now. I was working with a student the other day with simple nonparametric tests using SPSS, and under the nonpara

[R] Help with setting the y-axis in a plot.

2011-03-17 Thread Caitlin
Hi all. I'm working on an assignment for a psychology class and I am not sure how to adjust the y-axis so it displays the range: 0, 5, 10, 15 The code below is almost ideal: lsd = c(3, 5, 13) mar = c(1, 2, 3) g_range <- range(0, lsd, mar) plot(lsd, type="o", col="blue", ylim=g_range, axes=FALSE,

Re: [R] Help with setting the y-axis in a plot.

2011-03-17 Thread Joshua Wiley
Hi Caitlin, You just need to set the ylim argument to plot (you already did but g_range is 0, 13 not 0, 15. So... On Thu, Mar 17, 2011 at 10:49 PM, Caitlin wrote: > Hi all. > > I'm working on an assignment for a psychology class and I am not sure how to > adjust the y-axis so it displays the ra

<    1   2