Homework 10
Jasmine Hart
3/30/2022
vect <- function() {
count1 <- sample(c(0,1), replace=TRUE, size=10)
counter <- 0
for (i in count1) {
ifelse(i==0, counter<-counter+1, counter<-counter+0)
}
return(counter)
}
vect()## [1] 5count2 <- sample(c(0,1), replace=TRUE, size=10)
count3 <- length(which(count2==0))
print(count3)## [1] 6mat <- matrix(runif(16),
nrow=4,
ncol=4)
print(mat)## [,1] [,2] [,3] [,4]
## [1,] 0.6800183 0.8703001 0.3258383 0.9753394
## [2,] 0.4563434 0.8178313 0.1030108 0.7285098
## [3,] 0.9123433 0.3431711 0.2840801 0.8335479
## [4,] 0.1663272 0.4095295 0.3502869 0.0382572for (i in 1:nrow(mat)) {
for (j in 1:ncol(mat)) {
mat[i,j] <- mat[i,j] + i*j
} # end of column j loop
} # end or row i loop
print(mat)## [,1] [,2] [,3] [,4]
## [1,] 1.680018 2.870300 3.325838 4.975339
## [2,] 2.456343 4.817831 6.103011 8.728510
## [3,] 3.912343 6.343171 9.284080 12.833548
## [4,] 4.166327 8.409529 12.350287 16.0382574a.
library(tidyverse)## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ──## ✓ ggplot2 3.3.5 ✓ purrr 0.3.4
## ✓ tibble 3.1.6 ✓ dplyr 1.0.7
## ✓ tidyr 1.1.4 ✓ stringr 1.4.0
## ✓ readr 2.1.1 ✓ forcats 0.5.1## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()dat1 <- rnorm(n=15,mean=12,sd=1)
dat2 <- rnorm(n=15,mean=15,sd=1)
dat3 <- rnorm(n=15,mean=18,sd=1)
wholedat <- data.frame(dat1,dat2,dat3)
print(wholedat)## dat1 dat2 dat3
## 1 11.66178 14.96880 16.93555
## 2 11.72835 15.57507 17.85253
## 3 12.20975 14.12457 15.60680
## 4 10.77547 15.41338 17.28258
## 5 10.74713 14.43650 16.79180
## 6 11.62909 14.83959 17.64103
## 7 12.51621 14.59855 17.68594
## 8 12.51599 14.84876 18.61008
## 9 12.41881 14.28683 18.38705
## 10 11.34216 16.48631 17.52474
## 11 12.34568 14.15348 18.82061
## 12 12.23963 15.02303 17.27339
## 13 11.67485 14.82310 18.25707
## 14 13.61367 14.69968 19.23834
## 15 12.03355 15.49446 18.69561wholedat1 <- pivot_longer(wholedat,cols=dat1:dat3, names_to="DataSet",values_to="Response")
print(wholedat1)## # A tibble: 45 × 2
## DataSet Response
## <chr> <dbl>
## 1 dat1 11.7
## 2 dat2 15.0
## 3 dat3 16.9
## 4 dat1 11.7
## 5 dat2 15.6
## 6 dat3 17.9
## 7 dat1 12.2
## 8 dat2 14.1
## 9 dat3 15.6
## 10 dat1 10.8
## # … with 35 more rows4b.
###################################
# FUNCTION shuff_dat
# purpose: shuffles and randomizes the response variables and calculates means for each group
# input: created data set, wholedat1
# output: shuffled datasets and means
# --------------------------------------
shuff_dat <- function(x) {
shu <- transform(x, shuffle=sample(Response))
dat1 <- mean(shu$shuffle[shu$DataSet=="dat1"])
dat2 <- mean(shu$shuffle[shu$DataSet=="dat2"])
dat3 <- mean(shu$shuffle[shu$DataSet=="dat3"])
allmeans <- data.frame(dat1,dat2,dat3)
return(allmeans)
}
shuff_dat(wholedat1)## dat1 dat2 dat3
## 1 14.47955 15.33457 14.841044c.
df <- data.frame()
rep <- c(1:1000)
for(i in 1:1000) {
meanz <- shuff_dat(wholedat1)
df <- rbind(df,meanz)
}
datfram <- data.frame(rep,df)
print(datfram)## rep dat1 dat2 dat3
## 1 1 15.22314 15.10388 14.32814
## 2 2 14.39137 14.94410 15.31968
## 3 3 14.33129 15.22896 15.09490
## 4 4 14.95440 14.42132 15.27943
## 5 5 16.08471 15.43324 13.13721
## 6 6 14.92777 14.70916 15.01822
## 7 7 14.55079 13.93987 16.16450
## 8 8 14.69430 15.06446 14.89640
## 9 9 15.11657 14.08682 15.45176
## 10 10 14.41320 14.38663 15.85532
## 11 11 14.88801 14.84863 14.91851
## 12 12 14.60329 14.47476 15.57710
## 13 13 14.57294 15.40022 14.68199
## 14 14 15.24378 14.23947 15.17191
## 15 15 14.17531 15.15151 15.32834
## 16 16 15.14114 15.11261 14.40141
## 17 17 14.74201 14.49602 15.41713
## 18 18 15.08123 14.35161 15.22232
## 19 19 15.11073 15.23229 14.31214
## 20 20 14.80709 15.70166 14.14641
## 21 21 15.11684 14.68598 14.85233
## 22 22 15.41839 14.47431 14.76246
## 23 23 15.40017 14.02004 15.23494
## 24 24 14.04710 15.53536 15.07270
## 25 25 14.93190 15.23323 14.49003
## 26 26 14.69635 14.84079 15.11801
## 27 27 15.46461 14.45911 14.73143
## 28 28 13.76175 15.82878 15.06463
## 29 29 14.68262 15.31291 14.65963
## 30 30 15.60989 14.63422 14.41105
## 31 31 15.42941 14.93486 14.29088
## 32 32 14.08705 15.07412 15.49398
## 33 33 15.18172 15.25970 14.21374
## 34 34 13.87437 16.20759 14.57319
## 35 35 15.11657 14.55519 14.98340
## 36 36 14.40204 15.55501 14.69811
## 37 37 14.37877 14.95429 15.32209
## 38 38 15.81274 13.68704 15.15538
## 39 39 14.20676 15.41179 15.03661
## 40 40 15.16923 14.63940 14.84653
## 41 41 14.10779 15.36122 15.18614
## 42 42 15.55436 14.90052 14.20028
## 43 43 14.08724 14.89428 15.67364
## 44 44 14.54743 15.08084 15.02688
## 45 45 14.66382 14.72287 15.26847
## 46 46 15.10506 14.52046 15.02963
## 47 47 15.59356 14.01332 15.04828
## 48 48 13.94113 15.11658 15.59744
## 49 49 15.27575 15.21525 14.16416
## 50 50 15.40256 15.24506 14.00753
## 51 51 14.67125 14.49808 15.48582
## 52 52 14.53597 15.13809 14.98109
## 53 53 14.85699 14.15412 15.64404
## 54 54 15.09928 14.11891 15.43697
## 55 55 15.11354 14.61726 14.92435
## 56 56 15.15052 16.16700 13.33763
## 57 57 15.69547 14.49565 14.46404
## 58 58 14.33584 14.71472 15.60460
## 59 59 15.19448 15.00770 14.45297
## 60 60 15.27391 14.41174 14.96951
## 61 61 15.22919 15.51658 13.90938
## 62 62 14.59158 15.01132 15.05226
## 63 63 14.99098 14.57553 15.08864
## 64 64 14.77139 15.64461 14.23916
## 65 65 14.21336 15.48445 14.95735
## 66 66 15.75281 14.23111 14.67123
## 67 67 14.59494 15.68076 14.37945
## 68 68 15.52282 14.59102 14.54131
## 69 69 14.47963 14.21443 15.96110
## 70 70 15.42396 14.94877 14.28242
## 71 71 15.22486 14.29399 15.13630
## 72 72 14.49355 15.22153 14.94008
## 73 73 15.32454 14.36455 14.96607
## 74 74 15.33631 14.21997 15.09888
## 75 75 14.31461 15.05942 15.28113
## 76 76 15.00950 14.54945 15.09620
## 77 77 15.56988 14.53913 14.54614
## 78 78 14.87432 15.29566 14.48518
## 79 79 15.30639 14.58379 14.76497
## 80 80 14.78210 14.78359 15.08946
## 81 81 15.25787 14.76547 14.63181
## 82 82 14.22547 15.08176 15.34793
## 83 83 14.75723 15.58683 14.31110
## 84 84 13.83635 15.08259 15.73621
## 85 85 13.35764 15.99001 15.30750
## 86 86 15.07581 13.97754 15.60181
## 87 87 14.95109 14.97448 14.72959
## 88 88 14.94757 15.14317 14.56441
## 89 89 14.95903 15.42779 14.26834
## 90 90 14.59399 14.45585 15.60531
## 91 91 15.72063 14.67991 14.25462
## 92 92 14.97968 14.89471 14.78077
## 93 93 15.63717 14.30660 14.71138
## 94 94 14.48755 14.92950 15.23811
## 95 95 14.16862 14.28775 16.19879
## 96 96 15.63435 13.84518 15.17563
## 97 97 14.82734 15.23537 14.59245
## 98 98 15.51848 14.55167 14.58501
## 99 99 15.15466 14.86166 14.63884
## 100 100 14.84345 15.31630 14.49540
## 101 101 14.47524 15.24416 14.93576
## 102 102 14.37803 14.69179 15.58534
## 103 103 14.78716 15.06156 14.80643
## 104 104 14.89751 15.02405 14.73360
## 105 105 15.04723 14.76326 14.84466
## 106 106 15.34404 13.76645 15.54467
## 107 107 14.99027 14.66923 14.99566
## 108 108 14.42036 15.02393 15.21087
## 109 109 14.64507 14.68903 15.32106
## 110 110 14.73410 15.37623 14.54482
## 111 111 14.18441 14.96181 15.50893
## 112 112 15.09007 15.83132 13.73377
## 113 113 14.65038 14.44886 15.55591
## 114 114 14.89529 14.14988 15.60999
## 115 115 14.71382 14.73511 15.20623
## 116 116 16.35799 14.57405 13.72312
## 117 117 15.81541 14.35928 14.48047
## 118 118 14.03927 15.95708 14.65880
## 119 119 14.69588 15.52639 14.43289
## 120 120 15.63729 13.67231 15.34555
## 121 121 14.63757 14.79299 15.22460
## 122 122 15.47371 14.75510 14.42635
## 123 123 14.79297 14.67738 15.18480
## 124 124 15.46768 14.97831 14.20917
## 125 125 14.51390 15.43391 14.70735
## 126 126 14.74966 14.54944 15.35605
## 127 127 15.16022 14.17431 15.32063
## 128 128 15.10858 14.71802 14.82855
## 129 129 14.43646 14.94758 15.27111
## 130 130 14.72922 14.94589 14.98004
## 131 131 15.47383 15.04308 14.13824
## 132 132 15.56931 14.48967 14.59618
## 133 133 14.97653 14.89297 14.78565
## 134 134 14.89306 14.82895 14.93314
## 135 135 14.50462 14.70143 15.44911
## 136 136 14.09405 14.98849 15.57261
## 137 137 14.65673 15.08809 14.91033
## 138 138 14.98724 14.06927 15.59864
## 139 139 14.15555 15.84400 14.65560
## 140 140 14.77231 14.30327 15.57957
## 141 141 15.26925 14.42898 14.95693
## 142 142 14.91029 14.41036 15.33450
## 143 143 15.25656 14.80557 14.59303
## 144 144 14.76633 14.88131 15.00751
## 145 145 15.18721 14.98164 14.48631
## 146 146 13.82650 14.58605 16.24261
## 147 147 14.82120 15.00274 14.83121
## 148 148 14.41096 15.55282 14.69138
## 149 149 14.75171 15.17846 14.72499
## 150 150 14.17935 15.81682 14.65898
## 151 151 15.91122 14.46618 14.27775
## 152 152 14.67790 14.99192 14.98533
## 153 153 15.58217 14.18978 14.88320
## 154 154 15.48022 15.06073 14.11421
## 155 155 14.82694 14.51215 15.31606
## 156 156 14.90880 15.36918 14.37718
## 157 157 15.72841 13.61870 15.30805
## 158 158 14.57798 15.34666 14.73052
## 159 159 15.37446 13.83484 15.44586
## 160 160 14.38568 15.25046 15.01902
## 161 161 15.32263 14.56505 14.76748
## 162 162 14.75920 15.15830 14.73766
## 163 163 14.70090 15.27991 14.67435
## 164 164 15.16076 14.32018 15.17422
## 165 165 14.53224 15.31159 14.81133
## 166 166 14.97300 14.14445 15.53771
## 167 167 14.87466 15.84911 13.93138
## 168 168 15.96542 14.29226 14.39747
## 169 169 14.93848 14.77518 14.94150
## 170 170 15.27871 14.57645 14.79999
## 171 171 14.98265 15.53053 14.14197
## 172 172 15.74967 14.78085 14.12463
## 173 173 15.34745 14.42416 14.88355
## 174 174 15.21986 14.64907 14.78623
## 175 175 15.45074 14.98399 14.22042
## 176 176 15.17262 13.96269 15.51985
## 177 177 13.79123 15.43075 15.43317
## 178 178 15.19346 14.63180 14.82990
## 179 179 14.46471 15.75993 14.43051
## 180 180 14.64149 15.20759 14.80608
## 181 181 14.68392 15.86842 14.10282
## 182 182 15.22180 14.61978 14.81357
## 183 183 14.44509 14.22405 15.98601
## 184 184 14.40363 14.47915 15.77238
## 185 185 15.50226 14.10465 15.04824
## 186 186 15.17174 14.69115 14.79227
## 187 187 14.29879 15.30504 15.05132
## 188 188 14.11688 15.35005 15.18822
## 189 189 14.79593 14.27639 15.58283
## 190 190 16.04579 15.13180 13.47757
## 191 191 15.40859 13.96816 15.27841
## 192 192 15.86465 14.87064 13.91986
## 193 193 14.70952 15.17217 14.77347
## 194 194 14.85208 14.86962 14.93345
## 195 195 15.02134 14.02542 15.60839
## 196 196 13.68818 15.21545 15.75153
## 197 197 14.12006 15.03910 15.49600
## 198 198 15.17622 15.02863 14.45030
## 199 199 14.84402 15.61204 14.19910
## 200 200 15.16965 14.81447 14.67104
## 201 201 15.03042 14.93115 14.69359
## 202 202 15.66540 15.07901 13.91074
## 203 203 14.70709 14.80604 15.14202
## 204 204 14.96864 15.13310 14.55341
## 205 205 15.28636 15.03244 14.33636
## 206 206 15.08240 15.05145 14.52131
## 207 207 15.44275 14.79282 14.41959
## 208 208 15.13112 14.26712 15.25691
## 209 209 13.73556 15.73144 15.18815
## 210 210 15.83350 14.32027 14.50139
## 211 211 15.49418 15.52039 13.64059
## 212 212 14.81295 15.06247 14.77973
## 213 213 15.28959 14.43551 14.93006
## 214 214 15.15539 14.74757 14.75219
## 215 215 14.33768 14.61378 15.70369
## 216 216 15.32187 14.46669 14.86660
## 217 217 13.80922 15.35186 15.49408
## 218 218 15.06494 14.86642 14.72380
## 219 219 15.09732 14.27818 15.27966
## 220 220 14.70521 15.20403 14.74591
## 221 221 14.20824 15.27625 15.17066
## 222 222 14.89579 15.02216 14.73721
## 223 223 14.17825 15.82485 14.65206
## 224 224 13.91247 15.33227 15.41041
## 225 225 14.78163 14.63742 15.23610
## 226 226 15.56560 15.00484 14.08471
## 227 227 14.16026 14.19230 16.30259
## 228 228 14.30786 14.97201 15.37529
## 229 229 14.94711 15.74108 13.96696
## 230 230 15.03502 15.12953 14.49061
## 231 231 14.24900 15.36733 15.03883
## 232 232 15.01537 14.99795 14.64183
## 233 233 13.77055 15.59645 15.28816
## 234 234 14.70283 14.54464 15.40768
## 235 235 15.88255 14.28866 14.48395
## 236 236 15.46756 15.24522 13.94237
## 237 237 15.28726 15.46038 13.90751
## 238 238 14.65028 14.42040 15.58447
## 239 239 14.77383 14.59754 15.28379
## 240 240 15.37972 14.81502 14.46041
## 241 241 14.39994 15.12713 15.12808
## 242 242 14.72118 14.53710 15.39687
## 243 243 15.27289 14.60182 14.78045
## 244 244 14.89607 15.20026 14.55883
## 245 245 14.18199 15.49238 14.98078
## 246 246 14.72165 15.22205 14.71145
## 247 247 15.10325 14.11786 15.43404
## 248 248 14.59832 15.12103 14.93581
## 249 249 16.44123 13.26868 14.94524
## 250 250 15.28520 15.12210 14.24785
## 251 251 14.43069 15.54804 14.67642
## 252 252 15.71576 14.63111 14.30829
## 253 253 14.84933 14.31607 15.48976
## 254 254 14.63385 15.75124 14.27006
## 255 255 15.28743 14.65388 14.71385
## 256 256 15.03977 15.14690 14.46849
## 257 257 15.13964 14.74811 14.76740
## 258 258 15.36844 15.27412 14.01260
## 259 259 14.73119 15.27290 14.65107
## 260 260 14.61357 14.80974 15.23184
## 261 261 15.48991 14.02083 15.14442
## 262 262 15.15820 14.64756 14.84939
## 263 263 15.11966 13.81173 15.72377
## 264 264 14.29353 15.10330 15.25833
## 265 265 15.18244 15.10492 14.36780
## 266 266 15.10650 14.49152 15.05713
## 267 267 14.55091 14.79059 15.31366
## 268 268 14.40633 15.77014 14.47869
## 269 269 13.72479 15.55191 15.37845
## 270 270 15.44414 14.05731 15.15370
## 271 271 13.78719 16.12171 14.74625
## 272 272 14.65161 14.88914 15.11441
## 273 273 15.35962 14.45371 14.84182
## 274 274 15.64585 15.25011 13.75920
## 275 275 14.38338 15.00538 15.26640
## 276 276 14.30169 15.06079 15.29267
## 277 277 14.08485 15.82393 14.74637
## 278 278 14.58350 15.19400 14.87765
## 279 279 15.26898 14.52483 14.86135
## 280 280 14.73621 15.42726 14.49168
## 281 281 14.58168 15.46160 14.61188
## 282 282 15.65644 13.74323 15.25549
## 283 283 16.01153 14.19023 14.45339
## 284 284 14.99004 14.91925 14.74586
## 285 285 15.51198 14.50696 14.63622
## 286 286 15.05923 15.09501 14.50092
## 287 287 15.13847 14.36447 15.15222
## 288 288 15.21914 14.21822 15.21780
## 289 289 13.90801 15.11567 15.63147
## 290 290 14.87875 14.03754 15.73887
## 291 291 14.87740 16.04724 13.73051
## 292 292 15.05735 15.00598 14.59183
## 293 293 15.46740 14.69685 14.49090
## 294 294 14.67830 16.19171 13.78514
## 295 295 14.88118 14.78403 14.98995
## 296 296 15.52180 14.95775 14.17561
## 297 297 13.97934 15.76049 14.91533
## 298 298 15.33926 14.13684 15.17905
## 299 299 15.35219 14.37611 14.92685
## 300 300 14.43857 15.16258 15.05400
## 301 301 15.22406 15.15895 14.27214
## 302 302 14.31588 14.83923 15.50004
## 303 303 16.15821 14.11646 14.38049
## 304 304 14.96574 15.00327 14.68614
## 305 305 14.96263 14.52106 15.17147
## 306 306 14.84843 15.46806 14.33867
## 307 307 15.68603 13.35246 15.61666
## 308 308 15.30691 14.51286 14.83538
## 309 309 15.08075 14.89397 14.68044
## 310 310 15.39421 15.22694 14.03401
## 311 311 13.58912 15.85757 15.20846
## 312 312 14.92892 15.07295 14.65328
## 313 313 14.18137 14.28582 16.18797
## 314 314 15.22057 14.86909 14.56550
## 315 315 14.48585 15.36661 14.80269
## 316 316 15.12774 13.67286 15.85456
## 317 317 14.72136 14.25503 15.67877
## 318 318 15.06501 14.56265 15.02750
## 319 319 15.17224 14.38181 15.10110
## 320 320 15.54603 14.15159 14.95753
## 321 321 14.91666 15.33313 14.40537
## 322 322 14.30878 15.49023 14.85614
## 323 323 13.96809 15.91975 14.76731
## 324 324 14.67004 15.44912 14.53600
## 325 325 14.60715 14.99420 15.05380
## 326 326 15.98355 13.95176 14.71984
## 327 327 15.15463 14.84634 14.65419
## 328 328 15.12857 14.06783 15.45876
## 329 329 15.24561 14.66479 14.74475
## 330 330 14.46549 14.95364 15.23602
## 331 331 15.53658 14.53938 14.57920
## 332 332 15.84681 14.58167 14.22667
## 333 333 14.96526 14.93556 14.75433
## 334 334 15.88560 14.42768 14.34187
## 335 335 14.75122 15.32885 14.57508
## 336 336 14.47156 14.39083 15.79277
## 337 337 13.94888 15.76402 14.94225
## 338 338 14.92352 15.17484 14.55679
## 339 339 15.86701 14.24306 14.54508
## 340 340 14.99138 15.08186 14.58192
## 341 341 15.06032 15.73932 13.85551
## 342 342 14.96853 14.74256 14.94407
## 343 343 15.46756 14.09425 15.09334
## 344 344 14.86824 14.33745 15.44946
## 345 345 14.57283 15.09412 14.98820
## 346 346 15.33799 14.62180 14.69537
## 347 347 14.54811 15.17833 14.92872
## 348 348 15.18988 15.46012 14.00516
## 349 349 14.25887 15.36559 15.03070
## 350 350 13.58156 15.31382 15.75977
## 351 351 14.48721 14.32494 15.84300
## 352 352 15.33138 13.84414 15.47963
## 353 353 15.22994 14.20797 15.21725
## 354 354 15.54824 13.81492 15.29200
## 355 355 15.00060 14.31495 15.33961
## 356 356 15.80814 14.02034 14.82668
## 357 357 14.86346 15.05727 14.73442
## 358 358 14.40956 14.72339 15.52221
## 359 359 13.19312 16.23516 15.22688
## 360 360 15.14573 14.08694 15.42249
## 361 361 15.44081 14.97052 14.24382
## 362 362 14.47106 14.92363 15.26046
## 363 363 14.56638 14.73952 15.34926
## 364 364 14.70578 14.87958 15.06979
## 365 365 14.86292 14.53638 15.25586
## 366 366 14.89839 14.66906 15.08771
## 367 367 14.20152 14.33614 16.11749
## 368 368 13.40672 15.37708 15.87135
## 369 369 15.20483 14.78084 14.66949
## 370 370 15.27864 13.99943 15.37708
## 371 371 15.75358 14.68453 14.21704
## 372 372 15.28410 14.56320 14.80786
## 373 373 15.39875 15.10320 14.15321
## 374 374 14.45534 15.27565 14.92417
## 375 375 15.33630 14.99825 14.32061
## 376 376 15.71312 14.42449 14.51755
## 377 377 15.85869 14.59252 14.20395
## 378 378 13.83580 15.10603 15.71332
## 379 379 15.99182 14.81260 13.85073
## 380 380 15.48583 15.57464 13.59468
## 381 381 14.53665 15.16677 14.95173
## 382 382 14.59332 14.53510 15.52674
## 383 383 14.81663 15.56354 14.27499
## 384 384 15.41533 14.38679 14.85304
## 385 385 15.70056 14.51271 14.44188
## 386 386 14.86964 15.71984 14.06567
## 387 387 15.56969 14.56398 14.52148
## 388 388 15.32773 14.82282 14.50461
## 389 389 14.74991 14.60075 15.30450
## 390 390 15.17435 14.61958 14.86122
## 391 391 14.90242 15.14691 14.60583
## 392 392 15.15515 15.03505 14.46496
## 393 393 14.05430 16.21038 14.39048
## 394 394 14.59166 14.74000 15.32349
## 395 395 14.93585 14.54807 15.17124
## 396 396 15.53947 14.58703 14.52866
## 397 397 14.89045 15.94099 13.82372
## 398 398 14.13435 15.71405 14.80675
## 399 399 15.35345 14.60618 14.69553
## 400 400 14.58645 15.77094 14.29777
## 401 401 14.00786 15.04481 15.60248
## 402 402 14.61889 14.12805 15.90822
## 403 403 14.60217 15.06462 14.98836
## 404 404 15.34242 15.44559 13.86715
## 405 405 14.92407 14.16672 15.56436
## 406 406 15.03548 15.69911 13.92057
## 407 407 13.84425 15.29706 15.51385
## 408 408 15.32632 14.79677 14.53207
## 409 409 14.48083 14.35764 15.81669
## 410 410 15.68208 14.63448 14.33860
## 411 411 15.17596 14.80688 14.67232
## 412 412 15.73077 14.63939 14.28499
## 413 413 14.95474 15.00787 14.69255
## 414 414 15.30927 14.47307 14.87281
## 415 415 15.73441 14.12191 14.79883
## 416 416 14.88723 15.31618 14.45175
## 417 417 16.50704 13.77246 14.37566
## 418 418 14.83000 14.92900 14.89616
## 419 419 16.03805 14.13862 14.47849
## 420 420 15.35767 14.53168 14.76581
## 421 421 14.69504 14.12364 15.83648
## 422 422 14.70811 15.35159 14.59546
## 423 423 13.96002 15.84002 14.85512
## 424 424 16.14066 13.69155 14.82295
## 425 425 15.14188 14.40152 15.11176
## 426 426 15.38924 13.99350 15.27242
## 427 427 15.14735 15.21389 14.29391
## 428 428 14.37866 14.56228 15.71421
## 429 429 15.12599 14.50291 15.02626
## 430 430 14.58826 16.12276 13.94413
## 431 431 14.67137 15.26350 14.72029
## 432 432 14.54013 14.11363 16.00139
## 433 433 15.53047 15.68521 13.43948
## 434 434 14.44734 15.30554 14.90228
## 435 435 14.21094 15.53866 14.90556
## 436 436 15.71998 14.25926 14.67592
## 437 437 14.23861 14.96553 15.45102
## 438 438 15.69832 14.46524 14.49160
## 439 439 14.67197 14.13252 15.85067
## 440 440 15.12649 14.44546 15.08320
## 441 441 14.80917 15.48315 14.36284
## 442 442 15.51453 14.72592 14.41470
## 443 443 15.27960 15.13985 14.23571
## 444 444 14.37966 15.04352 15.23198
## 445 445 15.11428 15.32110 14.21978
## 446 446 15.55918 15.23703 13.85894
## 447 447 15.47618 14.98685 14.19214
## 448 448 15.29860 15.26451 14.09204
## 449 449 15.15871 13.83821 15.65823
## 450 450 14.37189 14.52618 15.75709
## 451 451 15.02602 13.85051 15.77862
## 452 452 14.56435 15.38741 14.70339
## 453 453 14.45902 15.23084 14.96530
## 454 454 15.61108 14.98249 14.06159
## 455 455 15.48134 14.51320 14.66062
## 456 456 14.60461 14.53247 15.51807
## 457 457 14.92351 15.33486 14.39679
## 458 458 16.04967 13.46847 15.13702
## 459 459 15.00119 14.53908 15.11488
## 460 460 14.25650 15.47306 14.92559
## 461 461 14.11499 15.65107 14.88910
## 462 462 14.93636 14.30382 15.41497
## 463 463 15.25999 14.42525 14.96992
## 464 464 14.94591 14.36836 15.34088
## 465 465 14.53740 15.02672 15.09104
## 466 466 15.33630 13.72024 15.59862
## 467 467 15.42211 14.95697 14.27607
## 468 468 13.90430 15.45686 15.29400
## 469 469 14.42421 15.68322 14.54773
## 470 470 15.28401 14.42057 14.95058
## 471 471 15.41827 15.14860 14.08829
## 472 472 15.45738 14.67135 14.52643
## 473 473 15.09721 14.90110 14.65685
## 474 474 14.53305 15.14569 14.97641
## 475 475 14.77546 15.75294 14.12675
## 476 476 14.94454 15.19090 14.51972
## 477 477 14.36002 14.86272 15.43241
## 478 478 14.98550 14.02457 15.64508
## 479 479 14.76799 15.24390 14.64327
## 480 480 14.78951 15.10410 14.76155
## 481 481 15.80823 14.00498 14.84195
## 482 482 15.46379 14.48484 14.70653
## 483 483 15.29761 14.30780 15.04974
## 484 484 15.13524 14.36970 15.15022
## 485 485 14.95795 14.90159 14.79561
## 486 486 14.99708 14.19131 15.46677
## 487 487 15.05329 15.05251 14.54936
## 488 488 14.92099 14.75240 14.98177
## 489 489 14.70515 15.14544 14.80456
## 490 490 14.27661 15.95823 14.42031
## 491 491 14.90289 15.43538 14.31688
## 492 492 14.78358 15.49405 14.37753
## 493 493 15.03676 15.20220 14.41620
## 494 494 15.30741 15.54982 13.79792
## 495 495 15.48269 14.80916 14.36330
## 496 496 15.66683 15.37171 13.61661
## 497 497 14.82307 15.03686 14.79522
## 498 498 14.04511 15.51012 15.09992
## 499 499 14.75047 14.99853 14.90615
## 500 500 14.94685 14.41890 15.28940
## 501 501 14.96404 14.98372 14.70740
## 502 502 14.47729 14.89178 15.28609
## 503 503 13.76693 14.60336 16.28487
## 504 504 14.34993 15.32560 14.97963
## 505 505 15.58753 14.50734 14.56028
## 506 506 15.34937 14.23710 15.06869
## 507 507 15.37892 14.65707 14.61917
## 508 508 14.89480 14.85686 14.90349
## 509 509 15.26671 14.56845 14.82000
## 510 510 14.91492 14.82693 14.91331
## 511 511 14.84892 13.97416 15.83207
## 512 512 15.61368 14.14639 14.89509
## 513 513 14.12678 15.54344 14.98493
## 514 514 14.42879 15.07841 15.14796
## 515 515 15.61944 13.89747 15.13824
## 516 516 15.03849 14.84127 14.77540
## 517 517 14.78523 15.53637 14.33355
## 518 518 15.32758 14.83478 14.49280
## 519 519 15.48362 14.81340 14.35814
## 520 520 15.20625 15.77319 13.67572
## 521 521 14.68146 15.61250 14.36120
## 522 522 14.34375 14.66914 15.64227
## 523 523 14.72666 14.55211 15.37638
## 524 524 14.79276 15.60384 14.25855
## 525 525 14.99991 14.92935 14.72590
## 526 526 14.91778 14.38789 15.34949
## 527 527 14.74820 14.92388 14.98308
## 528 528 14.77751 14.59820 15.27945
## 529 529 15.43910 13.49662 15.71944
## 530 530 14.04673 15.33111 15.27732
## 531 531 14.90588 15.13786 14.61142
## 532 532 15.35124 14.36626 14.93766
## 533 533 14.64740 14.27818 15.72957
## 534 534 15.34552 14.89638 14.41326
## 535 535 14.52086 15.64157 14.49272
## 536 536 15.07716 15.09197 14.48603
## 537 537 15.34719 14.37159 14.93637
## 538 538 13.55096 15.91529 15.18891
## 539 539 14.68828 15.01971 14.94717
## 540 540 15.39887 14.87182 14.38447
## 541 541 16.00113 14.25320 14.40082
## 542 542 15.56866 15.34474 13.74176
## 543 543 15.51226 14.75783 14.38506
## 544 544 15.45702 13.92252 15.27562
## 545 545 15.25925 13.96941 15.42650
## 546 546 13.71481 15.75208 15.18827
## 547 547 15.37434 14.97175 14.30906
## 548 548 14.39053 15.50388 14.76074
## 549 549 14.60397 15.80679 14.24440
## 550 550 15.09619 14.83567 14.72329
## 551 551 15.20109 15.72891 13.72515
## 552 552 15.37123 14.61144 14.67248
## 553 553 14.81476 13.75913 16.08127
## 554 554 15.07967 15.09122 14.48427
## 555 555 14.88765 14.19272 15.57478
## 556 556 15.50647 14.48484 14.66385
## 557 557 14.62768 15.59899 14.42849
## 558 558 14.49322 14.14562 16.01632
## 559 559 14.24694 15.96682 14.44139
## 560 560 15.22444 15.02077 14.40995
## 561 561 14.75879 15.58678 14.30958
## 562 562 15.52469 14.48673 14.64373
## 563 563 15.07269 13.75236 15.83011
## 564 564 14.25613 15.09900 15.30002
## 565 565 13.99256 15.62888 15.03372
## 566 566 14.36977 15.39742 14.88797
## 567 567 15.02913 15.82408 13.80195
## 568 568 14.32458 15.05681 15.27377
## 569 569 15.34670 15.23403 14.07442
## 570 570 14.76529 15.32150 14.56836
## 571 571 14.76531 13.94938 15.94047
## 572 572 14.15560 15.21643 15.28313
## 573 573 13.92767 15.34324 15.38425
## 574 574 15.62989 15.21258 13.81269
## 575 575 14.37227 16.03542 14.24747
## 576 576 15.69091 14.50831 14.45593
## 577 577 15.48645 14.17838 14.99032
## 578 578 14.81045 14.69004 15.15466
## 579 579 14.59565 15.13075 14.92876
## 580 580 13.97500 15.17837 15.50179
## 581 581 14.56731 15.06683 15.02102
## 582 582 15.57867 14.65368 14.42281
## 583 583 15.61291 15.01966 14.02258
## 584 584 14.42211 15.54914 14.68391
## 585 585 14.26722 15.69506 14.69288
## 586 586 14.94469 15.54692 14.16355
## 587 587 14.65304 14.93826 15.06385
## 588 588 14.86633 14.75244 15.03638
## 589 589 15.03833 14.69138 14.92544
## 590 590 15.10098 15.55501 13.99917
## 591 591 15.98943 14.65297 14.01275
## 592 592 15.10492 14.84649 14.70374
## 593 593 15.74611 14.54340 14.36565
## 594 594 15.17362 15.10119 14.38035
## 595 595 14.81585 15.30037 14.53893
## 596 596 15.52220 14.16865 14.96431
## 597 597 15.29004 14.55008 14.81503
## 598 598 15.08585 14.38587 15.18344
## 599 599 14.45643 14.91153 15.28720
## 600 600 13.86750 15.13788 15.64977
## 601 601 14.63964 14.88006 15.13545
## 602 602 14.97414 14.77011 14.91091
## 603 603 14.62251 14.72310 15.30955
## 604 604 15.47500 15.37265 13.80751
## 605 605 15.02315 14.57392 15.05809
## 606 606 14.53098 15.83944 14.28473
## 607 607 15.73690 14.65392 14.26433
## 608 608 15.42247 14.62696 14.60573
## 609 609 14.87791 14.78208 14.99516
## 610 610 15.45495 14.02647 15.17374
## 611 611 14.05170 15.35000 15.25345
## 612 612 14.30710 15.59493 14.75313
## 613 613 14.93762 15.18721 14.53033
## 614 614 14.44610 14.81838 15.39068
## 615 615 13.62801 15.65673 15.37041
## 616 616 14.39953 14.53529 15.72034
## 617 617 15.53345 14.80968 14.31203
## 618 618 14.13419 15.20710 15.31386
## 619 619 15.63810 14.71299 14.30407
## 620 620 15.54216 14.92752 14.18547
## 621 621 13.69377 15.44581 15.51558
## 622 622 14.59900 13.84009 16.21607
## 623 623 14.96633 14.62811 15.06071
## 624 624 14.91545 14.90573 14.83398
## 625 625 14.87030 14.48986 15.29500
## 626 626 14.45139 14.80029 15.40347
## 627 627 15.99676 14.78279 13.87560
## 628 628 14.24559 15.29222 15.11735
## 629 629 15.36794 14.25326 15.03395
## 630 630 13.64851 16.61588 14.39076
## 631 631 14.24783 15.07756 15.32976
## 632 632 15.91873 14.96169 13.77473
## 633 633 14.36867 15.68578 14.60071
## 634 634 15.99378 14.42679 14.23459
## 635 635 13.80072 15.20083 15.65361
## 636 636 14.66639 14.66848 15.32028
## 637 637 15.63320 13.86578 15.15617
## 638 638 15.14689 14.56741 14.94086
## 639 639 15.28616 14.79813 14.57086
## 640 640 14.59514 15.42701 14.63300
## 641 641 13.78426 15.39656 15.47434
## 642 642 14.87037 14.53940 15.24539
## 643 643 14.58007 15.22413 14.85095
## 644 644 14.84543 14.61195 15.19778
## 645 645 13.87686 14.65003 16.12827
## 646 646 15.23148 14.64586 14.77782
## 647 647 14.53459 15.22593 14.89463
## 648 648 14.30121 15.28121 15.07273
## 649 649 14.42781 14.97140 15.25594
## 650 650 14.89943 15.69024 14.06549
## 651 651 15.06262 14.92171 14.67083
## 652 652 14.74018 15.36652 14.54846
## 653 653 14.39248 14.76073 15.50195
## 654 654 14.81716 15.40910 14.42890
## 655 655 15.01581 15.03014 14.60921
## 656 656 14.77141 14.34691 15.53683
## 657 657 13.96937 15.16186 15.52392
## 658 658 15.26164 14.69255 14.70096
## 659 659 14.31060 14.93244 15.41212
## 660 660 14.65338 15.13647 14.86530
## 661 661 14.82728 14.72353 15.10434
## 662 662 14.95943 14.64908 15.04664
## 663 663 14.64880 15.20127 14.80508
## 664 664 14.86593 15.06926 14.71997
## 665 665 14.59316 15.32126 14.74074
## 666 666 14.31722 15.59096 14.74697
## 667 667 14.69245 15.29883 14.66387
## 668 668 15.01330 14.78987 14.85198
## 669 669 15.40725 14.28075 14.96715
## 670 670 14.37383 15.63883 14.64250
## 671 671 14.37201 14.98689 15.29626
## 672 672 15.35631 13.84425 15.45460
## 673 673 14.54380 15.63434 14.47702
## 674 674 14.69525 14.61066 15.34925
## 675 675 15.43965 14.10166 15.11385
## 676 676 14.29085 15.73958 14.62472
## 677 677 14.33705 14.58435 15.73376
## 678 678 14.55400 14.51297 15.58819
## 679 679 14.14991 15.06839 15.43686
## 680 680 15.38647 14.68772 14.58096
## 681 681 14.45853 14.51937 15.67725
## 682 682 14.54814 15.84195 14.26507
## 683 683 14.99683 15.22291 14.43541
## 684 684 15.17982 13.98496 15.49037
## 685 685 14.54418 14.70982 15.40115
## 686 686 14.97791 15.12150 14.55575
## 687 687 13.54710 15.99561 15.11245
## 688 688 14.26089 15.03190 15.36236
## 689 689 14.96945 15.37358 14.31213
## 690 690 14.19795 14.79468 15.66253
## 691 691 15.06593 14.52387 15.06536
## 692 692 14.25932 15.98622 14.40962
## 693 693 14.62882 15.26936 14.75698
## 694 694 14.55049 15.86121 14.24346
## 695 695 15.17942 15.18890 14.28684
## 696 696 15.43513 14.53953 14.68050
## 697 697 14.78918 15.53396 14.33202
## 698 698 15.15646 14.03430 15.46439
## 699 699 14.50190 14.27422 15.87904
## 700 700 14.27854 15.07383 15.30278
## 701 701 14.71686 14.45029 15.48801
## 702 702 14.47435 15.51456 14.66624
## 703 703 15.46951 14.46539 14.72025
## 704 704 15.37610 14.97256 14.30649
## 705 705 14.84564 14.67166 15.13786
## 706 706 13.84854 14.86494 15.94168
## 707 707 15.34317 15.12132 14.19066
## 708 708 14.13504 14.69496 15.82516
## 709 709 14.07038 15.46465 15.12012
## 710 710 15.21243 15.06822 14.37451
## 711 711 14.62597 14.51756 15.51162
## 712 712 14.93849 14.25597 15.46070
## 713 713 15.55536 14.84468 14.25512
## 714 714 14.45049 15.51186 14.69281
## 715 715 15.71470 13.65976 15.28070
## 716 716 14.51317 15.23290 14.90908
## 717 717 14.52740 15.10185 15.02591
## 718 718 14.55954 15.32806 14.76755
## 719 719 15.26386 13.94296 15.44834
## 720 720 14.73853 14.59252 15.32411
## 721 721 14.40948 15.21310 15.03258
## 722 722 14.82559 14.70572 15.12384
## 723 723 15.03993 14.64439 14.97084
## 724 724 14.86315 14.46224 15.32976
## 725 725 15.21833 14.28081 15.15602
## 726 726 15.26722 14.24413 15.14381
## 727 727 14.86378 14.67133 15.12005
## 728 728 15.14691 14.43154 15.07671
## 729 729 15.89098 13.54043 15.22374
## 730 730 14.26913 15.00885 15.37718
## 731 731 15.36107 13.86537 15.42872
## 732 732 13.82595 15.24085 15.58836
## 733 733 15.17277 14.45473 15.02766
## 734 734 13.90885 15.54700 15.19931
## 735 735 14.23632 14.37114 16.04770
## 736 736 14.34548 15.14250 15.16717
## 737 737 14.76884 15.17070 14.71562
## 738 738 14.58329 16.25862 13.81325
## 739 739 14.32233 15.64446 14.68836
## 740 740 14.53420 15.24514 14.87582
## 741 741 15.61894 13.69541 15.34080
## 742 742 14.52935 14.67608 15.44973
## 743 743 14.16794 15.45813 15.02909
## 744 744 14.97670 14.94249 14.73597
## 745 745 15.40867 13.88851 15.35798
## 746 746 14.97099 15.31914 14.36502
## 747 747 15.11908 14.39911 15.13697
## 748 748 15.07072 14.71989 14.86455
## 749 749 16.51432 14.05347 14.08736
## 750 750 14.15476 15.40924 15.09115
## 751 751 14.63339 15.86741 14.15435
## 752 752 14.50439 15.49551 14.65525
## 753 753 15.29169 14.54908 14.81439
## 754 754 14.73428 14.73956 15.18131
## 755 755 16.05860 14.52790 14.06866
## 756 756 15.06252 15.14322 14.44942
## 757 757 14.66567 15.45303 14.53646
## 758 758 15.66330 15.47143 13.52042
## 759 759 13.40917 15.07935 16.16664
## 760 760 14.83561 15.70506 14.11449
## 761 761 14.73247 15.18887 14.73382
## 762 762 14.34698 15.43824 14.86994
## 763 763 15.30109 14.36804 14.98602
## 764 764 14.74210 15.64009 14.27296
## 765 765 15.70998 13.78856 15.15662
## 766 766 14.67605 15.34224 14.63687
## 767 767 14.85438 15.03094 14.76984
## 768 768 15.12369 14.09215 15.43931
## 769 769 14.91253 13.97771 15.76491
## 770 770 15.53828 14.29356 14.82332
## 771 771 15.50546 15.30275 13.84694
## 772 772 14.19665 15.27753 15.18098
## 773 773 15.12081 15.19803 14.33631
## 774 774 15.60433 14.02850 15.02232
## 775 775 14.63522 15.49144 14.52849
## 776 776 15.35497 13.64276 15.65743
## 777 777 14.83232 14.41826 15.40458
## 778 778 14.82694 14.19831 15.62991
## 779 779 15.54098 15.01593 14.09825
## 780 780 15.96270 14.22852 14.46394
## 781 781 15.21856 14.05838 15.37821
## 782 782 15.32097 15.22727 14.10692
## 783 783 14.89754 14.41695 15.34067
## 784 784 14.55360 15.50180 14.59975
## 785 785 14.72202 15.36540 14.56774
## 786 786 15.41167 15.19842 14.04506
## 787 787 14.90304 14.14386 15.60826
## 788 788 13.91622 15.17927 15.55967
## 789 789 14.41126 15.87325 14.37064
## 790 790 14.81736 14.56697 15.27083
## 791 791 15.02216 15.25369 14.37931
## 792 792 15.21637 14.91981 14.51897
## 793 793 13.38453 15.60841 15.66222
## 794 794 15.43324 14.36109 14.86083
## 795 795 14.43183 15.77462 14.44870
## 796 796 14.40464 14.86883 15.38168
## 797 797 14.80464 14.59927 15.25124
## 798 798 14.66213 14.57226 15.42077
## 799 799 14.54551 14.94788 15.16176
## 800 800 14.28542 14.82244 15.54729
## 801 801 14.09203 16.40480 14.15832
## 802 802 14.72511 15.13794 14.79211
## 803 803 15.28392 14.75544 14.61579
## 804 804 14.14671 15.45372 15.05473
## 805 805 14.32425 14.61404 15.71687
## 806 806 15.03491 15.63179 13.98845
## 807 807 14.84108 15.14863 14.66544
## 808 808 15.42277 13.96181 15.27058
## 809 809 15.31286 14.38750 14.95479
## 810 810 15.49362 14.19061 14.97093
## 811 811 14.86769 14.91920 14.86826
## 812 812 15.47144 14.44350 14.74021
## 813 813 14.24974 15.06096 15.34446
## 814 814 15.04797 15.17485 14.43234
## 815 815 14.73165 14.32088 15.60262
## 816 816 16.48979 14.31643 13.84893
## 817 817 15.11991 14.70361 14.83163
## 818 818 15.70613 14.06816 14.88087
## 819 819 15.08592 14.46645 15.10278
## 820 820 14.39139 15.58582 14.67795
## 821 821 14.69416 15.04231 14.91869
## 822 822 14.75655 15.31957 14.57903
## 823 823 14.40859 14.81968 15.42689
## 824 824 15.51150 14.08631 15.05734
## 825 825 14.58465 15.37325 14.69725
## 826 826 15.16009 14.43661 15.05846
## 827 827 15.33682 14.09442 15.22392
## 828 828 14.89282 13.89500 15.86734
## 829 829 14.53528 15.46081 14.65906
## 830 830 15.07928 14.78887 14.78701
## 831 831 15.04849 15.00835 14.59832
## 832 832 15.08706 14.78770 14.78039
## 833 833 14.37449 14.45882 15.82185
## 834 834 14.85490 15.02471 14.77554
## 835 835 15.55515 13.77253 15.32748
## 836 836 14.85650 15.85414 13.94451
## 837 837 15.00038 14.42210 15.23267
## 838 838 15.12256 14.23548 15.29711
## 839 839 14.58180 15.47670 14.59666
## 840 840 14.52714 14.27011 15.85791
## 841 841 14.39634 15.17036 15.08845
## 842 842 14.56146 15.01171 15.08199
## 843 843 15.08304 15.18862 14.38349
## 844 844 15.08891 14.36683 15.19942
## 845 845 14.79734 14.68416 15.17365
## 846 846 15.24610 14.90969 14.49937
## 847 847 16.41600 14.33517 13.90399
## 848 848 14.84999 14.53340 15.27177
## 849 849 14.77896 14.85716 15.01904
## 850 850 15.09682 15.27727 14.28107
## 851 851 15.47145 14.40281 14.78090
## 852 852 14.84667 15.40025 14.40824
## 853 853 14.28533 15.09530 15.27453
## 854 854 14.52959 14.40539 15.72018
## 855 855 14.94056 15.08507 14.62952
## 856 856 13.82435 15.50685 15.32396
## 857 857 15.07402 14.35913 15.22200
## 858 858 14.48670 15.10154 15.06691
## 859 859 15.06548 14.63245 14.95723
## 860 860 14.41940 15.34909 14.88666
## 861 861 14.78217 15.55958 14.31341
## 862 862 14.28914 14.81968 15.54633
## 863 863 14.16993 15.56888 14.91634
## 864 864 14.70457 14.28464 15.66595
## 865 865 14.11722 15.29892 15.23901
## 866 866 14.85301 14.69090 15.11125
## 867 867 14.66715 15.47742 14.51058
## 868 868 15.15018 14.54202 14.96296
## 869 869 13.87841 15.41686 15.35989
## 870 870 14.60576 15.20477 14.84462
## 871 871 14.58963 15.08089 14.98464
## 872 872 14.64229 14.76306 15.24981
## 873 873 15.93766 14.27830 14.43920
## 874 874 15.45725 15.24174 13.95616
## 875 875 14.48414 14.94411 15.22691
## 876 876 15.41989 14.19358 15.04169
## 877 877 15.27198 14.89560 14.48758
## 878 878 14.35528 14.41162 15.88825
## 879 879 15.07360 14.54156 15.03999
## 880 880 14.45531 14.67080 15.52905
## 881 881 14.43893 14.59191 15.62432
## 882 882 15.47561 13.38273 15.79682
## 883 883 15.09683 14.85313 14.70520
## 884 884 14.90432 15.05799 14.69285
## 885 885 15.83668 14.21899 14.59949
## 886 886 14.53942 15.23960 14.87613
## 887 887 14.92693 14.69446 15.03376
## 888 888 14.78380 15.77631 14.09505
## 889 889 14.59562 14.40242 15.65712
## 890 890 14.41795 14.82418 15.41302
## 891 891 14.49369 15.17474 14.98673
## 892 892 14.96200 15.04039 14.65277
## 893 893 15.57525 14.06807 15.01183
## 894 894 15.46208 14.35629 14.83678
## 895 895 14.71988 15.34901 14.58626
## 896 896 15.14992 15.46924 14.03600
## 897 897 14.25021 16.68531 13.71964
## 898 898 15.59876 13.93521 15.12119
## 899 899 14.66070 15.33574 14.65872
## 900 900 15.27007 15.01461 14.37047
## 901 901 14.69146 14.37381 15.58989
## 902 902 15.09278 14.55116 15.01121
## 903 903 14.78352 15.39263 14.47900
## 904 904 15.48042 14.15937 15.01536
## 905 905 14.41790 14.42907 15.80818
## 906 906 14.84226 15.59582 14.21707
## 907 907 14.63461 15.00737 15.01318
## 908 908 14.29803 15.12725 15.22987
## 909 909 15.18353 15.28208 14.18954
## 910 910 15.38710 14.41651 14.85155
## 911 911 14.29164 14.86074 15.50277
## 912 912 14.45035 16.06432 14.14049
## 913 913 16.14389 14.15804 14.35323
## 914 914 15.21987 14.48827 14.94702
## 915 915 15.22444 14.64281 14.78791
## 916 916 15.38843 14.90264 14.36409
## 917 917 15.52968 13.79116 15.33431
## 918 918 15.49539 15.13029 14.02948
## 919 919 15.16541 14.08585 15.40390
## 920 920 15.36954 14.16131 15.12431
## 921 921 14.35781 14.83818 15.45916
## 922 922 15.28554 14.98351 14.38610
## 923 923 14.84114 14.80474 15.00928
## 924 924 15.78336 14.61786 14.25394
## 925 925 14.90823 15.57434 14.17258
## 926 926 14.77821 15.10249 14.77445
## 927 927 14.92792 14.28204 15.44520
## 928 928 15.34403 14.11637 15.19476
## 929 929 15.31560 14.09674 15.24282
## 930 930 16.04352 14.17156 14.44007
## 931 931 15.01692 14.07694 15.56130
## 932 932 15.27463 15.43662 13.94391
## 933 933 14.45530 15.35193 14.84793
## 934 934 14.93089 15.68534 14.03893
## 935 935 14.85710 15.10087 14.69719
## 936 936 14.44988 15.15535 15.04993
## 937 937 15.08801 15.27545 14.29169
## 938 938 14.60167 14.76959 15.28390
## 939 939 14.49104 15.06609 15.09803
## 940 940 13.98452 14.73217 15.93847
## 941 941 15.20472 14.64580 14.80463
## 942 942 15.50968 14.81373 14.33175
## 943 943 14.78897 16.41872 13.44747
## 944 944 14.48202 15.31693 14.85621
## 945 945 15.59043 14.41038 14.65435
## 946 946 13.79162 14.46890 16.39463
## 947 947 14.72150 15.84658 14.08707
## 948 948 15.78897 15.04325 13.82293
## 949 949 14.92287 15.19864 14.53364
## 950 950 16.17611 14.08375 14.39529
## 951 951 16.67485 14.19851 13.78180
## 952 952 14.90976 15.90029 13.84510
## 953 953 15.69236 14.82319 14.13960
## 954 954 15.11742 14.90847 14.62927
## 955 955 13.91093 15.36934 15.37489
## 956 956 14.62569 15.96960 14.05986
## 957 957 14.77177 15.28121 14.60218
## 958 958 14.65862 15.53241 14.46413
## 959 959 14.78038 15.04280 14.83197
## 960 960 14.66802 14.72134 15.26579
## 961 961 14.95084 15.17902 14.52529
## 962 962 14.78780 15.47134 14.39602
## 963 963 15.05512 15.27124 14.32879
## 964 964 14.40222 15.15079 15.10215
## 965 965 15.06599 14.90122 14.68794
## 966 966 15.26088 15.08221 14.31207
## 967 967 14.83044 15.01217 14.81254
## 968 968 15.35461 14.98916 14.31139
## 969 969 14.76996 15.60773 14.27747
## 970 970 13.47719 16.10270 15.07526
## 971 971 14.37157 14.99318 15.29040
## 972 972 14.17022 15.25547 15.22946
## 973 973 14.70670 14.43499 15.51346
## 974 974 15.32331 13.66435 15.66750
## 975 975 14.05951 14.77046 15.82519
## 976 976 13.86667 14.70480 16.08368
## 977 977 14.68106 15.63342 14.34067
## 978 978 14.57759 15.36327 14.71429
## 979 979 14.59671 15.28147 14.77697
## 980 980 14.59467 15.47968 14.58081
## 981 981 14.07588 15.57832 15.00096
## 982 982 15.36666 15.25251 14.03598
## 983 983 15.07122 14.88426 14.69968
## 984 984 15.28447 14.10094 15.26975
## 985 985 15.33808 14.27429 15.04279
## 986 986 15.23040 15.66553 13.75922
## 987 987 15.19895 15.10161 14.35459
## 988 988 15.84521 15.02219 13.78776
## 989 989 15.28409 14.73329 14.63778
## 990 990 14.79527 15.13676 14.72313
## 991 991 15.19079 14.10066 15.36371
## 992 992 14.77586 14.69211 15.18718
## 993 993 14.77490 14.46901 15.41125
## 994 994 14.97101 15.19966 14.48448
## 995 995 14.54174 14.34128 15.77214
## 996 996 15.39353 14.26670 14.99493
## 997 997 15.41949 13.93582 15.29985
## 998 998 14.66559 14.87694 15.11263
## 999 999 15.82584 14.78991 14.03940
## 1000 1000 15.51446 15.28130 13.859394d.
library(ggplot2)
library(ggpubr)
library(reshape2)##
## Attaching package: 'reshape2'## The following object is masked from 'package:tidyr':
##
## smithsp1 <- ggplot(data=datfram)+
aes(x=dat1)+
geom_histogram(aes(fill=I("mediumorchid"),
color=I("black")))
p2 <- ggplot(data=datfram)+
aes(x=dat2)+
geom_histogram(aes(fill=I("pink"),
color=I("black")))
p3 <- ggplot(data=datfram)+
aes(x=dat3)+
geom_histogram(aes(fill=I("olivedrab4"),
color=I("black")))
graphx3 <- ggarrange(p1,p2,p3, nrow=1)## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.print(graphx3)
# the reshuffled means are more similar than the original unshuffled means.Add slope or pvalue histogram at the end
####### ###########################################
# function: read_data
# read in (or generate) data set for analysis
# input: file name (or nothing, for this demo)
# output: 3 column data frame of observed data (x,y)
#-------------------------------------------------
read_data <- function(z=NULL) {
if(is.null(z)){
x_obs <- 1:20
y_obs <- x_obs + 10*rnorm(20)
df <- data.frame(x_obs,
y_obs)} else {
# set up data frame
df <-read.table(file=z,
header=TRUE,
sep=",")}
return(df)
}
get_metric <- function(z=NULL) {
if(is.null(z)){
x_obs <- 1:20
y_obs <- x_obs + 10*rnorm(20)
z <- data.frame(x_obs,y_obs)} # set up data frame
. <- lm(z[,2]~z[,1])
. <- summary(.)
. <- .$coefficients[2,1]
slope <- .
return(slope)
}
##################################################
# function: shuffle_data
# randomize data for regression analysis
# input: 2-column data frame (x_var,y_var)
# output: 2-column data frame (x_var,y_var)
#-------------------------------------------------
shuffle_data <- function(z=NULL) {
if(is.null(z)){
x_obs <- 1:20
y_obs <- x_obs + 3*rnorm(20)
z <- data.frame(x_obs,y_obs)} # set up data frame
z[,2] <- sample(z[,2]) # use sample function with defaults to reshuffle column
return(z)
}
##################################################
# function: get_pval
# calculate p value from simulation
# input: list of observed metric, and vector of simulated metrics
# output: lower, upper tail probability values
#-------------------------------------------------
get_pval <- function(z=NULL) {
if(is.null(z)){
z <- list(x_obs=runif(1),x_sim=runif(1000))}
p_lower <- mean(z[[2]]<=z[[1]])
p_upper <- mean(z[[2]]>=z[[1]])
return(c(pl=p_lower,pu=p_upper))
}
get_pval()## pl pu
## 0.131 0.869#summary(aov(Response~DataSet, data=wholedat1))
##################################################
# function: plot_ran_test
# create ggplot of histogram of simulated values
# input: list of observed metric and vector of simulated metrics
# output: saved ggplot graph
#-------------------------------------------------
plot_ran_test <- function(z=NULL) {
if(is.null(z)){
z <- list(rnorm(1),rnorm(1000)) }
df <- data.frame(ID=seq_along(z[[2]]),sim_x=z[[2]])
p1 <- ggplot(data=df) +
aes(x=sim_x)
p1 + geom_histogram(aes(fill=I("olivedrab4"),
color=I("black"))) +
geom_vline(aes(xintercept=z[[1]],
col="blue"))
}
plot_ran_test1 <- function(z=NULL) {
if(is.null(z)){
z <- list(rnorm(1),rnorm(1000)) }
df <- data.frame(ID=seq_along(z[[2]]),sim_x=z[[2]])
p1 <- ggplot(data=df) +
aes(x=sim_x)
p1 + geom_histogram(aes(fill=I("olivedrab4"),
color=I("black")))
}
n_sim <- 1000
x_sim <- rep(NA,n_sim) # vector of simulated slopes
df <- read_data()
x_obs <- get_metric(df)
for (i in seq_len(n_sim)) {
x_sim[i] <- get_metric(shuffle_data(df))
}
slopes <- list(x_obs,x_sim)
get_pval(slopes)## pl pu
## 0.964 0.036get_pval(datfram)## pl pu
## 0.986 0.014plot_ran_test(slopes)## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
plot_ran_test()## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
plot_ran_test1(datfram)## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.