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:
library(vcfR)
##
## ***** *** vcfR *** *****
## This is vcfR 1.13.0
## browseVignettes('vcfR') # Documentation
## citation('vcfR') # Citation
## ***** ***** ***** *****
getwd()
## [1] "/Users/nithinkodali/Downloads"
list.files()
## [1] "_Deck.java"
## [2] "_tmp__lease_documents_20221003060021.pdf"
## [3] "_tmp__lease_documents_20221003060050.pdf"
## [4] "10.11-cluster-analysis.pdf"
## [5] "1666239815.868967.jpg"
## [6] "2022 April Grad App.docx"
## [7] "2022 Summer & Fall Landmark Fundraising Schedule_as of 5.26.22.pdf"
## [8] "22.1059792-1299792.ALL.chr22_GRCh38.genotypes.20170504.vcf"
## [9] "7.11.pdf"
## [10] "770d265a-3c0b-4234-9faa-18212e9f5318.pdf"
## [11] "990-EZ_DesiDanceNetworkInc 2020.pdf"
## [12] "abc123_Assignment1_CMPINF0401.zip"
## [13] "abc123_Assignment2_CMPINF0401.zip"
## [14] "abc123_CommunityFoodVolunteeringManager"
## [15] "abc123_CommunityFoodVolunteeringManager 2"
## [16] "Act-GAL4_ Shroom Presentation.pptx"
## [17] "admin_20861-ch-4-nucleic-acids-and-the-rna-world.pdf"
## [18] "admin_22090-ch-12-solutions.pdf"
## [19] "admin_22092-ch-14-chemical-equilibrium.pdf"
## [20] "admin_22097-ch-18-thermodynamics-and-equilibrium.pdf"
## [21] "admin2_biology-12-biological-science-freeman-93-ch-15-dna-and-the-gene-synthesis-and-repair-3467.pdf"
## [22] "all_loci-1.vcf"
## [23] "all_loci.vcf"
## [24] "ALL.chr22_GRCh38.genotypes.20170504.vcf"
## [25] "ALL.chr22_GRCh38.genotypes.20170504.vcf.gz-2.download"
## [26] "allomtery_3_scatterplot3d (1).Rmd"
## [27] "â\u0080\u0094Pngtreeâ\u0080\u0094graduate class of 2022 with_6310583.png"
## [28] "Application-2022-Hybrid-Program.docx"
## [29] "Bench_E_Ashley_Adynn_Nithin_Mar_22_2022.jpg"
## [30] "Bench_E_Ashley_Adynn_Nithin_Mar_29_2022.jpg"
## [31] "bestcnn.sdf"
## [32] "BIOSC0160_R1.docx"
## [33] "bolinger JC2 part B_edit.pptx"
## [34] "c0120_expt3_datasheets.pdf"
## [35] "c0120_expt4_datasheets.pdf"
## [36] "CA Offer Letter - Kodali.pdf"
## [37] "center_function.R"
## [38] "Chapter2.pptx"
## [39] "CHEM0110 Cheat Sheet.pdf"
## [40] "clown-39861.png"
## [41] "cluster_analysis_portfolio.Rmd"
## [42] "cluster_analysis_with_Higgs_aa.R"
## [43] "code_checkpoint_vcfR.html"
## [44] "code_checkpoint_vcfR.Rmd"
## [45] "CODE_CHECKPOINT-first_rstudio_script.R"
## [46] "combinepdf-2.pdf"
## [47] "combinepdf.pdf"
## [48] "compile_actgal4_viability_11_6_22-2.xlsx"
## [49] "Copy of TaazaLogoSmall-04.png"
## [50] "DDN 501c3 Determination Letter.pdf"
## [51] "Dead_Pupae_ShroomB.8.jpg"
## [52] "Dead_Pupae_ShroomB.8x2.jpg"
## [53] "eclipse-inst-jre-mac64.dmg"
## [54] "Exam 2 Review.pdf"
## [55] "Exam 2.pdf"
## [56] "Exam Review 2 CHEM120 Spring 2022.docx"
## [57] "Exam_1_Practice (3).pdf"
## [58] "Exam_1_Practice_KEY.pdf"
## [59] "F48B4F30-FFA2-465D-BB4C-99EA57D44EC3.HEIC"
## [60] "feature_engineering.Rmd"
## [61] "Final 2.pdf"
## [62] "FInal Bio Lab 1 Presentation.pptx"
## [63] "Final.pdf"
## [64] "Flyer Eyewear Discount.pdf"
## [65] "flyer.pdf"
## [66] "Giant Eagle Proposal Management System.pdf"
## [67] "gnina"
## [68] "GPA-CALCULATOR.xls"
## [69] "Graph.java"
## [70] "Hardwig(1).pdf"
## [71] "HashCodeTester.java"
## [72] "hildebrand soriano assertion evidence templates f22.pdf"
## [73] "hildebrand soriano assertion evidence templates f22.pptx"
## [74] "Hoffman"
## [75] "Hoffman.zip"
## [76] "HPS 0612 essay topics 2.23.22.docx"
## [77] "I-9-Acceptable-Documents-List.pdf"
## [78] "IMG_0013.PNG"
## [79] "IMG_3721.jpg"
## [80] "IMG_4537.HEIC"
## [81] "IMG_4543.HEIC"
## [82] "IMG_4544.HEIC"
## [83] "IMG_4545.HEIC"
## [84] "IMG_4546.HEIC"
## [85] "IMG_4547.HEIC"
## [86] "IMG_4554 2.jpg"
## [87] "IMG_4554.jpg"
## [88] "IMG_4555.PNG"
## [89] "IMG_4758.HEIC"
## [90] "IMG_4807.HEIC"
## [91] "IMG_4809.HEIC"
## [92] "IMG_4811.HEIC"
## [93] "IMG_6020.heic"
## [94] "IMG_6081.HEIC"
## [95] "IMG_6145.HEIC"
## [96] "IMG_6149.HEIC"
## [97] "IMG_6156.JPG"
## [98] "IMG_6162.JPG"
## [99] "IMG_6165.JPG"
## [100] "IMG_6189.JPG"
## [101] "IMG_6460.JPG"
## [102] "IMG_6477 2.JPG"
## [103] "IMG_6477.JPG"
## [104] "IMG_6612.JPG"
## [105] "IMG_AE7E91570F83-1.jpeg"
## [106] "IMG_C50E5F1B5F81-1.jpeg"
## [107] "Intro (1).pptx"
## [108] "Intro (2).pptx"
## [109] "Intro-2.pptx"
## [110] "Intro.pptx"
## [111] "ja4052153_si_001.pdf"
## [112] "jab464_CommunityFoodVolunteeringManager"
## [113] "jab464_CommunityFoodVolunteeringManager.zip"
## [114] "jc2 bolinger et al f22_Martin edit.pptx"
## [115] "jm100112j.pdf"
## [116] "Journal Club1_DWS_SP22-2.pptx"
## [117] "Journal Club1_DWS_SP22.pptx"
## [118] "Journal club2_Hitsman_sp22_post.pptx"
## [119] "last slide lecture 3.docx"
## [120] "Late Withdrawal Form.pdf"
## [121] "LinkedList-2.java"
## [122] "LinkedList.java"
## [123] "LUX_2022_AE_Guide_PT_09_28_21.pdf"
## [124] "LuxotticaDigitalEmployeeGuide_09102018.pdf"
## [125] "Macro final.pdf"
## [126] "MALWARE_BYTES-setup-1.80.2.1012.exe"
## [127] "MATLAB_Book_2E_Solution_Manual.doc"
## [128] "micrometer1.25.jpg"
## [129] "micrometer2.5.jpg"
## [130] "Mid_term_2.pdf"
## [131] "Midterm Questions - Mind and Medicine Spring 2022.docx"
## [132] "midterm_exam_1.pdf"
## [133] "Mind and Med First Paper HW.pdf"
## [134] "Mini exam 1.pdf"
## [135] "Movie revisions.docx"
## [136] "MyHashCode.java"
## [137] "nk fine slip1.pdf"
## [138] "Note Jun 21, 2022.pdf"
## [139] "Oct25.docx"
## [140] "Oct25.Rmd"
## [141] "Online Survey Software _ Qualtrics Survey Solutions.pdf"
## [142] "Patel (A.4) & Kodali (B.8) Pyramid Poster Planning.pptx"
## [143] "Patel + Kodali Slides.pptx"
## [144] "PCR Prelab Homework_3-10-22 EDITED.pptx"
## [145] "Peer_Presentation_Evaluation_Duckweed_SP22_nk.docx"
## [146] "PHOTO-2022-05-24-14-29-48.jpg"
## [147] "Pollutant Research.pptx"
## [148] "pollutant websearch_1.7.22-2.docx"
## [149] "pollutant websearch_1.7.22.docx"
## [150] "portfolio_ggpubr_intro-2.Rmd"
## [151] "portfolio_ggpubr_log_transformation.Rmd"
## [152] "Presentation+planning+DWS+sp22%281%29 copy.docx"
## [153] "Quantification_graph.png"
## [154] "Question 1.pdf"
## [155] "Question 2.pdf"
## [156] "R2_DNARep.docx"
## [157] "readingVcf.Rmd"
## [158] "Rec 1 micro.pdf"
## [159] "Rec 2 micro.pdf"
## [160] "Rec 3.pdf"
## [161] "Rec 6.pdf"
## [162] "Recitation feb 21.pdf"
## [163] "removing_fixed_alleles.html"
## [164] "removing_fixed_alleles.Rmd"
## [165] "Resume Template.docx"
## [166] "rsconnect"
## [167] "RStudio-2022.07.1-554.dmg"
## [168] "rstudio-ide.pdf"
## [169] "Scanned Document 20.pdf"
## [170] "Screen Shot 2022-01-19 at 5.18.58 PM.png"
## [171] "Screen%20Shot%202022-01-10%20at%208.47.26%20PM.png"
## [172] "Screen%20Shot%202022-03-22%20at%2011.25.02%20PM.png"
## [173] "sept 19 2022 final slide.docx"
## [174] "Sponsorship Packet 4.0.docx.pdf"
## [175] "SS Econ Exam 1.pdf"
## [176] "Student_SP22anaylzing_course_population_2_20.pptx"
## [177] "Subject.png"
## [178] "T185_WIWNINLM66_rbcLF_IS-rbcL_forward.ab1"
## [179] "Taaza Fiscal Sponsorship Agreement .docx.pdf"
## [180] "Target_GiftCard-Request-Form.pdf"
## [181] "Teamwork semester evaluation_Duckweed_SP22_nk.docx"
## [182] "traj_0.pdb"
## [183] "transaction_history-2.csv"
## [184] "transaction_history.csv"
## [185] "transpose_VCF_data.html"
## [186] "transpose_VCF_data.Rmd"
## [187] "Tree.pdf"
## [188] "Untitled design-2.jpg"
## [189] "Untitled design-3.jpg"
## [190] "Untitled design.jpg"
## [191] "UPMC CFEI Sponsorship Application 2021.docx"
## [192] "vcfR_test.vcf"
## [193] "vcfR_test.vcf.gz"
## [194] "vegan_PCA_amino_acids-STUDENT.html"
## [195] "vegan_PCA_amino_acids-STUDENT.Rmd"
## [196] "vegan_pca_with_msleep-STUDENT.Rmd"
## [197] "viability_shroom_f22.xlsx"
## [198] "Vial1_parents.JPG"
## [199] "VitalSource-Bookshelf_10.1.1.1411.dmg"
## [200] "Warehouse_Donation_Request_090214.pdf"
## [201] "WDTS Offers Accepted_SULI - Fall FY22.xlsx"
## [202] "week08_cluster_analysis-1.pdf"
## [203] "WIWNINLM66_2_NaCl_D14_021522_NSK.jpg"
## [204] "WIWNINLM66_Al2_D14_021522_EW.jpg"
## [205] "WIWNINLM66_CuSO4_D14_021522_SG_1.jpg"
## [206] "WIWNINLM66_Fe1_D14_021522_RV.jpg"
## [207] "WIWNINLM66_Mn_D14_021522_NK_2.jpg"
## [208] "WIWNINLM66_Mn_D14_021522_NK_Cont.jpg"
## [209] "XC1.java"
## [210] "Zoom.pkg.download"
list.files(pattern = "vcf")
## [1] "22.1059792-1299792.ALL.chr22_GRCh38.genotypes.20170504.vcf"
## [2] "all_loci-1.vcf"
## [3] "all_loci.vcf"
## [4] "ALL.chr22_GRCh38.genotypes.20170504.vcf"
## [5] "ALL.chr22_GRCh38.genotypes.20170504.vcf.gz-2.download"
## [6] "code_checkpoint_vcfR.html"
## [7] "code_checkpoint_vcfR.Rmd"
## [8] "vcfR_test.vcf"
## [9] "vcfR_test.vcf.gz"
rd <- read.vcfR("22.1059792-1299792.ALL.chr22_GRCh38.genotypes.20170504.vcf",
convertNA = T)
## Scanning file to determine attributes.
## File attributes:
## meta lines: 130
## header_line: 131
## variant count: 0
## column count: 2513
##
Meta line 130 read in.
## All meta lines processed.
## gt matrix initialized.
## Character matrix gt created.
## Character matrix gt rows: 0
## Character matrix gt cols: 2513
## skip: 0
## nrows: 0
## row_num: 0
##
Processed variant: 0
## All variants processed
library(vcfR)
getwd()
list.files()
list.files(pattern = "vcf")
rd <- read.vcfR("22.1059792-1299792.ALL.chr22_GRCh38.genotypes.20170504.vcf",
convertNA = T)
You can also embed plots, for example:
Note that the echo = FALSE parameter was added to the
code chunk to prevent printing of the R code that generated the
plot.