OptimalRAR/Data_Generator.R at main · lukaspinpin/OptimalRAR · GitHub

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
generate_data <- function(N=200, dist="bern", parV = 0.5){
  if(dist=="bern"){
    return(rbinom(N,1, parV))
  }
  if(dist=="norm"){
    return(rnorm(n=N,mean=parV[1],sd=parV[2]))
  }
  if(dist=="Lickert"){
    return(round(rbeta(N, 5, parV)*5)+1)
  }
  if(dist=="expon"){
    return(rexp(n = N, rate = parV))
  }
  if(dist=="beta"){
    return(rbeta(N, 5, parV))
  }
}

superior <- function(par1=c(0,0), par2=FALSE, dist="bern"){
  if(dist=="bern"){
    if(par1[1]<par1[2]){
      return(1)
    }
    if(par1[1]>par1[2]){
      return(0)
    }
    if(par1[1]==par1[2]){
      return(2)
    }
  }
  if(dist=="norm"){
    if(par1[1]<par1[2]){
      return(1)
    }
    if(par1[1]>par1[2]){
      return(0)
    }
    if(par1[1]==par1[2]){
      return(2)
    }
  }
  if(dist=="Lickert" || dist=="beta"){
    exp1 <- 5/(5+par1[1])
    exp2 <- 5/(5+par1[2])
    if(exp1<exp2){
      return(1)
    }
    if(exp1>exp2){
      return(0)
    }
    if(exp1==exp2){
      return(2)
    }
  }
  if(dist=="expon"){
    if(par1[1]>par1[2]){
      return(1)
    }
    if(par1[1]<par1[2]){
      return(0)
    }
    if(par1[1]==par1[2]){
      return(2)
    }
  }
}