similar to: XML and odfWeave

Dear list, does anyone know of a R-package that has implemented the increasingly popular inclusion of the number of patients at risk below Kaplan-Meier curves like in http://bloodjournal.hematologylibrary.org/content/vol116/issue19/images/large/zh89991058760001.jpeg any hint (or negative answer) is much appreciated. Thanks Thorsten -- Thorsten Raff 2nd Medical Department, University

Dear R-helpers I am wondering if there is an option to the survplot function in the design package that allows for drawing Kaplan-Meier plots starting from 0 instead of 1, similar like fun = 'event' in the standard plotting function used on a survfit object. I apologize in advance for having missed any obvious informational sources but I really didn't find anything in the

Hello I have the problem that I want to transform a dataframe as generated by diagnosis <- rep(diagnosis[1:3], 3) marker <- gl(3,3) values <- rnorm(9) dataframe <- cbind(diagnosis, marker, values) dataframe <- dataframe[c(1:5, 7:9), ] into a matrix where levels(diagnosis) is indicating the rows of the matrix, levels(marker) the columns and values are the actual content of the

Dear R-helpers, I am trying to plot a Bland-Altman-Plot using the BA.plot function from the package MethComp. While there is a function to transform the values for analysis as shown in the snippet below, I would like to have logarithmic axes for display as well. The usual log = 'xy' does not work because of the properties of the y-asxis (positive and negative values). I am sure that

Hello, I am looking for a package implementing unconditional tests for comparing two binomials (Barnard's test or derivatives). As I am looking for rare genotypes within my groups, my number of subjects in one of the cells of the 2x2 table is often 2-4 resulting in a too conservative result by Fisher's exact test due to the discrete sample space of the test statistic. As I want to use

Hello R-experts, I have a simple question: how do I aggregate data using the median function with na.rm = TRUE, i.e. aggregate(x, list(x$y, x$z), FUN = '????') I have tried 'median(x, na.rm = TRUE)' but that doesn't help. Any suggestions are very much appreciated, Thorsten

Hi, I have a vector, which contains both strings and numbers, e.g. > foo <- c("str1",1234,"str2",0.9876) I want to know if a distinct element of the vector is a string or a number and took "is.numeric", but > is.numeric(foo[2]) [1] FALSE because R treats the numbers in a mixed vectors as strings: > foo [1] "str1" "1234"

Dear R Users: The functions of our "micEcon" package [1,2] can be subdivided into three categories: - microeconomic demand and firm models - sample selection models (mainly selection()) - routines for (likelihood) maximisation (e.g. maxLik(), maxNR(), maxBHHH()) (mainly used for ML estimation of sample selection models) Although sample selection models are often used in

Dear All, I would like to publish a function for 'heckit' estimations together with two examples from Greene's and Wooldridge's econometric textbooks. These examples use the dataset of Mroz (1987) that is also available in John Fox' "car" package. However, not all variables that are used in my examples are available in the "car" package. Therefore, I

Hi, I was really happy when I saw that in R version 2.3.0 "R CMD check" works for packages whose package name is different from the directory name in which it is located (see http://cran.r-project.org/src/base/NEWS). Now, I can have branches of my packages in directories like "[...]/branches/<version number>/", while I had to use "[...]/branches/<version

Hi, I have a dataframe (say myData) and want to get a list (say myList) that contains a matrix for each row of the dataframe myData. These matrices are calculated based on the corresponding row of myData. Using a for()-loop to do this is very slow. Thus, I tried to use apply(). However, afaik apply() does only return a list if the matrices have different dimensions, while my matrices have

Hi, I added a new function "predict.systemfit" to our package "systemfit" to make it closer to other packages (e.g. lm). Now "R CMD check" complains that the generic function "predict" has only the argument "object", while our function "predict.systemfit" has more arguments. However, the function "predict.lm" has also more

Dear R helpers, For an allometric analysis (allometric equation y = a*x^b) I would like to apply a non-linear regression instead of using log-log transformations of the measured parameters x and y and a Model II linear regression. Since both of the variables x and y are random, I would like to apply a Model II non-linear analog of either Reduced Major Axis or Major Axis Regression. The

?Hi: I always use subset the same way but now is returning 0 rows. What's wrong with the way I am subsetting? library(ggplot2) structure(list(first = c(38.2086, 43.1768, 43.146, 41.8044, 42.4232, 46.3646, 38.0813, 40.0745, 40.4889, 38.6246, 40.2826, 41.6056, 34.5353, 40.0768), second = c(43.3295, 42.4326, 38.8994, 37.0894, 42.3218, 46.1726, 39.1206, 41.2072, 42.4874, 40.2657, 38.7766,

Hi, I wrote a package for linear programming and want to submit it to CRAN. Since the package 'quadprog' has a function with the name 'solve.QP' to perform Quadratic Programming, I named my (main) function 'solve.LP'. However 'R CMD check' gives one warning: * checking generic/method consistency ... WARNING solve: function(a, b, ...) solve.LP:

Dr. Michael Breuer 22.04.05 ?kologiezentrum der Universit?t Kiel Olshausenstra?e 75 24118 Kiel Dear Ladies and Sirs, After updating the R-Windows-program (binary) by the latest version (2010), the R-Scripts that I want to execute are not shown in the File-Window anymore. In the former version it worked

Hello world, I am actually transferring a course in data management for students in biology, geography and agriculture from statistica to R - it works surprisingly well. If anyone is interested in my scratch/notepad (in German language), please see www.hydrology.uni-kiel.de/~schorsch/statistik/statistik_datenauswertung.pdf (pages 40-52) The dataset is:

Hi, I did not find any web page about using R in economics and econometrics so far. However, this does not mean that there is none (searching with google for "R" and "economics" gives many pages about economics and a name like Firstname R. Lastname on it ;-)). Does anybody in the list does know such a web page? If not, I will be happy if you, Ajay, could build and

The function 'pchisq' from the 'stats' library gives a wrong result if the argument equals exactly zero: # Upper tail of central 1-df chi^2 distribution > pchisq(1 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1] 0.3173105 > pchisq(0.5 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1] 0.4795001 > pchisq(0.01 , 1, ncp=0, lower.tail = F, log.p = FALSE) [1]

Dear All, we (Ott Toomet and I) would like to add functions for maximum likelihood (ML) estimations of generalized tobit models of type 2 and type 5 (*see below) in my R package for microeconomic analysis "micEcon". So far we have called these functions "tobit2( )" and "tobit5( )". Are these classifications well known? How are these functions called in other