This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:
Upload the Data
```r
library(vcfR)##
## ***** *** vcfR *** *****
## This is vcfR 1.13.0
## browseVignettes('vcfR') # Documentation
## citation('vcfR') # Citation
## ***** ***** ***** *****getwd()## [1] "/Users/oliviabrull/Downloads"setwd("/Users/oliviabrull/Downloads")
list.files(pattern= "vcf")## [1] "1.159051856-159301856.ALL.chr1_GRCh38.genotypes.20170504.vcf.gz"
## [2] "15.27865604-28105604.ALL.chr15_GRCh38.genotypes.20170504.vcf.gz"
## [3] "2.136483646-136733646.ALL.chr2_GRCh38.genotypes.20170504.vcf.gz"
## [4] "all_loci-1.vcf"
## [5] "all_loci.vcf"
## [6] "ALL.chr15_GRCh38.genotypes.20170504.vcf.gz"
## [7] "code_checkpoint_vcfR.docx"
## [8] "code_checkpoint_vcfR.html"
## [9] "code_checkpoint_vcfR.Rmd"
## [10] "vcfR_test.vcf"
## [11] "vcfR_test.vcf.gz"snps <- vcfR::read.vcfR(file = "/Users/oliviabrull/Downloads/15.27865604-28105604.ALL.chr15_GRCh38.genotypes.20170504.vcf.gz")## Scanning file to determine attributes.
## File attributes:
## meta lines: 130
## header_line: 131
## variant count: 6946
## column count: 2513
##
Meta line 130 read in.
## All meta lines processed.
## gt matrix initialized.
## Character matrix gt created.
## Character matrix gt rows: 6946
## Character matrix gt cols: 2513
## skip: 0
## nrows: 6946
## row_num: 0
##
Processed variant 1000
Processed variant 2000
Processed variant 3000
Processed variant 4000
Processed variant 5000
Processed variant 6000
Processed variant: 6946
## All variants processedExtract SNPs
my_snps <- vcfR::extract.gt(snps,
element = "GT",
IDtoRowNames = F,
as.numeric = T,
convertNA = T,
return.alleles = F)Note that the echo = FALSE parameter was added to the
code chunk to prevent printing of the R code that generated the
plot.
Remove make dataframe and transform data
snps_df <- t(my_snps)
snps_df <- data.frame(snps_df)
dim(snps_df)## [1] 2504 6946invar_omit <- function(x){
cat("Dataframe of dim", dim(x), "processed...\n")
sds <- apply(x, 2, sd, na.rm = TRUE)
i_var0 <- which(sds == 0)
cat(length(i_var0), "columns removed\n")
if(length(i_var0) > 0){
x <- x[, -i_var0] }
return(x) }
names(snps_df) #no invars## [1] "X1" "X2" "X3" "X4" "X5" "X6" "X7" "X8" "X9"
## [10] "X10" "X11" "X12" "X13" "X14" "X15" "X16" "X17" "X18"
## [19] "X19" "X20" "X21" "X22" "X23" "X24" "X25" "X26" "X27"
## [28] "X28" "X29" "X30" "X31" "X32" "X33" "X34" "X35" "X36"
## [37] "X37" "X38" "X39" "X40" "X41" "X42" "X43" "X44" "X45"
## [46] "X46" "X47" "X48" "X49" "X50" "X51" "X52" "X53" "X54"
## [55] "X55" "X56" "X57" "X58" "X59" "X60" "X61" "X62" "X63"
## [64] "X64" "X65" "X66" "X67" "X68" "X69" "X70" "X71" "X72"
## [73] "X73" "X74" "X75" "X76" "X77" "X78" "X79" "X80" "X81"
## [82] "X82" "X83" "X84" "X85" "X86" "X87" "X88" "X89" "X90"
## [91] "X91" "X92" "X93" "X94" "X95" "X96" "X97" "X98" "X99"
## [100] "X100" "X101" "X102" "X103" "X104" "X105" "X106" "X107" "X108"
## [109] "X109" "X110" "X111" "X112" "X113" "X114" "X115" "X116" "X117"
## [118] "X118" "X119" "X120" "X121" "X122" "X123" "X124" "X125" "X126"
## [127] "X127" "X128" "X129" "X130" "X131" "X132" "X133" "X134" "X135"
## [136] "X136" "X137" "X138" "X139" "X140" "X141" "X142" "X143" "X144"
## [145] "X145" "X146" "X147" "X148" "X149" "X150" "X151" "X152" "X153"
## [154] "X154" "X155" "X156" "X157" "X158" "X159" "X160" "X161" "X162"
## [163] "X163" "X164" "X165" "X166" "X167" "X168" "X169" "X170" "X171"
## [172] "X172" "X173" "X174" "X175" "X176" "X177" "X178" "X179" "X180"
## [181] "X181" "X182" "X183" "X184" "X185" "X186" "X187" "X188" "X189"
## [190] "X190" "X191" "X192" "X193" "X194" "X195" "X196" "X197" "X198"
## [199] "X199" "X200" "X201" "X202" "X203" "X204" "X205" "X206" "X207"
## [208] "X208" "X209" "X210" "X211" "X212" "X213" "X214" "X215" "X216"
## [217] "X217" "X218" "X219" "X220" "X221" "X222" "X223" "X224" "X225"
## [226] "X226" "X227" "X228" "X229" "X230" "X231" "X232" "X233" "X234"
## [235] "X235" "X236" "X237" "X238" "X239" "X240" "X241" "X242" "X243"
## [244] "X244" "X245" "X246" "X247" "X248" "X249" "X250" "X251" "X252"
## [253] "X253" "X254" "X255" "X256" "X257" "X258" "X259" "X260" "X261"
## [262] "X262" "X263" "X264" "X265" "X266" "X267" "X268" "X269" "X270"
## [271] "X271" "X272" "X273" "X274" "X275" "X276" "X277" "X278" "X279"
## [280] "X280" "X281" "X282" "X283" "X284" "X285" "X286" "X287" "X288"
## [289] "X289" "X290" "X291" "X292" "X293" "X294" "X295" "X296" "X297"
## [298] "X298" "X299" "X300" "X301" "X302" "X303" "X304" "X305" "X306"
## [307] "X307" "X308" "X309" "X310" "X311" "X312" "X313" "X314" "X315"
## [316] "X316" "X317" "X318" "X319" "X320" "X321" "X322" "X323" "X324"
## [325] "X325" "X326" "X327" "X328" "X329" "X330" "X331" "X332" "X333"
## [334] "X334" "X335" "X336" "X337" "X338" "X339" "X340" "X341" "X342"
## [343] "X343" "X344" "X345" "X346" "X347" "X348" "X349" "X350" "X351"
## [352] "X352" "X353" "X354" "X355" "X356" "X357" "X358" "X359" "X360"
## [361] "X361" "X362" "X363" "X364" "X365" "X366" "X367" "X368" "X369"
## [370] "X370" "X371" "X372" "X373" "X374" "X375" "X376" "X377" "X378"
## [379] "X379" "X380" "X381" "X382" "X383" "X384" "X385" "X386" "X387"
## [388] "X388" "X389" "X390" "X391" "X392" "X393" "X394" "X395" "X396"
## [397] "X397" "X398" "X399" "X400" "X401" "X402" "X403" "X404" "X405"
## [406] "X406" "X407" "X408" "X409" "X410" "X411" "X412" "X413" "X414"
## [415] "X415" "X416" "X417" "X418" "X419" "X420" "X421" "X422" "X423"
## [424] "X424" "X425" "X426" "X427" "X428" "X429" "X430" "X431" "X432"
## [433] "X433" "X434" "X435" "X436" "X437" "X438" "X439" "X440" "X441"
## [442] "X442" "X443" "X444" "X445" "X446" "X447" "X448" "X449" "X450"
## [451] "X451" "X452" "X453" "X454" "X455" "X456" "X457" "X458" "X459"
## [460] "X460" "X461" "X462" "X463" "X464" "X465" "X466" "X467" "X468"
## [469] "X469" "X470" "X471" "X472" "X473" "X474" "X475" "X476" "X477"
## [478] "X478" "X479" "X480" "X481" "X482" "X483" "X484" "X485" "X486"
## [487] "X487" "X488" "X489" "X490" "X491" "X492" "X493" "X494" "X495"
## [496] "X496" "X497" "X498" "X499" "X500" "X501" "X502" "X503" "X504"
## [505] "X505" "X506" "X507" "X508" "X509" "X510" "X511" "X512" "X513"
## [514] "X514" "X515" "X516" "X517" "X518" "X519" "X520" "X521" "X522"
## [523] "X523" "X524" "X525" "X526" "X527" "X528" "X529" "X530" "X531"
## [532] "X532" "X533" "X534" "X535" "X536" "X537" "X538" "X539" "X540"
## [541] "X541" "X542" "X543" "X544" "X545" "X546" "X547" "X548" "X549"
## [550] "X550" "X551" "X552" "X553" "X554" "X555" "X556" "X557" "X558"
## [559] "X559" "X560" "X561" "X562" "X563" "X564" "X565" "X566" "X567"
## [568] "X568" "X569" "X570" "X571" "X572" "X573" "X574" "X575" "X576"
## [577] "X577" "X578" "X579" "X580" "X581" "X582" "X583" "X584" "X585"
## [586] "X586" "X587" "X588" "X589" "X590" "X591" "X592" "X593" "X594"
## [595] "X595" "X596" "X597" "X598" "X599" "X600" "X601" "X602" "X603"
## [604] "X604" "X605" "X606" "X607" "X608" "X609" "X610" "X611" "X612"
## [613] "X613" "X614" "X615" "X616" "X617" "X618" "X619" "X620" "X621"
## [622] "X622" "X623" "X624" "X625" "X626" "X627" "X628" "X629" "X630"
## [631] "X631" "X632" "X633" "X634" "X635" "X636" "X637" "X638" "X639"
## [640] "X640" "X641" "X642" "X643" "X644" "X645" "X646" "X647" "X648"
## [649] "X649" "X650" "X651" "X652" "X653" "X654" "X655" "X656" "X657"
## [658] "X658" "X659" "X660" "X661" "X662" "X663" "X664" "X665" "X666"
## [667] "X667" "X668" "X669" "X670" "X671" "X672" "X673" "X674" "X675"
## [676] "X676" "X677" "X678" "X679" "X680" "X681" "X682" "X683" "X684"
## [685] "X685" "X686" "X687" "X688" "X689" "X690" "X691" "X692" "X693"
## [694] "X694" "X695" "X696" "X697" "X698" "X699" "X700" "X701" "X702"
## [703] "X703" "X704" "X705" "X706" "X707" "X708" "X709" "X710" "X711"
## [712] "X712" "X713" "X714" "X715" "X716" "X717" "X718" "X719" "X720"
## [721] "X721" "X722" "X723" "X724" "X725" "X726" "X727" "X728" "X729"
## [730] "X730" "X731" "X732" "X733" "X734" "X735" "X736" "X737" "X738"
## [739] "X739" "X740" "X741" "X742" "X743" "X744" "X745" "X746" "X747"
## [748] "X748" "X749" "X750" "X751" "X752" "X753" "X754" "X755" "X756"
## [757] "X757" "X758" "X759" "X760" "X761" "X762" "X763" "X764" "X765"
## [766] "X766" "X767" "X768" "X769" "X770" "X771" "X772" "X773" "X774"
## [775] "X775" "X776" "X777" "X778" "X779" "X780" "X781" "X782" "X783"
## [784] "X784" "X785" "X786" "X787" "X788" "X789" "X790" "X791" "X792"
## [793] "X793" "X794" "X795" "X796" "X797" "X798" "X799" "X800" "X801"
## [802] "X802" "X803" "X804" "X805" "X806" "X807" "X808" "X809" "X810"
## [811] "X811" "X812" "X813" "X814" "X815" "X816" "X817" "X818" "X819"
## [820] "X820" "X821" "X822" "X823" "X824" "X825" "X826" "X827" "X828"
## [829] "X829" "X830" "X831" "X832" "X833" "X834" "X835" "X836" "X837"
## [838] "X838" "X839" "X840" "X841" "X842" "X843" "X844" "X845" "X846"
## [847] "X847" "X848" "X849" "X850" "X851" "X852" "X853" "X854" "X855"
## [856] "X856" "X857" "X858" "X859" "X860" "X861" "X862" "X863" "X864"
## [865] "X865" "X866" "X867" "X868" "X869" "X870" "X871" "X872" "X873"
## [874] "X874" "X875" "X876" "X877" "X878" "X879" "X880" "X881" "X882"
## [883] "X883" "X884" "X885" "X886" "X887" "X888" "X889" "X890" "X891"
## [892] "X892" "X893" "X894" "X895" "X896" "X897" "X898" "X899" "X900"
## [901] "X901" "X902" "X903" "X904" "X905" "X906" "X907" "X908" "X909"
## [910] "X910" "X911" "X912" "X913" "X914" "X915" "X916" "X917" "X918"
## [919] "X919" "X920" "X921" "X922" "X923" "X924" "X925" "X926" "X927"
## [928] "X928" "X929" "X930" "X931" "X932" "X933" "X934" "X935" "X936"
## [937] "X937" "X938" "X939" "X940" "X941" "X942" "X943" "X944" "X945"
## [946] "X946" "X947" "X948" "X949" "X950" "X951" "X952" "X953" "X954"
## [955] "X955" "X956" "X957" "X958" "X959" "X960" "X961" "X962" "X963"
## [964] "X964" "X965" "X966" "X967" "X968" "X969" "X970" "X971" "X972"
## [973] "X973" "X974" "X975" "X976" "X977" "X978" "X979" "X980" "X981"
## [982] "X982" "X983" "X984" "X985" "X986" "X987" "X988" "X989" "X990"
## [991] "X991" "X992" "X993" "X994" "X995" "X996" "X997" "X998" "X999"
## [1000] "X1000" "X1001" "X1002" "X1003" "X1004" "X1005" "X1006" "X1007" "X1008"
## [1009] "X1009" "X1010" "X1011" "X1012" "X1013" "X1014" "X1015" "X1016" "X1017"
## [1018] "X1018" "X1019" "X1020" "X1021" "X1022" "X1023" "X1024" "X1025" "X1026"
## [1027] "X1027" "X1028" "X1029" "X1030" "X1031" "X1032" "X1033" "X1034" "X1035"
## [1036] "X1036" "X1037" "X1038" "X1039" "X1040" "X1041" "X1042" "X1043" "X1044"
## [1045] "X1045" "X1046" "X1047" "X1048" "X1049" "X1050" "X1051" "X1052" "X1053"
## [1054] "X1054" "X1055" "X1056" "X1057" "X1058" "X1059" "X1060" "X1061" "X1062"
## [1063] "X1063" "X1064" "X1065" "X1066" "X1067" "X1068" "X1069" "X1070" "X1071"
## [1072] "X1072" "X1073" "X1074" "X1075" "X1076" "X1077" "X1078" "X1079" "X1080"
## [1081] "X1081" "X1082" "X1083" "X1084" "X1085" "X1086" "X1087" "X1088" "X1089"
## [1090] "X1090" "X1091" "X1092" "X1093" "X1094" "X1095" "X1096" "X1097" "X1098"
## [1099] "X1099" "X1100" "X1101" "X1102" "X1103" "X1104" "X1105" "X1106" "X1107"
## [1108] "X1108" "X1109" "X1110" "X1111" "X1112" "X1113" "X1114" "X1115" "X1116"
## [1117] "X1117" "X1118" "X1119" "X1120" "X1121" "X1122" "X1123" "X1124" "X1125"
## [1126] "X1126" "X1127" "X1128" "X1129" "X1130" "X1131" "X1132" "X1133" "X1134"
## [1135] "X1135" "X1136" "X1137" "X1138" "X1139" "X1140" "X1141" "X1142" "X1143"
## [1144] "X1144" "X1145" "X1146" "X1147" "X1148" "X1149" "X1150" "X1151" "X1152"
## [1153] "X1153" "X1154" "X1155" "X1156" "X1157" "X1158" "X1159" "X1160" "X1161"
## [1162] "X1162" "X1163" "X1164" "X1165" "X1166" "X1167" "X1168" "X1169" "X1170"
## [1171] "X1171" "X1172" "X1173" "X1174" "X1175" "X1176" "X1177" "X1178" "X1179"
## [1180] "X1180" "X1181" "X1182" "X1183" "X1184" "X1185" "X1186" "X1187" "X1188"
## [1189] "X1189" "X1190" "X1191" "X1192" "X1193" "X1194" "X1195" "X1196" "X1197"
## [1198] "X1198" "X1199" "X1200" "X1201" "X1202" "X1203" "X1204" "X1205" "X1206"
## [1207] "X1207" "X1208" "X1209" "X1210" "X1211" "X1212" "X1213" "X1214" "X1215"
## [1216] "X1216" "X1217" "X1218" "X1219" "X1220" "X1221" "X1222" "X1223" "X1224"
## [1225] "X1225" "X1226" "X1227" "X1228" "X1229" "X1230" "X1231" "X1232" "X1233"
## [1234] "X1234" "X1235" "X1236" "X1237" "X1238" "X1239" "X1240" "X1241" "X1242"
## [1243] "X1243" "X1244" "X1245" "X1246" "X1247" "X1248" "X1249" "X1250" "X1251"
## [1252] "X1252" "X1253" "X1254" "X1255" "X1256" "X1257" "X1258" "X1259" "X1260"
## [1261] "X1261" "X1262" "X1263" "X1264" "X1265" "X1266" "X1267" "X1268" "X1269"
## [1270] "X1270" "X1271" "X1272" "X1273" "X1274" "X1275" "X1276" "X1277" "X1278"
## [1279] "X1279" "X1280" "X1281" "X1282" "X1283" "X1284" "X1285" "X1286" "X1287"
## [1288] "X1288" "X1289" "X1290" "X1291" "X1292" "X1293" "X1294" "X1295" "X1296"
## [1297] "X1297" "X1298" "X1299" "X1300" "X1301" "X1302" "X1303" "X1304" "X1305"
## [1306] "X1306" "X1307" "X1308" "X1309" "X1310" "X1311" "X1312" "X1313" "X1314"
## [1315] "X1315" "X1316" "X1317" "X1318" "X1319" "X1320" "X1321" "X1322" "X1323"
## [1324] "X1324" "X1325" "X1326" "X1327" "X1328" "X1329" "X1330" "X1331" "X1332"
## [1333] "X1333" "X1334" "X1335" "X1336" "X1337" "X1338" "X1339" "X1340" "X1341"
## [1342] "X1342" "X1343" "X1344" "X1345" "X1346" "X1347" "X1348" "X1349" "X1350"
## [1351] "X1351" "X1352" "X1353" "X1354" "X1355" "X1356" "X1357" "X1358" "X1359"
## [1360] "X1360" "X1361" "X1362" "X1363" "X1364" "X1365" "X1366" "X1367" "X1368"
## [1369] "X1369" "X1370" "X1371" "X1372" "X1373" "X1374" "X1375" "X1376" "X1377"
## [1378] "X1378" "X1379" "X1380" "X1381" "X1382" "X1383" "X1384" "X1385" "X1386"
## [1387] "X1387" "X1388" "X1389" "X1390" "X1391" "X1392" "X1393" "X1394" "X1395"
## [1396] "X1396" "X1397" "X1398" "X1399" "X1400" "X1401" "X1402" "X1403" "X1404"
## [1405] "X1405" "X1406" "X1407" "X1408" "X1409" "X1410" "X1411" "X1412" "X1413"
## [1414] "X1414" "X1415" "X1416" "X1417" "X1418" "X1419" "X1420" "X1421" "X1422"
## [1423] "X1423" "X1424" "X1425" "X1426" "X1427" "X1428" "X1429" "X1430" "X1431"
## [1432] "X1432" "X1433" "X1434" "X1435" "X1436" "X1437" "X1438" "X1439" "X1440"
## [1441] "X1441" "X1442" "X1443" "X1444" "X1445" "X1446" "X1447" "X1448" "X1449"
## [1450] "X1450" "X1451" "X1452" "X1453" "X1454" "X1455" "X1456" "X1457" "X1458"
## [1459] "X1459" "X1460" "X1461" "X1462" "X1463" "X1464" "X1465" "X1466" "X1467"
## [1468] "X1468" "X1469" "X1470" "X1471" "X1472" "X1473" "X1474" "X1475" "X1476"
## [1477] "X1477" "X1478" "X1479" "X1480" "X1481" "X1482" "X1483" "X1484" "X1485"
## [1486] "X1486" "X1487" "X1488" "X1489" "X1490" "X1491" "X1492" "X1493" "X1494"
## [1495] "X1495" "X1496" "X1497" "X1498" "X1499" "X1500" "X1501" "X1502" "X1503"
## [1504] "X1504" "X1505" "X1506" "X1507" "X1508" "X1509" "X1510" "X1511" "X1512"
## [1513] "X1513" "X1514" "X1515" "X1516" "X1517" "X1518" "X1519" "X1520" "X1521"
## [1522] "X1522" "X1523" "X1524" "X1525" "X1526" "X1527" "X1528" "X1529" "X1530"
## [1531] "X1531" "X1532" "X1533" "X1534" "X1535" "X1536" "X1537" "X1538" "X1539"
## [1540] "X1540" "X1541" "X1542" "X1543" "X1544" "X1545" "X1546" "X1547" "X1548"
## [1549] "X1549" "X1550" "X1551" "X1552" "X1553" "X1554" "X1555" "X1556" "X1557"
## [1558] "X1558" "X1559" "X1560" "X1561" "X1562" "X1563" "X1564" "X1565" "X1566"
## [1567] "X1567" "X1568" "X1569" "X1570" "X1571" "X1572" "X1573" "X1574" "X1575"
## [1576] "X1576" "X1577" "X1578" "X1579" "X1580" "X1581" "X1582" "X1583" "X1584"
## [1585] "X1585" "X1586" "X1587" "X1588" "X1589" "X1590" "X1591" "X1592" "X1593"
## [1594] "X1594" "X1595" "X1596" "X1597" "X1598" "X1599" "X1600" "X1601" "X1602"
## [1603] "X1603" "X1604" "X1605" "X1606" "X1607" "X1608" "X1609" "X1610" "X1611"
## [1612] "X1612" "X1613" "X1614" "X1615" "X1616" "X1617" "X1618" "X1619" "X1620"
## [1621] "X1621" "X1622" "X1623" "X1624" "X1625" "X1626" "X1627" "X1628" "X1629"
## [1630] "X1630" "X1631" "X1632" "X1633" "X1634" "X1635" "X1636" "X1637" "X1638"
## [1639] "X1639" "X1640" "X1641" "X1642" "X1643" "X1644" "X1645" "X1646" "X1647"
## [1648] "X1648" "X1649" "X1650" "X1651" "X1652" "X1653" "X1654" "X1655" "X1656"
## [1657] "X1657" "X1658" "X1659" "X1660" "X1661" "X1662" "X1663" "X1664" "X1665"
## [1666] "X1666" "X1667" "X1668" "X1669" "X1670" "X1671" "X1672" "X1673" "X1674"
## [1675] "X1675" "X1676" "X1677" "X1678" "X1679" "X1680" "X1681" "X1682" "X1683"
## [1684] "X1684" "X1685" "X1686" "X1687" "X1688" "X1689" "X1690" "X1691" "X1692"
## [1693] "X1693" "X1694" "X1695" "X1696" "X1697" "X1698" "X1699" "X1700" "X1701"
## [1702] "X1702" "X1703" "X1704" "X1705" "X1706" "X1707" "X1708" "X1709" "X1710"
## [1711] "X1711" "X1712" "X1713" "X1714" "X1715" "X1716" "X1717" "X1718" "X1719"
## [1720] "X1720" "X1721" "X1722" "X1723" "X1724" "X1725" "X1726" "X1727" "X1728"
## [1729] "X1729" "X1730" "X1731" "X1732" "X1733" "X1734" "X1735" "X1736" "X1737"
## [1738] "X1738" "X1739" "X1740" "X1741" "X1742" "X1743" "X1744" "X1745" "X1746"
## [1747] "X1747" "X1748" "X1749" "X1750" "X1751" "X1752" "X1753" "X1754" "X1755"
## [1756] "X1756" "X1757" "X1758" "X1759" "X1760" "X1761" "X1762" "X1763" "X1764"
## [1765] "X1765" "X1766" "X1767" "X1768" "X1769" "X1770" "X1771" "X1772" "X1773"
## [1774] "X1774" "X1775" "X1776" "X1777" "X1778" "X1779" "X1780" "X1781" "X1782"
## [1783] "X1783" "X1784" "X1785" "X1786" "X1787" "X1788" "X1789" "X1790" "X1791"
## [1792] "X1792" "X1793" "X1794" "X1795" "X1796" "X1797" "X1798" "X1799" "X1800"
## [1801] "X1801" "X1802" "X1803" "X1804" "X1805" "X1806" "X1807" "X1808" "X1809"
## [1810] "X1810" "X1811" "X1812" "X1813" "X1814" "X1815" "X1816" "X1817" "X1818"
## [1819] "X1819" "X1820" "X1821" "X1822" "X1823" "X1824" "X1825" "X1826" "X1827"
## [1828] "X1828" "X1829" "X1830" "X1831" "X1832" "X1833" "X1834" "X1835" "X1836"
## [1837] "X1837" "X1838" "X1839" "X1840" "X1841" "X1842" "X1843" "X1844" "X1845"
## [1846] "X1846" "X1847" "X1848" "X1849" "X1850" "X1851" "X1852" "X1853" "X1854"
## [1855] "X1855" "X1856" "X1857" "X1858" "X1859" "X1860" "X1861" "X1862" "X1863"
## [1864] "X1864" "X1865" "X1866" "X1867" "X1868" "X1869" "X1870" "X1871" "X1872"
## [1873] "X1873" "X1874" "X1875" "X1876" "X1877" "X1878" "X1879" "X1880" "X1881"
## [1882] "X1882" "X1883" "X1884" "X1885" "X1886" "X1887" "X1888" "X1889" "X1890"
## [1891] "X1891" "X1892" "X1893" "X1894" "X1895" "X1896" "X1897" "X1898" "X1899"
## [1900] "X1900" "X1901" "X1902" "X1903" "X1904" "X1905" "X1906" "X1907" "X1908"
## [1909] "X1909" "X1910" "X1911" "X1912" "X1913" "X1914" "X1915" "X1916" "X1917"
## [1918] "X1918" "X1919" "X1920" "X1921" "X1922" "X1923" "X1924" "X1925" "X1926"
## [1927] "X1927" "X1928" "X1929" "X1930" "X1931" "X1932" "X1933" "X1934" "X1935"
## [1936] "X1936" "X1937" "X1938" "X1939" "X1940" "X1941" "X1942" "X1943" "X1944"
## [1945] "X1945" "X1946" "X1947" "X1948" "X1949" "X1950" "X1951" "X1952" "X1953"
## [1954] "X1954" "X1955" "X1956" "X1957" "X1958" "X1959" "X1960" "X1961" "X1962"
## [1963] "X1963" "X1964" "X1965" "X1966" "X1967" "X1968" "X1969" "X1970" "X1971"
## [1972] "X1972" "X1973" "X1974" "X1975" "X1976" "X1977" "X1978" "X1979" "X1980"
## [1981] "X1981" "X1982" "X1983" "X1984" "X1985" "X1986" "X1987" "X1988" "X1989"
## [1990] "X1990" "X1991" "X1992" "X1993" "X1994" "X1995" "X1996" "X1997" "X1998"
## [1999] "X1999" "X2000" "X2001" "X2002" "X2003" "X2004" "X2005" "X2006" "X2007"
## [2008] "X2008" "X2009" "X2010" "X2011" "X2012" "X2013" "X2014" "X2015" "X2016"
## [2017] "X2017" "X2018" "X2019" "X2020" "X2021" "X2022" "X2023" "X2024" "X2025"
## [2026] "X2026" "X2027" "X2028" "X2029" "X2030" "X2031" "X2032" "X2033" "X2034"
## [2035] "X2035" "X2036" "X2037" "X2038" "X2039" "X2040" "X2041" "X2042" "X2043"
## [2044] "X2044" "X2045" "X2046" "X2047" "X2048" "X2049" "X2050" "X2051" "X2052"
## [2053] "X2053" "X2054" "X2055" "X2056" "X2057" "X2058" "X2059" "X2060" "X2061"
## [2062] "X2062" "X2063" "X2064" "X2065" "X2066" "X2067" "X2068" "X2069" "X2070"
## [2071] "X2071" "X2072" "X2073" "X2074" "X2075" "X2076" "X2077" "X2078" "X2079"
## [2080] "X2080" "X2081" "X2082" "X2083" "X2084" "X2085" "X2086" "X2087" "X2088"
## [2089] "X2089" "X2090" "X2091" "X2092" "X2093" "X2094" "X2095" "X2096" "X2097"
## [2098] "X2098" "X2099" "X2100" "X2101" "X2102" "X2103" "X2104" "X2105" "X2106"
## [2107] "X2107" "X2108" "X2109" "X2110" "X2111" "X2112" "X2113" "X2114" "X2115"
## [2116] "X2116" "X2117" "X2118" "X2119" "X2120" "X2121" "X2122" "X2123" "X2124"
## [2125] "X2125" "X2126" "X2127" "X2128" "X2129" "X2130" "X2131" "X2132" "X2133"
## [2134] "X2134" "X2135" "X2136" "X2137" "X2138" "X2139" "X2140" "X2141" "X2142"
## [2143] "X2143" "X2144" "X2145" "X2146" "X2147" "X2148" "X2149" "X2150" "X2151"
## [2152] "X2152" "X2153" "X2154" "X2155" "X2156" "X2157" "X2158" "X2159" "X2160"
## [2161] "X2161" "X2162" "X2163" "X2164" "X2165" "X2166" "X2167" "X2168" "X2169"
## [2170] "X2170" "X2171" "X2172" "X2173" "X2174" "X2175" "X2176" "X2177" "X2178"
## [2179] "X2179" "X2180" "X2181" "X2182" "X2183" "X2184" "X2185" "X2186" "X2187"
## [2188] "X2188" "X2189" "X2190" "X2191" "X2192" "X2193" "X2194" "X2195" "X2196"
## [2197] "X2197" "X2198" "X2199" "X2200" "X2201" "X2202" "X2203" "X2204" "X2205"
## [2206] "X2206" "X2207" "X2208" "X2209" "X2210" "X2211" "X2212" "X2213" "X2214"
## [2215] "X2215" "X2216" "X2217" "X2218" "X2219" "X2220" "X2221" "X2222" "X2223"
## [2224] "X2224" "X2225" "X2226" "X2227" "X2228" "X2229" "X2230" "X2231" "X2232"
## [2233] "X2233" "X2234" "X2235" "X2236" "X2237" "X2238" "X2239" "X2240" "X2241"
## [2242] "X2242" "X2243" "X2244" "X2245" "X2246" "X2247" "X2248" "X2249" "X2250"
## [2251] "X2251" "X2252" "X2253" "X2254" "X2255" "X2256" "X2257" "X2258" "X2259"
## [2260] "X2260" "X2261" "X2262" "X2263" "X2264" "X2265" "X2266" "X2267" "X2268"
## [2269] "X2269" "X2270" "X2271" "X2272" "X2273" "X2274" "X2275" "X2276" "X2277"
## [2278] "X2278" "X2279" "X2280" "X2281" "X2282" "X2283" "X2284" "X2285" "X2286"
## [2287] "X2287" "X2288" "X2289" "X2290" "X2291" "X2292" "X2293" "X2294" "X2295"
## [2296] "X2296" "X2297" "X2298" "X2299" "X2300" "X2301" "X2302" "X2303" "X2304"
## [2305] "X2305" "X2306" "X2307" "X2308" "X2309" "X2310" "X2311" "X2312" "X2313"
## [2314] "X2314" "X2315" "X2316" "X2317" "X2318" "X2319" "X2320" "X2321" "X2322"
## [2323] "X2323" "X2324" "X2325" "X2326" "X2327" "X2328" "X2329" "X2330" "X2331"
## [2332] "X2332" "X2333" "X2334" "X2335" "X2336" "X2337" "X2338" "X2339" "X2340"
## [2341] "X2341" "X2342" "X2343" "X2344" "X2345" "X2346" "X2347" "X2348" "X2349"
## [2350] "X2350" "X2351" "X2352" "X2353" "X2354" "X2355" "X2356" "X2357" "X2358"
## [2359] "X2359" "X2360" "X2361" "X2362" "X2363" "X2364" "X2365" "X2366" "X2367"
## [2368] "X2368" "X2369" "X2370" "X2371" "X2372" "X2373" "X2374" "X2375" "X2376"
## [2377] "X2377" "X2378" "X2379" "X2380" "X2381" "X2382" "X2383" "X2384" "X2385"
## [2386] "X2386" "X2387" "X2388" "X2389" "X2390" "X2391" "X2392" "X2393" "X2394"
## [2395] "X2395" "X2396" "X2397" "X2398" "X2399" "X2400" "X2401" "X2402" "X2403"
## [2404] "X2404" "X2405" "X2406" "X2407" "X2408" "X2409" "X2410" "X2411" "X2412"
## [2413] "X2413" "X2414" "X2415" "X2416" "X2417" "X2418" "X2419" "X2420" "X2421"
## [2422] "X2422" "X2423" "X2424" "X2425" "X2426" "X2427" "X2428" "X2429" "X2430"
## [2431] "X2431" "X2432" "X2433" "X2434" "X2435" "X2436" "X2437" "X2438" "X2439"
## [2440] "X2440" "X2441" "X2442" "X2443" "X2444" "X2445" "X2446" "X2447" "X2448"
## [2449] "X2449" "X2450" "X2451" "X2452" "X2453" "X2454" "X2455" "X2456" "X2457"
## [2458] "X2458" "X2459" "X2460" "X2461" "X2462" "X2463" "X2464" "X2465" "X2466"
## [2467] "X2467" "X2468" "X2469" "X2470" "X2471" "X2472" "X2473" "X2474" "X2475"
## [2476] "X2476" "X2477" "X2478" "X2479" "X2480" "X2481" "X2482" "X2483" "X2484"
## [2485] "X2485" "X2486" "X2487" "X2488" "X2489" "X2490" "X2491" "X2492" "X2493"
## [2494] "X2494" "X2495" "X2496" "X2497" "X2498" "X2499" "X2500" "X2501" "X2502"
## [2503] "X2503" "X2504" "X2505" "X2506" "X2507" "X2508" "X2509" "X2510" "X2511"
## [2512] "X2512" "X2513" "X2514" "X2515" "X2516" "X2517" "X2518" "X2519" "X2520"
## [2521] "X2521" "X2522" "X2523" "X2524" "X2525" "X2526" "X2527" "X2528" "X2529"
## [2530] "X2530" "X2531" "X2532" "X2533" "X2534" "X2535" "X2536" "X2537" "X2538"
## [2539] "X2539" "X2540" "X2541" "X2542" "X2543" "X2544" "X2545" "X2546" "X2547"
## [2548] "X2548" "X2549" "X2550" "X2551" "X2552" "X2553" "X2554" "X2555" "X2556"
## [2557] "X2557" "X2558" "X2559" "X2560" "X2561" "X2562" "X2563" "X2564" "X2565"
## [2566] "X2566" "X2567" "X2568" "X2569" "X2570" "X2571" "X2572" "X2573" "X2574"
## [2575] "X2575" "X2576" "X2577" "X2578" "X2579" "X2580" "X2581" "X2582" "X2583"
## [2584] "X2584" "X2585" "X2586" "X2587" "X2588" "X2589" "X2590" "X2591" "X2592"
## [2593] "X2593" "X2594" "X2595" "X2596" "X2597" "X2598" "X2599" "X2600" "X2601"
## [2602] "X2602" "X2603" "X2604" "X2605" "X2606" "X2607" "X2608" "X2609" "X2610"
## [2611] "X2611" "X2612" "X2613" "X2614" "X2615" "X2616" "X2617" "X2618" "X2619"
## [2620] "X2620" "X2621" "X2622" "X2623" "X2624" "X2625" "X2626" "X2627" "X2628"
## [2629] "X2629" "X2630" "X2631" "X2632" "X2633" "X2634" "X2635" "X2636" "X2637"
## [2638] "X2638" "X2639" "X2640" "X2641" "X2642" "X2643" "X2644" "X2645" "X2646"
## [2647] "X2647" "X2648" "X2649" "X2650" "X2651" "X2652" "X2653" "X2654" "X2655"
## [2656] "X2656" "X2657" "X2658" "X2659" "X2660" "X2661" "X2662" "X2663" "X2664"
## [2665] "X2665" "X2666" "X2667" "X2668" "X2669" "X2670" "X2671" "X2672" "X2673"
## [2674] "X2674" "X2675" "X2676" "X2677" "X2678" "X2679" "X2680" "X2681" "X2682"
## [2683] "X2683" "X2684" "X2685" "X2686" "X2687" "X2688" "X2689" "X2690" "X2691"
## [2692] "X2692" "X2693" "X2694" "X2695" "X2696" "X2697" "X2698" "X2699" "X2700"
## [2701] "X2701" "X2702" "X2703" "X2704" "X2705" "X2706" "X2707" "X2708" "X2709"
## [2710] "X2710" "X2711" "X2712" "X2713" "X2714" "X2715" "X2716" "X2717" "X2718"
## [2719] "X2719" "X2720" "X2721" "X2722" "X2723" "X2724" "X2725" "X2726" "X2727"
## [2728] "X2728" "X2729" "X2730" "X2731" "X2732" "X2733" "X2734" "X2735" "X2736"
## [2737] "X2737" "X2738" "X2739" "X2740" "X2741" "X2742" "X2743" "X2744" "X2745"
## [2746] "X2746" "X2747" "X2748" "X2749" "X2750" "X2751" "X2752" "X2753" "X2754"
## [2755] "X2755" "X2756" "X2757" "X2758" "X2759" "X2760" "X2761" "X2762" "X2763"
## [2764] "X2764" "X2765" "X2766" "X2767" "X2768" "X2769" "X2770" "X2771" "X2772"
## [2773] "X2773" "X2774" "X2775" "X2776" "X2777" "X2778" "X2779" "X2780" "X2781"
## [2782] "X2782" "X2783" "X2784" "X2785" "X2786" "X2787" "X2788" "X2789" "X2790"
## [2791] "X2791" "X2792" "X2793" "X2794" "X2795" "X2796" "X2797" "X2798" "X2799"
## [2800] "X2800" "X2801" "X2802" "X2803" "X2804" "X2805" "X2806" "X2807" "X2808"
## [2809] "X2809" "X2810" "X2811" "X2812" "X2813" "X2814" "X2815" "X2816" "X2817"
## [2818] "X2818" "X2819" "X2820" "X2821" "X2822" "X2823" "X2824" "X2825" "X2826"
## [2827] "X2827" "X2828" "X2829" "X2830" "X2831" "X2832" "X2833" "X2834" "X2835"
## [2836] "X2836" "X2837" "X2838" "X2839" "X2840" "X2841" "X2842" "X2843" "X2844"
## [2845] "X2845" "X2846" "X2847" "X2848" "X2849" "X2850" "X2851" "X2852" "X2853"
## [2854] "X2854" "X2855" "X2856" "X2857" "X2858" "X2859" "X2860" "X2861" "X2862"
## [2863] "X2863" "X2864" "X2865" "X2866" "X2867" "X2868" "X2869" "X2870" "X2871"
## [2872] "X2872" "X2873" "X2874" "X2875" "X2876" "X2877" "X2878" "X2879" "X2880"
## [2881] "X2881" "X2882" "X2883" "X2884" "X2885" "X2886" "X2887" "X2888" "X2889"
## [2890] "X2890" "X2891" "X2892" "X2893" "X2894" "X2895" "X2896" "X2897" "X2898"
## [2899] "X2899" "X2900" "X2901" "X2902" "X2903" "X2904" "X2905" "X2906" "X2907"
## [2908] "X2908" "X2909" "X2910" "X2911" "X2912" "X2913" "X2914" "X2915" "X2916"
## [2917] "X2917" "X2918" "X2919" "X2920" "X2921" "X2922" "X2923" "X2924" "X2925"
## [2926] "X2926" "X2927" "X2928" "X2929" "X2930" "X2931" "X2932" "X2933" "X2934"
## [2935] "X2935" "X2936" "X2937" "X2938" "X2939" "X2940" "X2941" "X2942" "X2943"
## [2944] "X2944" "X2945" "X2946" "X2947" "X2948" "X2949" "X2950" "X2951" "X2952"
## [2953] "X2953" "X2954" "X2955" "X2956" "X2957" "X2958" "X2959" "X2960" "X2961"
## [2962] "X2962" "X2963" "X2964" "X2965" "X2966" "X2967" "X2968" "X2969" "X2970"
## [2971] "X2971" "X2972" "X2973" "X2974" "X2975" "X2976" "X2977" "X2978" "X2979"
## [2980] "X2980" "X2981" "X2982" "X2983" "X2984" "X2985" "X2986" "X2987" "X2988"
## [2989] "X2989" "X2990" "X2991" "X2992" "X2993" "X2994" "X2995" "X2996" "X2997"
## [2998] "X2998" "X2999" "X3000" "X3001" "X3002" "X3003" "X3004" "X3005" "X3006"
## [3007] "X3007" "X3008" "X3009" "X3010" "X3011" "X3012" "X3013" "X3014" "X3015"
## [3016] "X3016" "X3017" "X3018" "X3019" "X3020" "X3021" "X3022" "X3023" "X3024"
## [3025] "X3025" "X3026" "X3027" "X3028" "X3029" "X3030" "X3031" "X3032" "X3033"
## [3034] "X3034" "X3035" "X3036" "X3037" "X3038" "X3039" "X3040" "X3041" "X3042"
## [3043] "X3043" "X3044" "X3045" "X3046" "X3047" "X3048" "X3049" "X3050" "X3051"
## [3052] "X3052" "X3053" "X3054" "X3055" "X3056" "X3057" "X3058" "X3059" "X3060"
## [3061] "X3061" "X3062" "X3063" "X3064" "X3065" "X3066" "X3067" "X3068" "X3069"
## [3070] "X3070" "X3071" "X3072" "X3073" "X3074" "X3075" "X3076" "X3077" "X3078"
## [3079] "X3079" "X3080" "X3081" "X3082" "X3083" "X3084" "X3085" "X3086" "X3087"
## [3088] "X3088" "X3089" "X3090" "X3091" "X3092" "X3093" "X3094" "X3095" "X3096"
## [3097] "X3097" "X3098" "X3099" "X3100" "X3101" "X3102" "X3103" "X3104" "X3105"
## [3106] "X3106" "X3107" "X3108" "X3109" "X3110" "X3111" "X3112" "X3113" "X3114"
## [3115] "X3115" "X3116" "X3117" "X3118" "X3119" "X3120" "X3121" "X3122" "X3123"
## [3124] "X3124" "X3125" "X3126" "X3127" "X3128" "X3129" "X3130" "X3131" "X3132"
## [3133] "X3133" "X3134" "X3135" "X3136" "X3137" "X3138" "X3139" "X3140" "X3141"
## [3142] "X3142" "X3143" "X3144" "X3145" "X3146" "X3147" "X3148" "X3149" "X3150"
## [3151] "X3151" "X3152" "X3153" "X3154" "X3155" "X3156" "X3157" "X3158" "X3159"
## [3160] "X3160" "X3161" "X3162" "X3163" "X3164" "X3165" "X3166" "X3167" "X3168"
## [3169] "X3169" "X3170" "X3171" "X3172" "X3173" "X3174" "X3175" "X3176" "X3177"
## [3178] "X3178" "X3179" "X3180" "X3181" "X3182" "X3183" "X3184" "X3185" "X3186"
## [3187] "X3187" "X3188" "X3189" "X3190" "X3191" "X3192" "X3193" "X3194" "X3195"
## [3196] "X3196" "X3197" "X3198" "X3199" "X3200" "X3201" "X3202" "X3203" "X3204"
## [3205] "X3205" "X3206" "X3207" "X3208" "X3209" "X3210" "X3211" "X3212" "X3213"
## [3214] "X3214" "X3215" "X3216" "X3217" "X3218" "X3219" "X3220" "X3221" "X3222"
## [3223] "X3223" "X3224" "X3225" "X3226" "X3227" "X3228" "X3229" "X3230" "X3231"
## [3232] "X3232" "X3233" "X3234" "X3235" "X3236" "X3237" "X3238" "X3239" "X3240"
## [3241] "X3241" "X3242" "X3243" "X3244" "X3245" "X3246" "X3247" "X3248" "X3249"
## [3250] "X3250" "X3251" "X3252" "X3253" "X3254" "X3255" "X3256" "X3257" "X3258"
## [3259] "X3259" "X3260" "X3261" "X3262" "X3263" "X3264" "X3265" "X3266" "X3267"
## [3268] "X3268" "X3269" "X3270" "X3271" "X3272" "X3273" "X3274" "X3275" "X3276"
## [3277] "X3277" "X3278" "X3279" "X3280" "X3281" "X3282" "X3283" "X3284" "X3285"
## [3286] "X3286" "X3287" "X3288" "X3289" "X3290" "X3291" "X3292" "X3293" "X3294"
## [3295] "X3295" "X3296" "X3297" "X3298" "X3299" "X3300" "X3301" "X3302" "X3303"
## [3304] "X3304" "X3305" "X3306" "X3307" "X3308" "X3309" "X3310" "X3311" "X3312"
## [3313] "X3313" "X3314" "X3315" "X3316" "X3317" "X3318" "X3319" "X3320" "X3321"
## [3322] "X3322" "X3323" "X3324" "X3325" "X3326" "X3327" "X3328" "X3329" "X3330"
## [3331] "X3331" "X3332" "X3333" "X3334" "X3335" "X3336" "X3337" "X3338" "X3339"
## [3340] "X3340" "X3341" "X3342" "X3343" "X3344" "X3345" "X3346" "X3347" "X3348"
## [3349] "X3349" "X3350" "X3351" "X3352" "X3353" "X3354" "X3355" "X3356" "X3357"
## [3358] "X3358" "X3359" "X3360" "X3361" "X3362" "X3363" "X3364" "X3365" "X3366"
## [3367] "X3367" "X3368" "X3369" "X3370" "X3371" "X3372" "X3373" "X3374" "X3375"
## [3376] "X3376" "X3377" "X3378" "X3379" "X3380" "X3381" "X3382" "X3383" "X3384"
## [3385] "X3385" "X3386" "X3387" "X3388" "X3389" "X3390" "X3391" "X3392" "X3393"
## [3394] "X3394" "X3395" "X3396" "X3397" "X3398" "X3399" "X3400" "X3401" "X3402"
## [3403] "X3403" "X3404" "X3405" "X3406" "X3407" "X3408" "X3409" "X3410" "X3411"
## [3412] "X3412" "X3413" "X3414" "X3415" "X3416" "X3417" "X3418" "X3419" "X3420"
## [3421] "X3421" "X3422" "X3423" "X3424" "X3425" "X3426" "X3427" "X3428" "X3429"
## [3430] "X3430" "X3431" "X3432" "X3433" "X3434" "X3435" "X3436" "X3437" "X3438"
## [3439] "X3439" "X3440" "X3441" "X3442" "X3443" "X3444" "X3445" "X3446" "X3447"
## [3448] "X3448" "X3449" "X3450" "X3451" "X3452" "X3453" "X3454" "X3455" "X3456"
## [3457] "X3457" "X3458" "X3459" "X3460" "X3461" "X3462" "X3463" "X3464" "X3465"
## [3466] "X3466" "X3467" "X3468" "X3469" "X3470" "X3471" "X3472" "X3473" "X3474"
## [3475] "X3475" "X3476" "X3477" "X3478" "X3479" "X3480" "X3481" "X3482" "X3483"
## [3484] "X3484" "X3485" "X3486" "X3487" "X3488" "X3489" "X3490" "X3491" "X3492"
## [3493] "X3493" "X3494" "X3495" "X3496" "X3497" "X3498" "X3499" "X3500" "X3501"
## [3502] "X3502" "X3503" "X3504" "X3505" "X3506" "X3507" "X3508" "X3509" "X3510"
## [3511] "X3511" "X3512" "X3513" "X3514" "X3515" "X3516" "X3517" "X3518" "X3519"
## [3520] "X3520" "X3521" "X3522" "X3523" "X3524" "X3525" "X3526" "X3527" "X3528"
## [3529] "X3529" "X3530" "X3531" "X3532" "X3533" "X3534" "X3535" "X3536" "X3537"
## [3538] "X3538" "X3539" "X3540" "X3541" "X3542" "X3543" "X3544" "X3545" "X3546"
## [3547] "X3547" "X3548" "X3549" "X3550" "X3551" "X3552" "X3553" "X3554" "X3555"
## [3556] "X3556" "X3557" "X3558" "X3559" "X3560" "X3561" "X3562" "X3563" "X3564"
## [3565] "X3565" "X3566" "X3567" "X3568" "X3569" "X3570" "X3571" "X3572" "X3573"
## [3574] "X3574" "X3575" "X3576" "X3577" "X3578" "X3579" "X3580" "X3581" "X3582"
## [3583] "X3583" "X3584" "X3585" "X3586" "X3587" "X3588" "X3589" "X3590" "X3591"
## [3592] "X3592" "X3593" "X3594" "X3595" "X3596" "X3597" "X3598" "X3599" "X3600"
## [3601] "X3601" "X3602" "X3603" "X3604" "X3605" "X3606" "X3607" "X3608" "X3609"
## [3610] "X3610" "X3611" "X3612" "X3613" "X3614" "X3615" "X3616" "X3617" "X3618"
## [3619] "X3619" "X3620" "X3621" "X3622" "X3623" "X3624" "X3625" "X3626" "X3627"
## [3628] "X3628" "X3629" "X3630" "X3631" "X3632" "X3633" "X3634" "X3635" "X3636"
## [3637] "X3637" "X3638" "X3639" "X3640" "X3641" "X3642" "X3643" "X3644" "X3645"
## [3646] "X3646" "X3647" "X3648" "X3649" "X3650" "X3651" "X3652" "X3653" "X3654"
## [3655] "X3655" "X3656" "X3657" "X3658" "X3659" "X3660" "X3661" "X3662" "X3663"
## [3664] "X3664" "X3665" "X3666" "X3667" "X3668" "X3669" "X3670" "X3671" "X3672"
## [3673] "X3673" "X3674" "X3675" "X3676" "X3677" "X3678" "X3679" "X3680" "X3681"
## [3682] "X3682" "X3683" "X3684" "X3685" "X3686" "X3687" "X3688" "X3689" "X3690"
## [3691] "X3691" "X3692" "X3693" "X3694" "X3695" "X3696" "X3697" "X3698" "X3699"
## [3700] "X3700" "X3701" "X3702" "X3703" "X3704" "X3705" "X3706" "X3707" "X3708"
## [3709] "X3709" "X3710" "X3711" "X3712" "X3713" "X3714" "X3715" "X3716" "X3717"
## [3718] "X3718" "X3719" "X3720" "X3721" "X3722" "X3723" "X3724" "X3725" "X3726"
## [3727] "X3727" "X3728" "X3729" "X3730" "X3731" "X3732" "X3733" "X3734" "X3735"
## [3736] "X3736" "X3737" "X3738" "X3739" "X3740" "X3741" "X3742" "X3743" "X3744"
## [3745] "X3745" "X3746" "X3747" "X3748" "X3749" "X3750" "X3751" "X3752" "X3753"
## [3754] "X3754" "X3755" "X3756" "X3757" "X3758" "X3759" "X3760" "X3761" "X3762"
## [3763] "X3763" "X3764" "X3765" "X3766" "X3767" "X3768" "X3769" "X3770" "X3771"
## [3772] "X3772" "X3773" "X3774" "X3775" "X3776" "X3777" "X3778" "X3779" "X3780"
## [3781] "X3781" "X3782" "X3783" "X3784" "X3785" "X3786" "X3787" "X3788" "X3789"
## [3790] "X3790" "X3791" "X3792" "X3793" "X3794" "X3795" "X3796" "X3797" "X3798"
## [3799] "X3799" "X3800" "X3801" "X3802" "X3803" "X3804" "X3805" "X3806" "X3807"
## [3808] "X3808" "X3809" "X3810" "X3811" "X3812" "X3813" "X3814" "X3815" "X3816"
## [3817] "X3817" "X3818" "X3819" "X3820" "X3821" "X3822" "X3823" "X3824" "X3825"
## [3826] "X3826" "X3827" "X3828" "X3829" "X3830" "X3831" "X3832" "X3833" "X3834"
## [3835] "X3835" "X3836" "X3837" "X3838" "X3839" "X3840" "X3841" "X3842" "X3843"
## [3844] "X3844" "X3845" "X3846" "X3847" "X3848" "X3849" "X3850" "X3851" "X3852"
## [3853] "X3853" "X3854" "X3855" "X3856" "X3857" "X3858" "X3859" "X3860" "X3861"
## [3862] "X3862" "X3863" "X3864" "X3865" "X3866" "X3867" "X3868" "X3869" "X3870"
## [3871] "X3871" "X3872" "X3873" "X3874" "X3875" "X3876" "X3877" "X3878" "X3879"
## [3880] "X3880" "X3881" "X3882" "X3883" "X3884" "X3885" "X3886" "X3887" "X3888"
## [3889] "X3889" "X3890" "X3891" "X3892" "X3893" "X3894" "X3895" "X3896" "X3897"
## [3898] "X3898" "X3899" "X3900" "X3901" "X3902" "X3903" "X3904" "X3905" "X3906"
## [3907] "X3907" "X3908" "X3909" "X3910" "X3911" "X3912" "X3913" "X3914" "X3915"
## [3916] "X3916" "X3917" "X3918" "X3919" "X3920" "X3921" "X3922" "X3923" "X3924"
## [3925] "X3925" "X3926" "X3927" "X3928" "X3929" "X3930" "X3931" "X3932" "X3933"
## [3934] "X3934" "X3935" "X3936" "X3937" "X3938" "X3939" "X3940" "X3941" "X3942"
## [3943] "X3943" "X3944" "X3945" "X3946" "X3947" "X3948" "X3949" "X3950" "X3951"
## [3952] "X3952" "X3953" "X3954" "X3955" "X3956" "X3957" "X3958" "X3959" "X3960"
## [3961] "X3961" "X3962" "X3963" "X3964" "X3965" "X3966" "X3967" "X3968" "X3969"
## [3970] "X3970" "X3971" "X3972" "X3973" "X3974" "X3975" "X3976" "X3977" "X3978"
## [3979] "X3979" "X3980" "X3981" "X3982" "X3983" "X3984" "X3985" "X3986" "X3987"
## [3988] "X3988" "X3989" "X3990" "X3991" "X3992" "X3993" "X3994" "X3995" "X3996"
## [3997] "X3997" "X3998" "X3999" "X4000" "X4001" "X4002" "X4003" "X4004" "X4005"
## [4006] "X4006" "X4007" "X4008" "X4009" "X4010" "X4011" "X4012" "X4013" "X4014"
## [4015] "X4015" "X4016" "X4017" "X4018" "X4019" "X4020" "X4021" "X4022" "X4023"
## [4024] "X4024" "X4025" "X4026" "X4027" "X4028" "X4029" "X4030" "X4031" "X4032"
## [4033] "X4033" "X4034" "X4035" "X4036" "X4037" "X4038" "X4039" "X4040" "X4041"
## [4042] "X4042" "X4043" "X4044" "X4045" "X4046" "X4047" "X4048" "X4049" "X4050"
## [4051] "X4051" "X4052" "X4053" "X4054" "X4055" "X4056" "X4057" "X4058" "X4059"
## [4060] "X4060" "X4061" "X4062" "X4063" "X4064" "X4065" "X4066" "X4067" "X4068"
## [4069] "X4069" "X4070" "X4071" "X4072" "X4073" "X4074" "X4075" "X4076" "X4077"
## [4078] "X4078" "X4079" "X4080" "X4081" "X4082" "X4083" "X4084" "X4085" "X4086"
## [4087] "X4087" "X4088" "X4089" "X4090" "X4091" "X4092" "X4093" "X4094" "X4095"
## [4096] "X4096" "X4097" "X4098" "X4099" "X4100" "X4101" "X4102" "X4103" "X4104"
## [4105] "X4105" "X4106" "X4107" "X4108" "X4109" "X4110" "X4111" "X4112" "X4113"
## [4114] "X4114" "X4115" "X4116" "X4117" "X4118" "X4119" "X4120" "X4121" "X4122"
## [4123] "X4123" "X4124" "X4125" "X4126" "X4127" "X4128" "X4129" "X4130" "X4131"
## [4132] "X4132" "X4133" "X4134" "X4135" "X4136" "X4137" "X4138" "X4139" "X4140"
## [4141] "X4141" "X4142" "X4143" "X4144" "X4145" "X4146" "X4147" "X4148" "X4149"
## [4150] "X4150" "X4151" "X4152" "X4153" "X4154" "X4155" "X4156" "X4157" "X4158"
## [4159] "X4159" "X4160" "X4161" "X4162" "X4163" "X4164" "X4165" "X4166" "X4167"
## [4168] "X4168" "X4169" "X4170" "X4171" "X4172" "X4173" "X4174" "X4175" "X4176"
## [4177] "X4177" "X4178" "X4179" "X4180" "X4181" "X4182" "X4183" "X4184" "X4185"
## [4186] "X4186" "X4187" "X4188" "X4189" "X4190" "X4191" "X4192" "X4193" "X4194"
## [4195] "X4195" "X4196" "X4197" "X4198" "X4199" "X4200" "X4201" "X4202" "X4203"
## [4204] "X4204" "X4205" "X4206" "X4207" "X4208" "X4209" "X4210" "X4211" "X4212"
## [4213] "X4213" "X4214" "X4215" "X4216" "X4217" "X4218" "X4219" "X4220" "X4221"
## [4222] "X4222" "X4223" "X4224" "X4225" "X4226" "X4227" "X4228" "X4229" "X4230"
## [4231] "X4231" "X4232" "X4233" "X4234" "X4235" "X4236" "X4237" "X4238" "X4239"
## [4240] "X4240" "X4241" "X4242" "X4243" "X4244" "X4245" "X4246" "X4247" "X4248"
## [4249] "X4249" "X4250" "X4251" "X4252" "X4253" "X4254" "X4255" "X4256" "X4257"
## [4258] "X4258" "X4259" "X4260" "X4261" "X4262" "X4263" "X4264" "X4265" "X4266"
## [4267] "X4267" "X4268" "X4269" "X4270" "X4271" "X4272" "X4273" "X4274" "X4275"
## [4276] "X4276" "X4277" "X4278" "X4279" "X4280" "X4281" "X4282" "X4283" "X4284"
## [4285] "X4285" "X4286" "X4287" "X4288" "X4289" "X4290" "X4291" "X4292" "X4293"
## [4294] "X4294" "X4295" "X4296" "X4297" "X4298" "X4299" "X4300" "X4301" "X4302"
## [4303] "X4303" "X4304" "X4305" "X4306" "X4307" "X4308" "X4309" "X4310" "X4311"
## [4312] "X4312" "X4313" "X4314" "X4315" "X4316" "X4317" "X4318" "X4319" "X4320"
## [4321] "X4321" "X4322" "X4323" "X4324" "X4325" "X4326" "X4327" "X4328" "X4329"
## [4330] "X4330" "X4331" "X4332" "X4333" "X4334" "X4335" "X4336" "X4337" "X4338"
## [4339] "X4339" "X4340" "X4341" "X4342" "X4343" "X4344" "X4345" "X4346" "X4347"
## [4348] "X4348" "X4349" "X4350" "X4351" "X4352" "X4353" "X4354" "X4355" "X4356"
## [4357] "X4357" "X4358" "X4359" "X4360" "X4361" "X4362" "X4363" "X4364" "X4365"
## [4366] "X4366" "X4367" "X4368" "X4369" "X4370" "X4371" "X4372" "X4373" "X4374"
## [4375] "X4375" "X4376" "X4377" "X4378" "X4379" "X4380" "X4381" "X4382" "X4383"
## [4384] "X4384" "X4385" "X4386" "X4387" "X4388" "X4389" "X4390" "X4391" "X4392"
## [4393] "X4393" "X4394" "X4395" "X4396" "X4397" "X4398" "X4399" "X4400" "X4401"
## [4402] "X4402" "X4403" "X4404" "X4405" "X4406" "X4407" "X4408" "X4409" "X4410"
## [4411] "X4411" "X4412" "X4413" "X4414" "X4415" "X4416" "X4417" "X4418" "X4419"
## [4420] "X4420" "X4421" "X4422" "X4423" "X4424" "X4425" "X4426" "X4427" "X4428"
## [4429] "X4429" "X4430" "X4431" "X4432" "X4433" "X4434" "X4435" "X4436" "X4437"
## [4438] "X4438" "X4439" "X4440" "X4441" "X4442" "X4443" "X4444" "X4445" "X4446"
## [4447] "X4447" "X4448" "X4449" "X4450" "X4451" "X4452" "X4453" "X4454" "X4455"
## [4456] "X4456" "X4457" "X4458" "X4459" "X4460" "X4461" "X4462" "X4463" "X4464"
## [4465] "X4465" "X4466" "X4467" "X4468" "X4469" "X4470" "X4471" "X4472" "X4473"
## [4474] "X4474" "X4475" "X4476" "X4477" "X4478" "X4479" "X4480" "X4481" "X4482"
## [4483] "X4483" "X4484" "X4485" "X4486" "X4487" "X4488" "X4489" "X4490" "X4491"
## [4492] "X4492" "X4493" "X4494" "X4495" "X4496" "X4497" "X4498" "X4499" "X4500"
## [4501] "X4501" "X4502" "X4503" "X4504" "X4505" "X4506" "X4507" "X4508" "X4509"
## [4510] "X4510" "X4511" "X4512" "X4513" "X4514" "X4515" "X4516" "X4517" "X4518"
## [4519] "X4519" "X4520" "X4521" "X4522" "X4523" "X4524" "X4525" "X4526" "X4527"
## [4528] "X4528" "X4529" "X4530" "X4531" "X4532" "X4533" "X4534" "X4535" "X4536"
## [4537] "X4537" "X4538" "X4539" "X4540" "X4541" "X4542" "X4543" "X4544" "X4545"
## [4546] "X4546" "X4547" "X4548" "X4549" "X4550" "X4551" "X4552" "X4553" "X4554"
## [4555] "X4555" "X4556" "X4557" "X4558" "X4559" "X4560" "X4561" "X4562" "X4563"
## [4564] "X4564" "X4565" "X4566" "X4567" "X4568" "X4569" "X4570" "X4571" "X4572"
## [4573] "X4573" "X4574" "X4575" "X4576" "X4577" "X4578" "X4579" "X4580" "X4581"
## [4582] "X4582" "X4583" "X4584" "X4585" "X4586" "X4587" "X4588" "X4589" "X4590"
## [4591] "X4591" "X4592" "X4593" "X4594" "X4595" "X4596" "X4597" "X4598" "X4599"
## [4600] "X4600" "X4601" "X4602" "X4603" "X4604" "X4605" "X4606" "X4607" "X4608"
## [4609] "X4609" "X4610" "X4611" "X4612" "X4613" "X4614" "X4615" "X4616" "X4617"
## [4618] "X4618" "X4619" "X4620" "X4621" "X4622" "X4623" "X4624" "X4625" "X4626"
## [4627] "X4627" "X4628" "X4629" "X4630" "X4631" "X4632" "X4633" "X4634" "X4635"
## [4636] "X4636" "X4637" "X4638" "X4639" "X4640" "X4641" "X4642" "X4643" "X4644"
## [4645] "X4645" "X4646" "X4647" "X4648" "X4649" "X4650" "X4651" "X4652" "X4653"
## [4654] "X4654" "X4655" "X4656" "X4657" "X4658" "X4659" "X4660" "X4661" "X4662"
## [4663] "X4663" "X4664" "X4665" "X4666" "X4667" "X4668" "X4669" "X4670" "X4671"
## [4672] "X4672" "X4673" "X4674" "X4675" "X4676" "X4677" "X4678" "X4679" "X4680"
## [4681] "X4681" "X4682" "X4683" "X4684" "X4685" "X4686" "X4687" "X4688" "X4689"
## [4690] "X4690" "X4691" "X4692" "X4693" "X4694" "X4695" "X4696" "X4697" "X4698"
## [4699] "X4699" "X4700" "X4701" "X4702" "X4703" "X4704" "X4705" "X4706" "X4707"
## [4708] "X4708" "X4709" "X4710" "X4711" "X4712" "X4713" "X4714" "X4715" "X4716"
## [4717] "X4717" "X4718" "X4719" "X4720" "X4721" "X4722" "X4723" "X4724" "X4725"
## [4726] "X4726" "X4727" "X4728" "X4729" "X4730" "X4731" "X4732" "X4733" "X4734"
## [4735] "X4735" "X4736" "X4737" "X4738" "X4739" "X4740" "X4741" "X4742" "X4743"
## [4744] "X4744" "X4745" "X4746" "X4747" "X4748" "X4749" "X4750" "X4751" "X4752"
## [4753] "X4753" "X4754" "X4755" "X4756" "X4757" "X4758" "X4759" "X4760" "X4761"
## [4762] "X4762" "X4763" "X4764" "X4765" "X4766" "X4767" "X4768" "X4769" "X4770"
## [4771] "X4771" "X4772" "X4773" "X4774" "X4775" "X4776" "X4777" "X4778" "X4779"
## [4780] "X4780" "X4781" "X4782" "X4783" "X4784" "X4785" "X4786" "X4787" "X4788"
## [4789] "X4789" "X4790" "X4791" "X4792" "X4793" "X4794" "X4795" "X4796" "X4797"
## [4798] "X4798" "X4799" "X4800" "X4801" "X4802" "X4803" "X4804" "X4805" "X4806"
## [4807] "X4807" "X4808" "X4809" "X4810" "X4811" "X4812" "X4813" "X4814" "X4815"
## [4816] "X4816" "X4817" "X4818" "X4819" "X4820" "X4821" "X4822" "X4823" "X4824"
## [4825] "X4825" "X4826" "X4827" "X4828" "X4829" "X4830" "X4831" "X4832" "X4833"
## [4834] "X4834" "X4835" "X4836" "X4837" "X4838" "X4839" "X4840" "X4841" "X4842"
## [4843] "X4843" "X4844" "X4845" "X4846" "X4847" "X4848" "X4849" "X4850" "X4851"
## [4852] "X4852" "X4853" "X4854" "X4855" "X4856" "X4857" "X4858" "X4859" "X4860"
## [4861] "X4861" "X4862" "X4863" "X4864" "X4865" "X4866" "X4867" "X4868" "X4869"
## [4870] "X4870" "X4871" "X4872" "X4873" "X4874" "X4875" "X4876" "X4877" "X4878"
## [4879] "X4879" "X4880" "X4881" "X4882" "X4883" "X4884" "X4885" "X4886" "X4887"
## [4888] "X4888" "X4889" "X4890" "X4891" "X4892" "X4893" "X4894" "X4895" "X4896"
## [4897] "X4897" "X4898" "X4899" "X4900" "X4901" "X4902" "X4903" "X4904" "X4905"
## [4906] "X4906" "X4907" "X4908" "X4909" "X4910" "X4911" "X4912" "X4913" "X4914"
## [4915] "X4915" "X4916" "X4917" "X4918" "X4919" "X4920" "X4921" "X4922" "X4923"
## [4924] "X4924" "X4925" "X4926" "X4927" "X4928" "X4929" "X4930" "X4931" "X4932"
## [4933] "X4933" "X4934" "X4935" "X4936" "X4937" "X4938" "X4939" "X4940" "X4941"
## [4942] "X4942" "X4943" "X4944" "X4945" "X4946" "X4947" "X4948" "X4949" "X4950"
## [4951] "X4951" "X4952" "X4953" "X4954" "X4955" "X4956" "X4957" "X4958" "X4959"
## [4960] "X4960" "X4961" "X4962" "X4963" "X4964" "X4965" "X4966" "X4967" "X4968"
## [4969] "X4969" "X4970" "X4971" "X4972" "X4973" "X4974" "X4975" "X4976" "X4977"
## [4978] "X4978" "X4979" "X4980" "X4981" "X4982" "X4983" "X4984" "X4985" "X4986"
## [4987] "X4987" "X4988" "X4989" "X4990" "X4991" "X4992" "X4993" "X4994" "X4995"
## [4996] "X4996" "X4997" "X4998" "X4999" "X5000" "X5001" "X5002" "X5003" "X5004"
## [5005] "X5005" "X5006" "X5007" "X5008" "X5009" "X5010" "X5011" "X5012" "X5013"
## [5014] "X5014" "X5015" "X5016" "X5017" "X5018" "X5019" "X5020" "X5021" "X5022"
## [5023] "X5023" "X5024" "X5025" "X5026" "X5027" "X5028" "X5029" "X5030" "X5031"
## [5032] "X5032" "X5033" "X5034" "X5035" "X5036" "X5037" "X5038" "X5039" "X5040"
## [5041] "X5041" "X5042" "X5043" "X5044" "X5045" "X5046" "X5047" "X5048" "X5049"
## [5050] "X5050" "X5051" "X5052" "X5053" "X5054" "X5055" "X5056" "X5057" "X5058"
## [5059] "X5059" "X5060" "X5061" "X5062" "X5063" "X5064" "X5065" "X5066" "X5067"
## [5068] "X5068" "X5069" "X5070" "X5071" "X5072" "X5073" "X5074" "X5075" "X5076"
## [5077] "X5077" "X5078" "X5079" "X5080" "X5081" "X5082" "X5083" "X5084" "X5085"
## [5086] "X5086" "X5087" "X5088" "X5089" "X5090" "X5091" "X5092" "X5093" "X5094"
## [5095] "X5095" "X5096" "X5097" "X5098" "X5099" "X5100" "X5101" "X5102" "X5103"
## [5104] "X5104" "X5105" "X5106" "X5107" "X5108" "X5109" "X5110" "X5111" "X5112"
## [5113] "X5113" "X5114" "X5115" "X5116" "X5117" "X5118" "X5119" "X5120" "X5121"
## [5122] "X5122" "X5123" "X5124" "X5125" "X5126" "X5127" "X5128" "X5129" "X5130"
## [5131] "X5131" "X5132" "X5133" "X5134" "X5135" "X5136" "X5137" "X5138" "X5139"
## [5140] "X5140" "X5141" "X5142" "X5143" "X5144" "X5145" "X5146" "X5147" "X5148"
## [5149] "X5149" "X5150" "X5151" "X5152" "X5153" "X5154" "X5155" "X5156" "X5157"
## [5158] "X5158" "X5159" "X5160" "X5161" "X5162" "X5163" "X5164" "X5165" "X5166"
## [5167] "X5167" "X5168" "X5169" "X5170" "X5171" "X5172" "X5173" "X5174" "X5175"
## [5176] "X5176" "X5177" "X5178" "X5179" "X5180" "X5181" "X5182" "X5183" "X5184"
## [5185] "X5185" "X5186" "X5187" "X5188" "X5189" "X5190" "X5191" "X5192" "X5193"
## [5194] "X5194" "X5195" "X5196" "X5197" "X5198" "X5199" "X5200" "X5201" "X5202"
## [5203] "X5203" "X5204" "X5205" "X5206" "X5207" "X5208" "X5209" "X5210" "X5211"
## [5212] "X5212" "X5213" "X5214" "X5215" "X5216" "X5217" "X5218" "X5219" "X5220"
## [5221] "X5221" "X5222" "X5223" "X5224" "X5225" "X5226" "X5227" "X5228" "X5229"
## [5230] "X5230" "X5231" "X5232" "X5233" "X5234" "X5235" "X5236" "X5237" "X5238"
## [5239] "X5239" "X5240" "X5241" "X5242" "X5243" "X5244" "X5245" "X5246" "X5247"
## [5248] "X5248" "X5249" "X5250" "X5251" "X5252" "X5253" "X5254" "X5255" "X5256"
## [5257] "X5257" "X5258" "X5259" "X5260" "X5261" "X5262" "X5263" "X5264" "X5265"
## [5266] "X5266" "X5267" "X5268" "X5269" "X5270" "X5271" "X5272" "X5273" "X5274"
## [5275] "X5275" "X5276" "X5277" "X5278" "X5279" "X5280" "X5281" "X5282" "X5283"
## [5284] "X5284" "X5285" "X5286" "X5287" "X5288" "X5289" "X5290" "X5291" "X5292"
## [5293] "X5293" "X5294" "X5295" "X5296" "X5297" "X5298" "X5299" "X5300" "X5301"
## [5302] "X5302" "X5303" "X5304" "X5305" "X5306" "X5307" "X5308" "X5309" "X5310"
## [5311] "X5311" "X5312" "X5313" "X5314" "X5315" "X5316" "X5317" "X5318" "X5319"
## [5320] "X5320" "X5321" "X5322" "X5323" "X5324" "X5325" "X5326" "X5327" "X5328"
## [5329] "X5329" "X5330" "X5331" "X5332" "X5333" "X5334" "X5335" "X5336" "X5337"
## [5338] "X5338" "X5339" "X5340" "X5341" "X5342" "X5343" "X5344" "X5345" "X5346"
## [5347] "X5347" "X5348" "X5349" "X5350" "X5351" "X5352" "X5353" "X5354" "X5355"
## [5356] "X5356" "X5357" "X5358" "X5359" "X5360" "X5361" "X5362" "X5363" "X5364"
## [5365] "X5365" "X5366" "X5367" "X5368" "X5369" "X5370" "X5371" "X5372" "X5373"
## [5374] "X5374" "X5375" "X5376" "X5377" "X5378" "X5379" "X5380" "X5381" "X5382"
## [5383] "X5383" "X5384" "X5385" "X5386" "X5387" "X5388" "X5389" "X5390" "X5391"
## [5392] "X5392" "X5393" "X5394" "X5395" "X5396" "X5397" "X5398" "X5399" "X5400"
## [5401] "X5401" "X5402" "X5403" "X5404" "X5405" "X5406" "X5407" "X5408" "X5409"
## [5410] "X5410" "X5411" "X5412" "X5413" "X5414" "X5415" "X5416" "X5417" "X5418"
## [5419] "X5419" "X5420" "X5421" "X5422" "X5423" "X5424" "X5425" "X5426" "X5427"
## [5428] "X5428" "X5429" "X5430" "X5431" "X5432" "X5433" "X5434" "X5435" "X5436"
## [5437] "X5437" "X5438" "X5439" "X5440" "X5441" "X5442" "X5443" "X5444" "X5445"
## [5446] "X5446" "X5447" "X5448" "X5449" "X5450" "X5451" "X5452" "X5453" "X5454"
## [5455] "X5455" "X5456" "X5457" "X5458" "X5459" "X5460" "X5461" "X5462" "X5463"
## [5464] "X5464" "X5465" "X5466" "X5467" "X5468" "X5469" "X5470" "X5471" "X5472"
## [5473] "X5473" "X5474" "X5475" "X5476" "X5477" "X5478" "X5479" "X5480" "X5481"
## [5482] "X5482" "X5483" "X5484" "X5485" "X5486" "X5487" "X5488" "X5489" "X5490"
## [5491] "X5491" "X5492" "X5493" "X5494" "X5495" "X5496" "X5497" "X5498" "X5499"
## [5500] "X5500" "X5501" "X5502" "X5503" "X5504" "X5505" "X5506" "X5507" "X5508"
## [5509] "X5509" "X5510" "X5511" "X5512" "X5513" "X5514" "X5515" "X5516" "X5517"
## [5518] "X5518" "X5519" "X5520" "X5521" "X5522" "X5523" "X5524" "X5525" "X5526"
## [5527] "X5527" "X5528" "X5529" "X5530" "X5531" "X5532" "X5533" "X5534" "X5535"
## [5536] "X5536" "X5537" "X5538" "X5539" "X5540" "X5541" "X5542" "X5543" "X5544"
## [5545] "X5545" "X5546" "X5547" "X5548" "X5549" "X5550" "X5551" "X5552" "X5553"
## [5554] "X5554" "X5555" "X5556" "X5557" "X5558" "X5559" "X5560" "X5561" "X5562"
## [5563] "X5563" "X5564" "X5565" "X5566" "X5567" "X5568" "X5569" "X5570" "X5571"
## [5572] "X5572" "X5573" "X5574" "X5575" "X5576" "X5577" "X5578" "X5579" "X5580"
## [5581] "X5581" "X5582" "X5583" "X5584" "X5585" "X5586" "X5587" "X5588" "X5589"
## [5590] "X5590" "X5591" "X5592" "X5593" "X5594" "X5595" "X5596" "X5597" "X5598"
## [5599] "X5599" "X5600" "X5601" "X5602" "X5603" "X5604" "X5605" "X5606" "X5607"
## [5608] "X5608" "X5609" "X5610" "X5611" "X5612" "X5613" "X5614" "X5615" "X5616"
## [5617] "X5617" "X5618" "X5619" "X5620" "X5621" "X5622" "X5623" "X5624" "X5625"
## [5626] "X5626" "X5627" "X5628" "X5629" "X5630" "X5631" "X5632" "X5633" "X5634"
## [5635] "X5635" "X5636" "X5637" "X5638" "X5639" "X5640" "X5641" "X5642" "X5643"
## [5644] "X5644" "X5645" "X5646" "X5647" "X5648" "X5649" "X5650" "X5651" "X5652"
## [5653] "X5653" "X5654" "X5655" "X5656" "X5657" "X5658" "X5659" "X5660" "X5661"
## [5662] "X5662" "X5663" "X5664" "X5665" "X5666" "X5667" "X5668" "X5669" "X5670"
## [5671] "X5671" "X5672" "X5673" "X5674" "X5675" "X5676" "X5677" "X5678" "X5679"
## [5680] "X5680" "X5681" "X5682" "X5683" "X5684" "X5685" "X5686" "X5687" "X5688"
## [5689] "X5689" "X5690" "X5691" "X5692" "X5693" "X5694" "X5695" "X5696" "X5697"
## [5698] "X5698" "X5699" "X5700" "X5701" "X5702" "X5703" "X5704" "X5705" "X5706"
## [5707] "X5707" "X5708" "X5709" "X5710" "X5711" "X5712" "X5713" "X5714" "X5715"
## [5716] "X5716" "X5717" "X5718" "X5719" "X5720" "X5721" "X5722" "X5723" "X5724"
## [5725] "X5725" "X5726" "X5727" "X5728" "X5729" "X5730" "X5731" "X5732" "X5733"
## [5734] "X5734" "X5735" "X5736" "X5737" "X5738" "X5739" "X5740" "X5741" "X5742"
## [5743] "X5743" "X5744" "X5745" "X5746" "X5747" "X5748" "X5749" "X5750" "X5751"
## [5752] "X5752" "X5753" "X5754" "X5755" "X5756" "X5757" "X5758" "X5759" "X5760"
## [5761] "X5761" "X5762" "X5763" "X5764" "X5765" "X5766" "X5767" "X5768" "X5769"
## [5770] "X5770" "X5771" "X5772" "X5773" "X5774" "X5775" "X5776" "X5777" "X5778"
## [5779] "X5779" "X5780" "X5781" "X5782" "X5783" "X5784" "X5785" "X5786" "X5787"
## [5788] "X5788" "X5789" "X5790" "X5791" "X5792" "X5793" "X5794" "X5795" "X5796"
## [5797] "X5797" "X5798" "X5799" "X5800" "X5801" "X5802" "X5803" "X5804" "X5805"
## [5806] "X5806" "X5807" "X5808" "X5809" "X5810" "X5811" "X5812" "X5813" "X5814"
## [5815] "X5815" "X5816" "X5817" "X5818" "X5819" "X5820" "X5821" "X5822" "X5823"
## [5824] "X5824" "X5825" "X5826" "X5827" "X5828" "X5829" "X5830" "X5831" "X5832"
## [5833] "X5833" "X5834" "X5835" "X5836" "X5837" "X5838" "X5839" "X5840" "X5841"
## [5842] "X5842" "X5843" "X5844" "X5845" "X5846" "X5847" "X5848" "X5849" "X5850"
## [5851] "X5851" "X5852" "X5853" "X5854" "X5855" "X5856" "X5857" "X5858" "X5859"
## [5860] "X5860" "X5861" "X5862" "X5863" "X5864" "X5865" "X5866" "X5867" "X5868"
## [5869] "X5869" "X5870" "X5871" "X5872" "X5873" "X5874" "X5875" "X5876" "X5877"
## [5878] "X5878" "X5879" "X5880" "X5881" "X5882" "X5883" "X5884" "X5885" "X5886"
## [5887] "X5887" "X5888" "X5889" "X5890" "X5891" "X5892" "X5893" "X5894" "X5895"
## [5896] "X5896" "X5897" "X5898" "X5899" "X5900" "X5901" "X5902" "X5903" "X5904"
## [5905] "X5905" "X5906" "X5907" "X5908" "X5909" "X5910" "X5911" "X5912" "X5913"
## [5914] "X5914" "X5915" "X5916" "X5917" "X5918" "X5919" "X5920" "X5921" "X5922"
## [5923] "X5923" "X5924" "X5925" "X5926" "X5927" "X5928" "X5929" "X5930" "X5931"
## [5932] "X5932" "X5933" "X5934" "X5935" "X5936" "X5937" "X5938" "X5939" "X5940"
## [5941] "X5941" "X5942" "X5943" "X5944" "X5945" "X5946" "X5947" "X5948" "X5949"
## [5950] "X5950" "X5951" "X5952" "X5953" "X5954" "X5955" "X5956" "X5957" "X5958"
## [5959] "X5959" "X5960" "X5961" "X5962" "X5963" "X5964" "X5965" "X5966" "X5967"
## [5968] "X5968" "X5969" "X5970" "X5971" "X5972" "X5973" "X5974" "X5975" "X5976"
## [5977] "X5977" "X5978" "X5979" "X5980" "X5981" "X5982" "X5983" "X5984" "X5985"
## [5986] "X5986" "X5987" "X5988" "X5989" "X5990" "X5991" "X5992" "X5993" "X5994"
## [5995] "X5995" "X5996" "X5997" "X5998" "X5999" "X6000" "X6001" "X6002" "X6003"
## [6004] "X6004" "X6005" "X6006" "X6007" "X6008" "X6009" "X6010" "X6011" "X6012"
## [6013] "X6013" "X6014" "X6015" "X6016" "X6017" "X6018" "X6019" "X6020" "X6021"
## [6022] "X6022" "X6023" "X6024" "X6025" "X6026" "X6027" "X6028" "X6029" "X6030"
## [6031] "X6031" "X6032" "X6033" "X6034" "X6035" "X6036" "X6037" "X6038" "X6039"
## [6040] "X6040" "X6041" "X6042" "X6043" "X6044" "X6045" "X6046" "X6047" "X6048"
## [6049] "X6049" "X6050" "X6051" "X6052" "X6053" "X6054" "X6055" "X6056" "X6057"
## [6058] "X6058" "X6059" "X6060" "X6061" "X6062" "X6063" "X6064" "X6065" "X6066"
## [6067] "X6067" "X6068" "X6069" "X6070" "X6071" "X6072" "X6073" "X6074" "X6075"
## [6076] "X6076" "X6077" "X6078" "X6079" "X6080" "X6081" "X6082" "X6083" "X6084"
## [6085] "X6085" "X6086" "X6087" "X6088" "X6089" "X6090" "X6091" "X6092" "X6093"
## [6094] "X6094" "X6095" "X6096" "X6097" "X6098" "X6099" "X6100" "X6101" "X6102"
## [6103] "X6103" "X6104" "X6105" "X6106" "X6107" "X6108" "X6109" "X6110" "X6111"
## [6112] "X6112" "X6113" "X6114" "X6115" "X6116" "X6117" "X6118" "X6119" "X6120"
## [6121] "X6121" "X6122" "X6123" "X6124" "X6125" "X6126" "X6127" "X6128" "X6129"
## [6130] "X6130" "X6131" "X6132" "X6133" "X6134" "X6135" "X6136" "X6137" "X6138"
## [6139] "X6139" "X6140" "X6141" "X6142" "X6143" "X6144" "X6145" "X6146" "X6147"
## [6148] "X6148" "X6149" "X6150" "X6151" "X6152" "X6153" "X6154" "X6155" "X6156"
## [6157] "X6157" "X6158" "X6159" "X6160" "X6161" "X6162" "X6163" "X6164" "X6165"
## [6166] "X6166" "X6167" "X6168" "X6169" "X6170" "X6171" "X6172" "X6173" "X6174"
## [6175] "X6175" "X6176" "X6177" "X6178" "X6179" "X6180" "X6181" "X6182" "X6183"
## [6184] "X6184" "X6185" "X6186" "X6187" "X6188" "X6189" "X6190" "X6191" "X6192"
## [6193] "X6193" "X6194" "X6195" "X6196" "X6197" "X6198" "X6199" "X6200" "X6201"
## [6202] "X6202" "X6203" "X6204" "X6205" "X6206" "X6207" "X6208" "X6209" "X6210"
## [6211] "X6211" "X6212" "X6213" "X6214" "X6215" "X6216" "X6217" "X6218" "X6219"
## [6220] "X6220" "X6221" "X6222" "X6223" "X6224" "X6225" "X6226" "X6227" "X6228"
## [6229] "X6229" "X6230" "X6231" "X6232" "X6233" "X6234" "X6235" "X6236" "X6237"
## [6238] "X6238" "X6239" "X6240" "X6241" "X6242" "X6243" "X6244" "X6245" "X6246"
## [6247] "X6247" "X6248" "X6249" "X6250" "X6251" "X6252" "X6253" "X6254" "X6255"
## [6256] "X6256" "X6257" "X6258" "X6259" "X6260" "X6261" "X6262" "X6263" "X6264"
## [6265] "X6265" "X6266" "X6267" "X6268" "X6269" "X6270" "X6271" "X6272" "X6273"
## [6274] "X6274" "X6275" "X6276" "X6277" "X6278" "X6279" "X6280" "X6281" "X6282"
## [6283] "X6283" "X6284" "X6285" "X6286" "X6287" "X6288" "X6289" "X6290" "X6291"
## [6292] "X6292" "X6293" "X6294" "X6295" "X6296" "X6297" "X6298" "X6299" "X6300"
## [6301] "X6301" "X6302" "X6303" "X6304" "X6305" "X6306" "X6307" "X6308" "X6309"
## [6310] "X6310" "X6311" "X6312" "X6313" "X6314" "X6315" "X6316" "X6317" "X6318"
## [6319] "X6319" "X6320" "X6321" "X6322" "X6323" "X6324" "X6325" "X6326" "X6327"
## [6328] "X6328" "X6329" "X6330" "X6331" "X6332" "X6333" "X6334" "X6335" "X6336"
## [6337] "X6337" "X6338" "X6339" "X6340" "X6341" "X6342" "X6343" "X6344" "X6345"
## [6346] "X6346" "X6347" "X6348" "X6349" "X6350" "X6351" "X6352" "X6353" "X6354"
## [6355] "X6355" "X6356" "X6357" "X6358" "X6359" "X6360" "X6361" "X6362" "X6363"
## [6364] "X6364" "X6365" "X6366" "X6367" "X6368" "X6369" "X6370" "X6371" "X6372"
## [6373] "X6373" "X6374" "X6375" "X6376" "X6377" "X6378" "X6379" "X6380" "X6381"
## [6382] "X6382" "X6383" "X6384" "X6385" "X6386" "X6387" "X6388" "X6389" "X6390"
## [6391] "X6391" "X6392" "X6393" "X6394" "X6395" "X6396" "X6397" "X6398" "X6399"
## [6400] "X6400" "X6401" "X6402" "X6403" "X6404" "X6405" "X6406" "X6407" "X6408"
## [6409] "X6409" "X6410" "X6411" "X6412" "X6413" "X6414" "X6415" "X6416" "X6417"
## [6418] "X6418" "X6419" "X6420" "X6421" "X6422" "X6423" "X6424" "X6425" "X6426"
## [6427] "X6427" "X6428" "X6429" "X6430" "X6431" "X6432" "X6433" "X6434" "X6435"
## [6436] "X6436" "X6437" "X6438" "X6439" "X6440" "X6441" "X6442" "X6443" "X6444"
## [6445] "X6445" "X6446" "X6447" "X6448" "X6449" "X6450" "X6451" "X6452" "X6453"
## [6454] "X6454" "X6455" "X6456" "X6457" "X6458" "X6459" "X6460" "X6461" "X6462"
## [6463] "X6463" "X6464" "X6465" "X6466" "X6467" "X6468" "X6469" "X6470" "X6471"
## [6472] "X6472" "X6473" "X6474" "X6475" "X6476" "X6477" "X6478" "X6479" "X6480"
## [6481] "X6481" "X6482" "X6483" "X6484" "X6485" "X6486" "X6487" "X6488" "X6489"
## [6490] "X6490" "X6491" "X6492" "X6493" "X6494" "X6495" "X6496" "X6497" "X6498"
## [6499] "X6499" "X6500" "X6501" "X6502" "X6503" "X6504" "X6505" "X6506" "X6507"
## [6508] "X6508" "X6509" "X6510" "X6511" "X6512" "X6513" "X6514" "X6515" "X6516"
## [6517] "X6517" "X6518" "X6519" "X6520" "X6521" "X6522" "X6523" "X6524" "X6525"
## [6526] "X6526" "X6527" "X6528" "X6529" "X6530" "X6531" "X6532" "X6533" "X6534"
## [6535] "X6535" "X6536" "X6537" "X6538" "X6539" "X6540" "X6541" "X6542" "X6543"
## [6544] "X6544" "X6545" "X6546" "X6547" "X6548" "X6549" "X6550" "X6551" "X6552"
## [6553] "X6553" "X6554" "X6555" "X6556" "X6557" "X6558" "X6559" "X6560" "X6561"
## [6562] "X6562" "X6563" "X6564" "X6565" "X6566" "X6567" "X6568" "X6569" "X6570"
## [6571] "X6571" "X6572" "X6573" "X6574" "X6575" "X6576" "X6577" "X6578" "X6579"
## [6580] "X6580" "X6581" "X6582" "X6583" "X6584" "X6585" "X6586" "X6587" "X6588"
## [6589] "X6589" "X6590" "X6591" "X6592" "X6593" "X6594" "X6595" "X6596" "X6597"
## [6598] "X6598" "X6599" "X6600" "X6601" "X6602" "X6603" "X6604" "X6605" "X6606"
## [6607] "X6607" "X6608" "X6609" "X6610" "X6611" "X6612" "X6613" "X6614" "X6615"
## [6616] "X6616" "X6617" "X6618" "X6619" "X6620" "X6621" "X6622" "X6623" "X6624"
## [6625] "X6625" "X6626" "X6627" "X6628" "X6629" "X6630" "X6631" "X6632" "X6633"
## [6634] "X6634" "X6635" "X6636" "X6637" "X6638" "X6639" "X6640" "X6641" "X6642"
## [6643] "X6643" "X6644" "X6645" "X6646" "X6647" "X6648" "X6649" "X6650" "X6651"
## [6652] "X6652" "X6653" "X6654" "X6655" "X6656" "X6657" "X6658" "X6659" "X6660"
## [6661] "X6661" "X6662" "X6663" "X6664" "X6665" "X6666" "X6667" "X6668" "X6669"
## [6670] "X6670" "X6671" "X6672" "X6673" "X6674" "X6675" "X6676" "X6677" "X6678"
## [6679] "X6679" "X6680" "X6681" "X6682" "X6683" "X6684" "X6685" "X6686" "X6687"
## [6688] "X6688" "X6689" "X6690" "X6691" "X6692" "X6693" "X6694" "X6695" "X6696"
## [6697] "X6697" "X6698" "X6699" "X6700" "X6701" "X6702" "X6703" "X6704" "X6705"
## [6706] "X6706" "X6707" "X6708" "X6709" "X6710" "X6711" "X6712" "X6713" "X6714"
## [6715] "X6715" "X6716" "X6717" "X6718" "X6719" "X6720" "X6721" "X6722" "X6723"
## [6724] "X6724" "X6725" "X6726" "X6727" "X6728" "X6729" "X6730" "X6731" "X6732"
## [6733] "X6733" "X6734" "X6735" "X6736" "X6737" "X6738" "X6739" "X6740" "X6741"
## [6742] "X6742" "X6743" "X6744" "X6745" "X6746" "X6747" "X6748" "X6749" "X6750"
## [6751] "X6751" "X6752" "X6753" "X6754" "X6755" "X6756" "X6757" "X6758" "X6759"
## [6760] "X6760" "X6761" "X6762" "X6763" "X6764" "X6765" "X6766" "X6767" "X6768"
## [6769] "X6769" "X6770" "X6771" "X6772" "X6773" "X6774" "X6775" "X6776" "X6777"
## [6778] "X6778" "X6779" "X6780" "X6781" "X6782" "X6783" "X6784" "X6785" "X6786"
## [6787] "X6787" "X6788" "X6789" "X6790" "X6791" "X6792" "X6793" "X6794" "X6795"
## [6796] "X6796" "X6797" "X6798" "X6799" "X6800" "X6801" "X6802" "X6803" "X6804"
## [6805] "X6805" "X6806" "X6807" "X6808" "X6809" "X6810" "X6811" "X6812" "X6813"
## [6814] "X6814" "X6815" "X6816" "X6817" "X6818" "X6819" "X6820" "X6821" "X6822"
## [6823] "X6823" "X6824" "X6825" "X6826" "X6827" "X6828" "X6829" "X6830" "X6831"
## [6832] "X6832" "X6833" "X6834" "X6835" "X6836" "X6837" "X6838" "X6839" "X6840"
## [6841] "X6841" "X6842" "X6843" "X6844" "X6845" "X6846" "X6847" "X6848" "X6849"
## [6850] "X6850" "X6851" "X6852" "X6853" "X6854" "X6855" "X6856" "X6857" "X6858"
## [6859] "X6859" "X6860" "X6861" "X6862" "X6863" "X6864" "X6865" "X6866" "X6867"
## [6868] "X6868" "X6869" "X6870" "X6871" "X6872" "X6873" "X6874" "X6875" "X6876"
## [6877] "X6877" "X6878" "X6879" "X6880" "X6881" "X6882" "X6883" "X6884" "X6885"
## [6886] "X6886" "X6887" "X6888" "X6889" "X6890" "X6891" "X6892" "X6893" "X6894"
## [6895] "X6895" "X6896" "X6897" "X6898" "X6899" "X6900" "X6901" "X6902" "X6903"
## [6904] "X6904" "X6905" "X6906" "X6907" "X6908" "X6909" "X6910" "X6911" "X6912"
## [6913] "X6913" "X6914" "X6915" "X6916" "X6917" "X6918" "X6919" "X6920" "X6921"
## [6922] "X6922" "X6923" "X6924" "X6925" "X6926" "X6927" "X6928" "X6929" "X6930"
## [6931] "X6931" "X6932" "X6933" "X6934" "X6935" "X6936" "X6937" "X6938" "X6939"
## [6940] "X6940" "X6941" "X6942" "X6943" "X6944" "X6945" "X6946"snps_noinvar <- snps_df
find_NAs <- function(x){
NAs_TF <- is.na(x)
i_NA <- which(NAs_TF ==T)
N_NA <- length(i_NA)
return(i_NA) }
N_rows <- nrow(snps_noinvar)
N_NA <- rep(x = 0, times = N_rows)
N_SNPs <- ncol(snps_noinvar)
for( i in 1:N_rows){
i_NA <- find_NAs(snps_noinvar[i,])
N_NA_i <- length(i_NA)
N_NA[i] <- N_NA_i }
cutoff50 <- N_SNPs * 0.5
percent_NA <- N_NA / N_SNPs * 100
any(percent_NA > 50)## [1] FALSEmean(percent_NA)## [1] 0N_meanNAs_per_row <- mean(percent_NA)mean imputation for NAs
mean_imputation <- function(df){
n_cols <- ncol(df)
for(i in 1:n_cols){
column_i <- df[, i]
mean_i <- mean(column_i, na.rm = TRUE)
NAs_i <- which(is.na(column_i))
N_NAs <- length(NAs_i)
column_i[NAs_i] <- mean_i
df[, i] <- column_i }
return(df) }
names(snps_noinvar)[1:10] ## [1] "X1" "X2" "X3" "X4" "X5" "X6" "X7" "X8" "X9" "X10"snps_noNA <- mean_imputation(snps_noinvar)
head(snps_noNA[,1:10])## X1 X2 X3 X4 X5 X6 X7 X8 X9 X10
## HG00096 0 0 0 0 0 0 0 0 0 0
## HG00097 0 0 0 0 0 0 0 0 0 0
## HG00099 0 0 0 0 0 0 0 0 0 0
## HG00100 0 0 0 0 0 0 0 0 0 0
## HG00101 0 0 0 0 0 0 0 0 0 0
## HG00102 0 0 0 0 0 0 0 0 0 0#scale function not working -> coercing NAs
##snps_noNA <- scale(snps_noNA)write csv file
write.csv(snps_noNA,file = "/Users/oliviabrull/Downloads/my_snps_noNA.csv")