similar to: Survreg Weibull lambda and p

Hi I have been asked to provide Weibull parameters from a paper using Kaplan Meir survival analysis. This is something I am not familiar with. The survival analysis in R works nicely and is the same as commercial software (only the graphs are superior in R). The Weibull does not and produces an error (see below). Any ideas why this error should occur? My approach may be spurious.

Hi All, I am currently conducting a mixed model. I have 7 repeated measures on a simulated clinical trial. If I understand the model correctly, the outcome is the measure (as a factor) the predictors are clinical group and trial (1-7). The fixed factors are the measure and group. The random factors are the intercept and id and group. I tried using 2 functions to calculate mixed effects.

Hi All., Following on from all the help I received earlier .. Just so as you know a well known text book (T&F) uses "repeated observations" - I think perhaps they meant repeated measures in retrospect I am a total beginer to ARIMA, but have read a a basic text. I read that ARIMA requires a minimum of 50 repeated measures in time (repeated measures)- Unfortunately, I have 8

Matching regular expressions Dear useRs! I have the following problem. I would like to find objects in my environment that have two strings in it. For example, I might want to find objects that have in their names "MY" and "TARGET". I do not care about the ordering of these two substrings in the name, neither what is in front, behind or between them, the only thing important

Dear useRs! I have written a function that should pass argument "m" to the next function, however it does not! Please have a look at the function below that shows a problem and tell me what I am missing. As you can see, the "blocks" argument is passed corectly, while "m" is not. Best, Ales Ziberna opt.par.new<-function( #function for optimizig

Hello! I am currently in the process of creating (my first) package, which (when ready) I intend to publish to CRAN. In the process of creating this package I have taken some code form existing packages. I have actually copied parts of functions in to new functions. This code is usually something very basic such as Rand index. What is the proper procedure for this? Since most of R (and also the

Hello! I would like to put an object in to a function that calls the current function. I thought the answer will be clear to me after reading the help files: ?assign ?sys.parent However it is not. Here is an example I thought should work, however it dose not exactly: f<-function(){s();print(a)} s<-function()assign(x="a",value="ok",pos=sys.parent()) f() #I want to get

Dear useRs? I have the following problem! I have a function that calls one or more functions, depending on the input parameters. I am searching for the fastest way to select and execute the selected functions and return their results in a list. The number of possible functions is 10, however usually only 2 are selected (although sometimes more, even all). For examples, if I have function

Dear useRs! I would like to generate a list of all possible (unique) partitions of N units into k classes. For example, all possible partitions of 4 units into 2 classes are (I hope I have not missed anyone): 1,1,1,2 (this can be read as {1,2,3},{4}) 1,1,2,1 1,2,1,1 2,1,1,1 1,1,2,2 1,2,1,2 1,2,2,1 The partitions 1,1,2,2 and 2,2,1,1 are the same and are therefore not two unique

Several times my packages have fallen foul of the latex check at the uploading-to-CRAN stage, creating extra work for Kurt Hornik. I've never bothered trying to get latex working on my own packages (happy with vanilla help) but am told that "latex is the closest thing to a syntax checker for Rd files" and that I should set it up. Now, there is already a latex on my (Windows XP)

Does anybody know how to produce a East-European character "??" - "c with a v-shaped hat " in R (in text or plot)? I do know how to get "??,??" - "s,z, with a v-shaped hat", but not this one. Thanks in advance for any suggestions, Ales Ziberna P.S.: I am using Windows XP and R version 1.9.1 (details below) $platform [1]

Hello, it seems that the main results of survival analysis with package survival are shown only as side effects of the print method. If I compute e.g. a Kaplan-Meier estimate by > km.survdur<-survfit(s.survdur) then I can simply print the results by > km.survdur Call: survfit(formula = s.survdur) n events median 0.95LCL 0.95UCL 100.0 58.0 46.8 41.0 79.3 Is

Hello R-devel, Best wishes in the new year. I am writing to kindly request new R function so NA_real_ can be more easily detected. Currently if one wants to test for NA_real_ (but not NaN) then extra work has to be done: `is.na(x) & !is.nan(x)` Required functionality is already at C level so to address my request there is not that much to do. Kevin Ushey made a nice summary of current R C api

Hi, I am a newbie with R, so I hope my question isn't too stupid. I am trying to generate dose-response curves using the "drfit" package. I have formatted my CSV files to the correct format, and have no trouble running drfit to get a summary of my data. The problem is that when I try to use "drplot" to graph my data I get an error. The message is: Error in

Hello! I would like to convert rows of a matrix to a elements of a list. #For example, if I have mat<-matrix(1:100,ncol=5, nrow=20) #I can do: list<-apply(mat,1,list) list #however this is not quite what I want. To get what I want, I have to do: list<-lapply(list,function(x)x[[1]]) list Is there a faster way? Best regards, Ales Ziberna

Hello! I have created a plot and would like to put a legend in the top left corner (above and left of plotting region), but non of the functions can plot there (I tried 'legend', 'text', 'points', 'lines', 'mtext', 'axes',...). I would be very grateful if someone could help me put something in that corrner, with any fuction. Thanks, Ales Ziberna

Hello! I have a problem checking the package. Firstly, I do not know how to specify the package to check. I tied specify it by supplying the path and by runing the R CMD check in the directory of the package. In addition to that, I get an error bellow. Any suggestions on how to set TMPDIR would be greatly appriciated! C:\Ales\Statistika\Blocno modeliranje\dr\blockmodeling>R CMD check *

"nana" is meant to express "NA, really NA". Your suggestion sounds good. On Thu 2 Jan, 2020, 3:38 AM Pages, Herve, <hpages at fredhutch.org> wrote: > Happy New Year everybody! > > The name (is.nana) doesn't make much sense to me. Can you explain it? > > One alternative would be to add an extra argument (e.g. 'strict') to > is.na(). FALSE by

Hello! I have writen some subrutines in Free-form Fortran. I would like to includ them in a package, which I would like to build on WinXP. I have all suggested tools/programs for bulding R packages on Windows (except latex). What is the best way of using these subrutines? Does sombody mybe know any translation tools for converting Free-form to fixed-form Fortran? Thanks for any suggestions,

Dear CSAG R users, I will be glad if someone can point out what I am doing wrong or not doing at all in this. I am trying to write out netcdf file in R. I have 26 time step but only the first time step is written. For example: >library(ncdf) >path <- '/home/work/' >forecast <- open.ncdf(paste(path,'cam.1980.2005.nc',sep="")) > fore <-