Displaying 1 result from an estimated 1 matches for "ziegelmann".
2010 Nov 18
0
On efficiency, Vectorize and loops
...that 'inner.product' operation 95*95*5 times on each of my B
bootstrap loops, where B should be 1000 but is more likely to be 500 or 200
if I'm willing to graduate before March!
# 'code0' is used for simulation and test purposes only;
# The model is the same as in Bathia, Yao, Ziegelmann p.13
# The code used for estimation is 'code1'
rm(list=ls(all=TRUE))
n<-100 # Number of observations
d<-4 # Dimension of X
a<-0 # Lower bound
b<-1 # Upper bound
p<-5 # [See Bathia, Yao, Ziegelmann p. 5)
#######################################
# Building the error...