similar to: Another packaging question

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

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

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.

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:

Dear expRts, I just realized that the "..." argument in a function cannot be used without taking precautions sometimes. The following behaviour is what I stumbled upon: > myrepl <- function(length, fun, ...) { + replicate(length, fun(...))} > myrepl(20, sample, 1:5) [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 > myotherrepl <- function(length, fun, ...) { + args <-

Dear expRts, when I try to build a package by myself, the process fails writing the following messages to 00install.out: * Installing *source* package 'dynalc' ... ** libs WARNING: no source files found chmod: Zugriff auf ?/home/hannes/texte/forschung/clot/programme/dynalc.Rcheck/dynalc/libs/*? nicht m?glich: No such file or directory ** R ** data ** inst ** help Note: removing empty

This is very possibly not a question on R. I was under the impression that the argument that gives rise to Fisher's LSD method in ANOVA works in other situations with three-way comparisons too, given that formal logic works the same ("if the omnibus test rejects, only two of the three groups may be equal, and therefore only one hypothesis can be rejected falsely"). However, when I

I have just installed version 2.6.0 of R using the Ubuntu Dapper Drake packages from CRAN. I picked the installation option of overwriting the existing config files. Since then, R starts with an error message which I don't quite get because the function it claims not to have found can be called nicely. It wouldn't bug me at all, but R CMD Sweave is a bit more picky than me and dies after

Dear all, I am trying to use the clipboard when writing a table. Typing: write.table(object, file="clipboard") leads to the message: Fehler in file(file, ifelse(append, "a", "w")) : 'mode' f?r die Zwischenablage muss unter Unix 'r' sein My interpretation is that I am not allowed to write into the clipboard from a program called from R. Is there a

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

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