search for: clot

...;optim" function to replicate the results from the "glm" using an example from the help page of "glm", but I could not get the "optim" function to work. Would you please point out where I did wrong? Thanks a lot. The following is the code: # Step 1: fit the glm clotting <- data.frame( u = c(5,10,15,20,30,40,60,80,100), lot1 = c(118,58,42,35,27,25,21,19,18), lot2 = c(69,35,26,21,18,16,13,12,12)) fit1 <- glm(lot1 ~ log(u), data=clotting, family=Gamma) # Step 2: use optim # define loglikelihood function to be maximized over # theta is a vector...

...he probabilities from a glm.predict() result: > str(res) Named num [1:9] 0.00814 0.01877 0.025 0.02941 0.03563 ... - attr(*, "names")= chr [1:9] "1" "2" "3" "4" ... I got from: # A Gamma example, from McCullagh & Nelder (1989, pp. 300-2) clotting <- data.frame( u = c(5,10,15,20,30,40,60,80,100), lot1 = c(118,58,42,35,27,25,21,19,18), lot2 = c(69,35,26,21,18,16,13,12,12)) model <- glm(lot1 ~ log(u), data=clotting, family=Gamma) res <- predict(model, clotting) I want to transfer the probabilities "0.00814 0.0187...

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 section \value Note: removing empty section \details Note: removing empty section \description Note: removing empty section \note Note: removing empty section \reference...

...xpecting. The y-axis of the plots is tau and I would like that axis to be the likelihood so that I have a curve that maxes at the parameter estimate. I am not sure where I find those likelihood values? I may just be misinterpreting the theory behind this. Thanks for any help you can give. Max clotting <- data.frame( + u = c(5,10,15,20,30,40,60,80,100), + lot1 = c(118,58,42,35,27,25,21,19,18), + lot2 = c(69,35,26,21,18,16,13,12,12)) glm(lot2 ~ log(u), data=clotting, family=Gamma) prof<-profile(glm2) plot(prof) [[alternative HTML version deleted]]

Hi R users, I want to use dglm Package. I run the examples and it give me an error: Error en dglm(lot1 ~ log(u), ~1, data = clotting, family = Gamma) : no se pudo encontrar la funci?n "as.family" dglm can't find "as.family" function why ? Thank you for your help

...-value p-value (Intercept) 52.79943 4.467854 21 11.817625 0.0000 time -0.44162 0.231605 21 -1.906793 0.0703 Correlation: (Intr) time -0.242 Number of Observations: 37 Number of Groups: 15 > > # Plot data > xyplot(value~time,group=subject,type="b",main="Clot size vs. time.") > # Add summary regression line > abline(fitxx) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : (list) object cannot be coerced to type 'double' > > John David Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics University of M...

Hi, I have a data set (weight) that does not follow the Gaussian (Normal) distribution. However, I have to transform the data before applying the Gaussian distribution. I used this syntax and used log(weight) as: posJy.model<-glm(log(weight) ~ factor(pos), family=gaussian(link='identity'), subset=Soil=="Jy"). This syntax COULD NOT transform the data. But if I transform the

...1 05-05-03 15:32 faraway/data/breaking.rda 2729 05-05-03 15:32 faraway/data/cathedral.rda 5909 05-05-03 15:32 faraway/data/chicago.rda 2181 05-05-03 15:32 faraway/data/chiczip.rda 5449 05-05-03 15:32 faraway/data/chmiss.rda 2133 05-05-03 15:32 faraway/data/clot.rda 2341 05-05-03 15:32 faraway/data/coagulation.rda 1389 05-05-03 15:32 faraway/data/corrosion.rda 24960 05-05-03 15:32 faraway/data/ctsib.rda 2061 05-05-03 15:32 faraway/data/death.rda 2065 05-05-03 15:32 faraway/data/drugpsy.rda 2989 05-05-03 15:32...

Anyone conceive of any reason why I should have a stream that RealPlayer and iTunes have no problem with, but MP9 refuses to play: just sits there dumbly having loaded the .m3u that points to it, ignoring the play button?

I am experiencing a curious event, and wondering if others have seen this... As well, I have a question related to it. Today, I noticed my puppet summary report from Foreman this morning, that 60 of my 160 hosts all stopped reporting at nearly the exact same time, and have not since restarted. Investigating, it appears that my puppetmaster temporarily ran out of disk space on the /var

...uot; function To: r-help at stat.math.ethz.ch Message-ID: <46C84066.6040108 at sinu.unicordoba.edu.co> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Hi R users, I want to use dglm Package. I run the examples and it give me an error: Error en dglm(lot1 ~ log(u), ~1, data = clotting, family = Gamma) : no se pudo encontrar la funci?n "as.family" dglm can't find "as.family" function why ? Thank you for your help ------------------------------ Message: 4 Date: Sun, 19 Aug 2007 16:00:36 +0100 (BST) From: Prof Brian Ripley <ripley a...