[R] How to it a "loess curve" and obtain the equation in R?

2011-10-21 Thread ali_protocol
Hi! How can I fit a loess curve to an array (384 x 2). How can I obtain the equation for thi fi? Thanks in advance. -- View this message in context: http://r.789695.n4.nabble.com/How-to-it-a-loess-curve-and-obtain-the-equation-in-R-tp3927622p3927622.html Sent from the R help mailing list archiv

Re: [R] Running R with browser without installing anything

2011-10-21 Thread Jörg Reuter
Thank for this hint. I installed RStudio and it in the night and test it few hour. It is a great tool. Jörg Am 20.10.2011 22:43 schrieb "Liviu Andronic" : > > On Thu, Oct 20, 2011 at 7:29 PM, Bogaso Christofer > wrote: > > Dear all, the company I work for has Matlab installed for > > statistica

Re: [R] Calculating difference between values in data frame based on separate column

2011-10-21 Thread Dennis Murphy
Here's another way, using the reshape2 package: library(reshape2) d <- dcast(df, vial ~ measure, value_var = 'value') d$diff <- with(d, B - A) > d vial A B diff 11 12 26 14 22 30 45 15 33 27 325 44 6 34 28 HTH, Dennis On Fri, Oct 21, 2011 at 3:31 PM, Nathan Miller

Re: [R] "Plotting" text?

2011-10-21 Thread baptiste auguie
Try this, library(grid) grid.newpage() grid.text("text") HTH, baptiste On 22 October 2011 13:26, wrote: > > I noticed that the text() command adds text to a plot. Is there a way to > either make the plot blank or add text to a "blank sheet". I would like > to "plot" a page that contains just

Re: [R] lattice::xyplot/ggplot2: plotting weighted data frames with lmline and smooth

2011-10-21 Thread Dennis Murphy
Hi Michael: The necessary argument to geom_smooth() is weight, not weights (my fault, sorry), so try this instead: ggplot(PearsonLee, aes(x = parent, y = child)) + geom_point(size = 1.5, position = position_jitter(width = 0.2)) + geom_smooth(method = lm, aes(weight = frequency,

Re: [R] Creating affybatch objects from matrix (data from qPCR array)

2011-10-21 Thread ali_protocol
Thanks Martin. Actually I want to do a "cyclic loess", so I chose Affy. How may I find other packages that perform Cyclic loess? -- View this message in context: http://r.789695.n4.nabble.com/Creating-affybatch-objects-from-matrix-data-from-qPCR-array-tp3921559p3927553.html Sent from the R help

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 8:14 PM, Rich Shepard wrote: On Fri, 21 Oct 2011, David Winsemius wrote: The only variable in that dataframe with what appears to be a continuous value (which is how I would expect "total dissolved solids" to be measured) is "quant" Are you saying that the value of quant

[R] Segfault and bad output with fOptions::rnorm.sobol

2011-10-21 Thread Robert McDonald
I have had the following problem with R 2.10, 2.13.1, and 2.13.2, running on Ubuntu linux 10.04, xubuntu 11.10, and a version of Redhat (I think 5). rnorm.sobol is producing impossible random values, and occasionally the routine crashes. Here are samples of the output and the crash message. librar

Re: [R] "Plotting" text?

2011-10-21 Thread Sarah Goslee
Sure. The simplest way is to call plot.new() directly, with no arguments. That will give you a square plot with range of 0-1 on each axis (roughly). If you want a more complex underlying plot, that's possible too. If you take a look at ?par, you'll see options for all the different components of

Re: [R] R v2.13.2 - Cannot find Rcmd on path?

2011-10-21 Thread Joshua Wiley
Hi Hardy, I would double check the PATH variable in Windows. You can check it through the control panel, command prompt, or PowerShell. If the correct path to R is not included, you will need to edit it (the PATH environment variable). I do not know a way to permanently edit it from the comma

Re: [R] Converting data frame into multidimensional array

2011-10-21 Thread Luk Arbuckle
Thank-you for the replies. I would say that the function xtabs() is the simplest method. And thanks for pointing out that a table is also an array, as I clearly wasn't aware of that. On 2011-10-21, at 19:07, Gabor Grothendieck wrote: > On Fri, Oct 21, 2011 at 7:04 PM, Gabor Grothendieck > wrot

[R] "Plotting" text?

2011-10-21 Thread rkevinburton
I noticed that the text() command adds text to a plot. Is there a way to either make the plot blank or add text to a "blank sheet". I would like to "plot" a page that contains just text, no plot lines, labels, etc. Suggestions? Kevin [[alternative HTML version deleted]] _

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, David Winsemius wrote: The only variable in that dataframe with what appears to be a continuous value (which is how I would expect "total dissolved solids" to be measured) is "quant" Are you saying that the value of quant is measuring something with different units depending

[R] estimation of ode from data

2011-10-21 Thread Vassily Shvets
Hello, I'm used to working with R in estimation of distributions and statistical models. I'd like to estimate a simple model of an output y(t) based on a concentration as independent variable, using a timed series of experimental data; I've seen this I think somewhere in R? I'd appreciate to h

[R] R v2.13.2 - Cannot find Rcmd on path?

2011-10-21 Thread Hardy Griesbauer
Hello, I recently updated to R version 2.13.2. With R version 2.10.0, I often created and installed R packages, however, I cannot do this since updating. In other words, when I type in Rcmd build -binary PACKAGENAME I receive an error message: "Rcmd is not recognized as an internal...". I sus

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 6:17 PM, Rich Shepard wrote: On Fri, 21 Oct 2011, David Winsemius wrote: What problem are you trying to solve? What I need now is to compare TDS (total dissolved solids) with specific conductivity and the ions that are normally comprise TDS. Before running any regr

Re: [R] Converting data frame into multidimensional array

2011-10-21 Thread Gabor Grothendieck
On Fri, Oct 21, 2011 at 7:04 PM, Gabor Grothendieck wrote: > On Fri, Oct 21, 2011 at 4:15 PM, Luk Arbuckle wrote: >> Consider the following data frame >> >> X <- data.frame(Titanic) >> >> Does anyone know of an easy way to convert X into a multidimensional >> array?  Example that doesn't work >>

Re: [R] rgl device on web

2011-10-21 Thread Ben qant
Thanks Duncan! I wish I had the time to work on something like that, but I have to stay focused on research. Thanks again for your extensive help! Have a good weekend everyone! Ben On Fri, Oct 21, 2011 at 4:24 PM, Duncan Murdoch wrote: > On 11-10-21 3:38 PM, Ben qant wrote: > >> Hello, >> >> I

Re: [R] Converting data frame into multidimensional array

2011-10-21 Thread Gabor Grothendieck
On Fri, Oct 21, 2011 at 4:15 PM, Luk Arbuckle wrote: > Consider the following data frame > > X <- data.frame(Titanic) > > Does anyone know of an easy way to convert X into a multidimensional > array?  Example that doesn't work > > X <- as.array(X, dim=c(4,2,2,2)) > > To do what I need, X needs to

Re: [R] Converting data frame into multidimensional array

2011-10-21 Thread Richard M. Heiberger
You are probably looking for the structable and related functions in the vcd package. Be sure to read the vignettes library(vcd) X <- data.frame(Titanic) XS <- structable(~Class+Sex+Age+Survived, data=X) as.table(XS) is.array(XS) Rich On Fri, Oct 21, 2011 at 5:56 PM, Luk Arbuckle wrote: > I nee

Re: [R] Calculating difference between values in data frame based on separate column

2011-10-21 Thread Sarah Goslee
Hi, It shouldn't be so complicated. What about simply: > td vial measure value 21 A12 11 B26 42 A30 32 B45 63 A27 53 B32 84 A 6 74 B34 > td <- td[order(td$vial, td$measure),] # make s

[R] Calculating difference between values in data frame based on separate column

2011-10-21 Thread Nathan Miller
Hi all, Say I have a data frame something like the one below with different sample vials, measured before(B) and after(A) some process, with a value recorded at each measurement point vialmeasurevalue 1 B26 1 A12 2 B 45 2

Re: [R] rgl device on web

2011-10-21 Thread Duncan Murdoch
On 11-10-21 3:38 PM, Ben qant wrote: Hello, I'm looking for help putting an interactive rgl package 3d device on the web so that it maintains full functionality. Where should I start? Is it possible? Is there an example I can see? (Note: I'm also looking at putting other normal plots on the web.

Re: [R] subplot strange behavoir

2011-10-21 Thread emorway
Hello Dr. Snow, With regard to your response from earlier this month: When I copy and paste your code I get what is expected, the 2 subplots line up on the same y-value. What version of R are you using, which version of subplot? What platform? I'm still troubled by the fact that layout and sub

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, David Winsemius wrote: What problem are you trying to solve? What I need now is to compare TDS (total dissolved solids) with specific conductivity and the ions that are normally comprise TDS. Before running any regression models I need to look at these data from three p

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 4:38 PM, Rich Shepard wrote: On Fri, 21 Oct 2011, David Winsemius wrote: How are we to determine which lines contain information about the "relationships" of param=="TDS" with whatever cases or variable has values of "Cond" and "SO4"? Are you really trying to compare two

Re: [R] Converting data frame into multidimensional array

2011-10-21 Thread Luk Arbuckle
I need something that will work with any, possibly very large, data frame. This dataset is only a very small example. On 2011-10-21, at 17:31, David L Carlson wrote: > How about? > >> x <- array(Titanic, dim=c(4,2,2,2)) >> str(x> > num [1:4, 1:2, 1:2, 1:2] 0 0 35 0 0 0 17 0 118 154 ... > >

Re: [R] quantmod package

2011-10-21 Thread R. Michael Weylandt
You have to convert the long numbers to time objects; they are kept as POSIXct values so just apply as.POSIXct() to them. However, you may just want to store the time data as a string containing the same info: try this, library(quantmod) tck = "YHOO" filename = paste(tck, ".txt", sep="") while(T

Re: [R] Converting data frame into multidimensional array

2011-10-21 Thread David L Carlson
How about? > x <- array(Titanic, dim=c(4,2,2,2)) > str(x> num [1:4, 1:2, 1:2, 1:2] 0 0 35 0 0 0 17 0 118 154 ... -- David L Carlson Associate Professor of Anthropology Texas A&M University College Station, TX 77843-4352 -Original Message- From

[R] Reading in and modifying multiple datasets in a loop

2011-10-21 Thread Debs Majumdar
Hi,   I have been given a set of around 300 files where there are 5 files corresponding to each chunk. E.g. Chunk 1 for chr1 contains these 5 files:         chr1.one.phased.impute2.chunk1     chr1.one.phased.impute2.chunk1_info     chr1.one.phased.impute2.chunk1_info_by_sample    

Re: [R] plotting with a symbol on every nth point

2011-10-21 Thread R. Michael Weylandt
Try something like this: plot(x,type="o", pch = c(5,rep(NA,9))) for, e.g., every 10th point. Michael Weylandt On Fri, Oct 21, 2011 at 5:18 PM, zugi young wrote: > Hi, > > I would like to produce a plot with a symbol on every nth point in a time > series data, like the one in the following: > >

[R] plotting with a symbol on every nth point

2011-10-21 Thread zugi young
Hi, I would like to produce a plot with a symbol on every nth point in a time series data, like the one in the following: http://www.phon.ucl.ac.uk/home/yi/ProsodyPro/EnglishFocus.png x <- seq(-100,1000,25) plot(x,type="l") Could someone help me out with the above example? Thanks

Re: [R] glm-poisson fitting 400.000 records

2011-10-21 Thread Ken
Your memory shouldn't be capped there, try ?memory.size and ?memory.limit. Background less things. Good luck, Ken Hutchison On Oct 21, 2554 BE, at 11:57 AM, D_Tomas wrote: > My apologies for my vague comment. > > My data comprises 400.000 x 21 (17 explanatory variables, plus response

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, David Winsemius wrote: How are we to determine which lines contain information about the "relationships" of param=="TDS" with whatever cases or variable has values of "Cond" and "SO4"? Are you really trying to compare two disjoint groups on some statistic like the means and

[R] Converting data frame into multidimensional array

2011-10-21 Thread Luk Arbuckle
Consider the following data frame X <- data.frame(Titanic) Does anyone know of an easy way to convert X into a multidimensional array? Example that doesn't work X <- as.array(X, dim=c(4,2,2,2)) To do what I need, X needs to be converted into an array of dimensions c(4,2,2,2) in this case, not

Re: [R] Find a particular point on a curve

2011-10-21 Thread Joanie
The most important thing is the point termed "C" (on the image): http://r.789695.n4.nabble.com/file/n3926631/courbe_temp%C3%A9rature.png which is the first point (time, temperature) where temperature stabilizes after the temperature drop (end of feeding). The definition of that particular point i

Re: [R] glm-poisson fitting 400.000 records

2011-10-21 Thread D_Tomas
My apologies for my vague comment. My data comprises 400.000 x 21 (17 explanatory variables, plus response variable, plus two offsets). If I build the full model (only linear) I get: Error: cannot allocate vector of size 112.3 Mb I have a 4GB RAM laptop... Would i get any improvemnt on a 8G

[R] Kleinberg's burst detection algorithm

2011-10-21 Thread Stavros Macrakis
Has anyone here implemented Jon Kleinberg's burst detection algorithm ("Bursty and Hierarchical Structure in Streams" http://www.cs.cornell.edu/home/kleinber/bhs.pdf)? I'd rather not reimplement if there's already running code available Thanks, -s [[alternative HTML vers

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 3:02 PM, Rich Shepard wrote: On Fri, 21 Oct 2011, David Winsemius wrote: First you need to clarify whether "TDS" is the name of a column or a possible value in a column named "param". This whole painful multi-question process would be greatly accelerated if you offered str

[R] rgl device on web

2011-10-21 Thread Ben qant
Hello, I'm looking for help putting an interactive rgl package 3d device on the web so that it maintains full functionality. Where should I start? Is it possible? Is there an example I can see? (Note: I'm also looking at putting other normal plots on the web.) I'd like to stay within R as much as

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, David Winsemius wrote: First you need to clarify whether "TDS" is the name of a column or a possible value in a column named "param". This whole painful multi-question process would be greatly accelerated if you offered str(chemdata). Yes, I did on a different thread, bu

Re: [R] lattice::xyplot/ggplot2: plotting weighted data frames with lmline and smooth

2011-10-21 Thread Michael Friendly
Thanks very much, Dennis. See below for something I don't understand. On 10/21/2011 12:15 PM, Dennis Murphy wrote: Hi Michael: Here's one way to get it from ggplot2. To avoid possible overplotting, I jittered the points horizontally by ± 0.2. I also reduced the point size from the default 2 an

Re: [R] PCA and Regression with complex categorical variables

2011-10-21 Thread David Winsemius
Did you perhaps send an HTML message? As detailed in the Posting Guide, those get scrubbed by the mail-server. On Oct 21, 2011, at 10:48 AM, seanstcl...@verizon.net wrote: -- David Winsemius, MD West Hartford, CT __ R-help@r-project.org mailing li

Re: [R] Scatterplot with the 3rd dimension = color?

2011-10-21 Thread Kerry
Beautiful! It works perfectly, thanks! kb On Oct 21, 7:42 am, Jim Lemon wrote: > On 10/21/2011 06:25 AM, Kerry wrote: > > > Can someone please help me out with this? The ggplot2 suggestion works > > great but I've spent a few days trying to figure out how to plot 2 > > variables with it and I'm

Re: [R] quantmod package

2011-10-21 Thread ATANU
thanks for the help. but with that code it is possible to save the current quotes in a text file(only the date-time in the first columnis not preserved). when i used read.table and tried to convert it into an xts object it shows error as it cannot take the indices as time object. same case happens

Re: [R] How to use gev.fit (package ismev) under box constraints?

2011-10-21 Thread NoSkill ButStyle
Hallo, it seems as if something did not work with my first email I would like to estimate parameters of a general extreme value (GEV) distribution using maximum likelihood as implemented in the gev.fit function of package ismev. If I do the follwing: y.training<- c(22, 22, 18, 19, 18, 18, 22,

[R] PCA and Regression with complex categorical variables

2011-10-21 Thread seanstclair
__ 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, self-contained, reproducible code.

Re: [R] Scatterplot with the 3rd dimension = color?

2011-10-21 Thread Kerry
Awesome, thank you so much for this! I plan to play around with this more next week with my actual data, but it provides a lot more options than I had before I posted. The link will help too. kb On Oct 20, 8:18 pm, Dennis Murphy wrote: > AFAIK, you can't 'add' two ggplot2 graphs together; the pr

[R] Serialization help.

2011-10-21 Thread rkevinburton
I have the following code: c <- file("c:/temp/r/SkuSalesModel.br", "rb") s <- unserialize(c) close(c) rm(c) And it worked as late as yesterday. Today when I came in I get the following error: Error in .Call("R_unserialize", connection, refhook, PACKAGE = "base") : negative length vectors ar

Re: [R] Arima Models - Error and jump error

2011-10-21 Thread Ken
Perhaps: require(forecast) ?auto.arima # Or look into package fitAR. The first performs seasonal optimization so it is likely better for your application. Ken Hutchison On Oct 21, 2554 BE, at 1:59 PM, Flávio Fagundes wrote: > Hi people, > > I´m trying to development a simple routine to ru

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 2:09 PM, Rich Shepard wrote: On Fri, 21 Oct 2011, David Winsemius wrote: The last part ("in the same column") does not make sense, since I was interpreting the term "parameter" to mean a value in a particular column. Assuming these are R NA's then logical indexing: wit

Re: [R] stacked plot

2011-10-21 Thread Henri-Paul Indiogine
Hi Dennis! Fantastic, great, wonderful, beautiful. I slightly changed your code to adapt it to my situation: ggplot(DF.2, aes(x=file.name, y=value, fill=codes))+geom_histogram(position="stack", stat="identity") + labs(x="document", y="number of codings") ### file.name codes

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, David Winsemius wrote: The last part ("in the same column") does not make sense, since I was interpreting the term "parameter" to mean a value in a particular column. Assuming these are R NA's then logical indexing: with( chemdata, chemdata[!is.na(param1) & !is.na(param2)

[R] Arima Models - Error and jump error

2011-10-21 Thread Flávio Fagundes
Hi people, I´m trying to development a simple routine to run many Arima models result from some parâmeters combination. My data test have one year and daily level. A part of routine is: for ( d in 0:1 ) { for ( p in 0:3 ) { for ( q in 0:3 ) { for ( sd in 0:1 ) {

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, David Winsemius wrote: The last part ("in the same column") does not make sense, since I was interpreting the term "parameter" to mean a value in a particular column. David, That's what I meant: two values from the 'param' column. Assuming these are R NA's then logica

[R] nls making R "not responding"

2011-10-21 Thread Schatzi
Here is the code I am running: library(nls2) modeltest<- function(A,mu,l,b,thour){ out<-vector(length=length(thour)) for (i in 1:length(thour)) { out[i]<-b+A/(1+exp(4*mu/A*(l-thour[i])+2)) } return(out) } A=1.3 mu=.22 l = 15 b = .07 thour = 1:25 Yvals<-modeltest(A,mu,l,b,thour)-.125+runif(25)/4

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 1:04 PM, Rich Shepard wrote: On Fri, 21 Oct 2011, Weidong Gu wrote: No easy way out with missing data problems, all imputations are based on some strong and untestable assumptions. Thanks for the insights. Let me rephrase my question in a way that should work: is th

Re: [R] plotting average effects.

2011-10-21 Thread Dennis Murphy
Hi: Your approach to computing the means is not efficient; a better way would be to use the aggregate() function. I would start by combining the grouping variable and the three prediction variables into a data frame. To get the groupwise mean for all three prediction variables, you can use a formu

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, B77S wrote: I know in my experience "Cond" (conductivity??) doesn't vary much within a stream except for during high flow events, and I would imagine the same is true for TDS. This is generally true, but not in the streams with which we're working. TDS values, for exampl

Re: [R] (no subject)

2011-10-21 Thread Daniel Nordlund
I believe you could also set your subscription to NOMAIL and then read the posts from the R-help archive. This would also allow you to post to R-help since you are still subscribed. Hope this is helpful, Dan Daniel Nordlund Bothell, WA USA > -Original Message- > From: r-help-boun...@

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, Weidong Gu wrote: No easy way out with missing data problems, all imputations are based on some strong and untestable assumptions. Thanks for the insights. Let me rephrase my question in a way that should work: is there a way to subset my comprehensive data frame ('ch

Re: [R] Specifying Greek Character in Lattice Plot Label

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, Luke Miller wrote: The following produces something very similar to David's method: plot(1,1, xlab = expression(paste("Conductivity (", mu, "S / cm)"))) but with a slightly different slash character. I think David's method is more "correct", but I've used the above method in

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread B77S
I know in my experience "Cond" (conductivity??) doesn't vary much within a stream except for during high flow events, and I would imagine the same is true for TDS. If these are all low flow values, you could possibly determine a mean/median value to use for the missing data points. Obviously this

Re: [R] Working With Variables Having Different Lengths

2011-10-21 Thread Weidong Gu
Sounds like you are dealing with missing data problem. At default, lm or glm would only keep observations with complete records (complete case analysis). This can be problematic if you have many missing variables and missing values occur not completely at random (i.e., missing values are dependent

Re: [R] Specifying Greek Character in Lattice Plot Label

2011-10-21 Thread Luke Miller
The following produces something very similar to David's method: plot(1,1, xlab = expression(paste("Conductivity (", mu, "S / cm)"))) but with a slightly different slash character. I think David's method is more "correct", but I've used the above method in the past with some success. On Fri, Oct

Re: [R] (no subject)

2011-10-21 Thread Trevor Davies
Alternatively, since you are on gmail you can set up a folder and filter so all r-help emails bypass you inbox and go right to an r-help folder (or something). I find it very useful for just browsing during down time so I can offer my assistance or move to an 'r-keepers ' folder the for little gem

Re: [R] lattice::xyplot/ggplot2: plotting weighted data frames with lmline and smooth

2011-10-21 Thread Dennis Murphy
Hi Michael: Here's one way to get it from ggplot2. To avoid possible overplotting, I jittered the points horizontally by +/- 0.2. I also reduced the point size from the default 2 and increased the line thickness to 1.5 for both fitted curves. In ggplot2, the term faceting is synonymous with condit

[R] Working With Variables Having Different Lengths

2011-10-21 Thread Rich Shepard
Because of regulatory requirement changes over several decades and weather conditions preventing site access the variables in my data set have different lengths. I'd like guidance on how to perform linear regressions and other models with these variables. For example, there are 2206 rows for

Re: [R] replicating SAS's "proc rank" procedure

2011-10-21 Thread David L Carlson
You can get the same results with the cut() function in R: cut(cars$speed, breaks=quantile(cars$speed, probs=c(0:15/15)), labels=1:15, include.lowest=TRUE) -- David L Carlson Associate Professor of Anthropology Texas A&M University College Station, TX 7

Re: [R] Specifying Greek Character in Lattice Plot Label

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, David Winsemius wrote: plot(1,1, xlab=expression(Conductivity~"("*mu*S/cm*")") ) Thank you, David. It did not occur to me to look for a help page. I'll read that now that I looked and found it. Rich __ R-help@r-project.org ma

Re: [R] Multiple factorial comparison LSD

2011-10-21 Thread Richard M. Heiberger
Vera, The glht function in the multcomp package provides the capability you are looking for. The MMC functions in the HH package build on the ghlt function. There are examples in ?MMC on data with more than one factor. Rich On Fri, Oct 21, 2011 at 4:55 AM, Vera Marjorie E. Velasco < velas...@uni

Re: [R] plotting average effects.

2011-10-21 Thread gradstudent
i will include the data to read if if you so choose. dat <- read.dta("http://quantoid.net/hw1_2011.dta";) model in question: mod99 <- glm(democracy ~ popc100kpc + ngrevpc, data=dat, family=binomial) -- looking for average effects code, with error on mod99. popckpc is coded in 1

Re: [R] cph/nomogram Design/RMS package hazard ratio: interquartile vs per unit

2011-10-21 Thread David Winsemius
On Oct 20, 2011, at 8:22 PM, renee wrote: Hello, I am constructing a nomogram using cph and nomogram commands in Dr. Harrell's Design/RMS package. The HR that I obtain for dichotomous and categorical variables are identical to those that I obtain using STATA stcox. When posting to r-help it

Re: [R] Specifying Greek Character in Lattice Plot Label

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 11:27 AM, Rich Shepard wrote: For an axis label I want to include the Greek letter mu within the string. I've not found the proper way of including that expression within the string. What I want is "Conductivity (uS/cm)" with the 'u' replaced by mu. When I try "Conduc

Re: [R] stair-step plot

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 6:39 AM, knut-o wrote: Hello Is it possible to map a plot with horizontal lines like in the step- plot, but without the vertical lines? There is no function named 'step-plot'. If you are talking about the plot.stepfun function then look at the "verticals" argument.

[R] Specifying Greek Character in Lattice Plot Label

2011-10-21 Thread Rich Shepard
For an axis label I want to include the Greek letter mu within the string. I've not found the proper way of including that expression within the string. What I want is "Conductivity (uS/cm)" with the 'u' replaced by mu. When I try "Conductivity (" expression(paste(mu)) "S/cm)" I get an error.

[R] lattice::xyplot/ggplot2: plotting weighted data frames with lmline and smooth

2011-10-21 Thread Michael Friendly
In the HistData package, I have a data frame, PearsonLee, containing observations on heights of parent and child, in weighted form: library(HistData) > str(PearsonLee) 'data.frame': 746 obs. of 6 variables: $ child: num 59.5 59.5 59.5 60.5 60.5 61.5 61.5 61.5 61.5 61.5 ... $ parent

Re: [R] Change column/row-name

2011-10-21 Thread David Winsemius
On Oct 21, 2011, at 1:57 AM, Jörg Reuter wrote: Hi, I am very happy. My problems are solved without one little thing: (Iske <- matrix(c(1, 1, 1, 2, 2, 2, 1, 1, 1, 5, 1, 2, 2, 2, 1, 1, 1, 4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 5, 4, 4, 4, 4, 4, 4, 2, 2, 2, 2, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5,

[R] question about aggregate

2011-10-21 Thread Adel ESSAFI
Hello I am discovering R and I find it is really very powerful. However, I find some newbie difficulties. Here, I have a data frame with manu values that I want to calculate the frequency (the nomber of line) of the some criteria. For exemple here, I want it to print the number of occurence wher

[R] question about aggregate

2011-10-21 Thread Adel ESSAFI
Hi list I am discovering R, and -- PhD candidate in Computer Science Address 3 avenue lamine, cité ezzahra, Sousse 4000 Tunisia tel: +216 97 246 706 (+33640302046 jusqu'au 15/6) fax: +216 71 391 166 [[alternative HTML version deleted]] __ R-h

Re: [R] glm-poisson fitting 400.000 records

2011-10-21 Thread Ben Bolker
D_Tomas hotmail.com> writes: > > Hi, > > I am trying to fi a glm-poisson model to 400.000 records. I have tried biglm > and glmulti but i have problems... can it really be the case that 400.000 > are too many records??? > > I am thinking of using random samples of my dataset. > "I hav

Re: [R] (no subject)

2011-10-21 Thread Uwe Ligges
On 21.10.2011 13:02, Lisa Henault wrote: can i be taken off of this mailing list please? is there another way that you can access this without having to get all the emails?? [[alternative HTML version deleted]] __ R-help@r-project.org maili

Re: [R] stair-step plot

2011-10-21 Thread Duncan Murdoch
On 21/10/2011 6:39 AM, knut-o wrote: Hello Is it possible to map a plot with horizontal lines like in the step-plot, but without the vertical lines? Not in the basic plot function, but you can write your own fairly easily, using segments(). For example: x <- y <- 1:10 plot(x,y, type='n') s

Re: [R] How to create time series objects combining two vectors

2011-10-21 Thread Gabor Grothendieck
On Fri, Oct 21, 2011 at 6:24 AM, sarelseerower wrote: > I am new to R and trying to understand time series objects. > > I have 2 vectors, one containing rainfall values  (lets call the vector > "rain") and the other the time/date in seconds (lets call it "time"). Is > there a method to create a ti

Re: [R] POT package

2011-10-21 Thread Amina Shahzadi
> npy <- length(events1[, "obs"])/(diff(range(ardieres[, "time"], + na.rm = TRUE)) - diff(ardieres[c(20945, 20947), "time"])) This line is from the mannual "A user's Guide to POT Approach". I am just trying to ask why the values 20945, 20947 are used?? Regards On Fri, Oct 21, 2011 at 4:19 PM, R.

Re: [R] plotting average effects.

2011-10-21 Thread Bart Joosen
How about posting a reproducible sample, so that we can see what is going on? Read the posting guide!!! -- View this message in context: http://r.789695.n4.nabble.com/plotting-average-effects-tp3923982p3925324.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] 'Apply' giving me errors

2011-10-21 Thread kickout
Thanks for the tips/adviceI actually used a different solution to circumvent this, but Uwe's solutions would also work -- View this message in context: http://r.789695.n4.nabble.com/Apply-giving-me-errors-tp3923880p3925377.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] R-help Digest, Vol 104, Issue 21

2011-10-21 Thread mihalicza . peter
Október 19-től 21-ig irodán kívül vagyok, és az emailjeimet nem érem el. Sürgős esetben kérem forduljon Kárpáti Edithez (karpati.e...@gyemszi.hu). Üdvözlettel, Mihalicza Péter I will be out of the office from 19 till 21 October with no access to my emails. In urgent cases please contact Ms. Ed

[R] add=TRUE or similar in spplot?

2011-10-21 Thread infochat
Dear Helper, I have a spatial lines data frame object 'spRiverDf'. The data frame consists of numbers {0,1,...,5}. And I have a vector 'colorS' of length 6 with different colours. If I make a plot with spplot I get a plot of the lines - colours depending on there number in the data frame colum

Re: [R] R-help Digest, Vol 104, Issue 19

2011-10-21 Thread Dénes TÓTH
> > On Oct 21, 2011, at 09:01 , Martin Maechler wrote: > >>> "ARE" == Alex Ruiz Euler >>>on Wed, 19 Oct 2011 14:05:16 -0700 writes: >> >>ARE> Motion supported. Very. >> >>ARE> On Wed, 19 Oct 2011 15:40:14 +0200 >>ARE> peter dalgaard wrote: >> Argh! Someone

[R] glm-poisson fitting 400.000 records

2011-10-21 Thread D_Tomas
Hi, I am trying to fi a glm-poisson model to 400.000 records. I have tried biglm and glmulti but i have problems... can it really be the case that 400.000 are too many records??? I am thinking of using random samples of my dataset. Many thanks, -- View this message in context: http://r.78

[R] How to use gev.fit (package ismev) under box constraints?

2011-10-21 Thread NoSkill ButStyle
__ 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, self-contained, reproducible code.

[R] (no subject)

2011-10-21 Thread Lisa Henault
can i be taken off of this mailing list please? is there another way that you can access this without having to get all the emails?? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

[R] stair-step plot

2011-10-21 Thread knut-o
Hello Is it possible to map a plot with horizontal lines like in the step-plot, but without the vertical lines? Thanks, knut -- View this message in context: http://r.789695.n4.nabble.com/stair-step-plot-tp3924903p3924903.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] How to remove multiple outliers

2011-10-21 Thread aajit75
Hi Michael, Thanks for the help. Yes, I have gone through the document for ?outlier. As it removes one outlier at a time, being new to R, I was woondering is there any function available for removing multiple outliers whithout calling say rm.outlier for n number of time because n is not finite he

[R] How to create time series objects combining two vectors

2011-10-21 Thread sarelseerower
I am new to R and trying to understand time series objects. I have 2 vectors, one containing rainfall values (lets call the vector "rain") and the other the time/date in seconds (lets call it "time"). Is there a method to create a time series object simply by giving the "rain" and "time" vectors

[R] varying coefficients model

2011-10-21 Thread Soberon Velez, Alexandra Pilar
Dear members, I'm trying to estimate a varying coefficients model using the local polynomial estimation method in two case (univariate and bivariate) and with two smooth functions. In the univariate case I use: m1<-smooth.lf(x=lp(z1,by=x1,deg=1,h=0.7,ev=z1)+lp(z3,by=expl2,deg=1,h=0.7,ev=z1

Re: [R] windows limits

2011-10-21 Thread Duncan Murdoch
On 21/10/2011 9:36 AM, Ben qant wrote: Hello, Using the rgl package, I can set the device window to any dimension (that I have tested): par3d(windowRect=c(1,1,700,700)) With windows I can't get the window to span from the top to the bottom of the monitor. In the following, no matter how large t

Re: [R] xyplot() or splom()?: two factors from same data frame

2011-10-21 Thread Rich Shepard
On Fri, 21 Oct 2011, Duncan Mackay wrote: Without a dataset I am not sure what you need. Duncan, Part of the problems I'm trying to resolve come from changing priorities from my client and the regulators. I end up stopping one process and starting on a different one. But, that's life in the

[R] windows limits

2011-10-21 Thread Ben qant
Hello, Using the rgl package, I can set the device window to any dimension (that I have tested): par3d(windowRect=c(1,1,700,700)) With windows I can't get the window to span from the top to the bottom of the monitor. In the following, no matter how large the ypinch value gets it stops, leaving ab

  1   2   >