WGN Decile Analysis
noah castelo
November 19, 2017
Eric’s new analysis
setwd("~/Documents/Dropbox/Research/Eric")
combined2<-read.csv("old and new data.csv", header=T, sep=",")
combined2$bd<-combined2$beta*((1/(1+combined2$delta))^(2/52))*1.5
combined2$bd<-combined2$bd-mean(combined2$bd)
library(binr)
table(combined2$cond)##
## LL no SS
## 388 302 398LL<-subset(combined2, cond=="LL")
no<-subset(combined2, cond=="no")
SS<-subset(combined2, cond=="SS")
#this tells you how to split the data into equal groups of 10 (can replace 10 with whatever number you like)
bins(LL$bd, 10, max.breaks=10)## $binlo
## [1] 1 40 79 118 156 195 234 272 311 350
##
## $binhi
## [1] 39 78 117 155 194 233 271 310 349 388
##
## $binct
## [-1.11227845662551, -0.26547540127159]
## 39
## [-0.264652198896631, -0.0974938409475334]
## 39
## [-0.097363873962121, -0.0584618640757131]
## 39
## [-0.0577745058735397, -0.0198004518905934]
## 38
## [-0.0195612540987744, 0.0669931696446873]
## 39
## [0.067371217911983, 0.0793009744058693]
## 39
## [0.0807925559701905, 0.0891775871803935]
## 38
## [0.0896905601442171, 0.132743672691635]
## 39
## [0.132927718649512, 0.141290070916923]
## 39
## [0.1487291411856, 0.534114524707168]
## 39
##
## $xtbl
## x
## -1.11227845662551 -1.0741272439068 -0.980436176506176
## 1 1 1
## -0.955453742004612 -0.828228226663682 -0.794353573009068
## 1 1 1
## -0.759574641950244 -0.724161458728211 -0.718799414709572
## 1 1 1
## -0.709913548874529 -0.681537839797338 -0.675236830879661
## 1 1 1
## -0.652553685418721 -0.602596911843903 -0.588778611147252
## 1 1 1
## -0.586863223589377 -0.580965897418416 -0.578479307073299
## 1 1 1
## -0.576856016266728 -0.55502816998055 -0.535505445245196
## 1 1 1
## -0.532718519129385 -0.489990690788054 -0.478429995570203
## 1 1 1
## -0.47718703129689 -0.476729446167565 -0.459894591747333
## 1 1 1
## -0.43728150536141 -0.42991816974438 -0.32338199869776
## 1 1 1
## -0.319730414244534 -0.316855935392764 -0.308731539529093
## 1 1 1
## -0.285685505085471 -0.284337382108776 -0.278300959647696
## 1 1 1
## -0.275531791048855 -0.274919339781037 -0.26547540127159
## 1 1 1
## -0.264652198896631 -0.263739851283789 -0.262179084825662
## 1 1 1
## -0.260252248051553 -0.253445102188586 -0.233329342119634
## 1 1 1
## -0.223237303091576 -0.221758775405746 -0.219524651342838
## 1 1 1
## -0.214525518408301 -0.212038682137928 -0.209912514512976
## 1 1 1
## -0.204762250791223 -0.187065678943863 -0.186978471994227
## 1 1 1
## -0.185023112495876 -0.184565443989401 -0.183018911199105
## 1 1 1
## -0.174493514157916 -0.174311867825896 -0.173310723623601
## 1 1 1
## -0.171868054759016 -0.170666914553193 -0.170340507377603
## 1 1 1
## -0.164112379657869 -0.163064516911348 -0.153056184767794
## 1 1 1
## -0.150027169253655 -0.139042325618539 -0.126292549952758
## 1 1 1
## -0.120058438591758 -0.115708377568568 -0.11400644867216
## 1 1 1
## -0.113982341074229 -0.106793428833768 -0.105392949692959
## 1 1 1
## -0.0996284638052505 -0.0984357352922178 -0.0974938409475334
## 1 1 1
## -0.097363873962121 -0.0967612174875339 -0.096741945204156
## 1 1 1
## -0.0959424892698171 -0.095850973910258 -0.0956168218275868
## 1 1 1
## -0.0921890113859387 -0.0876606034113034 -0.0873414695886401
## 1 1 1
## -0.0852766371159834 -0.0849736769146685 -0.0849636867337802
## 1 1 1
## -0.0837179308995111 -0.0826252995773322 -0.0824005654632489
## 1 1 1
## -0.0820231405476624 -0.0817051103642554 -0.0807037242445743
## 1 1 1
## -0.0804128755199742 -0.0798794967367553 -0.079127231687286
## 1 1 1
## -0.0788466380551367 -0.0785187441845145 -0.0784143868285139
## 1 1 1
## -0.0782173448063883 -0.0780560506069177 -0.077863568997582
## 1 1 1
## -0.0771139848781961 -0.0769328081113276 -0.0767047041329536
## 1 1 1
## -0.0763142631091502 -0.0756049941463259 -0.0750881012251634
## 1 1 1
## -0.0740396432801091 -0.0739622103153135 -0.0731815196753469
## 1 1 1
## -0.0730649775103875 -0.0729889789549578 -0.0584618640757131
## 1 1 1
## -0.0577745058735397 -0.0576519889512255 -0.0569084469849892
## 1 1 1
## -0.056564608716573 -0.0473421149744593 -0.0454837743626881
## 1 1 1
## -0.0454691106577747 -0.0436050478991907 -0.0423914502732696
## 1 1 1
## -0.0409968410860539 -0.0408670372720663 -0.0408103085136282
## 1 1 1
## -0.0397748001683915 -0.0396100771723336 -0.0390295345139715
## 1 1 1
## -0.038791940565726 -0.03770118720717 -0.0373861082776978
## 1 1 1
## -0.0373114880259915 -0.0371751511903775 -0.0370930478417251
## 1 1 1
## -0.0370456676002411 -0.0364874631460022 -0.0359572965764885
## 1 1 1
## -0.0356130654584776 -0.0347000055805267 -0.0337086828626847
## 1 1 1
## -0.0329396699574125 -0.0324761111093612 -0.0318578271241441
## 1 1 1
## -0.0317228346649263 -0.031349663490235 -0.0295781443270395
## 1 1 1
## -0.0295517643786043 -0.0290673523279834 -0.0271743967547238
## 1 1 1
## -0.020603600673917 -0.0198004518905934 -0.0195612540987744
## 1 1 1
## -0.0195068144289756 -0.0194236428112753 -0.0190938958174862
## 1 1 1
## -0.0168489655452562 -0.0167361298938484 -0.0148243689889869
## 1 1 1
## -0.0139613738348381 -0.0131723958691725 -0.0129993653674716
## 1 1 1
## -0.0119177674336883 -0.00842287723526725 0.0258272703450102
## 1 1 1
## 0.0259711762882917 0.0268668666866483 0.0281227223564451
## 1 1 1
## 0.0296321944777436 0.0384639604018764 0.0421986935063774
## 1 1 1
## 0.0451959678754501 0.0540982454599974 0.0594286911616404
## 1 1 1
## 0.0632778740606375 0.0635266769154059 0.064567757548001
## 1 1 1
## 0.0650521587807291 0.0651123071451956 0.0654878309145071
## 1 1 1
## 0.0657686600714864 0.0658146589384658 0.0659779425015103
## 1 1 1
## 0.0663157456155097 0.0666468055075535 0.0668077358648169
## 1 1 1
## 0.0668479525809289 0.0668686585892369 0.0669131334208037
## 1 1 1
## 0.0669213885356661 0.0669931696446873 0.067371217911983
## 1 1 1
## 0.0674659489970559 0.0676335722771342 0.0676443148183914
## 1 1 1
## 0.0677680874138917 0.0681895057477693 0.0683616885173284
## 1 1 1
## 0.0684440359522043 0.0685400558344238 0.0685918446967273
## 1 1 1
## 0.0686175716640278 0.0686837325135272 0.0687775503764354
## 1 1 1
## 0.0688159235153498 0.0689413088795603 0.0690589827669914
## 1 1 1
## 0.0693084877599326 0.0693565090990957 0.0694215310401398
## 1 1 1
## 0.0694633221909451 0.0697900898925563 0.0698200007522987
## 1 1 1
## 0.0700572271401647 0.0701213510354042 0.0701618342742689
## 1 1 1
## 0.0702021060152371 0.0703863059649144 0.0704277180706674
## 1 1 1
## 0.070626315558949 0.0712231902752802 0.0714436503916964
## 1 1 1
## 0.0715509327059016 0.071760948310275 0.0720090706116472
## 1 1 1
## 0.072321225516248 0.0727900766157545 0.0731000627880065
## 1 1 1
## 0.0788930377374908 0.0793009744058693 0.0807925559701905
## 1 1 1
## 0.0818187894888678 0.0821455729300287 0.0822693229642033
## 1 1 1
## 0.0823956585091068 0.0824153824003777 0.0824368405490712
## 1 1 1
## 0.082660184794304 0.0827787074363662 0.0829260083554191
## 1 1 1
## 0.0832966291639352 0.0833015732820335 0.0834132696203587
## 1 1 1
## 0.0838508036014489 0.0840499781549884 0.0841932651812725
## 1 1 1
## 0.0844545467867492 0.0847915253346301 0.0851304992696647
## 1 1 1
## 0.0852173808391568 0.0857039901060639 0.0859375183805167
## 1 1 1
## 0.086365521775108 0.0864456329470189 0.0864667611090684
## 1 1 1
## 0.0864690317893759 0.0865286448553595 0.08676986830447
## 1 1 1
## 0.0870053531134736 0.0875053603854545 0.087759837241788
## 1 1 1
## 0.0878758877895383 0.0879235104973892 0.0881653704557777
## 1 1 1
## 0.0882089237088208 0.0883067536443791 0.0889916491431686
## 1 1 1
## 0.0891775871803935 0.0896905601442171 0.0903792600316073
## 1 1 1
## 0.0923460835101215 0.0935116099345767 0.106585480535538
## 1 1 1
## 0.107257312804952 0.107788232981788 0.107852973558389
## 1 1 1
## 0.108441991678047 0.10913637077884 0.109290774396495
## 1 1 1
## 0.110430029348209 0.110501844469719 0.110561445787604
## 1 1 1
## 0.111386240753809 0.111568995261981 0.111872506185578
## 1 1 1
## 0.112082422620578 0.112086353674442 0.112391408786168
## 1 1 1
## 0.112484212520253 0.113375509140526 0.113413347647526
## 1 1 1
## 0.113776105407962 0.114063145233244 0.114959390196214
## 1 1 1
## 0.116084315174759 0.116323458294273 0.116599962462022
## 1 1 1
## 0.117432513930646 0.129176042676793 0.129352897312813
## 1 1 1
## 0.129601408385247 0.130115331972626 0.13153499448079
## 1 1 1
## 0.13182373428516 0.132641345987397 0.132740945707321
## 1 1 1
## 0.132743672691635 0.132927718649512 0.133100219084926
## 1 1 1
## 0.133546338491266 0.13354994560267 0.133570926637283
## 1 1 1
## 0.133709112453237 0.134055922138938 0.134307083568273
## 1 1 1
## 0.134372050527348 0.134389042487766 0.134467940835287
## 1 1 1
## 0.134785667818912 0.134854490393745 0.134874684586914
## 1 1 1
## 0.135004396765791 0.13501037974824 0.135045169511866
## 1 1 1
## 0.135793397035989 0.13581786948993 0.135890358759999
## 1 1 1
## 0.136054290761541 0.136118870031577 0.136136100013238
## 1 1 1
## 0.136505504962887 0.136574772228083 0.136649939152269
## 1 1 1
## 0.137143517901988 0.137558796325784 0.13784508768546
## 1 1 1
## 0.137911465830862 0.137959451090259 0.137998557376306
## 1 1 1
## 0.138349794140666 0.138825794129616 0.139049480332712
## 1 1 1
## 0.139596567849236 0.139931375103231 0.139938824757353
## 1 1 1
## 0.141290070916923 0.1487291411856 0.150437057747395
## 1 1 1
## 0.151788389236001 0.152242154263942 0.176463659944973
## 1 1 1
## 0.189472076472383 0.189878026263939 0.191075947503616
## 1 1 1
## 0.192727652039144 0.193678880451525 0.202734915417667
## 1 1 1
## 0.20437463115617 0.205566348991616 0.206369341620749
## 1 1 1
## 0.209509146339571 0.210442817622622 0.211404736111103
## 1 1 1
## 0.211422828512046 0.212132616485895 0.212208658946841
## 1 1 1
## 0.212451580519801 0.212536196840912 0.227842112265946
## 1 1 1
## 0.229080215825358 0.23044994116796 0.23197707320252
## 1 1 1
## 0.291460664336377 0.402750071307765 0.446186680133073
## 1 1 1
## 0.474858092726109 0.482001366054306 0.482423851718019
## 1 1 1
## 0.482851502762573 0.483616490812146 0.531649905824617
## 1 1 1
## 0.533074014854971 0.533469660949924 0.533872388650733
## 1 1 1
## 0.534114524707168
## 1
##
## $xval
## [1] -1.112278457 -1.074127244 -0.980436177 -0.955453742 -0.828228227
## [6] -0.794353573 -0.759574642 -0.724161459 -0.718799415 -0.709913549
## [11] -0.681537840 -0.675236831 -0.652553685 -0.602596912 -0.588778611
## [16] -0.586863224 -0.580965897 -0.578479307 -0.576856016 -0.555028170
## [21] -0.535505445 -0.532718519 -0.489990691 -0.478429996 -0.477187031
## [26] -0.476729446 -0.459894592 -0.437281505 -0.429918170 -0.323381999
## [31] -0.319730414 -0.316855935 -0.308731540 -0.285685505 -0.284337382
## [36] -0.278300960 -0.275531791 -0.274919340 -0.265475401 -0.264652199
## [41] -0.263739851 -0.262179085 -0.260252248 -0.253445102 -0.233329342
## [46] -0.223237303 -0.221758775 -0.219524651 -0.214525518 -0.212038682
## [51] -0.209912515 -0.204762251 -0.187065679 -0.186978472 -0.185023112
## [56] -0.184565444 -0.183018911 -0.174493514 -0.174311868 -0.173310724
## [61] -0.171868055 -0.170666915 -0.170340507 -0.164112380 -0.163064517
## [66] -0.153056185 -0.150027169 -0.139042326 -0.126292550 -0.120058439
## [71] -0.115708378 -0.114006449 -0.113982341 -0.106793429 -0.105392950
## [76] -0.099628464 -0.098435735 -0.097493841 -0.097363874 -0.096761217
## [81] -0.096741945 -0.095942489 -0.095850974 -0.095616822 -0.092189011
## [86] -0.087660603 -0.087341470 -0.085276637 -0.084973677 -0.084963687
## [91] -0.083717931 -0.082625300 -0.082400565 -0.082023141 -0.081705110
## [96] -0.080703724 -0.080412876 -0.079879497 -0.079127232 -0.078846638
## [101] -0.078518744 -0.078414387 -0.078217345 -0.078056051 -0.077863569
## [106] -0.077113985 -0.076932808 -0.076704704 -0.076314263 -0.075604994
## [111] -0.075088101 -0.074039643 -0.073962210 -0.073181520 -0.073064978
## [116] -0.072988979 -0.058461864 -0.057774506 -0.057651989 -0.056908447
## [121] -0.056564609 -0.047342115 -0.045483774 -0.045469111 -0.043605048
## [126] -0.042391450 -0.040996841 -0.040867037 -0.040810309 -0.039774800
## [131] -0.039610077 -0.039029535 -0.038791941 -0.037701187 -0.037386108
## [136] -0.037311488 -0.037175151 -0.037093048 -0.037045668 -0.036487463
## [141] -0.035957297 -0.035613065 -0.034700006 -0.033708683 -0.032939670
## [146] -0.032476111 -0.031857827 -0.031722835 -0.031349663 -0.029578144
## [151] -0.029551764 -0.029067352 -0.027174397 -0.020603601 -0.019800452
## [156] -0.019561254 -0.019506814 -0.019423643 -0.019093896 -0.016848966
## [161] -0.016736130 -0.014824369 -0.013961374 -0.013172396 -0.012999365
## [166] -0.011917767 -0.008422877 0.025827270 0.025971176 0.026866867
## [171] 0.028122722 0.029632194 0.038463960 0.042198694 0.045195968
## [176] 0.054098245 0.059428691 0.063277874 0.063526677 0.064567758
## [181] 0.065052159 0.065112307 0.065487831 0.065768660 0.065814659
## [186] 0.065977943 0.066315746 0.066646806 0.066807736 0.066847953
## [191] 0.066868659 0.066913133 0.066921389 0.066993170 0.067371218
## [196] 0.067465949 0.067633572 0.067644315 0.067768087 0.068189506
## [201] 0.068361689 0.068444036 0.068540056 0.068591845 0.068617572
## [206] 0.068683733 0.068777550 0.068815924 0.068941309 0.069058983
## [211] 0.069308488 0.069356509 0.069421531 0.069463322 0.069790090
## [216] 0.069820001 0.070057227 0.070121351 0.070161834 0.070202106
## [221] 0.070386306 0.070427718 0.070626316 0.071223190 0.071443650
## [226] 0.071550933 0.071760948 0.072009071 0.072321226 0.072790077
## [231] 0.073100063 0.078893038 0.079300974 0.080792556 0.081818789
## [236] 0.082145573 0.082269323 0.082395659 0.082415382 0.082436841
## [241] 0.082660185 0.082778707 0.082926008 0.083296629 0.083301573
## [246] 0.083413270 0.083850804 0.084049978 0.084193265 0.084454547
## [251] 0.084791525 0.085130499 0.085217381 0.085703990 0.085937518
## [256] 0.086365522 0.086445633 0.086466761 0.086469032 0.086528645
## [261] 0.086769868 0.087005353 0.087505360 0.087759837 0.087875888
## [266] 0.087923510 0.088165370 0.088208924 0.088306754 0.088991649
## [271] 0.089177587 0.089690560 0.090379260 0.092346084 0.093511610
## [276] 0.106585481 0.107257313 0.107788233 0.107852974 0.108441992
## [281] 0.109136371 0.109290774 0.110430029 0.110501844 0.110561446
## [286] 0.111386241 0.111568995 0.111872506 0.112082423 0.112086354
## [291] 0.112391409 0.112484213 0.113375509 0.113413348 0.113776105
## [296] 0.114063145 0.114959390 0.116084315 0.116323458 0.116599962
## [301] 0.117432514 0.129176043 0.129352897 0.129601408 0.130115332
## [306] 0.131534994 0.131823734 0.132641346 0.132740946 0.132743673
## [311] 0.132927719 0.133100219 0.133546338 0.133549946 0.133570927
## [316] 0.133709112 0.134055922 0.134307084 0.134372051 0.134389042
## [321] 0.134467941 0.134785668 0.134854490 0.134874685 0.135004397
## [326] 0.135010380 0.135045170 0.135793397 0.135817869 0.135890359
## [331] 0.136054291 0.136118870 0.136136100 0.136505505 0.136574772
## [336] 0.136649939 0.137143518 0.137558796 0.137845088 0.137911466
## [341] 0.137959451 0.137998557 0.138349794 0.138825794 0.139049480
## [346] 0.139596568 0.139931375 0.139938825 0.141290071 0.148729141
## [351] 0.150437058 0.151788389 0.152242154 0.176463660 0.189472076
## [356] 0.189878026 0.191075948 0.192727652 0.193678880 0.202734915
## [361] 0.204374631 0.205566349 0.206369342 0.209509146 0.210442818
## [366] 0.211404736 0.211422829 0.212132616 0.212208659 0.212451581
## [371] 0.212536197 0.227842112 0.229080216 0.230449941 0.231977073
## [376] 0.291460664 0.402750071 0.446186680 0.474858093 0.482001366
## [381] 0.482423852 0.482851503 0.483616491 0.531649906 0.533074015
## [386] 0.533469661 0.533872389 0.534114525
##
## $err
## [1] 0.4#then I actually split the data according to the previous output
LL1<-subset(LL, bd <=-0.26547540127159)
LL2<-subset(LL, bd > -0.26547540127159 & bd <= -0.0974938409475334)
LL3<-subset(LL, bd >-0.0974938409475334 & bd <=-0.0584618640757131)
LL4<-subset(LL, bd >-0.0584618640757131 & bd <=-0.0198004518905934)
LL5<-subset(LL, bd >-0.0198004518905934 & bd <=0.0669931696446873)
LL6<-subset(LL, bd >0.0669931696446873 & bd <=0.0793009744058693)
LL7<-subset(LL, bd >0.0793009744058693 & bd <=0.0891775871803935)
LL8<-subset(LL, bd >0.0891775871803935 & bd <=0.132743672691635)
LL9<-subset(LL, bd >0.132743672691635 & bd <=0.141290070916923)
LL10<-subset(LL, bd >0.141290070916923 & bd <=0.534114524707168)
#preparing to recombine data into 1 set
LL1$decile<-1
LL2$decile<-2
LL3$decile<-3
LL4$decile<-4
LL5$decile<-5
LL6$decile<-6
LL7$decile<-7
LL8$decile<-8
LL9$decile<-9
LL10$decile<-10
#selecting only the variables we need
LL1<-subset(LL1, select= c(23, 24, 28:29))
LL2<-subset(LL2, select= c(23, 24, 28:29))
LL3<-subset(LL3, select= c(23, 24, 28:29))
LL4<-subset(LL4, select= c(23, 24, 28:29))
LL5<-subset(LL5, select= c(23, 24, 28:29))
LL6<-subset(LL6, select= c(23, 24, 28:29))
LL7<-subset(LL7, select= c(23, 24, 28:29))
LL8<-subset(LL8, select= c(23, 24, 28:29))
LL9<-subset(LL9, select= c(23, 24, 28:29))
LL10<-subset(LL10, select= c(23, 24, 28:29))
#recombining data
LLcond<-rbind(LL1, LL2, LL3, LL4, LL5, LL6, LL7, LL8, LL9, LL10)
head(LLcond)## cond SS bd decile
## 11 LL 0 -0.5784793 1
## 12 LL 1 -0.2783010 1
## 14 LL 0 -0.4372815 1
## 40 LL 0 -0.3087315 1
## 48 LL 1 -0.7099135 1
## 70 LL 1 -0.4299182 1LLcond<-subset(LLcond, select=c(1,2,4))
#doing the same thing for SS and no default conditions
bins(SS$bd, 10, max.breaks=10)## $binlo
## [1] 1 41 81 121 160 200 240 279 319 359
##
## $binhi
## [1] 40 80 120 159 199 239 278 318 358 398
##
## $binct
## [-1.36805637964413, -0.264375779976942]
## 40
## [-0.262788549780799, -0.112954465359158]
## 40
## [-0.0999049130865957, -0.0427567472388666]
## 40
## [-0.0427299010942992, -0.00962850641581015]
## 39
## [-0.00852062110239693, 0.0686410557667005]
## 40
## [0.0686670746798195, 0.0848576665096226]
## 40
## [0.0848982901993935, 0.0957727781552671]
## 39
## [0.104636895319711, 0.129296165495585]
## 40
## [0.129524881400457, 0.13816098077111]
## 40
## [0.138167193448738, 0.534108536275362]
## 40
##
## $xtbl
## x
## -1.36805637964413 -0.808819659670332 -0.80046443703434
## 1 1 1
## -0.777135799722571 -0.68516269979883 -0.673100313221977
## 1 1 1
## -0.648752223156282 -0.606398836836818 -0.603319855030218
## 1 1 1
## -0.597722441512497 -0.578773423832448 -0.57684391130148
## 1 1 1
## -0.529815593352933 -0.524086468677554 -0.5176498904133
## 1 1 1
## -0.516149715584571 -0.513876759726287 -0.451141236874604
## 1 1 1
## -0.436150466454698 -0.434954034054065 -0.434496306585217
## 1 1 1
## -0.428984184400305 -0.402059628015283 -0.358793812173929
## 1 1 1
## -0.351660663801291 -0.345726981036656 -0.328353805716596
## 1 1 1
## -0.326353958744406 -0.3248265238159 -0.322162091069332
## 1 1 1
## -0.321141591004627 -0.31518567205632 -0.311485967063435
## 1 1 1
## -0.310351547583572 -0.301299233911081 -0.269408088941477
## 1 1 1
## -0.267462441764412 -0.264569469179821 -0.264482978144487
## 1 1 1
## -0.264375779976942 -0.262788549780799 -0.262200133939682
## 1 1 1
## -0.262111599205768 -0.261951876902474 -0.259938904901168
## 1 1 1
## -0.258973236395966 -0.25777547809935 -0.253434857864135
## 1 1 1
## -0.218733603131825 -0.218006734174117 -0.217730444324147
## 1 1 1
## -0.214817323161161 -0.201894965826654 -0.188362968240248
## 1 1 1
## -0.188141851432206 -0.18757377003693 -0.187045280368567
## 1 1 1
## -0.186702565921112 -0.184858610502618 -0.183365054405671
## 1 1 1
## -0.181834005699018 -0.173421912319852 -0.171956418593441
## 1 1 1
## -0.171672822303532 -0.171085351446871 -0.17036045830443
## 1 1 1
## -0.17016875157196 -0.169498754373239 -0.169498566921955
## 1 1 1
## -0.168241608878119 -0.163096299657128 -0.159778016210773
## 1 1 1
## -0.15935240832811 -0.149359424260547 -0.139941661517663
## 1 1 1
## -0.123704478240967 -0.122454315158257 -0.12153819497135
## 1 1 1
## -0.115985560776651 -0.112954465359158 -0.0999049130865957
## 1 1 1
## -0.0997412388970236 -0.0962669245127241 -0.0946507229463136
## 1 1 1
## -0.092673770449863 -0.0907918990372376 -0.0855495522692611
## 1 1 1
## -0.0851169648669301 -0.0847502415001444 -0.0841057522906004
## 1 1 1
## -0.0831865580272759 -0.0831745779281694 -0.0829063628389761
## 1 1 1
## -0.0823705167059938 -0.0812513774441816 -0.0810170700754389
## 1 1 1
## -0.0803725349680662 -0.0803204139868083 -0.0791230464208725
## 1 1 1
## -0.0769994282512922 -0.0764458326498461 -0.0764377523807669
## 1 1 1
## -0.0760357534176204 -0.0731628856014828 -0.0731499754331892
## 1 1 1
## -0.0719156843990656 -0.0643501441587004 -0.0630524818842193
## 1 1 1
## -0.0575775003029728 -0.0559920002686367 -0.0559122382831991
## 1 1 1
## -0.0558186941206076 -0.0551151867304409 -0.0550880527477129
## 1 1 1
## -0.0543189659914127 -0.0537001706969789 -0.0478457205223473
## 1 1 1
## -0.0432408181813853 -0.042927045923324 -0.0427567472388666
## 1 1 1
## -0.0427299010942992 -0.0423079440850944 -0.0405504723500536
## 1 1 1
## -0.0402410884283393 -0.0398338712489887 -0.0396971746968864
## 1 1 1
## -0.038712988122126 -0.0386195437005843 -0.0385039167104866
## 1 1 1
## -0.0383571254139667 -0.0375537853314238 -0.0370864442248191
## 1 1 1
## -0.0369505845336593 -0.0359655864313184 -0.0355092378736044
## 1 1 1
## -0.0353387274380246 -0.0347521395077528 -0.0338070240400294
## 1 1 1
## -0.0327613878810002 -0.0325581007390241 -0.0324053844679395
## 1 1 1
## -0.0320353856037565 -0.0316961815792847 -0.0316043451541312
## 1 1 1
## -0.0312915480774434 -0.0312646063318958 -0.0312351344283814
## 1 1 1
## -0.0283826005107266 -0.0283393441407833 -0.0279553329298206
## 1 1 1
## -0.023776505602445 -0.0221353627571266 -0.0193813321542904
## 1 1 1
## -0.0190942502310627 -0.013817866052648 -0.0136186738940598
## 1 1 1
## -0.0116029190105722 -0.0110980765458752 -0.00962850641581015
## 1 1 1
## -0.00852062110239693 -0.00707972411716384 -0.00651513952366622
## 1 1 1
## 0.0225333823634137 0.0272516135510663 0.0291096314913415
## 1 1 1
## 0.0299436731991072 0.0543533107747738 0.0602694835288953
## 1 1 1
## 0.0608972550200699 0.0631159214541692 0.0637619155096749
## 1 1 1
## 0.0646982740789968 0.0648193591138588 0.0651147465488349
## 1 1 1
## 0.065219734170729 0.0656137873541891 0.0656200464737038
## 1 1 1
## 0.0657453441815816 0.0657642931475309 0.0663360172103018
## 1 1 1
## 0.0664304017147743 0.0667665595679354 0.0668998485775432
## 1 1 1
## 0.0672195740145627 0.0672567868711127 0.0672917767552303
## 1 1 1
## 0.0673426999533038 0.0674596101305109 0.0675775397283189
## 1 1 1
## 0.0676029429979694 0.067708603663714 0.0680188107692128
## 1 1 1
## 0.0681240453331848 0.0681538204727605 0.0682859929980262
## 1 1 1
## 0.0685085775065772 0.0685225637617075 0.0685438793923394
## 1 1 1
## 0.0686410557667005 0.0686670746798195 0.0688092479012092
## 1 1 1
## 0.0688223282911125 0.0688926833805135 0.0689227157043164
## 1 1 1
## 0.0689700418711938 0.0691211239618643 0.0696267865155267
## 1 1 1
## 0.0697437207189495 0.0698454426954602 0.0706700058205256
## 1 1 1
## 0.0708774020881626 0.0708856092127905 0.0709802940804112
## 1 1 1
## 0.071212343586226 0.0712308641343145 0.0713680904153995
## 1 1 1
## 0.0716286384925997 0.0719197728233376 0.0727599683272464
## 1 1 1
## 0.0731107392875421 0.077001177317765 0.0804687414652105
## 1 1 1
## 0.0805439595112429 0.0811414836188058 0.0814308112207009
## 1 1 1
## 0.0822767140908034 0.0824179029933099 0.0827828823049981
## 1 1 1
## 0.0831965551638221 0.0834721247306169 0.083519293477444
## 1 1 1
## 0.0838163281987099 0.0838557737895356 0.0840314717893647
## 1 1 1
## 0.0840594060855568 0.084113968493305 0.0844500194817503
## 1 1 1
## 0.0845302839515716 0.0848576665096226 0.0848982901993935
## 1 1 1
## 0.0849359506740721 0.0849384534394755 0.0850863784060765
## 1 1 1
## 0.0851644454460172 0.0852143148233719 0.0854969303924624
## 1 1 1
## 0.085635206493176 0.085704947821946 0.0858400295980739
## 1 1 1
## 0.0858929289117762 0.0858986735137215 0.0859559183051268
## 1 1 1
## 0.0860894617698231 0.0861007858039762 0.0861937651386626
## 1 1 1
## 0.0863040457180677 0.0863988094427595 0.0866981351728124
## 1 1 1
## 0.0868686737434243 0.0873152270869046 0.0875642746418299
## 1 1 1
## 0.0876406295170085 0.0876547396571024 0.0884971980249976
## 1 1 1
## 0.0885293855911877 0.0888328271009624 0.089048683910101
## 1 1 1
## 0.0894890025241279 0.0898654780947599 0.0901201821236106
## 1 1 1
## 0.0904795698922434 0.0907959459485361 0.0925381734044699
## 1 1 1
## 0.0931396982311734 0.093452796889931 0.0942543305004822
## 1 1 1
## 0.0944053991869915 0.0957727781552671 0.104636895319711
## 1 1 1
## 0.105100107160417 0.106769089870077 0.107165684330758
## 1 1 1
## 0.107401231170575 0.107688012920448 0.10770043893522
## 1 1 1
## 0.107718914937874 0.108825471173018 0.109197545372269
## 1 1 1
## 0.109316225731382 0.109678638710294 0.110310546095292
## 1 1 1
## 0.110435987876015 0.110602539440114 0.110722068449796
## 1 1 1
## 0.110897116539827 0.111077676573089 0.111323847586308
## 1 1 1
## 0.111383892294853 0.111505377991434 0.111747824065669
## 1 1 1
## 0.11206409030545 0.112328391974226 0.112411778290658
## 1 1 1
## 0.112565794528298 0.112709529184964 0.113311299099927
## 1 1 1
## 0.113923604099377 0.114000768834237 0.114030330549959
## 1 1 1
## 0.114042000081265 0.115169574129629 0.11555191707509
## 1 1 1
## 0.115640729448355 0.116127650709667 0.117197511937642
## 1 1 1
## 0.118382993143462 0.128542579253951 0.129296165495585
## 1 1 1
## 0.129524881400457 0.130171910755517 0.130861259568829
## 1 1 1
## 0.13086249446566 0.131182564543327 0.131407245112461
## 1 1 1
## 0.131582365059481 0.131680980289726 0.131951267670931
## 1 1 1
## 0.131983474739337 0.13212954256098 0.132201441455728
## 1 1 1
## 0.132242540858868 0.132262291737041 0.132530424338477
## 1 1 1
## 0.132897483804346 0.133151424678357 0.133586885048663
## 1 1 1
## 0.133714706675567 0.134313241540789 0.13458832930597
## 1 1 1
## 0.134779784026656 0.134875233112687 0.134954245238555
## 1 1 1
## 0.134957500732565 0.135015589953717 0.135092965440318
## 1 1 1
## 0.135131266024712 0.135170314249714 0.135629629468506
## 1 1 1
## 0.135971967107579 0.136046440066218 0.136363622811779
## 1 1 1
## 0.13643872309721 0.136465893963259 0.136889243957943
## 1 1 1
## 0.137038146187937 0.13786883933441 0.137947288147466
## 1 1 1
## 0.13816098077111 0.138167193448738 0.138418867931289
## 1 1 1
## 0.138537614880816 0.146046594738226 0.146504924323236
## 1 1 1
## 0.150109505288199 0.150166404030363 0.178173706875201
## 1 1 1
## 0.178760143127007 0.184609351603891 0.192222609460657
## 1 1 1
## 0.193030735970841 0.202231484081623 0.203478176278183
## 1 1 1
## 0.204117165832111 0.205452064281375 0.206168092550338
## 1 1 1
## 0.208316552967181 0.208569478982057 0.209772146723001
## 1 1 1
## 0.210847856054049 0.210890942737765 0.211424092717761
## 1 1 1
## 0.211984344408966 0.212988425999724 0.21535379035544
## 1 1 1
## 0.215756845523931 0.228520012472161 0.233297767120711
## 1 1 1
## 0.303147680896973 0.331742504431085 0.333147922950815
## 1 1 1
## 0.360496612014482 0.480414507725707 0.48361609154738
## 1 1 1
## 0.525792946145805 0.530370377920974 0.530410799508147
## 1 1 1
## 0.532008442877618 0.534108536275362
## 1 1
##
## $xval
## [1] -1.368056380 -0.808819660 -0.800464437 -0.777135800 -0.685162700
## [6] -0.673100313 -0.648752223 -0.606398837 -0.603319855 -0.597722442
## [11] -0.578773424 -0.576843911 -0.529815593 -0.524086469 -0.517649890
## [16] -0.516149716 -0.513876760 -0.451141237 -0.436150466 -0.434954034
## [21] -0.434496307 -0.428984184 -0.402059628 -0.358793812 -0.351660664
## [26] -0.345726981 -0.328353806 -0.326353959 -0.324826524 -0.322162091
## [31] -0.321141591 -0.315185672 -0.311485967 -0.310351548 -0.301299234
## [36] -0.269408089 -0.267462442 -0.264569469 -0.264482978 -0.264375780
## [41] -0.262788550 -0.262200134 -0.262111599 -0.261951877 -0.259938905
## [46] -0.258973236 -0.257775478 -0.253434858 -0.218733603 -0.218006734
## [51] -0.217730444 -0.214817323 -0.201894966 -0.188362968 -0.188141851
## [56] -0.187573770 -0.187045280 -0.186702566 -0.184858611 -0.183365054
## [61] -0.181834006 -0.173421912 -0.171956419 -0.171672822 -0.171085351
## [66] -0.170360458 -0.170168752 -0.169498754 -0.169498567 -0.168241609
## [71] -0.163096300 -0.159778016 -0.159352408 -0.149359424 -0.139941662
## [76] -0.123704478 -0.122454315 -0.121538195 -0.115985561 -0.112954465
## [81] -0.099904913 -0.099741239 -0.096266925 -0.094650723 -0.092673770
## [86] -0.090791899 -0.085549552 -0.085116965 -0.084750242 -0.084105752
## [91] -0.083186558 -0.083174578 -0.082906363 -0.082370517 -0.081251377
## [96] -0.081017070 -0.080372535 -0.080320414 -0.079123046 -0.076999428
## [101] -0.076445833 -0.076437752 -0.076035753 -0.073162886 -0.073149975
## [106] -0.071915684 -0.064350144 -0.063052482 -0.057577500 -0.055992000
## [111] -0.055912238 -0.055818694 -0.055115187 -0.055088053 -0.054318966
## [116] -0.053700171 -0.047845721 -0.043240818 -0.042927046 -0.042756747
## [121] -0.042729901 -0.042307944 -0.040550472 -0.040241088 -0.039833871
## [126] -0.039697175 -0.038712988 -0.038619544 -0.038503917 -0.038357125
## [131] -0.037553785 -0.037086444 -0.036950585 -0.035965586 -0.035509238
## [136] -0.035338727 -0.034752140 -0.033807024 -0.032761388 -0.032558101
## [141] -0.032405384 -0.032035386 -0.031696182 -0.031604345 -0.031291548
## [146] -0.031264606 -0.031235134 -0.028382601 -0.028339344 -0.027955333
## [151] -0.023776506 -0.022135363 -0.019381332 -0.019094250 -0.013817866
## [156] -0.013618674 -0.011602919 -0.011098077 -0.009628506 -0.008520621
## [161] -0.007079724 -0.006515140 0.022533382 0.027251614 0.029109631
## [166] 0.029943673 0.054353311 0.060269484 0.060897255 0.063115921
## [171] 0.063761916 0.064698274 0.064819359 0.065114747 0.065219734
## [176] 0.065613787 0.065620046 0.065745344 0.065764293 0.066336017
## [181] 0.066430402 0.066766560 0.066899849 0.067219574 0.067256787
## [186] 0.067291777 0.067342700 0.067459610 0.067577540 0.067602943
## [191] 0.067708604 0.068018811 0.068124045 0.068153820 0.068285993
## [196] 0.068508578 0.068522564 0.068543879 0.068641056 0.068667075
## [201] 0.068809248 0.068822328 0.068892683 0.068922716 0.068970042
## [206] 0.069121124 0.069626787 0.069743721 0.069845443 0.070670006
## [211] 0.070877402 0.070885609 0.070980294 0.071212344 0.071230864
## [216] 0.071368090 0.071628638 0.071919773 0.072759968 0.073110739
## [221] 0.077001177 0.080468741 0.080543960 0.081141484 0.081430811
## [226] 0.082276714 0.082417903 0.082782882 0.083196555 0.083472125
## [231] 0.083519293 0.083816328 0.083855774 0.084031472 0.084059406
## [236] 0.084113968 0.084450019 0.084530284 0.084857667 0.084898290
## [241] 0.084935951 0.084938453 0.085086378 0.085164445 0.085214315
## [246] 0.085496930 0.085635206 0.085704948 0.085840030 0.085892929
## [251] 0.085898674 0.085955918 0.086089462 0.086100786 0.086193765
## [256] 0.086304046 0.086398809 0.086698135 0.086868674 0.087315227
## [261] 0.087564275 0.087640630 0.087654740 0.088497198 0.088529386
## [266] 0.088832827 0.089048684 0.089489003 0.089865478 0.090120182
## [271] 0.090479570 0.090795946 0.092538173 0.093139698 0.093452797
## [276] 0.094254331 0.094405399 0.095772778 0.104636895 0.105100107
## [281] 0.106769090 0.107165684 0.107401231 0.107688013 0.107700439
## [286] 0.107718915 0.108825471 0.109197545 0.109316226 0.109678639
## [291] 0.110310546 0.110435988 0.110602539 0.110722068 0.110897117
## [296] 0.111077677 0.111323848 0.111383892 0.111505378 0.111747824
## [301] 0.112064090 0.112328392 0.112411778 0.112565795 0.112709529
## [306] 0.113311299 0.113923604 0.114000769 0.114030331 0.114042000
## [311] 0.115169574 0.115551917 0.115640729 0.116127651 0.117197512
## [316] 0.118382993 0.128542579 0.129296165 0.129524881 0.130171911
## [321] 0.130861260 0.130862494 0.131182565 0.131407245 0.131582365
## [326] 0.131680980 0.131951268 0.131983475 0.132129543 0.132201441
## [331] 0.132242541 0.132262292 0.132530424 0.132897484 0.133151425
## [336] 0.133586885 0.133714707 0.134313242 0.134588329 0.134779784
## [341] 0.134875233 0.134954245 0.134957501 0.135015590 0.135092965
## [346] 0.135131266 0.135170314 0.135629629 0.135971967 0.136046440
## [351] 0.136363623 0.136438723 0.136465894 0.136889244 0.137038146
## [356] 0.137868839 0.137947288 0.138160981 0.138167193 0.138418868
## [361] 0.138537615 0.146046595 0.146504924 0.150109505 0.150166404
## [366] 0.178173707 0.178760143 0.184609352 0.192222609 0.193030736
## [371] 0.202231484 0.203478176 0.204117166 0.205452064 0.206168093
## [376] 0.208316553 0.208569479 0.209772147 0.210847856 0.210890943
## [381] 0.211424093 0.211984344 0.212988426 0.215353790 0.215756846
## [386] 0.228520012 0.233297767 0.303147681 0.331742504 0.333147923
## [391] 0.360496612 0.480414508 0.483616092 0.525792946 0.530370378
## [396] 0.530410800 0.532008443 0.534108536
##
## $err
## [1] 0.4SS1<-subset(SS, bd <=-0.264375779976942)
SS2<-subset(SS, bd > -0.264375779976942 & bd <= -0.112954465359158)
SS3<-subset(SS, bd >-0.112954465359158 & bd <=-0.0427567472388666)
SS4<-subset(SS, bd >-0.0427567472388666 & bd <=-0.00962850641581015)
SS5<-subset(SS, bd >-0.00962850641581015 & bd <=0.0686410557667005)
SS6<-subset(SS, bd >0.0686410557667005 & bd <=0.0848576665096226)
SS7<-subset(SS, bd >0.0848576665096226 & bd <=0.0957727781552671)
SS8<-subset(SS, bd >0.0957727781552671 & bd <=0.129296165495585)
SS9<-subset(SS, bd >0.129296165495585 & bd <=0.13816098077111)
SS10<-subset(SS, bd >0.13816098077111 & bd <=0.534108536275362)
SS1$decile<-1
SS2$decile<-2
SS3$decile<-3
SS4$decile<-4
SS5$decile<-5
SS6$decile<-6
SS7$decile<-7
SS8$decile<-8
SS9$decile<-9
SS10$decile<-10
SS1<-subset(SS1, select= c(23, 24, 28:29))
SS2<-subset(SS2, select= c(23, 24, 28:29))
SS3<-subset(SS3, select= c(23, 24, 28:29))
SS4<-subset(SS4, select= c(23, 24, 28:29))
SS5<-subset(SS5, select= c(23, 24, 28:29))
SS6<-subset(SS6, select= c(23, 24, 28:29))
SS7<-subset(SS7, select= c(23, 24, 28:29))
SS8<-subset(SS8, select= c(23, 24, 28:29))
SS9<-subset(SS9, select= c(23, 24, 28:29))
SS10<-subset(SS10, select= c(23, 24, 28:29))
SScond<-rbind(SS1, SS2, SS3, SS4, SS5, SS6, SS7, SS8, SS9, SS10)
head(SScond)## cond SS bd decile
## 38 SS 1 -0.3457270 1
## 45 SS 1 -0.5787734 1
## 93 SS 1 -0.3211416 1
## 111 SS 1 -1.3680564 1
## 122 SS 1 -0.5138768 1
## 123 SS 1 -0.3221621 1SScond<-subset(SScond, select=c(1,2,4))
LLSS<-rbind(SScond, LLcond)
bins(no$bd, 10, max.breaks=10)## $binlo
## [1] 1 32 62 92 122 152 182 212 242 271
##
## $binhi
## [1] 31 61 91 121 151 181 211 241 270 301
##
## $binct
## [-1.38645212505623, -0.227223217775932]
## 31
## [-0.203101268954379, -0.08382568091853]
## 30
## [-0.0821786946466201, -0.00835404473681023]
## 30
## [-0.00767433993924627, 0.0812511890237588]
## 30
## [0.0815625498081314, 0.0853754962890645]
## 30
## [0.0855232835434692, 0.0929411583232294]
## 30
## [0.0936741800566516, 0.129971931536811]
## 30
## [0.130421457193488, 0.135192438875584]
## 30
## [0.135288314093837, 0.191413204826878]
## 30
## [0.206955878594119, 0.579165143053758]
## 31
##
## $xtbl
## x
## -1.38645212505623 -1.36957637438239 -1.14320986116188
## 1 1 1
## -0.939363256448824 -0.854213958214141 -0.796858925844065
## 1 1 1
## -0.650107331282998 -0.637416559131636 -0.635897660702729
## 1 1 1
## -0.602017569167831 -0.562838592463594 -0.542615055876405
## 1 1 1
## -0.474555979659809 -0.41243895023298 -0.383327837830421
## 1 1 1
## -0.350027663079732 -0.333746647588586 -0.332823213087264
## 1 1 1
## -0.329763669879146 -0.29118984474779 -0.287564811868721
## 1 1 1
## -0.285327269937615 -0.267761866115502 -0.266733362207131
## 1 1 1
## -0.266502827053659 -0.264551704823618 -0.26108614085649
## 1 1 1
## -0.259146112025253 -0.254112293139151 -0.232259009377734
## 1 1 1
## -0.227223217775932 -0.203101268954379 -0.188395760972889
## 1 1 1
## -0.183202951171981 -0.182478012939671 -0.175286160800199
## 1 1 1
## -0.175245380431386 -0.173885367845595 -0.173724398148916
## 1 1 1
## -0.171956413940737 -0.171400176992773 -0.171359175216715
## 1 1 1
## -0.170275391914284 -0.160439377878859 -0.152625107610848
## 1 1 1
## -0.147046338549776 -0.137707699889599 -0.136339874275679
## 1 1 1
## -0.135650152527443 -0.134148140828038 -0.131190949332481
## 1 1 1
## -0.128650770498482 -0.119162380697966 -0.115352644666026
## 1 1 1
## -0.113849823135898 -0.101894636565401 -0.0987900357509119
## 1 1 1
## -0.0873517492936116 -0.0861815962046477 -0.0844626068936045
## 1 1 1
## -0.08382568091853 -0.0821786946466201 -0.0787494390282097
## 1 1 1
## -0.0777665108857502 -0.0768850750172234 -0.076077717623019
## 1 1 1
## -0.0742414433454344 -0.0741492197611777 -0.0626555364891952
## 1 1 1
## -0.0602679541667697 -0.0587381408044969 -0.0582414487245
## 1 1 1
## -0.0573589449672194 -0.0573179348586697 -0.0548521808474418
## 1 1 1
## -0.0508498198751737 -0.0460714736329706 -0.038710459784943
## 1 1 1
## -0.0383809501588825 -0.0379223187420317 -0.0377394196766421
## 1 1 1
## -0.0363736710371652 -0.0281614862344557 -0.0198758915711927
## 1 1 1
## -0.0197349816691936 -0.0178467507729718 -0.0171929116230003
## 1 1 1
## -0.0163709335412074 -0.0116258383461387 -0.00993782284360445
## 1 1 1
## -0.00835404473681023 -0.00767433993924627 -0.00689815683517225
## 1 1 1
## 0.00499088765447375 0.00788609225893633 0.0155392897347042
## 1 1 1
## 0.0195642778391472 0.0197313226311082 0.0294411284502423
## 1 1 1
## 0.0408990981648811 0.045254943906222 0.0627776361281887
## 1 1 1
## 0.0678597591956758 0.0685797048954415 0.0717654009877755
## 1 1 1
## 0.0744388126792255 0.0749891398111888 0.0757242860028233
## 1 1 1
## 0.0759063947834497 0.0762294683422406 0.0763229083798742
## 1 1 1
## 0.0768685534782678 0.0769242602762308 0.0772120248211645
## 1 1 1
## 0.0778210210330996 0.0778551436069064 0.0797868413064744
## 1 1 1
## 0.0806447292836037 0.0808396988624152 0.0808643301794365
## 1 1 1
## 0.0812511890237588 0.0815625498081314 0.0817541174598819
## 1 1 1
## 0.0819517132075471 0.0822424028923581 0.0823897904080666
## 1 1 1
## 0.0824956610620418 0.0825864355724746 0.0826814362123076
## 1 1 1
## 0.083381319079221 0.0834788484771498 0.0836180951871615
## 1 1 1
## 0.0837819199983469 0.0838357272216932 0.0839944820529031
## 1 1 1
## 0.0840202497248503 0.0840499781549884 0.0841374019664027
## 1 1 1
## 0.0841476734279096 0.0841690591613629 0.0841924561388079
## 1 1 1
## 0.0844108850326226 0.084432975797665 0.0846651938903946
## 1 1 1
## 0.0847241984679381 0.0847349690048884 0.0850154211278404
## 1 1 1
## 0.0851066807003265 0.0853241480947728 0.0853621370658568
## 1 1 1
## 0.0853754962890645 0.0855232835434692 0.085554231044479
## 1 1 1
## 0.0855836724468579 0.085682338440894 0.0856928486475943
## 1 1 1
## 0.0858590849246039 0.0861941052972268 0.0863638113872347
## 1 1 1
## 0.086456088443998 0.0864891876595979 0.0866230910746035
## 1 1 1
## 0.0869534453472984 0.0871089031426762 0.087543129196394
## 1 1 1
## 0.0877734511702668 0.087833348284617 0.087859624555946
## 1 1 1
## 0.0881082504869668 0.0881155594289056 0.0888385761599426
## 1 1 1
## 0.0889758870388297 0.0889956316803522 0.0894155158560928
## 1 1 1
## 0.0895401414161703 0.0906263376514862 0.0912137272513505
## 1 1 1
## 0.0913892218995598 0.0914017666396092 0.0928490548094063
## 1 1 1
## 0.0929411583232294 0.0936741800566516 0.0960600115333641
## 1 1 1
## 0.09859051179178 0.102533086322563 0.106867824653706
## 1 1 1
## 0.108179185448947 0.109378551115839 0.109676061413151
## 1 1 1
## 0.110848341442624 0.111713332131165 0.11190867631987
## 1 1 1
## 0.112355396049886 0.112481015114839 0.113072496146609
## 1 1 1
## 0.113644425233889 0.113866816390387 0.115636396225774
## 1 1 1
## 0.116040458735209 0.116705210125594 0.11742530237678
## 1 1 1
## 0.118622500710069 0.118959022275395 0.119070729862743
## 1 1 1
## 0.119651078123019 0.123514501892562 0.123938727208264
## 1 1 1
## 0.124476836589756 0.125934774532678 0.128775995801547
## 1 1 1
## 0.129971931536811 0.130421457193488 0.130569021768124
## 1 1 1
## 0.130643525491348 0.130785776592459 0.130856902535819
## 1 1 1
## 0.131123117871284 0.131521684633244 0.131768329102695
## 1 1 1
## 0.131922161127486 0.132942250524498 0.133466827876439
## 1 1 1
## 0.133534814897237 0.133556483633424 0.133686304442375
## 1 1 1
## 0.133727716351774 0.134104914297631 0.134170754600798
## 1 1 1
## 0.134220920277434 0.134273390828604 0.134362180268855
## 1 1 1
## 0.134519741122491 0.134539248726989 0.13458526313855
## 1 1 1
## 0.134835970694704 0.134878431794984 0.134887675452722
## 1 1 1
## 0.134919632025332 0.13512459882739 0.135178370463376
## 1 1 1
## 0.135192438875584 0.135288314093837 0.135486433587855
## 1 1 1
## 0.1356065290131 0.135659079442498 0.135834248511602
## 1 1 1
## 0.135986571130569 0.136263377792931 0.136270820324624
## 1 1 1
## 0.136272286465488 0.136288685591655 0.136449680414525
## 1 1 1
## 0.136598819229656 0.136684225213288 0.137140830417861
## 1 1 1
## 0.13743087622145 0.137452535786205 0.138316184654101
## 1 1 1
## 0.139031581826237 0.139653415933545 0.139894524663986
## 2 1 1
## 0.140264467501282 0.145255034945966 0.146526482877685
## 1 1 1
## 0.146854102369194 0.149145166367166 0.151710428807325
## 1 1 1
## 0.151836411347817 0.153443863578322 0.191413204826878
## 1 1 1
## 0.206955878594119 0.20747611964525 0.208294362605784
## 1 1 1
## 0.208680220072718 0.209825610557652 0.210140288381574
## 1 1 1
## 0.210795531416506 0.227744450564301 0.22780850220279
## 1 1 1
## 0.234551942276499 0.236091728991389 0.237351786876024
## 1 1 1
## 0.237642717351108 0.253633036799568 0.258134158583926
## 1 1 1
## 0.259599950285197 0.259793599126469 0.260214355496027
## 1 1 1
## 0.260530209829675 0.325724495071878 0.325993219017629
## 1 1 1
## 0.326367585098578 0.329324991151206 0.33600365853628
## 1 1 1
## 0.529579756357625 0.53246972845501 0.534488371136465
## 1 1 1
## 0.534991895639789 0.537933153626412 0.575323343077641
## 1 1 1
## 0.579165143053758
## 1
##
## $xval
## [1] -1.386452125 -1.369576374 -1.143209861 -0.939363256 -0.854213958
## [6] -0.796858926 -0.650107331 -0.637416559 -0.635897661 -0.602017569
## [11] -0.562838592 -0.542615056 -0.474555980 -0.412438950 -0.383327838
## [16] -0.350027663 -0.333746648 -0.332823213 -0.329763670 -0.291189845
## [21] -0.287564812 -0.285327270 -0.267761866 -0.266733362 -0.266502827
## [26] -0.264551705 -0.261086141 -0.259146112 -0.254112293 -0.232259009
## [31] -0.227223218 -0.203101269 -0.188395761 -0.183202951 -0.182478013
## [36] -0.175286161 -0.175245380 -0.173885368 -0.173724398 -0.171956414
## [41] -0.171400177 -0.171359175 -0.170275392 -0.160439378 -0.152625108
## [46] -0.147046339 -0.137707700 -0.136339874 -0.135650153 -0.134148141
## [51] -0.131190949 -0.128650770 -0.119162381 -0.115352645 -0.113849823
## [56] -0.101894637 -0.098790036 -0.087351749 -0.086181596 -0.084462607
## [61] -0.083825681 -0.082178695 -0.078749439 -0.077766511 -0.076885075
## [66] -0.076077718 -0.074241443 -0.074149220 -0.062655536 -0.060267954
## [71] -0.058738141 -0.058241449 -0.057358945 -0.057317935 -0.054852181
## [76] -0.050849820 -0.046071474 -0.038710460 -0.038380950 -0.037922319
## [81] -0.037739420 -0.036373671 -0.028161486 -0.019875892 -0.019734982
## [86] -0.017846751 -0.017192912 -0.016370934 -0.011625838 -0.009937823
## [91] -0.008354045 -0.007674340 -0.006898157 0.004990888 0.007886092
## [96] 0.015539290 0.019564278 0.019731323 0.029441128 0.040899098
## [101] 0.045254944 0.062777636 0.067859759 0.068579705 0.071765401
## [106] 0.074438813 0.074989140 0.075724286 0.075906395 0.076229468
## [111] 0.076322908 0.076868553 0.076924260 0.077212025 0.077821021
## [116] 0.077855144 0.079786841 0.080644729 0.080839699 0.080864330
## [121] 0.081251189 0.081562550 0.081754117 0.081951713 0.082242403
## [126] 0.082389790 0.082495661 0.082586436 0.082681436 0.083381319
## [131] 0.083478848 0.083618095 0.083781920 0.083835727 0.083994482
## [136] 0.084020250 0.084049978 0.084137402 0.084147673 0.084169059
## [141] 0.084192456 0.084410885 0.084432976 0.084665194 0.084724198
## [146] 0.084734969 0.085015421 0.085106681 0.085324148 0.085362137
## [151] 0.085375496 0.085523284 0.085554231 0.085583672 0.085682338
## [156] 0.085692849 0.085859085 0.086194105 0.086363811 0.086456088
## [161] 0.086489188 0.086623091 0.086953445 0.087108903 0.087543129
## [166] 0.087773451 0.087833348 0.087859625 0.088108250 0.088115559
## [171] 0.088838576 0.088975887 0.088995632 0.089415516 0.089540141
## [176] 0.090626338 0.091213727 0.091389222 0.091401767 0.092849055
## [181] 0.092941158 0.093674180 0.096060012 0.098590512 0.102533086
## [186] 0.106867825 0.108179185 0.109378551 0.109676061 0.110848341
## [191] 0.111713332 0.111908676 0.112355396 0.112481015 0.113072496
## [196] 0.113644425 0.113866816 0.115636396 0.116040459 0.116705210
## [201] 0.117425302 0.118622501 0.118959022 0.119070730 0.119651078
## [206] 0.123514502 0.123938727 0.124476837 0.125934775 0.128775996
## [211] 0.129971932 0.130421457 0.130569022 0.130643525 0.130785777
## [216] 0.130856903 0.131123118 0.131521685 0.131768329 0.131922161
## [221] 0.132942251 0.133466828 0.133534815 0.133556484 0.133686304
## [226] 0.133727716 0.134104914 0.134170755 0.134220920 0.134273391
## [231] 0.134362180 0.134519741 0.134539249 0.134585263 0.134835971
## [236] 0.134878432 0.134887675 0.134919632 0.135124599 0.135178370
## [241] 0.135192439 0.135288314 0.135486434 0.135606529 0.135659079
## [246] 0.135834249 0.135986571 0.136263378 0.136270820 0.136272286
## [251] 0.136288686 0.136449680 0.136598819 0.136684225 0.137140830
## [256] 0.137430876 0.137452536 0.138316185 0.139031582 0.139653416
## [261] 0.139894525 0.140264468 0.145255035 0.146526483 0.146854102
## [266] 0.149145166 0.151710429 0.151836411 0.153443864 0.191413205
## [271] 0.206955879 0.207476120 0.208294363 0.208680220 0.209825611
## [276] 0.210140288 0.210795531 0.227744451 0.227808502 0.234551942
## [281] 0.236091729 0.237351787 0.237642717 0.253633037 0.258134159
## [286] 0.259599950 0.259793599 0.260214355 0.260530210 0.325724495
## [291] 0.325993219 0.326367585 0.329324991 0.336003659 0.529579756
## [296] 0.532469728 0.534488371 0.534991896 0.537933154 0.575323343
## [301] 0.579165143
##
## $err
## [1] 0.4no1<-subset(no, bd <=-0.227223217775932)
no2<-subset(no, bd > -0.227223217775932 & bd <= -0.08382568091853)
no3<-subset(no, bd >-0.08382568091853 & bd <=-0.00835404473681023)
no4<-subset(no, bd >-0.00835404473681023 & bd <=0.0812511890237588)
no5<-subset(no, bd >0.0812511890237588 & bd <=0.0853754962890645)
no6<-subset(no, bd >0.0853754962890645 & bd <=0.0929411583232294)
no7<-subset(no, bd >0.0929411583232294 & bd <=0.129971931536811)
no8<-subset(no, bd >0.129971931536811 & bd <=0.135192438875584)
no9<-subset(no, bd >0.135192438875584 & bd <=0.191413204826878)
no10<-subset(no, bd >0.191413204826878 & bd <=0.579165143053758)
no1$decile<-1
no2$decile<-2
no3$decile<-3
no4$decile<-4
no5$decile<-5
no6$decile<-6
no7$decile<-7
no8$decile<-8
no9$decile<-9
no10$decile<-10
names(no1)## [1] "X" "mturkID" "serial.x"
## [4] "beta" "delta" "id"
## [7] "Completed" "Wave" "IP.address"
## [10] "reader" "gender" "age"
## [13] "income" "education" "smoking"
## [16] "email" "X.studies.x" "X.studiesweek.x"
## [19] "defaultstudies.x" "serial.y" "id.x"
## [22] "interviewtime" "cond" "SS"
## [25] "intent" "intent2" "believe"
## [28] "bd" "decile"no1<-subset(no1, select= c(23, 24, 28:29))
no2<-subset(no2, select= c(23, 24, 28:29))
no3<-subset(no3, select= c(23, 24, 28:29))
no4<-subset(no4, select= c(23, 24, 28:29))
no5<-subset(no5, select= c(23, 24, 28:29))
no6<-subset(no6, select= c(23, 24, 28:29))
no7<-subset(no7, select= c(23, 24, 28:29))
no8<-subset(no8, select= c(23, 24, 28:29))
no9<-subset(no9, select= c(23, 24, 28:29))
no10<-subset(no10, select= c(23, 24, 28:29))
nocond<-rbind(no1, no2, no3, no4, no5, no6, no7, no8, no9, no10)
head(nocond)## cond SS bd decile
## 208 no 1 -0.2665028 1
## 234 no 1 -1.3864521 1
## 253 no 0 -0.6374166 1
## 255 no 1 -0.6020176 1
## 260 no 1 -0.8542140 1
## 274 no 1 -0.9393633 1nocond<-subset(nocond, select=c(1,2,4))
allconds<-rbind(LLSS, nocond)
library(ggplot2)
library(plyr)
cdata <- ddply(allconds, c("decile", "cond"), summarise,
N = length(SS),
mean = mean(SS),
sd = sd(SS),
se = sd / sqrt(N)
)
ggplot(cdata, aes(x=decile, y=mean, colour=cond)) +
geom_errorbar(aes(ymin=mean-se, ymax=mean+se), width=.1) +
geom_line() +
geom_point()