search for: betav

Hello, I am trying to construct two matrices, F and V, composed of partial derivatives and then find the eigenvalues of F*Inverse(V). I have the following equations in ryacas notation: > library(Ryacas) > FIh <- Expr("betah*Sh*Iv") > FIv <- Expr("betav*Sv*Ih") > VIh <- Expr("(muh + gamma)*Ih") > VIv <- Expr("muv*Iv") I successfully found the partial derivatives: > f11 <- deriv(FIh, "Ih") > f12 <- deriv(FIh, "Iv") > f21 <- deriv(FIv, "Ih") > f22 <- d...

...# MODEL ESTIMATION ### # build the big X matrix needed for GLS estimation: X <- kronecker(diag(1,3),x) Y <- c(y) # stack the y in a vector # residual covariance matrix for each observation covar <- kronecker(sigma,diag(1,N)) # GLS betas covariance matrix inv.sigma <- solve(covar) betav <- solve(t(X)%*%inv.sigma%*%X) # GLS mean parameter estimates betam <- betav%*%t(X)%*%inv.sigma%*%Y ----- end of code ----

Hi all!!! How can I create a legend to my plot that automatically reads the values from two vectors? betav<-c(0.78,0.94,0.88,0.41,0.59,4.68) etav<-c(235.6,59.5,31.2,8.7,3.2,1174) I want my legend to show 6 lines beta(in greeks)=0.78,eta(in greeks)=235.6 beta(in greeks)=0.94,eta(in greeks)=6.59 .... beta(in greeks)=4.68,eta(in greeks)=1174 I know how to do it "by hand" but I have to wri...

...lation growth rate using ?popbio?. I would like to perform a stochastic model in which my demographic rates are sampled from beta distribution with known mean and SD. For example, rather than 0.50 (4th term of my matrix), I want to sample this parameter from a beta function at each time-step {e.g., betaval(0.5, 0.05)}. Then I want to run the model for 1000 time steps (or more) and estimate the asymptotic stochastic growth rate. Is there a command that enable this kind of modelling? Thank you, Best wishes, Privalan -- View this message in context: http://www.nabble.com/stochastic-growth-rate-%2...