Displaying 3 results from an estimated 3 matches for "modl2".
Did you mean:
model2
2009 Jun 05
2
p-values from VGAM function vglm
Anyone know how to get p-values for the t-values from the coefficients
produced in vglm?
Attached is the code and output ? see comment added to output to show
where I need p-values
+ print(paste("********** Using VGAM function gamma2 **********"))
+ modl2<-
vglm(MidPoint~Count,gamma2,data=modl.subset,trace=TRUE,crit="c")
+ print(coef(modl2,matrix=TRUE))
+ print(summary(modl2))
[1] "********** Using VGAM function gamma2 **********"
VGLM linear loop 1 : coefficients =
0.408464609241, 3.255887520104, -0.00022...
2009 Jun 26
1
predicted values after fitting gamma2 function
...nts))
colnames(Max.cnt)<-c("Simulation","MidPoint","Count")
then I fit this to a gamma distribution function:
modl<-
vglm
(Count
~
MidPoint
,gamma2
,data
=subset(Max.cnt,select=(simulation,MidPoint,Count),trace=TRUE,crit="c")
print(coef(modl2,matrix=TRUE))
print(summary(modl2))
This produces the output:
VGLM linear loop 1 : coefficients =
3.231473495, 0.597085743, -0.195591168
...
VGLM linear loop 20 : coefficients =
3.663316903, 0.897355891, -0.620449146
log(mu) log(shape)
(Intercept) 3.6633169 0...
2009 Jun 16
0
Generation from COX PH with gamma frailty
...e know how to get p-values for the t-values from the coefficients
> produced in vglm?
> Attached is the code and output ? see comment added to output to show
> where I need p-values
>
>
> + print(paste("********** Using VGAM function gamma2 **********"))
> + modl2<-
> vglm(MidPoint~Count,gamma2,data=modl.subset,trace=TRUE,crit="c")
> + print(coef(modl2,matrix=TRUE))
> + print(summary(modl2))
>
>
> [1] "********** Using VGAM function gamma2 **********"
> VGLM linear loop 1 : coefficients =
> 0.40...