Displaying 6 results from an estimated 6 matches for "ggamma".
Did you mean:
gamma
2005 Jul 12
1
three par. fitting with fitdistr
Hello,
I want to fit a tree parameter distribution to given data. I tried it with
sample data using the "fitdistr" function.
Here my workflow that didn't had any result:
I started with the generalized gamma distr, which is:
r*dgamma(x^r,shape,rate)
The R-function is:
ggamma = function (x,r,shape,rate) r*dgamma(x^r,shape,rate=rate)
For the first step I assumed r = 1 and I generated random numbers with the
"standard" Gamma distr.
rn=rgamma(1000,10,5)
In the last step I want to reconstruct the parameters from the dataset:
library(MASS)
fitdistr(rn, ggam...
2020 Oct 09
1
Aide pour finaliser ce code
...]
to=w*0
for (j in 1:R[i]){
to[j,]=w[j,]-Xbar
}
TO.jh[[i]]=to
}
#calcul des Lamda J
Lamda=matrix(0,p,p)
for (i in 1:q){
to=TO.jh[[i]]
w=matrix(0,p,p)
for (j in 1:R[i]){
w=w+(P11[[i]][j]*(to[j,]%*%t(to[j,])))
}
Lamda=Lamda+w
}
tr1=n*sum(diag(Lamda))
# Calcul de Gamma
GGamma=matrix(0,p*sum(R),p*sum(R))
PGamma=kronecker(diag(P11[[1]]),diag(p))
Ifin=p*R[1]
GGamma[1:Ifin,1:Ifin]=PGamma
for (i in 2:q){
PGamma=kronecker(diag(P11[[i]]),diag(p))
Idebut=((p*sum(R[1:(i-1)]))+1)
Ifin=(p*sum(R[1:i]))
GGamma[Idebut:Ifin,Idebut:Ifin]=PGamma
}
#Calcul de Sigma
# Calcul d...
2020 Oct 10
3
Please need help to finalize my code
...w[j,]-Xbar
93 }
94 TO.jh[[i]]=to
95 }
96 #calculation of Lamda J
97 Lamda=matrix(0,p,p)
98 for (i in 1:q){
99 to=TO.jh[[i]]
100 w=matrix(0,p,p)
101 for (j in 1:R[i]){
102 w=w+(P11[[i]][j]*(to[j,]%*%t(to[j,])))
103 }
104 Lamda=Lamda+w
105 }
106 tr1=n*sum(diag(Lamda))
107 # Gamma Calculation
108 GGamma=matrix(0,p*sum(R),p*sum(R))
109 PGamma=kronecker(diag(P12[[1]]),diag(p))
110 Ifin=p*R [1]
111 GGamma[1:Ifin,1:Ifin]=PGamma
112 for (i in 2:q){
113 PGamma=kronecker(diag(P12[[i]]),diag(p))
114 Idebut=((p*sum(R[1:(i-1)]))+1)
115 Ifin=(p*sum(R[1:i]))
116 GGamma [Idebut:Ifin,Idebut:Ifin]=PGamma
117 }...
2020 Oct 13
1
Please need help to finalize my code
...Lamda=matrix(0,p,p)
> 98 for (i in 1:q){
> 99 to=TO.jh[[i]]
> 100 w=matrix(0,p,p)
> 101 for (j in 1:R[i]){
> 102 w=w+(P11[[i]][j]*(to[j,]%*%t(to[j,])))
> 103 }
> 104 Lamda=Lamda+w
> 105 }
>
> 106 tr1=n*sum(diag(Lamda))
>
> 107 # Gamma Calculation
>
> 108 GGamma=matrix(0,p*sum(R),p*sum(R))
> 109 PGamma=kronecker(diag(P12[[1]]),diag(p))
> 110 Ifin=p*R [1]
> 111 GGamma[1:Ifin,1:Ifin]=PGamma
> 112 for (i in 2:q){
> 113 PGamma=kronecker(diag(P12[[i]]),diag(p))
> 114 Idebut=((p*sum(R[1:(i-1)]))+1)
> 115 Ifin=(p*sum(R[1:i]))
> 116 GGamma...
2020 Oct 13
0
Please need help to finalize my code
....jh[[i]]
> > 100 w=matrix(0,p,p)
> > 101 for (j in 1:R[i]){
> > 102 w=w+(P11[[i]][j]*(to[j,]%*%t(to[j,])))
> > 103 }
> > 104 Lamda=Lamda+w
> > 105 }
> >
> > 106 tr1=n*sum(diag(Lamda))
> >
> > 107 # Gamma Calculation
> >
> > 108 GGamma=matrix(0,p*sum(R),p*sum(R))
> > 109 PGamma=kronecker(diag(P12[[1]]),diag(p))
> > 110 Ifin=p*R [1]
> > 111 GGamma[1:Ifin,1:Ifin]=PGamma
> > 112 for (i in 2:q){
> > 113 PGamma=kronecker(diag(P12[[i]]),diag(p))
> > 114 Idebut=((p*sum(R[1:(i-1)]))+1)
> > 115 If...
2006 Oct 27
0
VGAM package released on CRAN
...al) Family Function
genbetaII Generalized Beta Distribution of the Second Kind
genpoisson Generalized Poisson distribution
geometric Geometric Distribution
gev Generalized Extreme Value Distribution Family
Function
ggamma Generalized Gamma distribution family function
gpd Generalized Pareto Distribution Family Function
grc Fitting Goodman's RC Association Model
gumbel Gumbel Distribution Family Function
hzeta Haight's...