mbk0asis

library(gplots) 

dta<-mouse_epi_HTT


breaks=c(seq(-2,-0.5,length=100),seq(-0.5,0.5,length=100),seq(0.5,2,length=100))  # manual range setting
my_color <- colorRampPalette(c("blue","white","red"))(n = 299) 
#rc <- rainbow(nrow(dta), start=0, end=.3)
#cc <- rainbow(ncol(dta), start=0, end=.1)

mbk0asis

#source("http://bioconductor.org/biocLite.R")
#biocLite(c("DESeq2","gplots","pcaExplorer","bovine.db"))
#biocLite("calibrate")
#biocLite("AnnotationFuncs")

library(DESeq2)
#library(pcaExplorer)
#library(bovine.db)
library(ggplot2)
library(gplots)

mbk0asis

cntNorm_ER <- read.csv("/home/bio1/00-NGS/RNAseq/single_ICM_TE/DESeq/cntNorm_ER.csv",row.names = 1)
head(cntNorm_ER)

colnames(cntNorm_ER)=c("ii1","ii2","ii3","ii4","ii5","ii6","ii7",
                       "ii8","ii9","ii10","ii11","ii12","ii13","ii14",
                       "ii15","ii16","ii17","ii18","ii19","ii20","ii21",
                       "ti1","ti2","ti3","ti4","ti5","ti6","ti7",
                       "ti8","ti9","ti10","ti11","ti12","ti13","ti14",
                       "ti15","ti16","ti17","ti18","ti19","ti20","ti21",
                       "in1","in2","in3","in4","in5","in6","in7",

mbk0asis

#!/bin/bash

cat $1 $2 | muscle |fasta_formatter | paste - - - - | cut -f 1,2,4 | sed 's/>M.//g' 
cat $1 $2 | muscle -clw | grep "*" | sed 's/^ *//g;s/ /!/g' | tr '\n' ' ' | sed 's/ //g' | grep -bo "!" | tr '\n' ' '
printf "\n"

mbk0asis

#!/bin/bash

export LC_ALL=C

printf "*** SNVcounter\n
	*** 'samtools, bcftools, and bedtools' must be installed to run this script\n
	*** Prepare 'reference fasta' and 'SNV list' in BED format\n
	*** merge consecutive SNVs using\n
	*** 'sort -k1,1 -k2,2n SNV.bed | mergeBed -c 4 -o collapse'\n
	*** usage : ./SNVcounter Input.Bam Reference.Fasta mergedSNV.bed\n\n"

mbk0asis

17	7676265	COSM45241	AA	A	.	.	GENE=TP53;STRAND=-;CDS=c.103delT;AA=p.L35fs*9;CNT=1
17	7676266	COSM96588	A	AG	.	.	GENE=TP53;STRAND=-;CDS=c.102_103insC;AA=p.P36fs*7;CNT=4
17	7676266	COSM2745164	AG	A	.	.	GENE=TP53;STRAND=-;CDS=c.102delC;AA=p.L35fs*9;CNT=2
17	7676266	COSM48976	AGG	A	.	.	GENE=TP53;STRAND=-;CDS=c.101_102delCC;AA=p.P34fs*8;CNT=2
17	7676267	COSM5387432	G	C	.	.	GENE=TP53;STRAND=-;CDS=c.102C>G;AA=p.P34P;CNT=1
17	7676268	COSM43672	G	A	.	.	GENE=TP53;STRAND=-;CDS=c.101C>T;AA=p.P34L;CNT=4
17	7676270	COSM4765600	G	GA	.	.	GENE=TP53;STRAND=-;CDS=c.98_99insT;AA=p.P36fs*7;CNT=1
17	7676270	COSM1172527	GG	G	.	.	GENE=TP53;STRAND=-;CDS=c.98delC;AA=p.L35fs*9;CNT=2
17	7676270	COSM22928	GG	G	.	.	GENE=TP53;STRAND=-;CDS=c.98delC;AA=p.L35fs*9;CNT=2
17	7676271	COSM4272154	G	A	.	.	GENE=TP53;STRAND=-;CDS=c.98C>T;AA=p.S33F;CNT=2

mbk0asis

grep -Po 'c..*' COSMIC_TP53_MUT | sed 's/;/\t/g' | cut -f 1,3 | sed 's/CNT=//g'| grep -Po 'c.\d*(\_?\-?\+?)\d*' > position

grep -Po 'c..*' COSMIC_TP53_MUT | sed 's/;/\t/g' | cut -f 1,3 | sed 's/CNT=//g'| grep -Po '[a-zA-Z]*\>?[a-zA-Z]*\t' > sequence

grep -Po 'c..*' COSMIC_TP53_MUT | sed 's/;/\t/g' | cut -f 3 | sed 's/CNT=//g' > count

paste position sequence count | sed 's/c.//g' > COSMIC_TP53_MUT_all

mbk0asis

102_103	insC		4
102	delC		2
101_102	delCC		2
102	C>G		1
101	C>T		4
98_99	insT		1
98	delC		2
98	delC		2
98	C>T		2
97	delT		2

mbk0asis

library(limma)
library(gplots) 
library(FactoMineR)

##########################################
# Data importing
setwd("D:/00-LabData/Lab/Experiments/2017-01/EpiDriver/T_cell/")
dta<-read.csv("EpiDriver_Tcells_biasCorrection_Final_StdQt.csv", row.names = 1, header = T)

mbk0asis

library(WGCNA)
allowWGCNAThreads()

#################################################################################
# Loading expression data

setwd('/home/bio0/00-NGS/TCGA_met/data')

set1 <- read.csv("600samples.cosmic170k.csv", header = T, row.names = 1)
all.set <- read.csv("2400samples.27k.set2.2.csv", header = T, row.names = 1)