[R] forest plots in column

2008-09-23 Thread carol white
Hi, Is it possible to display different forest plots in the same graphical device (in the same fashion as par(mfrow())) or ideally, display the forest plot of different data in column, labels on the left, data points and CI on the right having the same zero? Any suggestion is welcome. Best, Caro

[R] Unordered multinomial logistic regression with repeated measures

2008-09-23 Thread Naomi Richmond
> Hello > > I would like to know if there is a command in R software for an > unordered multinomial logistic regression with repeated measures using > GEE? > > Any help will be appreciated. > > Kind regards > > Naomi Richmond. [[alternative HTML version deleted]] _

Re: [R] how to set rownames / colnames for matrices in a list

2008-09-23 Thread Antje
Thanks for the hint with the "dimnames". I found a rather similar problem in the mailing list solved like this: x <- matrix(1:4,2) y <- matrix(5:8,2) z <- list(x,y) nm <- c("a","b") nms <- list(nm,nm) z <- lapply(z,function(x){ dimnames(x)<-nms x }) Is there anything wrong using a lis

Re: [R] Help for SUR model

2008-09-23 Thread Arne Henningsen
Hi Bill! On Monday 22 September 2008 23:53, Xianchun Liao wrote: > I am an R beginner and trying to run a SUR model in R framework. > > > > subset(esasp500, Obs <=449 & Obs>=197, select = -Date) ->ev13sub > > c(Obs>=397) & c(Obs<=399) ->d13 > > c(Obs>=400) & c(Obs<=449) ->f13 > > SP500*f13 ->SP500

[R] Predictive Analytics event Oct 24-25 (DC) and Nov 6-7 (SF)

2008-09-23 Thread Elise Johnson
Hi, I wanted to make sure you were all aware of these upcoming events. There is a seminar in Predictive Analytics on Oct. 24-25 in DC, and in San Francisco Nov 6-7. This is intensive training for managers, marketers, and IT people who need to make sense of customer data to predict buying behavior

Re: [R] R-2.7.2 infected?

2008-09-23 Thread Peter Dalgaard
Peter Dalgaard wrote: > Dave DeBarr wrote: >>> Did you check the md5 checksum on it? >>> >> >> Yes; it matched: 540090dd892657804d1099c54d6f770d >> >> > And it is binary identical to the Austria CRAN one. >> >>> You're the first to report it, and 2.7.2 has been out for almost a >>> month,

Re: [R] Building binary package fails because of missing dependent package

2008-09-23 Thread Uwe Ligges
Hans-Peter Suter wrote: On an (Intel Leopard) Mac I try to build a package (mxFinance) which depends on another package (mxGraphics). The dependendy is 1) a 'Depends:' in DESCRIPTION and 2) an import in NAMESPACE. - The build fails if the dependent package (mxGraphics) is not installed in the

Re: [R] R-2.7.2 infected?

2008-09-23 Thread Ajay ohri
This is what it does. It seems like a false alarm because in case of actual infection it seems quite conspicious Ajay www.decisionstats.com http://www.spywareguide.com/product_show.php?id=2569 Full Name: Win32.AdClicker Websearch Read More Type:Trojan SG Index: 5 [Explain] Removal tools:Li

[R] Generalising to n-dimensions

2008-09-23 Thread Laura Bonnett
Hi R-helpers, I have two queries relating to generalising to n dimensions: What I want to do in the first one is generalise the following statement: expand<-expand.grid(1:x[1],1:x[2],...1:x[n]) where x is a vector of integers and expand.grid gives every combination of the set of numbers, so for

[R] odds ratio: how to create reference

2008-09-23 Thread Bunny, lautloscrew.com
HI there, i know this is a basic question, though i need some help because this is somewhat away from my current issue, but nevertheless interesting to me... Lets assume i have some estimated probabilities, say estimated by a logit model. i know i can also state them as an odds ratio. N

Re: [R] Error: subscript out of bounds.

2008-09-23 Thread Patrick Burns
There is a good reason for the error in the matrix case. Consider what would need to happen for: m[3, 5] <- 100 to work. Consistency might be a good thing here, but if so I think it should be an error for the 1-D case. This might be a useful difference between 1-D arrays and vectors. However

Re: [R] lme problems

2008-09-23 Thread Mark Difford
Hi Tommaso, >> I struggle to understand the discrepancy in df between the anova and lme, >> and the >> fact that the interaction term is not significant in the anova but >> significant in lme. To begin with, why try to compare things that are obviously quite different? Surely you can see that t

[R] read.table & readLines behaviour?

2008-09-23 Thread J . delasHeras
Hi, I have been using 'read.table' regularly to read tab-delimited text files with data. No problem, until now. Now I have a file that appeared to have read fine, and the data inside looks correct (structure etc), except I only had 15000+ rows out of the expected 24000. Using 'readLines'

[R] Installing R 2.7.2 on Windows with /LOADINF flag doesn't do anything

2008-09-23 Thread michael watson (IAH-C)
Hi I want to install R from a script, so am following http://cran.r-project.org/bin/windows/base/rw-FAQ.html#Can-I-customize-t he-installation_003f. I first installed R with the /SAVEINF="r_install_config.txt" flag and it created a nice text file for me - perfect. I then uninstalled R, and re-

Re: [R] Generalising to n-dimensions

2008-09-23 Thread Laura Bonnett
crosstable[,,expand[d,1],expand[d,2],expand[d,3],...expand[d,n]] crosstable is just a crosstabulation of an n+2-dimensional dataset and I am trying to pick out those that are in combination 'd' of expand. So for example, for 5-dimensional data using your example: Var1 Var2 Var3 1 11

Re: [R] odds ratio: how to create reference

2008-09-23 Thread Eik Vettorazzi
Hi, You can state a probability p as odds p/(1-p) and vice versa. To get an odds ratio you need actually two odds. Then you can get the odds ration of being/having "a" instead of "b" by odds(a)/odds(b), where "b" is the reference level. If you fit a logistic regression model (which means that y

[R] plot error

2008-09-23 Thread glaporta
HI there, why these lines of code are correct plot(count~spray, data = InsectSprays) plot(InsectSprays$count) but this return an error: plot(count, data = InsectSprays); "data" method is not implemented in plot?! Thanx, Gianandrea -- View this message in context: http://www.nabble.com/plot-er

Re: [R] Generalising to n-dimensions

2008-09-23 Thread Robin Hankin
First bit: > x <- c(3,2,2) > expand.grid(sapply(x,seq_len)) Var1 Var2 Var3 1 111 2 211 3 311 4 121 5 221 6 321 7 112 8 212 9 312 10122 11222 12322 >

Re: [R] plot error

2008-09-23 Thread Barry Rowlingson
2008/9/23 glaporta <[EMAIL PROTECTED]>: > > HI there, > > why these lines of code are correct > plot(count~spray, data = InsectSprays) > plot(InsectSprays$count) > > but this return an error: > plot(count, data = InsectSprays); > > "data" method is not implemented in plot?! The 'plot' function, l

Re: [R] Generalising to n-dimensions

2008-09-23 Thread Robin Hankin
Laura Bonnett wrote: crosstable[,,expand[d,1],expand[d,2],expand[d,3],...expand[d,n]] crosstable is just a crosstabulation of an n+2-dimensional dataset and I am trying to pick out those that are in combination 'd' of expand. So for example, for 5-dimensional data using your example: Var1

Re: [R] R Map using SAS data

2008-09-23 Thread Roger Bivand
Junjie Zhang hotmail.com> writes: > > I'd like to plot some maps. Is it possible for me to use SAS map data in R? It is always a good idea to look at the Task Views on CRAN. In the Spatial Task View, you'll find a description of what you can do using contributed packages, but - as you would

Re: [R] Need help creating spatial correlation for MC simulation

2008-09-23 Thread Roger Bivand
jjh21 gmail.com> writes: > > > Thank you for the input. > > Which command in the spatstat package am I looking for? The documentation is > unclear to me. > > milton ruser wrote: > > > > I think that on spatial stat you will find several ways of simulate > > spatial > > pattern that (point o

Re: [R] odds ratio: how to create reference

2008-09-23 Thread Eik Vettorazzi
You notice that your reference is for the SAS-procedure glm?! To play around and since you provided no data example I've created some data (by an idea from Frank Harrell jr. http://www.biostat.wustl.edu/archives/html/s-news/2002-04/msg00103.html) n <- 30 treat <- rep(c('a','b','c'), length.ou

Re: [R] Installing R 2.7.2 on Windows with /LOADINF flag doesn't do anything

2008-09-23 Thread Duncan Murdoch
On 23/09/2008 5:27 AM, michael watson (IAH-C) wrote: Hi I want to install R from a script, so am following http://cran.r-project.org/bin/windows/base/rw-FAQ.html#Can-I-customize-t he-installation_003f. I first installed R with the /SAVEINF="r_install_config.txt" flag and it created a nice te

Re: [R] read.table & readLines behaviour?

2008-09-23 Thread Peter Dalgaard
[EMAIL PROTECTED] wrote: > > Hi, > > > I have been using 'read.table' regularly to read tab-delimited text > files with data. No problem, until now. > Now I have a file that appeared to have read fine, and the data inside > looks correct (structure etc), except I only had 15000+ rows out of > the e

Re: [R] R-2.7.2 infected?

2008-09-23 Thread Rory.WINSTON
>could this be an intentional attack to compromise a very popular download, and >infect thousands of people.what could be the motivations...i hope its not >some corporate thug here No. False positives are relatively common. What exactly does the Win32/Adclicker.JO trojan do ??? Ajay www.

Re: [R] read.table & readLines behaviour?

2008-09-23 Thread Gabor Grothendieck
Try looking at the result of count.fields to diagnose it. On Tue, Sep 23, 2008 at 5:19 AM, <[EMAIL PROTECTED]> wrote: > > Hi, > > > I have been using 'read.table' regularly to read tab-delimited text files > with data. No problem, until now. > Now I have a file that appeared to have read fine, an

Re: [R] odds ratio: how to create reference

2008-09-23 Thread Philip Twumasi-Ankrah
To select another reference level from the default that R chooses you can use the "relevel" function example; Reorder Levels of Factorwarpbreaks$tension <- relevel(warpbreaks$tension, ref="M") summary(lm(breaks ~ wool + tension, data=warpbreaks)) function findlink(pkg, fn) { var Y, link; Y =

Re: [R] odds ratio: how to create reference

2008-09-23 Thread Frank E Harrell Jr
Philip Twumasi-Ankrah wrote: To select another reference level from the default that R chooses you can use the "relevel" function example; Reorder Levels of Factorwarpbreaks$tension <- relevel(warpbreaks$tension, ref="M") summary(lm(breaks ~ wool + tension, data=warpbreaks)) function findli

Re: [R] Coefficients, OR and 95% CL

2008-09-23 Thread Eik Vettorazzi
l.mod<-glm(menarche~age,binomial) you will get odds-ratios by exponentiating the coefficients of this model, so exp(coef(l.mod)) will do this job. You may notice that this will produce an "OR" for the intercept part as well - which is not interpretable. For the confidence intervals for this

Re: [R] Coefficients, OR and 95% CL

2008-09-23 Thread Nutter, Benjamin
You might also consider ?confint -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jorge Ivan Velez Sent: Monday, September 22, 2008 5:36 PM To: Luciano La Sala Cc: R mailing list Subject: Re: [R] Coefficients, OR and 95% CL Dear Luciano, See ?logistic.dis

Re: [R] Coefficients, OR and 95% CL

2008-09-23 Thread Frank E Harrell Jr
Eik Vettorazzi wrote: l.mod<-glm(menarche~age,binomial) you will get odds-ratios by exponentiating the coefficients of this model, so exp(coef(l.mod)) will do this job. You may notice that this will produce an "OR" for the intercept part as well - which is not interpretable. For the confide

[R] WG: Problem during porting R-2.7.2 on HP-UP 11.11 PA-Risc

2008-09-23 Thread claus-juergen . neumann
Claus-Juergen Neumann/BASF-AG/B ASF

Re: [R] WG: Problem during porting R-2.7.2 on HP-UP 11.11 PA-Risc

2008-09-23 Thread Peter Dalgaard
[EMAIL PROTECTED] wrote: > > Claus-Juergen > Neumann/BASF-AG/B > ASF

Re: [R] WG: Problem during porting R-2.7.2 on HP-UP 11.11 PA-Risc

2008-09-23 Thread Martyn Plummer
On Tue, 2008-09-23 at 15:35 +0200, Peter Dalgaard wrote: > [EMAIL PROTECTED] wrote: > > > > Claus-Juergen > > Neumann/BASF-AG/B

[R] Problem during porting R-2.7.2 on HP-UP 11.11 PA-Risc

2008-09-23 Thread claus-juergen . neumann
Dear "R"-Team, we had some problem to install R-2.7.2 on HP-UX 11.11 on PA-Risc platform. 1.) i load down the Software and made a extract of the tar-file in a specific Directory. 2.) i follow the instructions on the INSTALL file: ./configure 3.) during "configure" we got a error message se

[R] additional parameters in function called by tapply

2008-09-23 Thread Heidemeier Dr, Joachim
Der R-Gurus, first apologies if this is a FAQ, but I due to lack of R-knowledge and terminology I wasn't able to find it. I have the following problem in aggregating results of a model calculation: The results are yearly values of several parameters with several hierarchical spatial factors take

[R] xyplot problem

2008-09-23 Thread Pascal A. Niklaus
Hi all, I am trying to produce some panels with dots in an X/Y plane where the diameter of the dots indicates a Z value (like e.g. earthquake maps where dot sizes indicate magnitudes and X/Y the location). This works fine with xyplot, e.g.: xyplot(1:3~1:3,cex=1:3,pch=16) However, when I do th

Re: [R] gbm error

2008-09-23 Thread jtf2076
Darrin, what version of R are you using have the same question. I have been using the gbm package without problem with R 2.7.1. but I recently switched computers and downloaded 2.7.2. I used the following code, which works fine on my old computer with 2.7.1 but I get an error with the new versio

Re: [R] Statistical question re assessing fit of distribution functions.

2008-09-23 Thread Ted Byers
Thanks Timur While assessing whether or not the best option would be a normal distribution (it won't be, the data in this case LOOKS more poisson, or if I explude the first week of results, a negative exponential; and in my other case, cauchy is more likely), I really need a test that can be appl

Re: [R] xyplot problem

2008-09-23 Thread Henrique Dallazuanna
Is this what you want? xyplot(y~x | grp, panel = function(x, y, ...){ panel.xyplot(x, y, cex = 1:3, pch = 19, col = 2:5) }) On Tue, Sep 23, 2008 at 10:27 AM, Pascal A. Niklaus <[EMAIL PROTECTED]> wrote: > Hi all, > > I am trying to produce some panels with dots in an X/Y p

Re: [R] xyplot problem

2008-09-23 Thread Richard . Cotton
> I am trying to produce some panels with dots in an X/Y plane where the > diameter of the dots indicates a Z value (like e.g. earthquake maps where dot > sizes indicate magnitudes and X/Y the location). > > This works fine with xyplot, e.g.: > > xyplot(1:3~1:3,cex=1:3,pch=16) > > However, wh

Re: [R] R-2.7.2 infected?

2008-09-23 Thread Duncan Murdoch
On 22/09/2008 8:38 PM, Dave DeBarr wrote: I tried downloading R-2.7.2 (http://cran.cnr.berkeley.edu/bin/windows/base/R-2.7.2-win32.exe, both from Berkeley and cran) and both times I got a warning from Computer Associates eTrust Antivirus (version 7.1.710) that the Win32/Adclicker.JO trojan was

[R] Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness?

2008-09-23 Thread zhijie zhang
Dear R users, I hope to explain the concepts of skewness and kurtosis by generating series of distributions with same skewness and different kurtosis or with same kurtosis and different skewness, but it seems that i cannot find the right functions. I have searched the mailing list, but no answ

Re: [R] xyplot problem

2008-09-23 Thread Gabor Grothendieck
xyplot will pass a vector of indices through the subscripts argument to the panel so try this: xyplot(y ~ x | grp, panel = function(x, y, subscripts, ...) { panel.xyplot(x, y, cex = sz[subscripts], ...) } ) On Tue, Sep 23, 2008 at 9:27 AM, Pascal A. Niklaus <[EMAIL PROTECTED]> wrote: > H

[R] Summing up columns with NA elements in a data frame

2008-09-23 Thread ravi
Hi all, When I try to sum a few columns in a data frame (which has NA’s randomly distributed as elements), I am not able to get the summing operation to ignore the NA’s. Let me explain with the following code :   a<-structure(list(a1 = c(NA, 2, 3, NA, 5, 6, NA, 8,

Re: [R] Summing up columns with NA elements in a data frame

2008-09-23 Thread Henrique Dallazuanna
Try this: a$a4 <- rowSums(a, na.rm = T) On Tue, Sep 23, 2008 at 12:20 PM, ravi <[EMAIL PROTECTED]> wrote: > Hi all, > When I try to sum a few columns in a data frame (which has NA's randomly > distributed as elements), I am not able to get the summing operation to > ignore the NA's. > Let me ex

Re: [R] Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness?

2008-09-23 Thread Yihui Xie
Hi, Certainly it's possible. Use any distribution function as long as you can change its skewness and kurtosis, e.g. the Chi-square distribution. The corresponding R functions are p*, q*, d*, and r* - I think you know these functions already (e.g. rchisq()). The only thing that you should be clea

Re: [R] xyplot problem

2008-09-23 Thread Deepayan Sarkar
On 9/23/08, Pascal A. Niklaus <[EMAIL PROTECTED]> wrote: > Hi all, > > I am trying to produce some panels with dots in an X/Y plane where the > diameter of the dots indicates a Z value (like e.g. earthquake maps where dot > sizes indicate magnitudes and X/Y the location). > > This works fine wi

Re: [R] graphing netCDF files

2008-09-23 Thread Paul Hiemstra
Hi Steve, (I'm replying to your e-mail at the mailing list to get it included there also) I've written some stuff on graphs with sp at the R wiki: http://wiki.r-project.org/rwiki/doku.php?id=tips:spatial-data:spatial_data_visualization or have a look at the gallery containing plots with the

Re: [R] Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness?

2008-09-23 Thread hadley wickham
On Tue, Sep 23, 2008 at 9:59 AM, zhijie zhang <[EMAIL PROTECTED]> wrote: > Dear R users, > I hope to explain the concepts of skewness and kurtosis by generating > series of distributions with same skewness and different kurtosis or with > same kurtosis and different skewness, but it seems that i

[R] Modifying output to Google Docs

2008-09-23 Thread Ajay ohri
Hi List, Graphical output to PDF's ,RTF ,CSV is known through R. Can it be modified for outputting to Google Docs (which is basically uploaded files ,published to become html pages) Is there any package on this ? Regards, Ajay -- [[alternative HTML version deleted]] ___

Re: [R] xyplot problem

2008-09-23 Thread hadley wickham
On Tue, Sep 23, 2008 at 8:27 AM, Pascal A. Niklaus <[EMAIL PROTECTED]> wrote: > Hi all, > > I am trying to produce some panels with dots in an X/Y plane where the > diameter of the dots indicates a Z value (like e.g. earthquake maps where dot > sizes indicate magnitudes and X/Y the location). > > T

[R] How to draw the graph of f(x,y) = x * y ?

2008-09-23 Thread Paul Smith
Dear All, The function curve() draws the graph of functions from R to R. Is there some homologous function to curve() to draw functions from R^2 to R? Thanks in advance, Paul __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] superscript text style in axis label in graph

2008-09-23 Thread Benoit Boulinguiez
Hi all, I'd like to use superscript or subscript in the axis label of a graph. Is that possible in R? For instance in xlab="Ce (mmol/m3)", I'd need the "3" of cubic meter in superscript. Regards/Cordialement - Benoit Boulinguiez Ph.D Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 E

[R] Creating a Legend

2008-09-23 Thread Steve Murray
Dear all, I'm trying to create a legend for my graph. I hope to have the title as "Land Use Type" and the two elements being "Urban" and "Rural" with a red point and green point respectively. So far I have the following command, but obviously it isn't correct: > legend("topright", title="Land

Re: [R] superscript text style in axis label in graph

2008-09-23 Thread Jorge Ivan Velez
Dear Benoit, Perhaps plot(1:10,xlab=expression(Ce (mmol/m^3))) See also ?plotmath for more information. HTH, Jorge On Tue, Sep 23, 2008 at 12:57 PM, Benoit Boulinguiez < [EMAIL PROTECTED]> wrote: > Hi all, > > I'd like to use superscript or subscript in the axis label of a graph. Is > that

[R] R2HTML: output from for-loops

2008-09-23 Thread Werner Wernersen
Hi, I am trying to prepare a report with R2HTML using HTMLStart(outdir="./html", filename="report", echo=T, HTMLframe=F) then, for instance, I want to get the output of a loop: for (i in 1:20) print(summary(rnorm(1000))) but only the first of summaries really ends up in the html file. What am I

[R] gsub difficulty

2008-09-23 Thread Charles Annis, P.E.
Greetings R-ians: I know what doesn’t work but I don’t know why, nor how to remedy things. I have a character string containing "." which I want to replace with " " gsub(".", " ", file.label) replaces the every character with a blank. However gsub(".xls", " ", file.label) replaces ".xls" with a

[R] bandwidth selection for locpoly

2008-09-23 Thread Chinthaka Kuruwita
Hello All, Is there a local bandwidth selection routine for local polynomial regression (locpoly) ? Thanks Chinthaka Kuruwita __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www

[R] How to draw arrows or segment lines in 3 dimensional graph

2008-09-23 Thread yabado
Hi I have a dataset which is like this: ID X1 Y1 Z1 X2 Y2 Z2 1 32 127 85 31.7121 85 2 37.9153 88 34.6137 90 3 36.8134 88 36.1155 97 4 31.4136 87 29 113 70

Re: [R] Modifying output to Google Docs

2008-09-23 Thread Barry Rowlingson
2008/9/23 Ajay ohri <[EMAIL PROTECTED]>: > Hi List, > Graphical output to PDF's ,RTF ,CSV is known through R. > Can it be modified for outputting to Google Docs (which is basically > uploaded files ,published to become html pages) > > Is there any package on this ? Do you basically mean an R impl

[R] Oja Median will soon be available

2008-09-23 Thread Daniel Vogel
Hello, in reply to the recent discussion about the Oja median, I would like to announce that by the end of this year we will make a package available that offers a routine for calculating the Oja median. This routine will have both, an exact computation as well as a randomized search heuristic, w

[R] Nonparametric multiple comparison with 'npmc'.

2008-09-23 Thread zollin
Hi dear R users, I want to do multiple comparison with 'npmc' package. But I don't understand which tests use in it? What limitations they have? Than they differ from Tukey or Dunnett tests? Or at least where can I find the detailed information on the tests used in 'npmc' (manual from package is

[R] help with large gamm

2008-09-23 Thread Dave Depew
Hi All, I was wondering if anyone has had any experience running gamm(mgcv) with a large spatial data set? I know the gamm help says this is a bad idea, but I wondered if there was a way to use something like filehash to store the matrices that gamm creates? -- David Depew PhD Candidate Depar

Re: [R] How to draw the graph of f(x,y) = x * y ?

2008-09-23 Thread Duncan Murdoch
On 9/23/2008 12:54 PM, Paul Smith wrote: Dear All, The function curve() draws the graph of functions from R to R. Is there some homologous function to curve() to draw functions from R^2 to R? No, you would have to write it yourself. The basic idea is to use outer() to construct a matrix of f

Re: [R] gsub difficulty

2008-09-23 Thread Marc Schwartz
on 09/23/2008 12:16 PM Charles Annis, P.E. wrote: > Greetings R-ians: > > I know what doesn’t work but I don’t know why, nor how to remedy things. > > I have a character string containing "." which I want to replace with " " > > gsub(".", " ", file.label) replaces the every character with a blan

Re: [R] gsub difficulty

2008-09-23 Thread Charles Annis, P.E.
Thanks! Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com -Original Message- From: Phil Spector [mailto:[EMAIL PROTECTED] Sent: Tuesday, September 23, 2008 1:30 PM To: Charles Annis, P.E. Subject: Re: [R] gsub difficult

Re: [R] Modifying output to Google Docs

2008-09-23 Thread Ajay ohri
I am India based. Winters aren't long at all. I dont have that hot programming skills in C or R packages so need help. Maybe we can get Google on our side this way . They are supposed to have some good stats package too. Ajay On Tue, Sep 23, 2008 at 10:49 PM, Barry Rowlingson < [EMAIL PROTECTED

[R] Options for viewing / using results from lm

2008-09-23 Thread Michael Just
Hello, I would like to state what I am (trying) to do. I have data set. It has 5749 rows (including the header) and 23 columns. The data contains values related to spatial aspects of the 412 landscapes (over various years). I will be making 2 groups from the data based on spatial extent. I will th

Re: [R] Modifying output to Google Docs

2008-09-23 Thread Ajay ohri
Usage of Google Docs for this- PDFs especially Graphical PDFs can be a bit heavy. If sending by email to other users, they may find it difficult to open ,download etc (older computers or low bandwidth). By outputting PDF directly to a Google Doc, it becomes light and portable for people to view

Re: [R] Creating a Legend

2008-09-23 Thread John Fox
Dear Steve, Try legend("topright", title="Land Use Type", cex=0.75, pch=16, col=c("red", "green"), legend=c("Urban", "Rural"), ncol=2) I hope this helps, John -- John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada web: socser

[R] quantile regression: plotting coefficients on only one variable (rq)

2008-09-23 Thread Michael Faye
Dear all. I have a question on plotting the coefficients from a series of mutivariate quantile regressions. The following code plots the coefficients for each RHS variable x1 and x2. Is there a way to plot only the coefficients on x1? In the data I am using, I have a large number of fixed effects

Re: [R] gsub difficulty

2008-09-23 Thread Charles Annis, P.E.
Thanks, Jorge, for another alternative. Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com _ From: Jorge Ivan Velez [mailto:[EMAIL PROTECTED] Sent

Re: [R] gsub difficulty

2008-09-23 Thread Gabor Grothendieck
Any of these will replace every dot with a space: x <- "a.b.c" gsub("\\.", " ", x) gsub("[.]", " ", x) gsub(".", " ", x, fixed = TRUE) chartr(".", " ", x) See ?regex for more info on regular expressions. On Tue, Sep 23, 2008 at 1:16 PM, Charles Annis, P.E. <[EMAIL PROTECTED]> wrote: > Greetings

Re: [R] Options for viewing / using results from lm

2008-09-23 Thread Nordlund, Dan (DSHS/RDA)
I haven't seen a response to this yet, so I will give my $0.25US worth (which is not worth that much anymore ). :-) > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Michael Just > Sent: Tuesday, September 23, 2008 11:01 AM > To: r-help@r-project.or

Re: [R] quantile regression: plotting coefficients on only one variable (rq)

2008-09-23 Thread Mark Difford
Hi Michael, It's in the manual: ?plot.summary.rqs plot(summary(rq(..., tau=c(...)), parm = x1, ...) Regards, Mark. Michael Faye wrote: > > Dear all. > > I have a question on plotting the coefficients from a series of > mutivariate > quantile regressions. The following code plots the coeff

Re: [R] R2HTML: output from for-loops

2008-09-23 Thread Greg Snow
The problem is probably due to what is considered a top level expression to be evaluated and how R2HTML processes these. A simpler alternative is to use txtStart() or etxtStart() from the TeachingDemos package. These work similar to HTMLStart(), but don't do the html markup (the etxt version d

Re: [R] Modifying output to Google Docs

2008-09-23 Thread Barry Rowlingson
2008/9/23 Ajay ohri <[EMAIL PROTECTED]>: > By outputting PDF directly to a Google Doc, it becomes light and portable > for people to view and collaborate on projects in different parts of the > world. > If you extend this to Google spreadsheets this functionality is especially > useful for financ

Re: [R] How to draw the graph of f(x,y) = x * y ?

2008-09-23 Thread Paul Smith
On Tue, Sep 23, 2008 at 6:32 PM, Duncan Murdoch <[EMAIL PROTECTED]> wrote: >> The function curve() draws the graph of functions from R to R. Is >> there some homologous function to curve() to draw functions from R^2 >> to R? > > No, you would have to write it yourself. The basic idea is to use out

[R] How to convert backslash to slash?

2008-09-23 Thread Shengqiao Li
How to use sub, gsub, etc. to replace "\" in a string to "/"? For example,convert "C:\foo\bar" to "C:/foo/bar". Thanks, Shengqiao Li __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide h

Re: [R] How to convert backslash to slash?

2008-09-23 Thread Gabor Grothendieck
This was just posted today: https://stat.ethz.ch/pipermail/r-help/2008-September/174696.html On Tue, Sep 23, 2008 at 4:00 PM, Shengqiao Li <[EMAIL PROTECTED]> wrote: > > How to use sub, gsub, etc. to replace "\" in a string to "/"? > > For example,convert "C:\foo\bar" to "C:/foo/bar". > > Thanks,

Re: [R] Coloring spaces between lines in xyplot

2008-09-23 Thread Seth W Bigelow
R-listers: I am summarizing input I got on a recent query about creating stacked, colored line plots in xyplot. The only input I got regarding xyplot was to use the polygon() command, but this seemed to require some awkward data manipulation. In contrast, I got several snippets of functional code

Re: [R] Options for viewing / using results from lm

2008-09-23 Thread Michael Just
Hello, >Why do you think it is better done in another program? Keeping it in R saves you from the exporting, which you say you are having trouble with. I think it might be better for viewing capabilities. If I had columns of residuals right next to each other I could spot differences amongst loca

Re: [R] How to convert backslash to slash?

2008-09-23 Thread hadley wickham
Except it's a bit tricker because: s <- "C:\foo\bar" print(s) You'd need to work with encodeString( "C:\foo\bar") but that won't work in general: "c:\xjobs" "c:\" "c:\help" Hadley On Tue, Sep 23, 2008 at 3:00 PM, Shengqiao Li <[EMAIL PROTECTED]> wrote: > > How to use sub, gsub, etc. to repla

Re: [R] Options for viewing / using results from lm

2008-09-23 Thread Greg Snow
You can do things like this pseudocode: > fit1 <- lm( y ~ x1, data=mydata ) > fit2 <- lm( y ~ x2, data=mydata ) > fit3 <- lm( y ~ x1*x2, data=mydata ) ... > mydata$r1 <- resid(fit1) > mydata$r2 <- resid(fit2) > mydata$r3 <- resid(fit3) Or > mydata <- cbind(mydata, r1=resid(fit1), r2=resid(fit2)

[R] Proper power computation for one-sided binomial tests.

2008-09-23 Thread Collin Lynch
Hi, I trying to determine the best way to compute the power for a one-sample one-sided binomial test. Specifically I need to sample a population of individuals and ask whether a sample rate of 0% is compatable with a minimum threshold of 3% and how many samples are needed. I have made use of powe

[R] Contional

2008-09-23 Thread Michael Pearmain
Hi All, I'm having trouble selecting rows to delete, that i can't seem to overcome. Below is some sample data, i am trying to dedup the data based on each user, and simultaneously the timestamp (at the side i have highlighted expected row to be removed) I've looked at the lag function but can't

Re: [R] Options for viewing / using results from lm

2008-09-23 Thread Michael Just
Greg, Thanks for the pseudocode and cbind suggestions I think these should be quite helpful. What package is the function 'view' in? As for exporting elsewhere, basically I would to have the values of residuals saved in a spreadsheet in addition to any visualizations created within R. Thanks, M Ju

[R] how to model this

2008-09-23 Thread senne
Dear list, I met with a business case and feel frustrated to model it,scenario is like this 1. logistics department running 2. a lot of customers place orders and goods need to be send to their site (nationwide) 3. order received and then choose from a list of truck carriers to deliver goods,a ch

Re: [R] Statistical question re assessing fit of distribution functions.

2008-09-23 Thread David Scott
On Tue, 23 Sep 2008, Ted Byers wrote: Thanks Timur While assessing whether or not the best option would be a normal distribution (it won't be, the data in this case LOOKS more poisson, or if I explude the first week of results, a negative exponential; and in my other case, cauchy is more likel

[R] Can R load a saved workspace from the Internet?

2008-09-23 Thread Paul Smith
Dear All, Can load() load a saved workspace from an Internet address? I have tried it but I am getting errors: In addition: Warning message: In readChar(con, 5L, useBytes = TRUE) : cannot open compressed file 'http://X/workspace.RData', probable reason 'No such file or directory' Paul ___

Re: [R] Contional

2008-09-23 Thread jim holtman
Is this what you want: TRUE marks the ones to be removed > mark <- (head(x$timestamp, -1) == tail(x$timestamp, -1)) & + (head(x$user_id, -1) == tail(x$user_id, -1)) > x$flag <- c(FALSE, mark) > x Source_type timestamp user_id flag 75381 0 07-07-2008-21:03:55 84

Re: [R] How to convert backslash to slash?

2008-09-23 Thread Duncan Murdoch
On 23/09/2008 4:00 PM, Shengqiao Li wrote: How to use sub, gsub, etc. to replace "\" in a string to "/"? For example,convert "C:\foo\bar" to "C:/foo/bar". If those are R strings, there are no backslashes in the first one. It has a formfeed and a backspace in it. Duncan Murdoch ___

Re: [R] Proper power computation for one-sided binomial tests.

2008-09-23 Thread Peter Dalgaard
Collin Lynch wrote: Hi, I trying to determine the best way to compute the power for a one-sample one-sided binomial test. Specifically I need to sample a population of individuals and ask whether a sample rate of 0% is compatable with a minimum threshold of 3% and how many samples are needed. I

Re: [R] Can R load a saved workspace from the Internet?

2008-09-23 Thread Peter Dalgaard
Paul Smith wrote: Dear All, Can load() load a saved workspace from an Internet address? I have tried it but I am getting errors: In addition: Warning message: In readChar(con, 5L, useBytes = TRUE) : cannot open compressed file 'http://X/workspace.RData', probable reason 'No such file or d

[R] how to model this

2008-09-23 Thread senne
Dear list, I met with a business case and feel frustrated to model it,scenario is like this 1. logistics department running 2. a lot of customers place orders and goods need to be send to their site (nationwide) 3. order received and then choose from a list of truck carriers to deliver goods,a ch

Re: [R] Contional

2008-09-23 Thread Peter Alspach
Michael test[!duplicated(paste(test$timestamp, test$user_id)),] should remove the second (and subsequent) occurrences of duplicates. Your example suggests you don't always want to keep the first occurrence, but the rule which determines which occurrence you want to keep is not obvious to me. HT

Re: [R] Can R load a saved workspace from the Internet?

2008-09-23 Thread Paul Smith
On Tue, Sep 23, 2008 at 11:06 PM, Peter Dalgaard <[EMAIL PROTECTED]> wrote: >> Can load() load a saved workspace from an Internet address? I have >> tried it but I am getting errors: >> >> In addition: Warning message: >> In readChar(con, 5L, useBytes = TRUE) : >> cannot open compressed file 'http

Re: [R] how to model this

2008-09-23 Thread jim holtman
One way of modeling this is to use discrete event simulation since you want some randomness to the time it takes trucks to travel and therefore be available for the next load. I have used this to model various aspects of loading and routing trucks. On Tue, Sep 23, 2008 at 4:58 PM, senne <[EMAIL P

  1   2   >