On 9/17/19 2:08 PM, David Winsemius wrote:
On 9/17/19 1:35 PM, varin sacha wrote:
Many thanks David, it perfectly works.
Now, one last think.
If I want my R code here below to run let's say B=500 times and at
the end I want to get the average for the MSE_GAM and for the
MSE_MARS. How can I do that ?
The `replicate` function is designed for that purpose.
Although I also just noticed that you were separately computing
residuals. Many R regression functions return a residual vector. Your
code would be a lot faster over the course of 500 repeats if you used
the resid function:
> str( resid(gam_model))
num [1:2000] -0.1385 0.1848 -0.0567 0.0605 -0.3297 ...
> str( resid(mars_model))
num [1:2000, 1] -0.2181 0.294 -0.0773 0.1626 -0.3512 ...
- attr(*, "dimnames")=List of 2
..$ : chr [1:2000] "1" "2" "3" "4" ...
..$ : chr "y_obs"
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.