- for (i in 1:nf)
- {
- avg_param <- as.double(params_hat[[i]])
- std_param <- as.double(stdev[[i]])
- matplot(cbind(params[o],avg_param[o],avg_param[o]+std_param[o],avg_param[o]-std_param[o]),
- col=1, lty=c(1,5,2,2), type="l", lwd=2, xlab="param", ylab="value")
- }
+ avg_param <- as.double(params_hat)
+ std_param <- as.double(stdev)
+ matplot(cbind(params[o],avg_param[o],avg_param[o]+std_param[o],avg_param[o]-std_param[o]),
+ col=1, lty=c(1,5,3,3), type="l", lwd=2, xlab="Parameter", ylab="Value")