Displaying 1 result from an estimated 1 matches for "sigmax2".
Did you mean:
sigma2
2001 Nov 05
1
Why doesn't outer work?
...function
sx<-function(x) {
u<-p.vec[2]*x+p.vec[1];
if(u>=0.7) u<-0.7;
if(u<0) u<-0;
return(u);
}
#flowering functiom
fx<-function(x) {
u<-exp(p.vec[3]+p.vec[4]*x);
return(u/(1+u));
}
#growth function
gxy<-function(x,y) {
sigmax2<-p.vec[7];
sigmax<-sqrt(sigmax2);
mux<-p.vec[5]+p.vec[6]*x;
fac1<-sqrt(2*pi)*sigmax;
fac2<-((y-mux)^2)/(2*sigmax2);
return(exp(-fac2)/fac1);
}
#survival-growth function
pxy<-function(x,y) { return(sx(x)*(1-fx(x))*gxy(x,y)) }
#fecundity function
fxy<-function(x,y) {
nk...