Displaying 2 results from an estimated 2 matches for "stdr".
Did you mean:
std
2006 Sep 19
1
-Need help with function
....8
a<-B0/R0*(1-(z-0.2)/(0.8*z))
b<-(z-0.2)/(0.8*z*R0)
dat<-data.frame(years=seq(1983,2004),cobs=c(19032324,19032324,17531618,2
0533029,20298099,20793744,23519369,23131780,19922247,17274513,17034419,1
2448318,4551585,4226451,7183688,7407924,7538366,7336039,8869193,7902341,
6369089,6211886))
stdr<-runif(100,0,0.5)
stdc<-runif(100,0,0.5)
BC<-runif(1000,0,100)
#model calibration
calib<-function(x){
v<-sample(stdr,1)
cr<-sample(stdc,1)
N<-rnorm(1)
Bq<-sample(BC,1)
Rerr<-exp(N*v-(v^2/2))
Cerr<-exp(N*cr-(cr^2/2))
Bt<-vector();Bt[1]=B0;Bt[2]=B0
Rt<-v...
2006 Sep 20
1
help with function
....8
a<-B0/R0*(1-(z-0.2)/(0.8*z))
b<-(z-0.2)/(0.8*z*R0)
dat<-data.frame(years=seq(1983,2004),cobs=c(19032324,19032324,17531618,2
0533029,20298099,20793744,23519369,23131780,19922247,17274513,17034419,1
2448318,4551585,4226451,7183688,7407924,7538366,7336039,8869193,7902341,
6369089,6211886))
stdr<-runif(100,0,0.5)
stdc<-runif(100,0,0.5)
BC<-runif(1000,0,100)
#model calibration
calib<-function(x){
v<-sample(stdr,1)
cr<-sample(stdc,1)
N<-rnorm(1)
Bq<-sample(BC,1)
Rerr<-exp(N*v-(v^2/2))
Cerr<-exp(N*cr-(cr^2/2))
Bt<-vector();Bt[1]=B0;Bt[2]=B0
Rt<-v...