Re: [R] R-2.7.2 infected?

2008-09-23 Thread Dave DeBarr
Did you check the md5 checksum on it? Yes; it matched: 540090dd892657804d1099c54d6f770d You're the first to report it, and 2.7.2 has been out for almost a month, so I think it's likely that the CRAN copy is uninfected. Sounds promising. Perhaps it's a false positive from eTrust. If it

Re: [R] R-2.7.2 infected?

2008-09-23 Thread Peter Dalgaard
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, so I think it's likely that the CRAN copy is

Re: [R] Manage huge database

2008-09-23 Thread José E. Lozano
Hello: I've been reading all the replies, and I think I have some good ideas to work on. Right now the code I programmed is running, It has been running in a batch process 20h now, and It has imported 1750 rows out of 2000. I will read docs for the bioconductor package, and I will check the gawk

[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,

[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 list instead of

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 -SP500f13

[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

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, so I think it's likely

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

[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.

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.

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 the

[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: [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 your

[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:

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, like many

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 thujacky at 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

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

2008-09-23 Thread Roger Bivand
jjh21 jjharden at 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 or

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'),

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

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 expected

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

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, and the

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

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

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

[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 see

[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

[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

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

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

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 plane

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, when I do

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

[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

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: Hi

[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 explain

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

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 with

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). This

[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

[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

[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 Use

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 blank as

[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

[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

[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,

[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

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

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 blank. However

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

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

[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

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

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:

[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

Re: [R] gsub difficulty

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

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.org

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

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

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 financial

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 outer() to

[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

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, Shengqiao Li

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

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

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), ... To

[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

[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

[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

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

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

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.

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

[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

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.

Re: [R] How to put given values in lower triangle of splom-plot?

2008-09-23 Thread Deepayan Sarkar
On Sun, Sep 21, 2008 at 2:11 AM, Hofert Marius [EMAIL PROTECTED] wrote: Dear R-experts, I have found a splom-modification online which is given below. This works perfectly, but I would like to have a matrix of given correlation values to be used in the lower triangular part (lower.panel) of

[R] Counting character occurrences in data frame

2008-09-23 Thread Hutchinson,David [PYR]
Hi R-Users, I have a data frame containing year, month, day, and code columns. The code column is a unique character of set ('E','A','B') - I am trying to determine an efficient way of summarizing the count of each of these codes by month and year without having to use for...loops and subsets.

Re: [R] Counting character occurrences in data frame

2008-09-23 Thread Henrique Dallazuanna
Try this: with(DF, tapply(code, list(year, month, code), length)) On Tue, Sep 23, 2008 at 8:10 PM, Hutchinson,David [PYR] [EMAIL PROTECTED] wrote: Hi R-Users, I have a data frame containing year, month, day, and code columns. The code column is a unique character of set ('E','A','B') - I am

Re: [R] Counting character occurrences in data frame

2008-09-23 Thread Charles C. Berry
See ?ftable ?as.data.frame ?xtabs e.g. ftable( xtabs( ~code+year+month, your.df ), col.vars=1 ) as.data.frame( xtabs(~code+year+month, your.df ) ) HTH, Chuck On Tue, 23 Sep 2008, Hutchinson,David [PYR] wrote: Hi R-Users, I have a data frame

  1   2   >