similar to: This site may harm your computer - Google warning about cran website

Hi, I noticed whether some one could explain why "&" and "&&" behave differently in data frame transformations. Consider the following : a<-data.frame(r=c(0,0,2,3),g=c(0,2,0,2.1)) Then: > transform(a,R=ifelse(r>0 && g> 0,log(r/g),NA)) r g R 1 0 0.0 NA 2 0 2.0 NA 3 2 0.0 NA 4 3 2.1 NA but > transform(a,R=ifelse(r>0 & g>

Hi, I am using GAMs (package mgcv) to smooth event rates in a penalized regression setting and I was wondering if/how one can select the order of the derivative penalty. For my particular problem the order of the penalty (parameter "m" inside the "s" terms of the formula argument) appears to have a larger effect on the AIC/deviance of the estimated model than the

Hello, everyone, I have a question about integration of product of two densities. Here is the sample code; however the mean of first density is a function of another random variable, which is to be integrated. ## f=function(x) {dmvnorm(c(0.6, 0.8), mean=c(0.75, 0.75/x))*dnorm(x, mean=0.6, sd=0.15)} integrate(f, lower=-Inf, upper=Inf) ## error message Error in dmvnorm(c(0.6, 0.8), mean = c(0.75,

Goodmorning, This is a documentation related question about the B-spline function in R. In the help file it is stated that: "df degrees of freedom; one can specify df rather than knots; bs() then chooses df-degree-1 knots at suitable quantiles of x (which will ignore missing values)." So if one were to specify a spline with 6 degrees of freedom (and no intercept) then a basis

Dear R helpers I have some data and through some other software, it is understood that I can fit the Frechet Distribution to it. However, I need to fit the distribution using R code only. I have searched many R packages and one R helper has suggested some sites too, but unfortunately parameters couldn't be estimated. Using LMOM package, I know how to estimate the parameters of Generalized

Hi All, I've been fiddling around with various ways to estimate the popularity of R, SAS, SPSS, Stata, JMP, Minitab, Statistica, Systat, BMDP, S-PLUS, R-PLUS and Revolution R. It's not an easy task. You can see what I've come up with so far at http://r4stats.com/popularity . I'm sure people will have plenty of ideas on how to improve this, so please let me know what you think.

I can run one-sample t-test on an array, for example a matrix myData1, with the following apply(myData1, 2, t.test) Is there a similar fashion using apply() or something else to run 2-sample t-test with datasets from two groups, myData1 and myData2, without looping? TIA, Gang

Dear r-helpers, This is my first time to run survival analysis. Currently, I have a data set which contains two variables, the variable of time to event (or time to censoring) and the variable of censor indicator. For the indicator variable, it was coded as 0 and 1. 0 represents right censor, 1 means event of interest. Now I try to use "survfit" in the package of "survival". I

Hi, You may want to try the double exponential transformation on the numerator and the denominator on this one. The method is described in detail here: http://projecteuclid.org/DPubS?service=UI&version=1.0&verb=Display&handle=euclid.prims/1145474600 If you want to give it a shot outside R there are a couple of C/C++ implementations out there:

Hello again R users, I have a devilishly hard problem, which should be very simple. I hope someone out there will have the answer to this on the tip of their tongue. Please consider the following toy example: x <- read.table(textConnection("y x1 x2 indv.1 bagels 4 6 indv.2 donuts 5 1 indv.3 donuts 1 10 indv.4 donuts 10 9 indv.5 bagels 0 2 indv.6 bagels 2 9 indv.7 bagels 8 5 indv.8

Hi all, I have a large dataset that has >10k entries. The dataset is stored in a dataframe with the headers: SubID Condition1 Condition2 Result1 Result2 There are multiple entries for a given SubID(Subject ID). Condition 1 has 3 levels and condition2 has 2 levels (therefore there are 6 possible combinations all together e.g. Cond1 Level1 x Cond2 Level 1 etc.) and i need to compute for 1. The

Hello, I was hoping that someone could answer a few questions for me (the background is given below): 1) Can the coxph accept an interaction between a covariate and a frailty term 2) If so, is it possible to a) test the model in which the covariate and the frailty appear as main terms using the penalized likelihood (for gaussian/t frailties) b)augment model 1) by stratifying on the variable that

If you mean this paper by Fahrmeir: http://biomet.oxfordjournals.org/cgi/content/abstract/81/2/317 I would recommend BayesX: http://www.stat.uni-muenchen.de/~bayesx/. BayesX interfaces with R and estimates discrete (and continuous) time survival data with penalized regression methods. If you are looking for a bona fide Bayesian survival analysis method and do not wish to spend a lot of time

Hi, I am trying to "convince" lattice to use bold face in an expression containing Greek characters without any luck. Example code (without Greek characters, renders the expression in bold) library(lattice) data(Cars93,package=="MASS") splom(~Cars93[,5:8]|Origin,data=Cars93,panel=function(x,y,...) { panel.splom(x,y,...)

Sir I have a problem that from which country i can assess R package "utils" for data entry for loading it in R Please help me in this regard. -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. Email: amnakhan493@gmail.com amna_989@hotmail.com amna_989@yahoo.com [[alternative HTML version deleted]]

Full_Name: Christos Hatzis Version: 2.0.1 patched (2005-02-18) OS: WinXP SP2 Submission from: (NULL) (24.61.19.101) I had encountered a similar problem when changing the font in the R Console and then then open the graphics device by plot(), the console window whites out. this might be related to bug report 7271. I mentioned it sometime back and was told is was not reproducible. Now I have

yes I have tries, and nothing. It just shows the strings with the slashn, just like i typed it. -----Original Message----- From: Rashid Nassar [mailto:rnassar at duke.edu] Sent: ÐÝì?ôç, 16 Óå?ôåìâñßïõ 2004 15:44 To: Christos Rodopoulos Subject: Re: [R] How do I insert a newline in my title in a plot? Have you not tried what you have already suggested: title("this is a title\nIn 2

Running Version 2.6.0 GUI 1.21. Generally, I enjoy working with the OSX GUI, which has some very pleasant features. Kudos to the developers! There is a nice history panel which slides out of the Console window. The buttons on the panel work as expected. I am puzzled by the behavior of history() and savehistory() run from the Console window. I would think that these functions would view and

Hello Colleagues, I've recently become a fan of Mark Bravington's mvbutils package for organizing analysis projects in a tree. Using cd(), Save(), fixr(), mlazy() etcetera solves nicely some of the nuisances that have worried or annoyed me and sometimes caused big problems over the years. Well thought out. Now one feature that would be fabulous would be automatic time-stamping of

Hello, Could someone recommend a good book on regular expressions with focus on applications/use as it might relate to R. I remember there was a mention of such a reference book recently, but I could not locate that message on the archive. Thanks. -Christos Christos Hatzis, Ph.D. Nuvera Biosciences, Inc. 400 West Cummings Park Suite 5350 Woburn, MA 01801 Tel: 781-938-3830