error.R 1.31 KB
Newer Older
Jean-Marie Lepioufle's avatar
Jean-Marie Lepioufle committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76

#' error_input

#' error_input
#' @param x true input
#' @param phi_f1 phi_f1
#' @keywords aqoslogen
#' @export
#' @examples
#' \dontrun{
#' error_input()
#' }
#
error_input <- function(x,phi_f1=3/100){
    res <- x
    n <- length(res)
    res <- stats::rnorm(n,x,sd=abs(phi_f1)*x)
    return(res)
}


#' error_output

#' error_output
#' @param y true output
#' @param phi_g1 phi_g1
#' @keywords aqoslogen
#' @export
#' @examples
#' \dontrun{
#' error_output()
#' }
#
error_output <- function(y,phi_g1=3/100){
    res <- y
    n <- length(res)
    res <- stats::rnorm(n,y,sd=abs(phi_g1)*y)
    return(res)
}

#' error_structural

#' error_structural
#' @param m true model output
#' @param phi_h1 phi_h1
#' @keywords aqoslogen
#' @export
#' @examples
#' \dontrun{
#' error_structural()
#' }
#
error_structural <- function(m,phi_h1=8/100){
  n <- length(m)
  res <- stats::rnorm(n,m,sd=abs(phi_h1)*m)
  return(res)
}

#' error_remnant

#' error_remnant
#' @param y_hat true prediction
#' @param beta_r beta_r
#' @param sigma_r sigma_r
#' @keywords aqoslogen
#' @export
#' @examples
#' \dontrun{
#' error_remnant()
#' }
#
error_remnant <- function(y_hat,beta_r=0.8,sigma_r=0.5){
  n <- length(y_hat)
  res <- 1/beta_r*stats::rnorm(n,y_hat,sigma_r)
  return(res)
}