similar to: extract score vector and covariance matrix in glm package

Hello; Does anybody knows that R have a function for Generelized Negative Binomial model, something like "gnbreg" in "STATA" where dispersion parameter itself is a function of covaraites ? Thanks; [[alternative HTML version deleted]]

Hello ; I am doing a regression of count data (number of award and there are some covariates) I have estiamted the parameters of negative binomial distribuion (lambda is a function of covaraites, GLM model) by glm.nb function and training dataset. Now I want to predict the number of award (for example y=0, y=1, y=2,) or testing dataset. I dont know how to calculate this numbers? I would be very

  Hello list, i have a very simple question about matrix assignment. i did like this. res<-1:30 dim(res)<-c(5,6) ind<-1:6 now i want to assign the value of this variable ind to first coloumn in matrix res. like res[,1]<-ind but this code is giving error , Actualy i have a for loop and value of ind variable is changing every time ,,total 6 times,,i want to assign every vlaue

Dear All: I used optim() to minimise the loglikelihood function for fitting data to negative binomial distribution. But there initial value of log-likelihood and iteration 10 value are reasonable. for example: initial value 1451657.994524 iter 10 value 47297.534905 iter 20 value -623478636.8236478 Then the iter 20 vlaue suddelnly changes to a negative value and in the end the error mesage is

I want to use a function as an argument to ingtegrate it twice. See the following (senseless) example of a double integration: test<-function(sf,lo,up,rest) { innerFkn<-function(sf,lo) { inte=integrate(f=sf,lower=lo,upper=4) return( inte$value ) } integral=integrate(f=innerFkn,lower=1,upper=2,sf=sf,lo=lo,up=up) return( integral$vlaue+rest ) }

Dear R users: I encountered difficulties in michaelis-menten equation. I found that when I use right model definiens, I got wrong Km vlaue, and I got right Km value when i use wrong model definiens. The value of Vd and Vmax are correct in these two models. #-----right model definiens-------- PKindex<-data.frame(time=c(0,1,2,4,6,8,10,12,16,20,24),

Currently, if you are running CentOS-4.x on a Vmware box you end up with time skew problems. The main fix has been using a set of recompiled kernel that has a HZ that is more fitting with what Vmware expects (eg 100 hz). The 4.7 kernels look like they will not need this: http://jons-thoughts.blogspot.com/2008/04/rhel-and-vmware-time-skew-problems.html Well, I've gotten quite the number of

Hello, I watched your YouTube videos for explanation for R package. They are really helpful for new beginner for using R. Please I need your help for solving inequalities that contain beta functions using R studio. I attached the file that contains this inequalities. I appreciate your help and looking forward to hear from you. Akram -------------- next part -------------- A non-text attachment

Hi, I need to write two inequalities depend on cumulative distribution (CDF) of beta binomial distribution where alpha and beta are unknown and need to find them. CDF of betabinomial(2,10,alpha,beta) <0.3<=CDF of betabinomial(3,10,alpha,beta) and CDF of betabinomial(5,10,alpha,beta) <0.8<=CDF of betabinomial(6,10,alpha,beta) How I can do that using r studio package? I tried to do

Hi all, I''ve installed paperclips but paperclip is not saving the files my model has has_attached_file :attachment,:styles => { :medium => "300x300>", :thumb => "100x100>" } and i have db migration as class AddAttachmentToPolicies < ActiveRecord::Migration def self.up add_column :policies,

Hi All, I have a doubt regarding join tables I''m having 2 models 1)Fac 2)Cont and both models have " has and belong to many" relationships so there are 3 tables 1)facs 2)conts 3)conts_facs then i''m fetching the data in controller as @conts=Cont.find(:all]) @cfacs=Fac.all(:joins=>:conts, :select=>"facs.name") but i dont know how to get the

Hello everyone, I am doing the hypothesis test and I need the score vector of each iteration in glm (family=binomial); how could I do it. I tried trace option but it gives me just AIC of each iteration nothing more. thanks; [[alternative HTML version deleted]]

Is there a way to compute score residuals from glm() fits (my current example is Poisson regression)? I'm referring to the types of residuals used in computing the sandwich covariance matrix estimator. The sandwich package will compute the Huber-White sandwich estimator but I need to use the cluster sandwich estimator to account for intra-cluster correlation. Thanks Frank ----- Frank

I have just committed some code to the r-devel branch to implement the Rao efficient score test. This is asymptotically equivalent to the LRT, but there is some indication that it might have better properties in smaller samples since it is based more directly on the distribution of the sufficient sums under the null hypothesis (e.g., if you have a divergent fit to the model under the alternative,

I want to predict a set of proportions from a linear regression using glm. The model includes a logit link. However I want to extract the variance:covariance matrix among the predicted proportions rather than on the logit scale. Is there a way to do this using Vcov or a similar package in R? Thanks Chris

Hello, Is there a way to obtain the variance-covariance matrix of the estimated parameters from GLM? my.glm<-glm(mat ~X,family = binomial, data =myDATA) out1<-predict(my.glm,se.fit = TRUE) std<-out1$se.fit se.fit is for getting the standard errors of the estimated parameters (\betas). Is there a way to get the variance-covariance matrix of the estimated parameters? Many thanks,

can anybody guide me on how to create inner joins and left and right outer joins in RoR as both the tables has "has and belongs to many relationships" and there is ajoin table as well

Hi, probably there's some parameter I don't know how to set. Using R 0.64 on a Sun sparc, I use the following command axis(1, at=xspos(times), labels=result, cex.axis=0.7) where result is sth. like c("May", "Jun" ,"Jul" ,"Aug", "Sep", ...) Now it plots it fine.. but when two xspos are too (?) close together the label gets dropped.

Hi, I am writing a function that requires a specific package to be installed. Is there a way of checking if the package is installed and returning a TRUE / FALSE result so my function can return an appropriate error message and exit the function gracefully rather than just bombing out? I'm thinking along the following lines (but want code that works), f_checkpackage <- function() { if

Am 04.10.2016 um 15:43 schrieb Rowland Penny via samba: > On Tue, 4 Oct 2016 15:16:17 +0200 > Achim Gottinger via samba <samba at lists.samba.org> wrote: > >> >> Am 04.10.2016 um 10:21 schrieb Rowland Penny: >>> On Tue, 4 Oct 2016 02:35:21 +0200 >>> Achim Gottinger via samba <samba at lists.samba.org> wrote: >>> >>>> Am