[R] Odp: adding rows to table

2008-09-18 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 17.09.2008 11:08:29: Greetings everyone, I'm trying to add a specific table or a specific number of rows (e.g.44) to a table with no success. This is my basic table head(dataA) yearplot spp prop.BDCA1DCA2DCA3DCA4 1

Re: [R] ggplot2 - deprecated guide= argument in

2008-09-18 Thread Pedro Barros
Hi Hadley, As for now, I have not yet delved deep enough into the new version to be sure if I like it or not. It's just that I like the freedom to choose, and the previous form was great. Any chance you may put it back as an optional parameter? Thanks, Pedro hadley wrote: On Wed, Sep 17, 2008

Re: [R] inserting values for null

2008-09-18 Thread nalluri pratap
The group would be happy to answer your question if you can show an example data frame on what is there and what is required.   Cheers, Pratap --- On Wed, 17/9/08, Rajasekaramya [EMAIL PROTECTED] wrote: From: Rajasekaramya [EMAIL PROTECTED] Subject: [R] inserting values for null To:

Re: [R] unix-type commandline keystrokes in the windows RGUI

2008-09-18 Thread Peter Dalgaard
Adaikalavan Ramasamy wrote: ... Anyway, here is how to do what you want: 1) Install bash on your Windows machine - You can use cgywin. Or download and unzip http://www.steve.org.uk/Software/bash/ 2) Make the directory to bash.exe and R.exe are in your PATH variable. 3) Start - Run - cmd 4)

Re: [R] help on sampling from the truncated normal/gamma distribution on the far end (probability is very low)

2008-09-18 Thread Moshe Olshansky
Hi Sonia, If I did not make a mistake, the conditional distribution of X given that X 0 is very close to exponential distribution with parameter lambda = 40, so you can sample from this distribution. --- On Mon, 15/9/08, Daniel Davis [EMAIL PROTECTED] wrote: From: Daniel Davis [EMAIL

[R] axis text size in scatter plot

2008-09-18 Thread Benoit Boulinguiez
Hi, Quite simple question I guess. Where could we specify the size of the text (axis,legend) for a scatter plot graphic made with plot()? Regards/Cordialement - Benoit Boulinguiez Ph.D Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 Sciences Chimiques de

Re: [R] axis text size in scatter plot

2008-09-18 Thread Yihui Xie
Generally speaking, this is controlled by the argument cex (or cex.axis, cex.lab, ...). See ?par (and ?plot, ?text, ?axis, ?legend, ?text, ...) Regards, Yihui -- Yihui Xie [EMAIL PROTECTED] Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.name

[R] Tool for simple two class Fisher linear discriminant analysis?

2008-09-18 Thread Yip
Could anyone kindly let me know what tool(s) in R are available/convenient for Fisher linear discriminant analysis? I only need to handle the two class tasks for now. Thank you, Yiping -- View this message in context:

[R] Tool for simple two class Fisher linear discriminant analysis?

2008-09-18 Thread justin bem
Use lda in MASS package Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : Yip [EMAIL PROTECTED] À : r-help@r-project.org Envoyé le : Jeudi, 18 Septembre 2008, 8h04mn 23s Objet : [R] Tool for simple two class Fisher linear discriminant analysis?

[R] Joint distributions

2008-09-18 Thread Jorunn Slagstad
Dear R-help! I need to draw contour lines in a plot of wave heights (Hs) versus peak periods (Tp) showing the joint probabilities of 1-year wave heights~peak periods, 10-year wave heights~peak periods and 100-year wave heights~peak periods. I've used the contourplot() function in the plot I've

[R] Converting EnSeMBL Probe names into Gene Name

2008-09-18 Thread Gundala Viswanath
Dear all, Is there a way with Bioconductor in which I can convert such EnSemBL probe names into the standard gene names? AFFX-M27830_5_at AFFX-M27830_M_at ENSG003_at ENSG005_at ENSG419_at - Gundala Viswanath Jakarta - Indonesia

[R] character position

2008-09-18 Thread Paulo Cardoso
Hi, I'm not being able to capture a position of a 'string' in a character string. In this example: 'There are 20 species in this grid' I would like to capture the string (number) after 'are' and before 'species'. Consider they do not change. I wouldn't like to use substr because stop position

Re: [R] character position

2008-09-18 Thread jim holtman
If it always occurs after 'are', this will work: x - 'There are 20 species in this grid' y - sub(.*are (\\d+).*, \\1, x, perl=TRUE) y [1] 20 On Thu, Sep 18, 2008 at 4:51 AM, Paulo Cardoso [EMAIL PROTECTED] wrote: Hi, I'm not being able to capture a position of a 'string' in a character

Re: [R] character position

2008-09-18 Thread Dimitris Rizopoulos
you could use gsub(), e.g., strg - There are 20 species in this grid gsub([^0-9], , strg) I hope it helps. Best, Dimitris Paulo Cardoso wrote: Hi, I'm not being able to capture a position of a 'string' in a character string. In this example: 'There are 20 species in this grid' I would

[R] Oja median

2008-09-18 Thread Rahul-A.Agarwal
Hi, Can we get the code for calculating Oja median for multivariate data Thanks and Regards Rahul Agarwal Analyst Equities Quantitative Research UBS_ISC, Hyderabad On Net: 19 533 6363 [[alternative HTML version deleted]] __

Re: [R] Oja median

2008-09-18 Thread Richard . Cotton
Can we get the code for calculating Oja median for multivariate data RSiteSearch(oja median) returns a link to this R-help post with code http://finzi.psych.upenn.edu/R/Rhelp02a/archive/12781.html Regards, Richie. Mathematical Sciences Unit HSL

[R] Writing greek letters and subscripts in graph

2008-09-18 Thread Delphine COURVOISIER
Dear R-users, I would like to write Greek letters followed by subscripts in a graph (on the X-axis and in a legend). I would appreciate any help. thank you, delphine courvoisier Delphine Courvoisier Clinical Epidemiology Division University of Geneva

[R] PNG file don't run on mac's?

2008-09-18 Thread Tom Willems
Dear R users, I 'm having problems with creating PNG graphic outputs. Usualy i create reports in HTML format, containing PNG graphics, so they can ealsaly be exported to word and xl and so on. On a windows pc that i use at work all works fine, but it never works on my mac. The HTML's i create

Re: [R] Oja median

2008-09-18 Thread Rahul-A.Agarwal
Hi Richard, Thanks for the code.but I already have a code for bivariate data...I want it for data when it is multivariate or its dimension is greater than two...i mean to say if it is trivariate or higher dimensioncan I get a more generalised form of oja median code Thanks and Regards

Re: [R] Writing greek letters and subscripts in graph

2008-09-18 Thread Richard . Cotton
I would like to write Greek letters followed by subscripts in a graph (on the X-axis and in a legend). I would appreciate any help. See FAQ 7.13 and ?plotmath for examples. Regards, Richie. Mathematical Sciences Unit HSL

[R] R2HTML, HTMLStart(): are plots send to HTML?

2008-09-18 Thread Werner Wernersen
Hi, I just want to make sure that I really misunderstood the documentation: Is HTMLStart redirecting plots as well so that they are included in the HTML report? I tried the below and the plot is not included in the HTML file. Or did I miss some switch or so? Many thanks, Werner

Re: [R] using for variable as rowname

2008-09-18 Thread Ralikwen
Thanks for your help, I didn't know where you can or can't use a for variable. Seems that I have to use rownames to do this anyway, so finally I went for a more procedural solution: v=1:6 for (a in 1:3){ for (b in 4:5) { v-rbind(v,1+a+b:6+a+b) rownames(v)[nrow(v)]=

Re: [R] character position

2008-09-18 Thread Paulo Cardoso
Perfect but when one wants to get something different, like x - 'INPUT FILE record.dat' sub(.*FILE (\\w+).*, \\1, x, perl=TRUE) will return 'record' and not 'record.dat' the [:punct:] is not retained. Paulo -Original Message- From: jim holtman [mailto:[EMAIL PROTECTED] Sent:

Re: [R] character position

2008-09-18 Thread Henrique Dallazuanna
Try: sub(.*FILE (\\w+)|\\.*, \\1, x, perl=TRUE) On Thu, Sep 18, 2008 at 7:39 AM, Paulo Cardoso [EMAIL PROTECTED] wrote: Perfect but when one wants to get something different, like x - 'INPUT FILE record.dat' sub(.*FILE (\\w+).*, \\1, x, perl=TRUE) will return 'record' and not 'record.dat'

Re: [R] matrix to list

2008-09-18 Thread Henrique Dallazuanna
Try this: m - matrix(rep(1:3, 3), 3, dimnames = list(letters[1:3], letters[1:3])) as.data.frame.table(t(m)) On Wed, Sep 17, 2008 at 10:52 PM, Pinder, Adrian [EMAIL PROTECTED] wrote: Hi, can someone please tell me how to convert a square matrix to a list in R? i.e. I want to convert from:

Re: [R] Car.proper C[] matrix

2008-09-18 Thread Roger Bivand
jgasper Jason.Gasper at noaa.gov writes: I have been using the nb2WB() package for the car.normal function in WinBUGS, but it will not create the C[] matrix; it only creates adj[], num[], and weights[]. I was planning on using this function to create the C[] matrix (by using the num

Re: [R] rgl: How to position a window during open3d call

2008-09-18 Thread Duncan Murdoch
Yihui Xie wrote: Thanks Duncan, I have also been wondering about this problem for a long time. If you suspect a bug in a package, it's good to contact the package maintainer. Often (as in this case) they are unaware of the bug and appreciate the report. Duncan Murdoch Regards, Yihui

Re: [R] Use of colour in plots

2008-09-18 Thread Steve Murray
Dear all, I've finally got round to plotting my data and trying to apply colour (had some problems with the data which I needed to rectify first!). I'm having trouble however getting the colour to work as I'd hoped, despite the help offered in previous messages. Just to recap, and with more

[R] Ability estimates for partial credit model

2008-09-18 Thread Andrej Kastrin
Dear all, I'm working on ability estimates using Rasch model. Using the ltm package, the procedure is quite simple: ## Factor Scores for the Rasch model fit - rasch(LSAT) factor.scores(fit) What about Partial Credit Model (PCM)? For PCM I use PCM function from eRm package. Is there any

Re: [R] Oja median

2008-09-18 Thread Martin Maechler
[EMAIL PROTECTED] on Thu, 18 Sep 2008 05:20:56 -0400 writes: Hi, Can we get the code for calculating Oja median for multivariate data Excuse me, but must you really? The Oja median has (finite) breakdown point 2/n, i.e., is not robust in any reasonable sense, and is

Re: [R] Stacked Area Plot

2008-09-18 Thread Jim Lemon
Josip Dasovic wrote: Hi: I've searched the archives and the Internet for hours but have yet to find a way to do stacked area plots (like the kind in Excel) in R. I think that polygon may work but it would require a bit of manipulation of my data. I was hoping for an easier way. An example

Re: [R] character position

2008-09-18 Thread Paulo Cardoso
Thank you all. Works fine. -Original Message- From: Henrique Dallazuanna [mailto:[EMAIL PROTECTED] Sent: quinta-feira, 18 de Setembro de 2008 12:14 To: Paulo Cardoso Cc: jim holtman; r Subject: Re: [R] character position Try: sub(.*FILE (\\w+)|\\.*, \\1, x, perl=TRUE) On

Re: [R] Unexpected behaviour when testing for independence with multiple factors

2008-09-18 Thread Michael Dewey
At 16:03 17/09/2008, Javier Acuña wrote: Hi, I'm a new user of R. My background is Electrical Engineering, so please bear with me if this is a silly question. For future reference you might find ?interaction helpful as another tool in your box. I'm trying to assess whether the results of an

Re: [R] creating rainbow gradients

2008-09-18 Thread Jim Lemon
Gillian Silver wrote: Hi, how can I create a rainbow gradient in R? For example, let's say I have a plot of y = x...and I want the plot to go from red - orange - yellow - green - blue - etc. Right now, I know how to do something like go from red to blue, using the plotrix library:

Re: [R] Coxph and loglik converged before variable X

2008-09-18 Thread Terry Therneau
Loglik converged before variable 1,2 ; beta may be infinite. I have two (related) questions: 1. Can I in general ignore this kind of warning and just pay attention to predictors which produce significant p values? 2. In the specific case of stepAIC, can I ignore this warning and trust

Re: [R] Use of colour in plots

2008-09-18 Thread ONKELINX, Thierry
Steve, Have a look at the ggplot2 package: library(ggplot2) ggplot(Jan, aes(x = PopDensity, y = Average.Burnt.Area.Fraction, colour = factor(Urban.Rural 1.25))) + geom_point() ir. Thierry Onkelinx Instituut voor

Re: [R] ggplot2 - deprecated guide= argument in

2008-09-18 Thread hadley wickham
As for now, I have not yet delved deep enough into the new version to be sure if I like it or not. It's just that I like the freedom to choose, and the previous form was great. Any chance you may put it back as an optional parameter? I doubt it, but you can override the defaults as follows:

Re: [R] character position

2008-09-18 Thread Gabor Grothendieck
Try this: library(gsubfn) strapply(There are 23 species, [0-9]+, as.numeric)[[1]] The gsubfn home pae is at: http://gsubfn.googlecode.com On Thu, Sep 18, 2008 at 4:51 AM, Paulo Cardoso [EMAIL PROTECTED] wrote: Hi, I'm not being able to capture a position of a 'string' in a character

[R] Fatal error .RData

2008-09-18 Thread Green, Paul
We are using R in our class. One student claims to be receiving the message Fatal error: unable to restore saved data in .RData The student is working in a Windows environment with Vista using the precompiled binary distribution (R-2.7.2). After searching the R site and reading FAQs, I

Re: [R] Fatal error .RData

2008-09-18 Thread Duncan Murdoch
On 18/09/2008 9:41 AM, Green, Paul wrote: We are using R in our class. One student claims to be receiving the message Fatal error: unable to restore saved data in .RData The student is working in a Windows environment with Vista using the precompiled binary distribution (R-2.7.2). After

Re: [R] Plotting curves in lattice panels

2008-09-18 Thread Deepayan Sarkar
On 9/18/08, David Scott [EMAIL PROTECTED] wrote: I have a data set concerning ferritin levels in blood. There are three relevant columns for this question, ferritin (continuous), score (ordered, from 0 to 8) and gender. There is a good linear relationship between log(ferritin) and score for

Re: [R] Plotting curves in lattice panels

2008-09-18 Thread David Scott
On Thu, 18 Sep 2008, Deepayan Sarkar wrote: On 9/18/08, David Scott [EMAIL PROTECTED] wrote: I have a data set concerning ferritin levels in blood. There are three relevant columns for this question, ferritin (continuous), score (ordered, from 0 to 8) and gender. There is a good linear

Re: [R] Use of colour in plots

2008-09-18 Thread Steve Murray
Jim and all, Maybe I've misunderstood ?color.scale (appologies if this is so), but I don't think this is what I need. I'm not looking to scale the colours of points, instead I simply want to assign each point a colour (either red or green) based on it's value in the Urban.Rural column. To

Re: [R] unix-type commandline keystrokes in the windows RGUI

2008-09-18 Thread Michael Frumin
wow, ok, thanks everyone lots for me to explore. I will take some time, and if none of these achieve what I'm trying to get at I will come back here. Thanks! Mike Adaikalavan Ramasamy wrote: Well, I don't see why you need the CTRL-R functionality when you can just as rapidly and

Re: [R] PNG file don't run on mac's?

2008-09-18 Thread Don MacQueen
Tom, First of all, please in the future follow the posting guide and include a reproducible example. Yours is not, because the objects named bxpdtrail and ResultDir (and maybe others) are not available to readers of this list. And you defined a dataframe named size2 but then used one named

[R] rcurl and cookies

2008-09-18 Thread zubin
Hello, need some help on using Rcurl to navigate a site and the use of session cookies. I suspect the issue i am having presently is I am not handling session cookies properly. At a high level, in need to create a dataset for some analysis, my background is in using R for statistical work, I

[R] Odp: Fatal error .RData

2008-09-18 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 18.09.2008 15:41:07: We are using R in our class. One student claims to be receiving the message Fatal error: unable to restore saved data in .RData I sometimes get such message when I save .RData with some nonstandard packages loaded and trying to open R

[R] Suggestion: Re: Fatal error .RData

2008-09-18 Thread Henrik Bengtsson
R-core, may I suggest that the error message returns the absolute path (or even the relative path) and not just the filename/basename, e.g. Fatal error: unable to restore saved data in C:/Users/foo/.RData /Henrik On Thu, Sep 18, 2008 at 6:41 AM, Green, Paul [EMAIL PROTECTED] wrote: We are

[R] Fwd: Reporting with down and across variables

2008-09-18 Thread Matthew Pettis
[Reposting with changed example and font spec removed] Hi, I want to take the dataframe df generated below and reshape the data with column names being w, x, y, and the different levels of z. The values under the different levels of z are the corresponding values of r. I've tried reshape and

[R] Problem installing packages in newer versions of R

2008-09-18 Thread Miha Staut
Dear all, I was wandering what could be wrong with my system (regularly updated Fedora core 8) so that installing packages does not succeed with almost every package. I follow the procedure specified in the help file R-admin section 6.3. This is not a feature new to the current version but

Re: [R] Reporting with down and across variables

2008-09-18 Thread hadley wickham
On Thu, Sep 18, 2008 at 10:03 AM, Matthew Pettis [EMAIL PROTECTED] wrote: [Reposting with changed example] Hi, I want to take the dataframe df generated below and reshape the data with column names being w, x, y, and the different levels of z. The values under the different levels of z are

[R] Reporting with down and across variables

2008-09-18 Thread Matthew Pettis
[Reposting with changed example] Hi, I want to take the dataframe df generated below and reshape the data with column names being w, x, y, and the different levels of z. The values under the different levels of z are the corresponding values of r. I've tried reshape and cast, and I can't seem

Re: [R] Problem installing packages in newer versions of R

2008-09-18 Thread hadley wickham
Hi Miha, One tip for diagnosing build problems is to ignore all warnings and look for the first error. In your case it's: init.c:2:15: error: R.h: No such file or directory init.c:3:24: error: Rinternals.h: No such file or directory Which possibly suggests that you don't have the necessary

Re: [R] Suggestion: Re: Fatal error .RData

2008-09-18 Thread hadley wickham
And does it need to be a fatal error? load(test.txt) doesn't produce a fatal error even though the file clearly isn't in the correct format. Or is there the possibility that loading a malformed rdata file introduces memory corruption? On Thu, Sep 18, 2008 at 9:55 AM, Henrik Bengtsson [EMAIL

Re: [R] Reporting with down and across variables

2008-09-18 Thread Matthew Pettis
It does, thanks! Did I miss in the documentation that the variable has to be named 'value'? On Thu, Sep 18, 2008 at 10:08 AM, hadley wickham [EMAIL PROTECTED] wrote: On Thu, Sep 18, 2008 at 10:03 AM, Matthew Pettis [EMAIL PROTECTED] wrote: [Reposting with changed example] Hi, I want to

Re: [R] Spectrum of a kernel

2008-09-18 Thread Rory.WINSTON
Maybe something like: k - kernel(...) spectrum(c(rev(k$coef),k$coef)) Since the kernel is symmetric? Rory Winston RBS Global Banking Markets Office: +44 20 7085 4476 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: 17 September

Re: [R] Suggestion: Re: Fatal error .RData

2008-09-18 Thread Duncan Murdoch
On 18/09/2008 10:55 AM, Henrik Bengtsson wrote: R-core, may I suggest that the error message returns the absolute path (or even the relative path) and not just the filename/basename, e.g. Fatal error: unable to restore saved data in C:/Users/foo/.RData That's a good suggestion, but

Re: [R] Reporting with down and across variables

2008-09-18 Thread hadley wickham
Possibly - in the next version, you'll get a helpful warning message and cast will take it's best guess at the column that it should use. Hadley On Thu, Sep 18, 2008 at 10:15 AM, Matthew Pettis [EMAIL PROTECTED] wrote: It does, thanks! Did I miss in the documentation that the variable has to

Re: [R] Fatal error .RData

2008-09-18 Thread Rodrigo Aluizio
If the students have the permission, maybe run R as administrator solves the question too. Rodrigo. -- From: Duncan Murdoch [EMAIL PROTECTED] Sent: Thursday, September 18, 2008 10:48 AM To: Green, Paul [EMAIL PROTECTED] Cc: r-help@r-project.org

Re: [R] Problem installing packages in newer versions of R

2008-09-18 Thread Peter Dalgaard
Miha Staut wrote: Dear all, I was wandering what could be wrong with my system (regularly updated Fedora core 8) so that installing packages does not succeed with almost every package. I follow the procedure specified in the help file R-admin section 6.3. This is not a feature new to the

[R] propensity score adjustment using R

2008-09-18 Thread Bunny, lautloscrew.com
Hi all, i am looking to built a simple example of a very basic propensity score adjustment, just using the estimated propensity scores as inverse probability weights (respectively 1-estimated weights for the non-treated). As far as i understood, MLE predictions of a logit model can

[R] Odp: Fatal error .RData

2008-09-18 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 18.09.2008 15:41:07: We are using R in our class. One student claims to be receiving the message Fatal error: unable to restore saved data in .RData Well, throw my suggestion away. I sometimes get such message when I save .RData with some

Re: [R] rgl: How to position a window during open3d call

2008-09-18 Thread Duncan Murdoch
On 18/09/2008 12:08 AM, Yihui Xie wrote: Thanks Duncan, I have also been wondering about this problem for a long time. This is fixed now on R-forge; eventually it will make it into the next rgl release on CRAN. You should be able to download a binary of the development version from R-forge

Re: [R] propensity score adjustment using R

2008-09-18 Thread Ben Domingue
Bunny, lautloscrew.com bunny at lautloscrew.com writes: ix of some covariates. I wonder right now if te glm respectively summary(glm(...)) puts out something comparable to ML estimates that can be used as the estimated pscores, in such a way that there is one value for every observation.

Re: [R] propensity score adjustment using R

2008-09-18 Thread Frank E Harrell Jr
Bunny, lautloscrew.com wrote: Hi all, i am looking to built a simple example of a very basic propensity score adjustment, just using the estimated propensity scores as inverse probability weights (respectively 1-estimated weights for the non-treated). As far as i understood, MLE predictions

Re: [R] Problem installing packages in newer versions of R

2008-09-18 Thread Miha Staut
Thanks for your help Hadley and Peter. Your suggestions solved the question of installing new packages. But what can be done to get rid of the other mentioned problem in html help files. I am repeating what the problem is. After installation of new packages R-html help files try to follow

Re: [R] Suggestion: Re: Fatal error .RData

2008-09-18 Thread Henrik Bengtsson
Hi. On Thu, Sep 18, 2008 at 8:16 AM, Duncan Murdoch [EMAIL PROTECTED] wrote: On 18/09/2008 10:55 AM, Henrik Bengtsson wrote: R-core, may I suggest that the error message returns the absolute path (or even the relative path) and not just the filename/basename, e.g. Fatal error: unable to

[R] PDF fonts problem

2008-09-18 Thread Mihalicza Péter
Dear List, I am writing a paper in Hungarian, that I Sweave and than pdfLaTeX. Everything is fine, except for two accented letters in the graphs that behave strange, though on the screen and in eps exports they look perfect. The problem is that I need pdf graphs, since I would like to have a

Re: [R] Suggestion: Re: Fatal error .RData

2008-09-18 Thread Duncan Murdoch
On 18/09/2008 12:13 PM, Henrik Bengtsson wrote: Hi. On Thu, Sep 18, 2008 at 8:16 AM, Duncan Murdoch [EMAIL PROTECTED] wrote: On 18/09/2008 10:55 AM, Henrik Bengtsson wrote: R-core, may I suggest that the error message returns the absolute path (or even the relative path) and not just the

Re: [R] Suggestion: Re: Fatal error .RData

2008-09-18 Thread Richard M. Heiberger
Isn't the full pathname just paste(getwd(), .RData, sep=/) Rich __ 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,

Re: [R] R2HTML, HTMLStart(): are plots send to HTML?

2008-09-18 Thread Greg Snow
You need to tell R when to include the plots (there is the possibility that you will add to the current plot, or redo it, etc., so R cannot tell when to include plots automatically). Look at the HTMLplot function, you will issue that when you have made a plot (see it on the screen) that you

Re: [R] propensity score adjustment using R

2008-09-18 Thread ran2
Frank E Harrell Jr wrote: That is a high variance procedure as compared with covariate adjustment using the propensity score, or stratification. Frank Harrell Yes, I guess the foo$fitted.values was the syntax i missed. I know this method is not optimal and that it yields high

[R] Adding 1 month to a dataframe column

2008-09-18 Thread ANGELO.LINARDI
Dear R experts, I have a problem in modifying one column of a dataframe with a datatime format using a datetime operator. Here is my dataframe A: DATACONT PROVINCIA VALORE 1 2007-12-31MI 1 2 2007-12-31PV 2 3 2007-12-31NA 3 4 2007-12-31MI

Re: [R] Spectrum of a kernel

2008-09-18 Thread stephen sefick
every value is the same 0.004975 besides you have to apply the smooth to a time series what you get when you plot plot(kernel(daniell, 100)) it isn't in the time domain because there is no association with a time just values, I think. On Thu, Sep 18, 2008 at 11:18 AM, [EMAIL PROTECTED] wrote:

Re: [R] Writing greek letters and subscripts in graph

2008-09-18 Thread Mike Prager
Delphine COURVOISIER [EMAIL PROTECTED] wrote: I would like to write Greek letters followed by subscripts in a graph (on the X-axis and in a legend). I would appreciate any help. As others have said, try demo(plotmath) and also ? plotmath Here's a brief example: plot(1:10, 1:10)

Re: [R] Coxph and loglik converged before variable X

2008-09-18 Thread Greg Snow
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] project.org] On Behalf Of Terry Therneau Sent: Thursday, September 18, 2008 6:49 AM To: [EMAIL PROTECTED] Cc: r-help@r-project.org Subject: Re: [R] Coxph and loglik converged before variable X [snip] b. The

Re: [R] Use of colour in plots

2008-09-18 Thread Greg Snow
Try something like: x - runif(25) y - rnorm(25) z - rnorm(25, 3*x) plot(x, y, col=ifelse( z 1.25, 'red', 'green') ) Does this help, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] 801.408.8111 -Original Message- From: [EMAIL

Re: [R] Suggestion: Re: Fatal error .RData

2008-09-18 Thread Duncan Murdoch
On 18/09/2008 12:23 PM, Richard M. Heiberger wrote: Isn't the full pathname just paste(getwd(), .RData, sep=/) Only on Unix-alikes, and the code that prints the error message is platform neutral. And at the point where we need this message, R isn't running, so we'd need this in C code.

Re: [R] propensity score adjustment using R

2008-09-18 Thread ran2
Frank E Harrell Jr wrote: That is a high variance procedure as compared with covariate adjustment using the propensity score, or stratification. Frank Harrell Ah, wait what if I got very high dimensional X ? Even with 20 binary covariates i would end up with more than 1 million

[R] Difficulty understanding sem errors / failed confirmatory factor analysis

2008-09-18 Thread Adam D. I. Kramer
Hello, I'm trying to fit a pretty simple confirmatory factor analysis using the sem package. There's a CFA example in the examples, which is helpful, but the output for my (failing) model is hard to understand. I'd be interested in any other ways to do a CFA in R, if this proves

Re: [R] Use of colour in plots

2008-09-18 Thread Steve Murray
Dear Thierry and all, I've tried out ggplot from the ggplot2 package and it seems to provide much more favourable results! Just a few questions I have after consulting the 'help' file for ggplot. Is there a way of preventing overplotting? Some of the red points are being obscured by the

[R] detecting null values in a CSV file

2008-09-18 Thread Jason Thibodeau
Hello all, I have a CSV file, that is 2411 columns wide. There are certain instances in teh file, where null values are located. That is: two commas together, without anything in the middle. In a certain section, the only possible values are NULL, 0,1,and 2. I need to be able to detect these

[R] about the whisker in boxplot

2008-09-18 Thread catherineLF
Hi, Dear R-users, I have a problem when I drawing a boxplot. I want to extend the whisker to the 5% and the 95% quantiles and only show the most extreme outlier, like 0.01% and 99.99% percentiles. What should I do? I saw something on boxplot.stat, but even I define the parameter in

[R] Removing subtitles in levelplot

2008-09-18 Thread Sebastian Weirich
Hello, probable a simple question, but I’m unaware with it. When plotting something like levelplot(z1+z1~x*y|g), each table has a subtitle containing variable name of z1 respectively z2 and the value of g. How can I remove the subtitles? Thanks a lot, Sebastian -- Ist Ihr Browser

Re: [R] Creating smooth color regions with panel.contourplot()

2008-09-18 Thread Waichler, Scott R
Thank you very much, Deepayan. There is just one more feature I'd like to get, the ability to add the contour lines. My revision to your code below prints too many lines. What needs to be changed? --Thanks, Scott Waichler [EMAIL PROTECTED] library(gridBase) library(lattice) data(volcano)

Re: [R] detecting null values in a CSV file

2008-09-18 Thread Henrik Bengtsson
What have you tried this far? Can't you parse them as missing values, i.e. NAs? See ?read.csv and arguments '...', i.e. the arguments '...' are passed to read.table() which takes argument 'na.strings' - a character *vector* of strings that you want to be interpreted as NAs. See ?read.table for

Re: [R] propensity score adjustment using R

2008-09-18 Thread Ben Domingue
I'm not quite sure what you mean. If all you need is propensity scores to run an IPW analysis, the fitted values should work. Having many binary covariates shouldn't be a problem, the whole point of the propensity score is boiling down many dimensions to a single one. I use matchit() for my psm

Re: [R] detecting null values in a CSV file

2008-09-18 Thread Jason Thibodeau
I haven't given it a shot yet. To complicate matters further, this file I have already passed through a filter, which has already converted my null values to NA. your insight might be of assistance. On Thu, Sep 18, 2008 at 1:30 PM, Sebastian Weirich [EMAIL PROTECTED]wrote: Hello, maybe you

Re: [R] Difficulty understanding sem errors / failed confirmatory factor analysis

2008-09-18 Thread Adam D. I. Kramer
No new info, but the model and correlation table are pasted at the end of this message. --Adam On Thu, 18 Sep 2008, Adam D. I. Kramer wrote: Hello, I'm trying to fit a pretty simple confirmatory factor analysis using the sem package. There's a CFA example in the examples, which is

Re: [R] detecting null values in a CSV file

2008-09-18 Thread Jason Thibodeau
Phil's suggestion worked like a charm. My NA's were counted in the frequency table. Thanks for the help, all! On Thu, Sep 18, 2008 at 1:28 PM, Henrik Bengtsson [EMAIL PROTECTED]wrote: What have you tried this far? Can't you parse them as missing values, i.e. NAs? See ?read.csv and

[R] Ability estimates for partial credit model

2008-09-18 Thread Andrej Kastrin
Dear all, I'm working on ability estimates using Rasch model. Using the ltm package, the procedure is quite simple: ## Factor Scores for the Rasch model fit - rasch(LSAT) factor.scores(fit) What about Partial Credit Model (PCM)? For PCM I use PCM function from eRm package. Is there any

Re: [R] propensity score adjustment using R

2008-09-18 Thread ran2
Ben Domingue wrote: I'm not quite sure what you mean. If all you need is propensity scores to run an IPW analysis, the fitted values should work. Having many binary covariates shouldn't be a problem, the whole point of the propensity score is boiling down many dimensions to a single

Re: [R] Difficulty understanding sem errors / failed confirmatory factor analysis

2008-09-18 Thread John Fox
Dear Adam, (1) Note that your input correlation matrix appears to be numerically singular: solve(R) Error in solve.default(R) : system is computationally singular: reciprocal condition number = 2.38183e-17 det(R) [1] -1.753523e-25 qr(R)$rank [1] 23 (2) In addition, you have specified what

Re: [R] Difficulty understanding sem errors / failed confirmatory factor analysis

2008-09-18 Thread Adam D. I. Kramer
Dear John, On Thu, 18 Sep 2008, John Fox wrote: I'm trying to fit a pretty simple confirmatory factor analysis using the sem package. There's a CFA example in the examples, which is helpful, but the output for my (failing) model is hard to understand. I'd be interested in any other

Re: [R] detecting null values in a CSV file

2008-09-18 Thread Hutchinson,David [PYR]
Try length(na.omit(the particular data column)) Here's an example: data - runif(100,0,10) data[runif(20,0,100)] - NA file.contents - matrix(data, ncol = 5, byrow = TRUE) for (i in 1:5) { print (length(na.omit(file.contents[,i]))) } -Original Message- From: [EMAIL PROTECTED]

[R] caret package: arguments passed to the classification or regression routine

2008-09-18 Thread Peter Tait
Hi, I am having problems passing arguments to method=gbm using the train() function. I would like to train gbm using the laplace distribution or the quantile distribution. here is the code I used and the error: gbm.test - train(x.enet, y.matrix[,7], method=gbm,

[R] outlier and whisker in boxplot

2008-09-18 Thread catherineLF
Hi, Dear R-users: Sorry for bothering your guys again. I think I should rewrite my question. I know how to extend whisker by using range. The question is that I will set the range=1.5, and at the same time, I only want to show the extreme outlier, like 0.01% and 99.99% percentile, so what

[R] graphing netCDF files

2008-09-18 Thread Steve_Friedman
Hello I'm working with a large hydrological data set stored in a netCDF format. The file stores x and y coordinates in the UTM projected coordinate system, yet when I use image to graphically display the z variable, the image is distorted in the sense that it does not plot the map in the correct

Re: [R] propensity score adjustment using R

2008-09-18 Thread Frank E Harrell Jr
ran2 wrote: Frank E Harrell Jr wrote: That is a high variance procedure as compared with covariate adjustment using the propensity score, or stratification. Frank Harrell Ah, wait what if I got very high dimensional X ? Even with 20 binary covariates i would end up with more than 1

Re: [R] Difficulty understanding sem errors / failed confirmatory factor analysis

2008-09-18 Thread John Fox
Dear Adam, I'm afraid that our emails have crossed. Please see my previous message. John -- John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada web: socserv.mcmaster.ca/jfox -Original Message- From: Adam D. I. Kramer

  1   2   >