Displaying 1 result from an estimated 1 matches for "kidsiz".
Did you mean:
kidsize
2001 Nov 05
1
Why doesn't outer work?
...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) {
nkids<-p.est*exp(p.vec[8]+p.vec[9]*x);
kidsize.mean<- p.vec[10];
kidsize.var<- p.vec[11];
fac1<-sqrt(2*pi)*sqrt(kidsize.var);
fac2<-((y-kidsize.mean)^2)/(2*kidsize.var);
f<-fx(x)*nkids*exp(-fac2)/fac1;
return(f);
}
# Part (III)
################################################################
# THE KERNEL K(y,x)
# note th...