similar to: AIC for tweedie glm

Dear Catarina, I prefer to leave the AIC value as NA for the tweedie GLM family because it takes extra time to compute and is only occasionally wanted. It's easy to compute the AIC yourself using the dtweedie() function of the tweedie package. Best wishes Gordon At 03:05 AM 14/08/2007, Catarina Miranda wrote: >Dear Gordon; > >I have also sent this email to R help mailing list,

Say we have a tab delimited file called bug.txt Part Rep Cage Hb pupae 1 1 S 32 1 M 34 L 42 2 S 36 M 28 L 36 read.delim("bug.txt") Part Rep Cage Hb.pupae 1 1 1 S 32 2 NA 1 M 34 3 NA NA L 42 4 NA NA NA 5 NA NA NA 6 NA 2 S 36 7 NA NA M 28 8 NA NA L 36 >

I am new to R and have two scripts written slightly different but should to relatively the same thing but my lack of experience with the program I can not figure out the what I need to do to correct it. The first script gives me a consistent mean fold change values with every run but can generate negative p values for some. For the second version of the script, the fold changes seem to be very

On Sat, Oct 15, 2022 at 12:21:41PM -0700, syzbot wrote: > Hello, > > syzbot found the following issue on: > > HEAD commit: a1b6b102df03 Merge branch 'phylink_set_mac_pm' > git tree: net > console output: https://syzkaller.appspot.com/x/log.txt?x=179af0c2880000 > kernel config: https://syzkaller.appspot.com/x/.config?x=85495c44a2c25446 > dashboard

Hello, I'm wondering if someone may be able to help me, and do apologize if there is a simple and obvious solution for this. I am somewhat new to R, and have been searching for a simple solution for a couple of days. I am interested in finding a tool that allows me to plot a stacked bar plot. My data set is in the following format:

Hi all I've been getting some emails from the R package maintainers that I need to update some code in a CRAN packge that uses FORTRAN, to comply with (not so recent) changes. I've been a little busy... I'm having trouble adjusting my code. I hope someone can help. The package was working fine, and a few R functions in my package had lines like this: tmp <- .Fortran(

Dear R users; I am modelling densities of some species of birds, so I have a problem with a great ammount of zeros. I have decided to try GLMs with the tweedie family, but in all the models I have tried I got an NA for the AIC value. Just to check the problem I've compared the a glm using the Gaussian family with the identity link and a glm using the tweedie family with var.power=0 and

Hello all, I have count data with about 36% of observations being zeros. I found in some of the examples of the r-help mail archives that a tweedie family of distributions could be used to fit a model with random effects. Upon installing the tweedie package and attempting to fit the following model: lmer(SUS ~ 1 + (1|

Post-hoc test for a zero inflated continuous data set with a tweedie distribution? I have a zero inflated continuous data set of aphid feeding duration on 10+ species of plant. I have fitted a glm model with a tweedie distribution and used anova() function to show that there is significance between the plant species. However, I would now like to perform of post-hoc test, ideally a Tukey-Kramer

Hi there, I've been wanting to fit a GLMM and I'm not completely sure I'm doing things right. As I said in a previous message my response variable is continuous with many zeros, so I was having a hard time finding an appropriate error distribution. I read some previous help mails given to other people advising them to use the Tweedie distribution. I'm still not sure if this

Good afternoon My name is Joanne Lenehan, I am a post grad at UNE using R version 2.9.0 I came across the Tweedie package in old R help posts and was interested in giving it a go for some data. The data is below and also attached as BaregroundLitterLogs Site Treatment Graze Dam Plot Time Bare Litter Logs 1 C remote yes 1 A 0 2 0 1 C remote yes 2 A 0 15 0 1 HE

Hi there, I'm using the tweedie distribution package and I cant figure out how to run a model using an identity link. I know I can use a log link by having link.power=0 and I think identity would be link.power=1, but I'm not sure. Furthermore when I try running it with link.power=1 it requires starting values which I cant manage to give appropriately so I'm not sure if its

Dear list, When running the script below, the optim function returns the message 'Error in function (par) : could not find function "fn"'. I've re-read the code a few times and just can't figure out why it's happening. Any help/suggestion is appreciated. Regards, Rubem ## R script library(tweedie) period<-factor(1:4) herd<-factor(1:10)

Hi I'm encountering an error/warning when doing multiple comparisons with the package multcomp on a coxme model. My data: I'm looking at the removal of brood from the nest according to three treatments I applied on the brood. The brood and the workers caring about the brood in the nest, belonged to different colonies. Factor: treatment (3 levels: tx,uv,meta) Random effect 1: origin of

Hi, i tried to insmod the imq.o module from http://pupa.da.ru/imq after a successful compile, but it thows this error: > insmod imq.o imq.o: unresolved symbol nf_unregister_hook imq.o: unresolved symbol nf_register_hook I think I have some kernel options disabled, does anyone know which one(s)? Thanks cord _______________________________________________ LARTC mailing list /

Hi, I'm not sure if this is apprriate or not, but here goes anyway. The patch below updates MAINTAIER address Individuals (Only Andrew :): osdl.org -> linux-foundation.org Lists: osdl.org -> lists.osdl.org I assume the latter will change at some stage, but at least with this change the osdl/linux-foundation lists are consistent. Signed-off-by: Simon Horman

Hello there, I was wondering if somebody could offer me some advice on which error distribution would be appropriate for the type of data I have. I'm studying what continuous predictor variables such as grooming received, rank, etc. affect the amount of grooming given. This response variable is continuous with many zeros, and so positively skewed. So I can't use gamma because of the

Hi R-users: New to R and I am trying to run a GLM with random effects. I have 3 replicates ('Replicate) of counts of parasites ('nor.tot.lep') before and after an experiment ('In.Out'). When I run lmer I get the error messages (16 of each) below... > lmer(nor.tot.lep ~ In.Out + (In.Out|Replicate),data=coho, family =tweedie(var.power = 1, + link.power = 1)) Generalized

SunOS fluke 5.8 Generic_108528-07 sun4u sparc SUNW,Ultra-Enterprise R-1.4.0 Running from the command line R CMD check chlib results in the log file 00check.log: * using log directory `/home/woodstock/hoffmann/R/Sources/chlib.Rcheck' * checking for file `chlib/DESCRIPTION' ... OK * checking package directory ... OK * checking for sufficient/correct file permissions ... WARNING *

Stephen, all, Thank u very much for your answer and i wish you an happy thanksgiving. Currently we tried to migrate our base on RAW devices already using e41smp, SecurePath (HP) 3C and Qlogic 7.00. The average response time of our sql request on ia32 is 0.32Sec very regularly. On our Itanium on OCFS that's very randomly between 1 and 15Sec, and on RAW between <1s and 3 sec. David is an