Re: [R] creating a variable using concatenation

On Mar 30, 2010, at 11:15 PM, Steve Chen wrote: You can try this: library(tutoR) That's an interesting set of goals: Package to mask common functions so that inputs in error are explained and able to be corrected, prior to execution. 'assist' offers step-by-step assistance to correctly

[R] Predict function in ARIMA/GARCH- Real time forecastinng

Hi, I m working on time series forecasting.I was checking ARIMA and GARCH model for forecasting.Is there any possibility that we can include test data as input for predict function.I am working on real time forecasting where i would be getting real time feeds using which i need

Re: [R] nlsList{nlme} control of min and max parameter bounds?

Steve Oswald wrote: I'm analysing growth rates using a gompertz (logistic) curve and am attempting to fit parameters for all of my study birds using nlsList. I've been looking for an option in nlsList to set min and max limits for parameter values during estimation, although I have

Re: [R] Install package of BRugs?

Hua Li wrote: I'm trying to install the package BRugs and could not find it on the install package list. Then I tried the suggestion to run the following command in R install.packages(BRugs) Use R2WinBUGS instead http://cran.at.r-project.org/web/packages/R2WinBUGS/index.html

Re: [R] Install package of BRugs?

On 30.03.2010 22:01, David Winsemius wrote: On Mar 30, 2010, at 3:42 PM, Hua Li wrote: Hi All, I'm trying to install the package BRugs and could not find it on the install package list. Then I tried the suggestion to run the following command in R install.packages(BRugs) but get the

Re: [R] Install package of BRugs?

On 31.03.2010 08:36, Dieter Menne wrote: Hua Li wrote: I'm trying to install the package BRugs and could not find it on the install package list. Then I tried the suggestion to run the following command in R install.packages(BRugs) Use R2WinBUGS instead

Re: [R] Install package of BRugs?

Uwe Ligges-3 wrote: Use R2WinBUGS instead Dieter, thanks for the pointer, but I highly recommend *not* to use R2WinBUGS if you are under Windows but use BRugs which has a much cleaner interface and . Ahem... yes. I am aware of this, and I had mixed things up! Dieter --

[R] decomposing the error term into time-specific indiviual specific effect + transitory term

hi all! i am a beginner in R so i need a command that i couldnt find on the net. i need to transform the error term and i should decompose it into time-specific, indiviual specific effect and transitory term like here (åit =ui + vt + wit) the individual effect in my model is called stno and

Re: [R] is and the story of a typo

On 30-Mar-10 23:23:09, Jim Lemon wrote: Hi all, The gurus may pour scorn on me for not knowing this, but I happened to mistype if as is in the heat of debugging a function. As I scanned the debugged function with some satisfaction, I noticed the error. How could this have worked? I assume

Re: [R] MySQL and RODBC - limitations

Thanks guys - that is better. You all rock! -- View this message in context: http://n4.nabble.com/MySQL-and-RODBC-limitations-tp1692743p1746457.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] Weird R behaviour?

Dear list, I have observed a weird behaviour from R --- apologies if I am missing something obvious! df3f826f28 df3f826f28 Say you type in R: c.preec - 10074 c.gd - 2200 p1 - .2 c.neo - p1*9451 + (1-p1)*3883 n.preec - 3710 n.gd - 2650 n.neo - 2120 n.pcos - 53000 unit.met - 94 cost.met -

Re: [R] Weird R behaviour?

FAQ 7.31: Why doesn't R think these numbers are equal? http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f b On Wed, Mar 31, 2010 at 10:10 AM, Gianluca Baio gianl...@statistica.it wrote: Dear list, I have observed a weird behaviour from R ---

Re: [R] Create a new variable

Hi r-help-boun...@r-project.org napsal dne 30.03.2010 18:46:35: Hello, Thomas Jensen wrote: Dear R-list, Sorry for spamming the list lately, I am just learning the more advanced aspects of R! I have some data that looks like this: Out Country1 Country 2 Country 3 ...

[R] Odp: (no subject)

Hi r-help-boun...@r-project.org napsal dne 31.03.2010 06:15:38: Hi r-users, I would like to display my legend with fill box and line symbol. The first color will be in the box and the second colour will just be a line. legend(topright, c(observed,gamma sum fit), col= c

[R] Truncated Distributions - Estimating parameters

Dear R users, I am working on the Value at Risk (VaR) for the Operational risk. For a given loss data, we try to fit some statistical distributions using Kolmogorov-Smirnov (KS) test and A-D test and then for fitted distribution using the estimated parameters, the losses are simulated and the

Re: [R] (no subject)

On 03/31/2010 03:15 PM, Roslina Zakaria wrote: Hi r-users, I would like to display my legend with fill box and line symbol. The first color will be in the box and the second colour will just be a line. legend(topright, c(observed,gamma sum fit), col= c(greenyellow,red),pch=c(15,

[R] Simplifying particular piece of code

Hello, everyone I have a piece of code that looks like this: mrets - merge(mrets, BMM.SR=apply(mrets, 1, MyFunc, ret=BMM.AV120, stdev=BMM.SD120)) mrets - merge(mrets, GM1.SR=apply(mrets, 1, MyFunc, ret=GM1.AV120, stdev=GM1.SD120)) mrets - merge(mrets, IYC.SR=apply(mrets, 1, MyFunc,

Re: [R] Simplifying particular piece of code

How about this (not tested, since you did not provide example data nor function code): --- SRnames - paste(colnames.mrets, .SR, sep=) AVnames - paste(colnames.mrets, .AV120, sep=) SDnames - paste(colnames.mrets, .SD120, sep=)

[R] time series datasets

Hello, I am searching large multivariate time series datasets to use with R. Can any one help me??? Regards, Vibha [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Generative Topographic Map

I tried to use R version of package I noticed the original MatLab Pckage is much better documented. I had a look at the R demo code gtm_demo and found that variable Y is used in advanced of being created: I wrote my own few lines as follows: inDir - C:/Documents and Settings/Monville/Alanine

Re: [R] Value-at-Risk Portfolio(both equity and option)

Thank you very much. It is very helpful. As far as I understand, not easy to have a function to combine both the equity part and the option part? -- View this message in context: http://n4.nabble.com/Value-at-Risk-Portfolio-both-equity-and-option-tp1745179p1746520.html Sent from the R help

[R] Error singular gradient matrix at initial parameter estimates in nls

If you have a perfect fit, you have zero residuals. But in the nls manual page we have: Warning: *Do not use ‘nls’ on artificial zero-residual data.* So this is a case of complaining that your diesel car is broken because you ignored the Diesel fuel only sign on the filler cap and put

[R] interpretation of p values for highly correlated logistic analysis

Dear list, I want to perform a logistic regression analysis with multiple categorical predictors (i.e., a logit) on some data where there is a very definite relationship between one predicator and the response/independent variable. The problem I have is that in such a case the p value goes very

[R] RCOM Save

Dear all, I use the package RCOM to manipulate Excel-Files (Windows). After opening the Excel-File and formatting the cells I want to save the file automatically without request. If I use y[[Save]] - TRUE it opens the saving dialog and I have to select the destination for saving. But I simply

Re: [R] RCOM Save

Try adding this: y[[DisplayAlerts]] - FALSE On Wed, Mar 31, 2010 at 8:28 AM, koj jens.k...@gmx.li wrote: Dear all, I use the package RCOM to manipulate Excel-Files (Windows). After opening the Excel-File and formatting the cells I want to save the file automatically without request. If I

[R] ggplot2: Adding points to a density plot

Hi, Consider something like library(ggplot2) movies$decade - round_any(movies$year, 10) m - qplot(rating,data=movies,colour=factor(decade),geom=density) m (modified from ?stat_density). I'd like to add on the line y=0 a dot for the median of each decade category (using the same colour

Re: [R] interpretation of p values for highly correlated lo gistic analysis

claus orourke claus.orourke at gmail.com writes: look up Hauck-Donner effect ... __ 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

[R] Printing the function t.test() in R

Dear R Users, I have a question, how does one print an r function like t.test in R (to get the whole function and not just a summary of the environment)? For example if I type the following: t.test function (x, ...) UseMethod(t.test) environment: namespace:stats print(t.test) function (x,

Re: [R] Value-at-Risk Portfolio(both equity and option)

One thing that comes to mind immediately is the 'blotter' package which I believe is designed to handle and track multiple instruments. A sloppy temporary solution I guess would be to do some type of lookup or a flag that denotes which instrument is an option and apply the VaR methodology

Re: [R] Printing the function t.test() in R

Hi Luwis, this task is solved in the R Help Desk, see R News 6/4 http://cran.r-project.org/doc/Rnews/Rnews_2006-4.pdf in short t.test is a S3-method and you can have a look at implemented methods with methods(t.test) and then use e.g. getAnywhere(t.test.default) hth. Luwis Tapiwa Diya

Re: [R] time series datasets

What kind of time series data? Finance? Econometrics? The 'quantmod' package allows for downloading of stock and FRED data online for analysis. I believe PerformanceAnalytics also has some datasets as well. -c On 3/31/10, vibha patel vibhapatel...@gmail.com wrote: Hello, I am searching large

Re: [R] Printing the function t.test() in R

On Mar 31, 2010, at 8:35 AM, Luwis Tapiwa Diya wrote: Dear R Users, I have a question, how does one print an r function like t.test in R (to get the whole function and not just a summary of the environment)? For example if I type the following: t.test function (x, ...)

Re: [R] Problem comparing hazard ratios

Thank you, gentlemen. I greatly appreciate your help. -- Michal J. Figurski, PhD HUP, Pathology Laboratory Medicine Biomarker Research Laboratory 3400 Spruce St. 7 Maloney Philadelphia, PA 19104 tel. (215) 662-3413 __ R-help@r-project.org mailing

[R] position of mismatches in all.equal()?

Dear R users, I would like to compare two dataframes, actually their categorical variables (as factors) only (there are 12, from column 1 to 12). The reason I do that is that I got 2 datasets from two different methods and I would like to be sure that each method used the same data (3D

Re: [R] Error singular gradient matrix at initial parameter estimates in nls

Dear JN, Bert, 1) It is not a perfect fit. I do not think I have ever said that. I said that an external algorithms fits the model without any problems: with ~ 500,000 data points and 19 paramters (ki in the original equation), it fits the model in less than 1 second. The data are not

Re: [R] ggplot2: Adding points to a density plot

Hi: Try this: library(ggplot2) movmed - ddply(movies, .(decade), summarise, med = median(rating)) m + geom_point(data = movmed, aes(x = med), y = 0, size = 2) HTH, Dennis On Wed, Mar 31, 2010 at 4:46 AM, Johannes Graumann johannes_graum...@web.de wrote: Hi, Consider something like

[R] regular expression help to extract specific strings from text

Dear all, Lets say I have the following: x - c(Eve: Going to try something new today..., Adam: Hey @Eve, how are you finding R? #rstats, Eve: @Adam, It's awesome, so much better at statistics that #Excel ever was! @Cain @Able disagree though :(, Adam: @Eve I'm sure they'll sort it out

Re: [R] for loop; lm() regressions; list of vectors - lapply - accolades and square brackets??

Hello and thank you both for your answers! Dennis, I tried to simply run lm(a ~ b) after re-importing a as a matrix, but I get the following error message: Error in model.frame.default(formula = a ~ b, drop.unused.levels = TRUE) : invalid type (list) for variable 'a' so maybe I have to

[R] Installing JGR and rJava

Hello, I am running: R version 2.10.0 (2009-10-26) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 on a 64-bit RedHat box. To encourage collegues to use R, I am trying to install JGR, which depends on rJava, on the above machine. At first I received and error

Re: [R] for loop; lm() regressions; list of vectors - lapply - accolades and square brackets??

On Mar 31, 2010, at 9:13 AM, Driss Agramelal wrote: Hello and thank you both for your answers! Dennis, I tried to simply run lm(a ~ b) after re-importing a as a matrix, but I get the following error message: Error in model.frame.default(formula = a ~ b, drop.unused.levels = TRUE)

Re: [R] Error singular gradient matrix at initial parameter estimates in nls

Try the function called `nls.lm' which is contained in the minpack.lm package. See if it solves your problem. Ravi. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Corrado Sent: Wednesday, March 31, 2010 9:13 AM Cc:

[R] (no subject)

The answer depends on what your goals are. If you are simply trying to put all the hazard ratios on a similar scale, then I often use the HR comparing risk of a subject at the 25th percentile to one at the 75th percentile. This has already been suggested by David W. But this is not perfect.

Re: [R] ggplot2: Adding points to a density plot

Thanks a lot! This got me started! Joh Dennis Murphy wrote: Hi: Try this: library(ggplot2) movmed - ddply(movies, .(decade), summarise, med = median(rating)) m + geom_point(data = movmed, aes(x = med), y = 0, size = 2) HTH, Dennis On Wed, Mar 31, 2010 at 4:46 AM, Johannes

Re: [R] Problem comparing hazard ratios

Frank, Is there an article that discusses this idea of bootstrapping the ranks of the likelihood ratio chi-square Statistics to assess relative importance of predictors in time-to-event data (specifically Cox PH model)? Thanks, Ravi. -Original Message- From:

Re: [R] Error grid must have equal distances in each direction

M Joshi, I don't know but I guess that some might have looked at your previous thread on 14 March (also about the geoR package). You received help and good advice then, but it doesn't appear that you are following it. It appears to be a similar problem this time. Also, this list is the wrong

Re: [R] regular expression help to extract specific strings from text

strapply in gsubfn can extract matches based on content which seems to be what you want: library(gsubfn) f - function(...) sapply(list(...), paste, collapse = , ) DF - data.frame(x, Source = strapply(x, ^(\\w+):, c, simplify = f), Mentions = strapply(x, @(\\w+), c, simplify =

[R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

When using generalized sets from sets library A EMPTY does not work, but gset_intersection(A,EMPTY) works: example code below Is it a bug? Best regards, Ryszard library(sets) support - universe - c('a','b','c','d') A - gset(support=support, memberships=c(0.2, 0.2, 0.9, 0), universe=universe)

Re: [R] regular expression help to extract specific strings from text

On Wed, Mar 31, 2010 at 8:20 AM, Tony B tony.bre...@googlemail.com wrote: Dear all, Lets say I have the following: x - c(Eve: Going to try something new today..., Adam: Hey @Eve, how are you finding R? #rstats, Eve: @Adam, It's awesome, so much better at statistics that #Excel ever was!

Re: [R] Question about 'logit' and 'mlogit' in Zelig

Abraham, This appears to be your 3rd unanswered post to r-help in March, all 3 have been about the Zelig package. Please read the posting guide and find out the correct place to send questions about packages. Then you might get an answer. HTH Matthew Mathew, Abraham T amat...@ku.edu wrote

[R] tm package- remove stowords failling

Hi, I just noticed that by inspecting the matrix term that no all stopwords are removed, does someone know how to fix that? library(tm) data(crude) d-tm_map(crude, removeWords, stopwords(language='english')) dt-DocumentTermMatrix(d,control=list(minWordLength=3, minDocFreq=2)) inspect( dt) I am

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

Based on your success I thought I would try again. I am not sure why I had more success this time but this time I got this far even though I did not change my path at all or make any system changes from what I had before. I did reinstall Rgraphviz but used the previously installed graphviz. I

Re: [R] zero standard errors with geeglm in geepack

You may not have got an answer because you posted to the wrong place. Its a question about a package. Please read the posting guide. miriza miri...@sfwmd.gov wrote in message news:1269886286228-1695430.p...@n4.nabble.com... Hi! I am using geeglm to fit a Poisson model to a timeseries of

Re: [R] GEE for a timeseries of count (one cluster)

Contact the authors of those packages ? miriza miri...@sfwmd.gov wrote in message news:1269981675252-1745896.p...@n4.nabble.com... Hi! I was wondering if there were any packages that would allow me to fit a GEE to a single timeseries of counts so that I could account for autocorrelation

Re: [R] for loop; lm() regressions; list of vectors - lapply - accolades and square brackets??

Driss, David is right - you should include data/code. Nevertheless, the error you get suggests that what you call a matrix is in fact a data.frame; that's usually a crucial distinction. There are two things you might try: Dennis' way: lm(as.matrix(a) ~ b) David's way: for(i in

Re: [R] mcmcglmm starting value example

Apparently not, since this your 3rd unanswered thread to r-help this month about this package. Please read the posting guide and find out where you should send questions about packages. Then you might get an answer. ping chen chen1984...@yahoo.com.cn wrote in message

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

By the way, just in case you did not read the entire message R crashed when I tried to run the code from the vignette. On Wed, Mar 31, 2010 at 10:52 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: Based on your success I thought I would try again. I am not sure why I had more success this

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

Unless I'm missing something, I don't see any method in pkg:sets for intersection other than gset_intersection. So you're using the base R function `` whose help page tells you that its arguments should be vectors. Yours aren't. -Peter Ehlers On 2010-03-31 8:50, Czerminski, Ryszard wrote:

Re: [R] Simplifying particular piece of code

On Wed, Mar 31, 2010 at 5:11 PM, Sergey Goriatchev serg...@gmail.com wrote: but data - merge(data,data.list) works. Neither data or data.list is a list, so do.call does not work. I am very weak on lists, never used them before Best, Sergey Hej Sergey, Ok; I was wondering if the apply

[R] predict.Arima: warnings from xreg magic

When I run predict.Arima in my code, I get warnings like: Warning message: In cbind(intercept = rep(1, n), xreg) : number of rows of result is not a multiple of vector length (arg 1) I think this is because I'm not running predict.Arima in the same environment that I did the fit, so the

Re: [R] Simplifying particular piece of code

Thanks for help, Gustaf! (Kan man säga man tackar oxa?) :-) On Wed, Mar 31, 2010 at 17:34, Gustaf Rydevik gustaf.ryde...@gmail.com wrote: On Wed, Mar 31, 2010 at 5:11 PM, Sergey Goriatchev serg...@gmail.com wrote: but data - merge(data,data.list) works. Neither data or data.list is a

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

On 2010-03-31 9:30, Peter Ehlers wrote: Unless I'm missing something, I don't see any method in pkg:sets for intersection other than gset_intersection. Whoops, a bit quick on the draw. There are of course also set_intersection and cset_intersection, but not AFAICS any method for ``. -Peter

[R] Boundary kernel for kernel density plot?

Hello, Is a kernel called the boundary kernel available for kernel density plots? In the help, I can only see the following kernels: gaussian, rectangular, triangular, epanechnikov, biweight, cosine or optcosine. The boundary kernel apparently improves the estimate in a neighborhood of zero.

Re: [R] position of mismatches in all.equal()?

I've actually found a solution to my problem. In case someone could be interested, the which() function seems to do the trick: test2$let2 - factor(test2$let2, levels=letters[1:10]) which(test1[,3]!=test2[,3]) [1] 7 which(test1[,1:3]!=test2[,1:3]) [1] 27 Le 3/31/2010 15:13, Ivan Calandra

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

It seems that A B works the same way as gset_intersection(A,B) as long as A and B are not empty... see below: fuzzy_logic(Yager, p=2) support - universe - c('a','b','c','d') A - gset(support=support, memberships=c(0.2, 0.2, 0.9, 0), universe=universe) B - gset(support=support,

Re: [R] Error singular gradient matrix at initial parameter estimates in nls

Corrado: Over parameterization/non-identifiability is not determined by the ratio of the number of data values to the number of parameters: if you try to fit a scatter of a zillion points that lie near a straight line to a model with curvature -- the 4p logistic function, say -- you're

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

On 03/31/2010 08:18 AM, Gabor Grothendieck wrote: By the way, just in case you did not read the entire message R crashed I think you mean that you got an error message, not that R became unusable? More below... when I tried to run the code from the vignette. On Wed, Mar 31, 2010 at 10:52

Re: [R] GLM / large dataset question

Geelman, This appears to be your first post to this list. Welcome to R. Nearly 2 days is quite a long time to wait though, so you are unlikely to get a reply now. Feedback : the question seems quite vague and imprecise. It depends on which R you mean (32bit/64bit) and how much ram you have.

Re: [R] dot plot

Hi kayj kayj schrieb: Also, I do nto know how to combine two plots into one plot? Have a look at ?par . There is some argument add or new (i am not sure) that you have to set to TRUE after making the first plot and then you can add a second plot with running plot() once more. thanks and

Re: [R] Off Topic: teenie weenie numbers -- Was: Precision level

I spoke with a theoretical physicist and he said he encountered 10^-120. Has something to do with attempts to describe/explain the universe... Dimitri On Fri, Mar 26, 2010 at 2:33 PM, Steve Lianoglou mailinglist.honey...@gmail.com wrote: On Fri, Mar 26, 2010 at 1:11 PM, Barry Rowlingson

[R] Comparing times in two time zones

I have data that is collected in two different time zones z$OnDateTime - as.POSIXct(runway$OnDateTime, tz = EST5EDT, format=%m/%d/%Y %H:%M) is in Eastern time, and zi$ActualOnLocal - as.POSIXct(oooi$ActualOnLocal, tz=MST7MDT, format=%m/%d/%Y %H:%M) is in Mountain. I converted the runway time to

Re: [R] library sets: A EMPTY does not work; gset_intersection(A,EMPTY) works

Ryszard, You've made me take a closer look and now I do think that you've found a bug. After a quick look at the package vignette, I see that the authors have indeed overloaded and so it should work for your example. The problem seems to be the order of the class attribute which is used to

[R] Bar plots with bars made of stacked text

I would like to make bar plots where the bars are composed of text like this: http://www.thomaslevine.com/lowres/text_bars.png Is there a package that will help me with this? Thanks Tom __ R-help@r-project.org mailing list

Re: [R] interpretation of p values for highly correlated logistic analysis

Hi Claus, welcome to the wonderful world of collinearity (or multicollinearity, as some call it)! You have a near linear relationship between some of your predictors, which can (and in your case does) lead to extreme parameter estimates, which in some cases almost cancel out (a coefficient of

[R] reshaping data

Hi All, Can someone help me reshape following data: Var1 Var2 Val A X 1 A Y 2 A Z 3 B X 4 B Y 5 B Z 6 to some kind of matrix/tabular format (preferably as a matrix), may be like Var1 X Y Z A 1 2 3 B 4 5 6 Any help would be greatly appreciated, Kim [[alternative HTML version deleted]]

Re: [R] reshaping data

Try this: xtabs(Val ~ Var1 + Var2, data = x) On Wed, Mar 31, 2010 at 3:23 PM, Kim Jung Hwa kimhwamaill...@gmail.com wrote: Hi All, Can someone help me reshape following data: Var1 Var2 Val A X 1 A Y 2 A Z 3 B X 4 B Y 5 B Z 6 to some kind of matrix/tabular format (preferably as a

Re: [R] reshaping data

Hi Kim, look at the reshape() command with direction=wide. Or at the reshape package. HTH, Stephan Kim Jung Hwa schrieb: Hi All, Can someone help me reshape following data: Var1 Var2 Val A X 1 A Y 2 A Z 3 B X 4 B Y 5 B Z 6 to some kind of matrix/tabular format (preferably as a matrix),

Re: [R] for loop; lm() regressions; list of vectors - lapply - accolades and square brackets??

Peter, it works!!! Thanks a lot to all of you! Yes I mistook a data.frame for a matrix Now both David's and Dennis' solutions work! I am really a beginner with R, and I think it's a wonderful tool, but sometimes a little difficult to use only learning from the manuals and help files... I

Re: [R] reshaping data

Thanks Henrique and Stephan for your reply! Henrique, I'm planning to do some arthitmetic operations on tranformed (matrix) data and then would like to convert it back to original format (as a data frame)... is there an equivalent easy command for this too? Thanks, On Wed, Mar 31, 2010 at 2:26

Re: [R] reshaping data

Actually, apart from melt() in reshape package. On Wed, Mar 31, 2010 at 2:37 PM, Kim Jung Hwa kimhwamaill...@gmail.comwrote: Thanks Henrique and Stephan for your reply! Henrique, I'm planning to do some arthitmetic operations on tranformed (matrix) data and then would like to convert it back

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

I got an error message *AND* R becomes unusable and had to be restarted. Using dependency walker it seems to be complaining about R.dll so I copied all of my ...\R\bin\R*.dll files to .. ...\win-library\2.10\Rgraphviz\libs. Then it complained about IESHIMS.dll so I copied \Program Files\Internet

Re: [R] Off Topic: teenie weenie numbers -- Was: Precision level

There is a story (apocryphal?) about Fred Hoyle, many years ago, having come to the close of a public lecture about his work in Cosmology (I seem to recall that it was to the British Astronomical Association, a society of amateur astronomers, and therefore knowledgeable). He was taking questions.

Re: [R] reshaping data

Try: as.data.frame(with(x, tapply(Val, list(Var1, Var2), sum))) On Wed, Mar 31, 2010 at 3:37 PM, Kim Jung Hwa kimhwamaill...@gmail.com wrote: Thanks Henrique and Stephan for your reply! Henrique, I'm planning to do some arthitmetic operations on tranformed (matrix) data and then would like

Re: [R] reshaping data

On Mar 31, 2010, at 2:37 PM, Kim Jung Hwa wrote: Thanks Henrique and Stephan for your reply! Henrique, I'm planning to do some arthitmetic operations on tranformed (matrix) data and then would like to convert it back to original format (as a data frame)... is there an equivalent easy

[R] AHRQ Patient Quality Indicators

Is anyone aware of R code that mimic's AHRQ's SAS code for their Prevention Quality Indicators (PQI)? Don't see it anywhere, but wanted to see if anyone else knew of anything. Many thanks -- View this message in context:

[R] lattice: how to add points to the plot generated by levelplot()?

I am using levelplot() function from pkg:lattice and I want to add some points to the plot generated by levelplot() similarly as in basic R graphic sequence: plot(...); points(...), but it does not work. I was reading documentation for lattice, but so far without much success. I would very much

Re: [R] Problem comparing hazard ratios

Ravi Varadhan wrote: Frank, Is there an article that discusses this idea of bootstrapping the ranks of the likelihood ratio chi-square Statistics to assess relative importance of predictors in time-to-event data (specifically Cox PH model)? Thanks, Ravi. Do require(rms); ?anova.rms and

[R] ODD and EVEN numbers

Hi, anyone here who knows how to determine if an integer is odd or even in R? Thanks. tj -- View this message in context: http://n4.nabble.com/ODD-and-EVEN-numbers-tp1747032p1747032.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

Hi Gabor, I just used your R code below and the code worked properly without any error messages. Attached is the output graph, which may be the one you expected. library(Rgraphviz) set.seed(123) V - letters[1:10] M - 1:4 g1 - randomGraph(V, M, 0.2) g1 plot(g1) Thanks. Howard On Wed

[R] change mbmdr odds ratio

I am using the package 'mbmdr' based on the model-base multifactor dimensionality reduction method by Calle et al. In the paper, mbmdr is used for a case-control study. My data is not case-control, i.e. I have a lot more controls than cases. Does anyone know if the odds ratio in the fuction

[R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

Hi Martin, It is really a 'PATH' problem. After adding C:\Program Files\Graphviz2.20\bin to the 'PATH' environment variable, the Rgraphviz package can be loaded without any error messages. Sorry that I ONLY set my 'PATH' environment variable correctly for R but not for Graphviz. Thank

[R] How to work with the result of JRip?

Hello, I want to work with the result of JRip, but I don't know how to access it. I produce a rule like below: rule = JRip(Cls ~., data=dd) The result is: JRIP rules: === (C7 = 1) = Cls=TRUE (384.0/0.0) = Cls=FALSE (1344.0/0.0) \par Number of Rules : 2 I want to work with the

[R] Summing data based on certain conditions

Dear all, I have a dataset of 1073 rows, the first 15 which look as follows: data[1:15,]     date year month day rammday thmmday 1   3/8/1988 1988 3   8    1.43    0.94 2  3/15/1988 1988 3  15    2.86    0.66 3  3/22/1988 1988 3  22    5.06    3.43 4  3/29/1988 1988 3  29  

Re: [R] You are right and the problem is solved. Re: about the possible errors in Rgraphviz Package

On 03/31/2010 11:44 AM, Gabor Grothendieck wrote: I got an error message *AND* R becomes unusable and had to be restarted. Using dependency walker it seems to be complaining about R.dll so I copied all of my ...\R\bin\R*.dll files to .. ...\win-library\2.10\Rgraphviz\libs. Then it

Re: [R] ODD and EVEN numbers

See %% under ?Arithmetic tj wrote: Hi, anyone here who knows how to determine if an integer is odd or even in R? Thanks. tj __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] lattice: how to add points to the plot generated by levelplot()?

Czerminski, Ryszard wrote: I am using levelplot() function from pkg:lattice and I want to add some points to the plot generated by levelplot() It's best to include a small example of what you have done so far, and what you want to further do. Actual R code with a small sample dataset is

Re: [R] ODD and EVEN numbers

On Mar 31, 2010, at 12:43 PM, tj wrote: Hi, anyone here who knows how to determine if an integer is odd or even in R? Thanks. 2 %% 2 == 0 [1] TRUE 3 %% 2 == 0 [1] FALSE is.even - function(x){ x %% 2 == 0 } is.even(2) [1] TRUE -- David Winsemius, MD West Hartford, CT

Re: [R] Summing data based on certain conditions

?ave or ?tapply Bert Gunter Genentech Nonclinical Biostatistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Steve Murray Sent: Wednesday, March 31, 2010 10:20 AM To: r-help@r-project.org Subject: [R] Summing data based on

[R] strip.custom with strip on left for three conditioning variables

I want to use a strip.custom as with useOuterStrips for three conditioning variables. useOuterStrips restricts this to only two conditioning variables, and I cannot figure out how to write strip.custom properly to do this. library(lattice) mtcars$HP - equal.count(mtcars$hp) #with two factors

Re: [R] ODD and EVEN numbers

(foo %% 2) == 0 See ?%% HTH Stephan tj schrieb: Hi, anyone here who knows how to determine if an integer is odd or even in R? Thanks. tj __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] lattice: how to add points to the plot generated bylevelplot()?

Tip 1: This sort of thing is done with panel functions. ?xyplot and ?panel.xyplot is where to start. Tip 2: grid functions are generally not necessary, as translations of the common ones are part of lattice. See ?panel.points Tip 3: Deepayan's book provides a gentler explanation of these matters

  1   2   >