Usage Arguments Details Value Author(s) Examples
1 2 |
ions |
List of dataframes corresponding to ion. |
prev |
Numeric vector of observed prevalence values. |
model |
String, "NTE" or "TE" for non-targeted or targeted effects. |
w |
Numeric vector of experimental weights. |
Weight vector elements should correspond to dataframe element order. i.e. w[n] = ions[length(ions[:, 1]) / n][length(ions[:, 1])
Tested for the two examples below.
Numeric weighted mean errors.
Edward Greg Huang <eghuang@berkeley.edu>
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | ##=================== Cross validation ====================#
# Seperate Data into 8 blocks, i.e. test/training sets:
HZE_data <- select(filter(ion_data, Z > 3), 1:length(ion_data))
data_len <- 1:length(HZE_data)
O_350 <- select(filter(HZE_data, LET == 20), data_len)
Ne_670 <- select(filter(HZE_data, LET == 25), data_len)
Si_260 <- select(filter(HZE_data, LET == 70), data_len)
Ti_1000 <- select(filter(HZE_data, LET == 100), data_len)
Fe_600 <- select(filter(HZE_data, LET == 193), data_len)
Fe_350 <- select(filter(HZE_data, LET == 253), data_len)
Nb_600 <- select(filter(HZE_data, LET == 464), data_len)
La_593 <- select(filter(HZE_data, LET == 953), data_len)
set_list <- list(O_350, Ne_670, Si_260, Ti_1000,
Fe_600, Fe_350, Nb_600, La_593)
actual_prev <- HZE_data$Prev
NTE_cv <- loocv(set_list,"NTE", HZE_data$Prev, HZE_data$NWeight)
TE_cv <- loocv(set_list, "TE", HZE_data$Prev, HZE_data$NWeight)
cv_table <- cbind(NTE_cv, TE_cv)
cv_table
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.