similar to: setting attributes

I am trying to model survival data with a Weibull distribution using survreg. Units are clustered two apiece, sometimes receiving the same treatment and sometimes opposing treatment.

What is the least surprising way of initializing a factor with predefined levels and with length 0? as.factor(c("eins", "zwei", "drei"))[FALSE] does the job but looks a bit weird. -- Johannes H?sing There is something fascinating about science. One gets such wholesale returns of conjecture mailto:johannes at

I am failing to uncover Sweave chunks step by step using the LaTeX beamer class. The following minimal example: \documentclass{beamer} \usepackage{Sweave} \begin{document} \begin{frame}[fragile] In the year \uncover<2->{25}\uncover<3->{\Sexpr{5*5}} \uncover<4->{ <<echo=TRUE, print=TRUE>>= 5*5*101 @ } \end{frame} \end{document} leads to an error message when

Hi, I was wondering if anybody might have a reference for me: My R code is growing and getting more and more confusing. Thus, I figure it's time to switch to object-oriented again. I have done oo programming in C++ and Java before but the first few tutorial on R oo were a bit confusing for me. Is there any brief tutorial on oo programming in R especially for people who have done oo in Java

I am about to give an introduction to R to some clinical data managers used to SAS. There is already a lot of material in printed form and on the web that paves the way. What I haven't found so far are text wrapping capabilities in setting tables in raw text as in SAS PROC REPORT. At the moment i would direct them at producing HTML output from R and pipe the result through lynx. Coming from

Whenever I try to load the Matrix package, I get the following error message: libRlapack.so: cannot open shared object file: No such file or directory A file with that name is indeed not on the hard disk. I am using the R version which comes with Ubuntu Hardy Heron LTS. Here is the output of R.Version(): R.Version() $platform [1] "i486-pc-linux-gnu" $arch [1] "i486" $os [1]

I'm guessing there's a more efficient way to do the following using the index features of R. Appreciate any thoughts.... for (i in 1:nrow(dbs1)){ if(dbs1$Payor[i] %in% Payor.Group.Medicaid) dbs1$Payor.Group[i] = "Medicaid" if(dbs1$Payor[i] %in% Payor.Group.Medicare) dbs1$Payor.Group[i] = "Medicare" if(dbs1$Payor[i] %in% Payor.Group.Commercial)

The other day I was looking into one of the classics in resampling, Eugene Edgington's "Randomization Tests". This type of test is simple to do in R with things like a simple correlation, the sample () function is perfect for the purpose. However, things are more complex if you have grouped data, like a one-way ANOVA. The reason is that you have to avoid the consideration of

Hi all, I can't come over a problem with formula. Suppose I have a coxmod model with the following formula: > somemod$formula Surv(lebzeit, tot == 1) ~ sex + (alter >= 65) + diff3k + zelltyp_k_c + q_nuc_3k + kar_k80_g80 + stadium and I want to drop the stadium explanatory variable from the model with > update(somemod, ". ~ . - stadium") I get the following messages:

I'm interested in using R's plotting capabilities to try to generate a graph showing the relationship/pairing frequency of words appearing in a block of unstructured text. I don't have a specific algorithm or approach in mind, just looking to portray text in an interesting fashion. The output I'm hoping for is something like the following called a "Phrase Network":

Hi all, Can anyone here please tell me whether is it possible to produce a chart displayed in http://www.datawolf.blogspot.com/ in R for visualizing multivariate time series? If possible how? Regards, --------------------------------- [[alternative HTML version deleted]]

I have strings of the form TICKER.GGG.XXXXXX.dat but GGG is not always three characters so I can't use substr to pull it out of the string. Could someone tell me how to use sub to pull out the GGG but more generally the string between the dot after the R in TICKER and the next dot. I still don't have a very good understanding of regular expressions but I'm trying to get that

Hey, I am having trouble graphing the following function √2Γ(n/2)/[√n − 1Γ((n − 1)/2 for the values of n between 2 and 50. i know that Γ(n) = (n-1)!, which in R is factorial(n-1) When i type that into R, using y <- function(n)..... and then plot(y,2,50), it doesnt give me anything meaningful, in fact, it comes up with a message saying something like "in gamma(n+1) ploted" or

I am looking for a mailing list for general statistical questions that are not R related. Do you have any suggestions for lists that are busy and helpful and/or lists that you use and recommend? Thanks in advance, Ralf

Hello, ? What I really want to do is to add a rejection region in the form of a long rectangle to a density plot I have drawn.? I am getting? 2 plots.? How can I add rectangle to first plot?? see code below. First section works fine.? It just is not quite what I want. # NORMAL DISTRIBUTION PLOT OF RAW DATA WITH UPPER CRITICAL LEVEL - ok xcrit=144.1 # *** single-sample Upper one-tailed

It seems that Sweave is supposed to be used from Latex and R is called during the LaTeX compilation process whenever R chunks appear. What about the other way round? I would like to run it triggered by R. Is this possible? I understand that this does not correspond to the idea of literate programming since it means that there is R code running outside the document, but for my practical approach, I

Dear List, I have a series of nested loops with the structure shown below, and I am struggling to figure out how to divert output to folders created with dir.create() within the loops. What I need is for the output to end up as topNameK/subNameL/objNameM.pdf; what I get instead is a series of directories topNameK/, directories subNameL/, and files objNameM.pdf, all in the working directory. Any

Error in lda.default(x, grouping, ...) : variables 1 3 5 8 10 15 17 20 27 29 34 appear to be constant within groups what this sentence means?? [[alternative HTML version deleted]]

Hello Can someone help me with the following, specifically in judging whether a matrix exists. I have trouble with the first line below. In this case, matrix obj$hessian exists and is 74 x 74. I receive the error message: Warning message: In all(w$hessian) : coercing argument of type 'double' to logical Thank you all. --- if (!all(obj$hessian)|OPG){ vb<-obj$gradientObs;

Dear sir, I want to know about R command for parameter estimation for generalized Poisson regression. yours faithfully Arafat Ud Zaman 4th year(Hons) Applied statistics University of Dhaka Mob- +8801715173176 E-mail: arafat.zaman@yahoo.com azaman1@isrt.ac.bd azaman1@quantummethod.org.bd [[alternative HTML version deleted]]