This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#To calculate Fisher's p and look at the genome distribution of hypermutated sites in your amplicon, please input an aligned fasta with reference sequence at the top | |
#Step 00. Setup R | |
rm(list=ls()) | |
closeAllConnections() | |
graphics.off() | |
cat("\014") #send CTRL+L to R, clean console | |
StartSysTime <- Sys.time() | |
MyWD <- getwd() | |
setwd(MyWD) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#Goal: To find out time-to-: VSupp401 VReb401 etc (days and event0/censor1) | |
#Keyboard shortcuts | |
#Run-all F5 | |
#Run-line Ctrl+R | |
#CSV date MUST be in yyyy-mm-dd or yyyy/mm/dd (otherwise error) | |
#format: A character string. If not specified, it will try '"%Y-%m-%d"' then '"%Y/%m/%d"' on the first non-'NA' element, and give an error if neither works. | |
#csv must be sorted already by PATID and coldate | |
#Step1: Settings: Input filename |