similar to: Model building ...

Un bon mot s'il vous plait. I have coded an R routine to decompose messy data following Murray Aitkin and Brian Francis's NPMLE GLIM macros for the normal distribution only but extended to incorporate censoring and variance heterogeneity. Essentially it is a wrapper for nlm() in the M-part while the E-part re-estimates the weights in the same way as the GLIM macros. The big problem is

> >>>>> "PD" == Peter Dalgaard BSA <p.dalgaard at biostat.ku.dk> writes: > > PD> "Venables, Bill (CMIS, Cleveland)" <Bill.Venables at cmis.CSIRO.AU> > PD> writes: > >> The fact that every elementary book on statistics does it this way > >> does not make it correct. To be helpful, a histogram really

Dear all, I want to compute the NPMLE for interval censored data using the EMICM function in Icens package, but I didn't find much detailed instructions on how to use this. I couldn't even figure out how to specify the status of right censored, interval censored, observed. (probably due to my level of understanding.) Does anyone have experience on using this and could you please help me?

Hi all, I have uploaded an update of the interval package, which calculates nonparametric maximum likelihood estimates (NPMLE) of distributions and exact or asymptotic weighted logrank-type tests for interval censored data. An important part of the update is that now the default initial estimate of the NPMLE uses the MLEcens package and is much faster. This increases the speed of the logrank-type

Hi all, I have uploaded an update of the interval package, which calculates nonparametric maximum likelihood estimates (NPMLE) of distributions and exact or asymptotic weighted logrank-type tests for interval censored data. An important part of the update is that now the default initial estimate of the NPMLE uses the MLEcens package and is much faster. This increases the speed of the logrank-type

Dear R community, I am using the R package `interval` in order to perform some modelling tests of the NPMLE convergence in the case of censoring. So all I am doing is drawing a sample from exponential distribution, making it a censored sample and computing the NPMLE of its distribution function. But when run on Linux Calculate 10.4 the program keeps crashing and reporting a segmentation fault

Digium is sending a couple of Asterisk representatives to ISPcon next week in Baltimore: Greg Boehnlein (author of AstWind) and Brian Capouch (professor at St. Joseph's College and author of an upcoming Asterisk book by Prentice-Hall). We will also have a 10 x 10' booth there dedicated to Asterisk. *Greg and Brian need some help manning the booth.* If you have a business centered

First off I'll admit that I'm a bit of a newbie to Samba and Linux so sorry if this is a dumb question, but I had a look and couldn't find the answer anywhere. How do I map a username on a Windows 98 box to a different username on a Linux box? For example, if my username on a Windows machine was 'Daniel Foster' and my username on the linux box was 'dan' how do I make

I have in my inbox a note saying books are expected in early August. However, for R users there is relatively little advantage in the third edition, except that almost everything works unchanged. The real advantages of the third edition are for users of S-PLUS 5.x and 4.5/2000. > From: Marc R. Feldesman > To: guido at hal.stat.unipd.it; 'Guido Masarotto '; 'Maria Wolters '

What I've got... Software: Linux: Slackware 9.1 Asterisk: out of CVS - so its new. isdn4k-utils: to test the ISDN Card Hardware: PII Pentium 400Mhz (Its a test of concept machine) with 320Kb RAM 1 x ISDN BRI Card - DIVA EICON (Installed + working) 2 x Grandstream (Barbie?) BT100 SIP Phones. What Works.. I can call from one phone to the other... get read voicemail... I can

Hello, I am interested in fitting a generalized nonlinear regression (gnlr) model with negative binomial errors. I have found Jim Lindsay's package that will do gnlr, but I have having trouble with the particular model I am interested in fitting. It is a threshhold model, where below a certain value of one of the parameters being fitted, the model changes. Here is a sample: Cones:

Is there a simple say of extracting the shape parameter from gnlr? The return given is the actual function rather than the value. I have looked at all the values returned by names(gnlr). ie gmod<-gnlr(data,....) gmod$shape returns the function definition rather than the value found. Simple case I am using it is for censored normal data (rather simpler than survival4 it seems). The value

Hi list, I'd like to fit generalized gamma and weibull distributions to a number of data sets. I've been searching around and found references to R and Jim Lindsey's GNLM package, which has the gnlr and gnlr3 procedures that can do this. Now, I'm completely new to R, and I'm working my way through the introduction... Nevertheless, I'd like to ask if someone could post

Dear users of the R mailing list, I am a ph.d. student in biology working on red deer in Norway, who would like to conduct an analysis with random factor where the response is binomially distributed. This cannot be conducted in S-plus, and I was told by others that it may be possible in R. However, I soon got into trouble which I hope you can help me to solve. My model is on this form:

It may be my lack of unerstanding, but round() seems to me to give inconsistent results when rounding 5s as in the following examples? > round(1.45, 1) [1] 1.4 # OK > round(2.45, 1) [1] 2.5 # shouldn't this be 2.4? > round(1.05, 1) [1] 1.1 # 1.0 ? and signif(): > signif(2.445, 3) [1] 2.44 # OK > signif(3.445, 3) [1]

Is there any documentation about replacement and extraction functions in R, these fuenctions don't work as in S. Example with funtion factor.names (of S+ 4.5r2) > fnombres <- list(tratamiento=LETTERS[1:4], veneno=c("I","II","III")) > veneno.diseno <- data.frame( + factor(rep(c(1,2,3,4), rep(12,4))), + factor(rep(rep(c(1,2,3), rep(4,3)),4)) + ) >

I just updated rpart to the latest version (3.1.0). There are a number of changes between this and previous versions, and some of the code I've been using with earlier versions (e.g. 3.0.2) no longer work. Here is a simple illustration of a problem I'm having with xpred.rpart. iris.test.rpart<-rpart(iris$Species~., data=iris[,1:4], parms=list(prior=c(0.5,0.25, 0.25))) + ) >

On 2009-02-18, at 14:53, DeLesley Hutchins wrote: > On 2009-02-18, at 08:06, Gordon Henriksen wrote: > >> Still, there are a large number of potential foibles here. For >> instance, passing an argument can require platform-specific >> contortions to conform to the platform ABI... > > Are those contortions done by the native code generator back-end, or > are

I've created a bunch of postscript trees (post.rpart), which I subsequently edit in Adobe Acrobat 4.05 (I need additional labelling that isn't easy to add with R). After editing the labels and annotations using Acrobat, I then export the file to Adobe's .eps format so I can insert the picture into a Microsoft Word 2000 document. The .pdf file distilled from the .ps file is

Are there any S-Plus or R libraries/packages that do Finite Mixture Analysis following the algorithms similar to those implemented in Geoffrey MacLachlan's EMMIX program? Thanks. Dr. Marc R. Feldesman email: feldesmanm at pdx.edu email: feldesman at attglobal.net fax: 503-725-3905 "Don't know where I'm going. Don't like where I've been. There may be no exit. But