search for: betavers

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

Dear helpers, I wrote a code which estimates a multi-equation model with generalized least squares (GLS). I can use GLS because I know the covariance matrix of the residuals a priori. However, it is a bit slow and I wonder if anybody would be able to point out a way to make it faster (it is part of a bigger code and needs to run several times). Any suggestion would be greatly appreciated. Carlo

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

Dear all, I am running matrix population models using package "popbio". In a deterministic model {i.e., transition matrix is defined as A <- matrix(c(0.70, 0.70,0.35,0.50), nrow=2,byrow=TRUE}, population growth rate can be estimated from the dominant eigenvalue {command "eigen.analysis"}. However, I cannot figure out the way to compute the asymptotic stochastic population