1 | wrap(a)
|
a |
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 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (a)
{
Intensity2 = Enjoyment2 = Activity2 = Fatigue2 = Drowsy2 = Sleep2 = Thinking2 = Constipation2 = Sharpness2 = Hotness2 = Sensitivity2 = Block2 = Day2 = rep(NA,
nrow(a$data))
for (i in 1:nrow(a$data)) {
Day2[i] = as.numeric(as.Date(a$data[[1]][i]) - as.Date(a$metadata[[3]]))
Block2[i] = a$data[[3]][i]
Intensity2[i] = a$data[[4]][i]
Enjoyment2[i] = a$data[[5]][i]
Activity2[i] = a$data[[6]][i]
Fatigue2[i] = a$data[[7]][i]
Drowsy2[i] = a$data[[8]][i]
Sleep2[i] = a$data[[9]][i]
Thinking2[i] = a$data[[10]][i]
Constipation2[i] = a$data[[11]][i]
Sharpness2[i] = a$data[[12]][i]
Hotness2[i] = a$data[[13]][i]
Sensitivity2[i] = a$data[[14]][i]
}
Day2 = Day2 + 1
Pain2 = Intensity2 + Enjoyment2 + Activity2
Neuropain2 = Sharpness2 + Hotness2 + Sensitivity2
Treat2 = rep(NA, nrow(a$data))
for (i in 1:nrow(a$data)) {
if (!is.na(a$data[[2]][i]) & a$data[[2]][i] == "A") {
Treat2[i] = 0
}
if (!is.na(a$data[[2]][i]) & a$data[[2]][i] == "B") {
Treat2[i] = 1
}
}
Covs2 = cbind(Day2, Block2)
for (i in 1:length(Pain2)) {
if (Pain2[i] < 1) {
Pain2[i] = 1
}
if (Pain2[i] > 30) {
Pain2[i] = 30
}
if (Fatigue2[i] < 1) {
Fatigue2[i] = 1
}
if (Fatigue2[i] > 4) {
Fatigue2[i] = 4
}
if (Drowsy2[i] < 1) {
Drowsy2[i] = 1
}
if (Drowsy2[i] > 5) {
Drowsy2[i] = 5
}
if (Sleep2[i] < 1) {
Sleep2[i] = 1
}
if (Sleep2[i] > 4) {
Sleep2[i] = 4
}
if (Thinking2[i] < 1) {
Thinking2[i] = 1
}
if (Thinking2[i] > 4) {
Thinking2[i] = 4
}
if (Constipation2[i] < 1) {
Constipation2[i] = 1
}
if (Constipation2[i] > 4) {
Constipation2[i] = 4
}
if (Neuropain2[i] < 1) {
Neuropain2[i] = 1
}
if (Neuropain2[i] > 30) {
Neuropain2[i] = 30
}
}
observations = cbind(Day2, Pain2, Fatigue2, Drowsy2, Sleep2,
Thinking2, Constipation2, Neuropain2, Treat2, Block2)
logit = function(x) log(x/(1 - x))
nof1 = analyze(Pain = Pain2, Fatigue = Fatigue2, Drowsy = Drowsy2,
Sleep = Sleep2, Thinking = Thinking2, Constipation = Constipation2,
Neuropain = Neuropain2, Treat = Treat2, score.range = c(30,
4, 5, 4, 4, 4, 30), Covs = NULL, slopeprior = list("norm",
0, 0.1), nChains = 3, conv.limit = 1.05, niters = 10000,
setsize = 1000, alphaprior = list("norm", 0, 1e-06),
beta.norm.prior = list("norm", 0, 1e-06), beta.ord.prior = list("norm",
0, 1e-06), dcprior = list("unif", 0, 20), c1prior = list("unif",
-20, 20), varprior = list("Sd", "unif"), varprior.params = c(0,
5), path = "")
P025 = t(cbind(nof1$Pain$interval$P025, nof1$Fatigue$interval$P025,
nof1$Drowsy$interval$P025, nof1$Sleep$interval$P025,
nof1$Thinking$interval$P025, nof1$Constipation$interval$P025,
nof1$"Neuropathic Pain"$interval$P025))
Median = t(cbind(nof1$Pain$interval$Median, nof1$Fatigue$interval$Median,
nof1$Drowsy$interval$Median, nof1$Sleep$interval$Median,
nof1$Thinking$interval$Median, nof1$Constipation$interval$Median,
nof1$"Neuropathic Pain"$interval$Median))
P975 = t(cbind(nof1$Pain$interval$P975, nof1$Fatigue$interval$P975,
nof1$Drowsy$interval$P975, nof1$Sleep$interval$P975,
nof1$Thinking$interval$P975, nof1$Constipation$interval$P975,
nof1$"Neuropathic Pain"$interval$P975))
P975 = t(cbind(nof1$Pain$interval$P975, nof1$Fatigue$interval$P975,
nof1$Drowsy$interval$P975, nof1$Sleep$interval$P975,
nof1$Thinking$interval$P975, nof1$Constipation$interval$P975,
nof1$"Neuropathic Pain"$interval$P975))
Prob1 = t(cbind(nof1$Pain$probs$"Proportion < -0.2", nof1$Fatigue$probs$"Proportion < -0.2",
nof1$Drowsy$probs$"Proportion < -0.2", nof1$Sleep$probs$"Proportion < -0.2",
nof1$Thinking$probs$"Proportion < -0.2", nof1$Constipation$probs$"Proportion < -0.2",
nof1$"Neuropathic Pain"$probs$"Proportion < -0.2"))
Prob2 = t(cbind(nof1$Pain$probs$"Proportion -0.2 - 0", nof1$Fatigue$probs$"Proportion -0.2 - 0",
nof1$Drowsy$probs$"Proportion -0.2 - 0", nof1$Sleep$probs$"Proportion -0.2 - 0",
nof1$Thinking$probs$"Proportion -0.2 - 0", nof1$Constipation$probs$"Proportion -0.2 - 0",
nof1$"Neuropathic Pain"$probs$"Proportion -0.2 - 0"))
Prob3 = t(cbind(nof1$Pain$probs$"Proportion 0 - 0.2", nof1$Fatigue$probs$"Proportion 0 - 0.2",
nof1$Drowsy$probs$"Proportion 0 - 0.2", nof1$Sleep$probs$"Proportion 0 - 0.2",
nof1$Thinking$probs$"Proportion 0 - 0.2", nof1$Constipation$probs$"Proportion 0 - 0.2",
nof1$"Neuropathic Pain"$probs$"Proportion 0 - 0.2"))
Prob4 = t(cbind(nof1$Pain$probs$"Proportion > 0.2", nof1$Fatigue$probs$"Proportion > 0.2",
nof1$Drowsy$probs$"Proportion > 0.2", nof1$Sleep$probs$"Proportion > 0.2",
nof1$Thinking$probs$"Proportion > 0.2", nof1$Constipation$probs$"Proportion > 0.2",
nof1$"Neuropathic Pain"$probs$"Proportion > 0.2"))
Results = cbind(P025, Median, P975, Prob1, Prob2, Prob3,
Prob4)
colnames(Results) <- c("P025", "Median", "P975", "P(< - 0.2)",
"P(-0.2 - 0)", "P(0 - 0.2)", "P(> 0.2)")
rownames(Results) <- c("Pain", "Fatigue", "Drowsy", "Sleep",
"Thinking", "Constipation", "Neuropain")
Results_mod = matrix(NA, ncol = 10, nrow = 7)
for (i in 1:nrow(Results_mod)) {
if (Results[i, 2] < 0) {
Results_mod[i, 1] = "B"
}
if (Results[i, 2] > 0) {
Results_mod[i, 1] = "A"
}
if (Results[i, 2] == 0) {
Results_mod[i, 1] = "Neither"
}
if (Results[i, 3] < 0) {
Results_mod[i, 4] = "B"
}
if (Results[i, 3] > 0) {
Results_mod[i, 4] = "A"
}
if (Results[i, 1] < 0) {
Results_mod[i, 6] = "B"
}
if (Results[i, 1] > 0) {
Results_mod[i, 6] = "A"
}
Results_mod[i, 2] = abs(Results[i, 2])
Results_mod[i, 3] = abs(Results[i, 3])
Results_mod[i, 5] = abs(Results[i, 1])
Results_mod[i, 7] = Results[i, 4]
Results_mod[i, 8] = Results[i, 5]
Results_mod[i, 9] = Results[i, 7]
Results_mod[i, 10] = Results[i, 6]
}
colnames(Results_mod) = c("more_effective_regimen", "median_effect",
"upper_bound", "upper_bound_regimen", "lower_bound",
"lower_bound_regimen", "b_clinically_better", "b_marginally_better",
"a_clinically_better", "a_marginally_better")
rownames(Results_mod) = c("pain", "fatigue", "drowsiness",
"sleep_problems", "thinking_problems", "constipation",
"neuropathic_pain")
graph_5 = list(more_effective_regimen = Results_mod[1, 1],
median_effect = as.numeric(Results_mod[1, 2]), upper_bound = as.numeric(Results_mod[1,
3]), upper_bound_regimen = Results_mod[1, 4], lower_bound = as.numeric(Results_mod[1,
5]), lower_bound_regimen = Results_mod[1, 6])
graph_6 = list(b_clinically_better = as.numeric(Results_mod[1,
7]), b_marginally_better = as.numeric(Results_mod[1,
8]), a_clinically_better = as.numeric(Results_mod[1,
9]), a_marginally_better = as.numeric(Results_mod[1,
10]))
pain = list(graph_5, graph_6)
names(pain) = c("graph_5", "graph_6")
graph_5 = list(more_effective_regimen = Results_mod[4, 1],
median_effect = as.numeric(Results_mod[4, 2]), upper_bound = as.numeric(Results_mod[4,
3]), upper_bound_regimen = Results_mod[4, 4], lower_bound = as.numeric(Results_mod[4,
5]), lower_bound_regimen = Results_mod[4, 6])
graph_6 = list(b_clinically_better = as.numeric(Results_mod[4,
7]), b_marginally_better = as.numeric(Results_mod[4,
8]), a_clinically_better = as.numeric(Results_mod[4,
9]), a_marginally_better = as.numeric(Results_mod[4,
10]))
sleep_problems = list(graph_5, graph_6)
names(sleep_problems) = c("graph_5", "graph_6")
graph_5 = list(more_effective_regimen = Results_mod[6, 1],
median_effect = as.numeric(Results_mod[6, 2]), upper_bound = as.numeric(Results_mod[6,
3]), upper_bound_regimen = Results_mod[6, 4], lower_bound = as.numeric(Results_mod[6,
5]), lower_bound_regimen = Results_mod[6, 6])
graph_6 = list(b_clinically_better = as.numeric(Results_mod[6,
7]), b_marginally_better = as.numeric(Results_mod[6,
8]), a_clinically_better = as.numeric(Results_mod[6,
9]), a_marginally_better = as.numeric(Results_mod[6,
10]))
constipation = list(graph_5, graph_6)
names(constipation) = c("graph_5", "graph_6")
graph_5 = list(more_effective_regimen = Results_mod[3, 1],
median_effect = as.numeric(Results_mod[3, 2]), upper_bound = as.numeric(Results_mod[3,
3]), upper_bound_regimen = Results_mod[3, 4], lower_bound = as.numeric(Results_mod[3,
5]), lower_bound_regimen = Results_mod[3, 6])
graph_6 = list(b_clinically_better = as.numeric(Results_mod[3,
7]), b_marginally_better = as.numeric(Results_mod[3,
8]), a_clinically_better = as.numeric(Results_mod[3,
9]), a_marginally_better = as.numeric(Results_mod[3,
10]))
drowsiness = list(graph_5, graph_6)
names(drowsiness) = c("graph_5", "graph_6")
graph_5 = list(more_effective_regimen = Results_mod[5, 1],
median_effect = as.numeric(Results_mod[5, 2]), upper_bound = as.numeric(Results_mod[5,
3]), upper_bound_regimen = Results_mod[5, 4], lower_bound = as.numeric(Results_mod[5,
5]), lower_bound_regimen = Results_mod[5, 6])
graph_6 = list(b_clinically_better = as.numeric(Results_mod[5,
7]), b_marginally_better = as.numeric(Results_mod[5,
8]), a_clinically_better = as.numeric(Results_mod[5,
9]), a_marginally_better = as.numeric(Results_mod[5,
10]))
thinking_problems = list(graph_5, graph_6)
names(thinking_problems) = c("graph_5", "graph_6")
graph_5 = list(more_effective_regimen = Results_mod[2, 1],
median_effect = as.numeric(Results_mod[2, 2]), upper_bound = as.numeric(Results_mod[2,
3]), upper_bound_regimen = Results_mod[2, 4], lower_bound = as.numeric(Results_mod[2,
5]), lower_bound_regimen = Results_mod[2, 6])
graph_6 = list(b_clinically_better = as.numeric(Results_mod[2,
7]), b_marginally_better = as.numeric(Results_mod[2,
8]), a_clinically_better = as.numeric(Results_mod[2,
9]), a_marginally_better = as.numeric(Results_mod[2,
10]))
fatigue = list(graph_5, graph_6)
names(fatigue) = c("graph_5", "graph_6")
graph_5 = list(more_effective_regimen = Results_mod[7, 1],
median_effect = as.numeric(Results_mod[7, 2]), upper_bound = as.numeric(Results_mod[7,
3]), upper_bound_regimen = Results_mod[7, 4], lower_bound = as.numeric(Results_mod[7,
5]), lower_bound_regimen = Results_mod[7, 6])
graph_6 = list(b_clinically_better = as.numeric(Results_mod[7,
7]), b_marginally_better = as.numeric(Results_mod[7,
8]), a_clinically_better = as.numeric(Results_mod[7,
9]), a_marginally_better = as.numeric(Results_mod[7,
10]))
neuropathic_pain = list(graph_5, graph_6)
names(neuropathic_pain) = c("graph_5", "graph_6")
y = list(pain, sleep_problems, constipation, drowsiness,
thinking_problems, fatigue, neuropathic_pain)
names(y) = c("pain", "sleep_problems", "constipation", "drowsiness",
"thinking_problems", "fatigue", "neuropathic_pain")
return(y)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.