Skip to content

Instantly share code, notes, and snippets.

View MarcinKosinski's full-sized avatar
💭
on a sabbatical

Marcin MarcinKosinski

💭
on a sabbatical
View GitHub Profile
extractSurvival <- function(cohorts){
survivalData <- list()
for(i in cohorts){
get(paste0(i, ".clinical"), envir = .GlobalEnv) %>%
select(patient.bcr_patient_barcode,
patient.vital_status,
patient.days_to_last_followup,
patient.days_to_death ) %>%
mutate(bcr_patient_barcode = toupper(patient.bcr_patient_barcode),
# install.packages(c('magrittr', 'FSelectorRcpp'))
library(magrittr)
library(FSelectorRcpp)
information_gain( # Calculate the score for each attribute
formula = Species ~ ., # that is on the right side of the formula.
data = iris, # Attributes must exist in the passed data.
type = "infogain", # Choose the type of a score to be calculated.
threads = 2 # Set number of threads in a parallel backend.
) %>%
cut_attrs( # Then take attributes with the highest rank.
devtools::install_github("kassambara/survminer", build_vignettes = TRUE)
library("survminer")
library("survival")
fit<- survfit(Surv(time, status) ~ sex, data = lung)
ggsurvplot(
fit, # survfit object with calculated statistics.
risk.table = TRUE, # show risk table.
pval = TRUE, # show p-value of log-rank test.
conf.int = TRUE, # show confidence intervals for
# point estimaes of survival curves.
library(shiny)
library(shinydashboard)
ui <- dashboardPage(
dashboardHeader(),
dashboardSidebar(
sidebarMenu(
menuItem(
"Comments",
tabName = "disqus_here",
icon = icon("info")
library(stringi)
library(stringr)
library(rvest)
aktualne_oferty <- function(link) {
linki<-read_html(link) %>%
html_nodes('.href-link') %>%
library(VennDiagram)
venn.diagram(list(first = c("a", "b", "c"),
second = c("a", "c", "d"),
third = c("a", "d"))
imagetype = "png",
col = "transparent",
filename = "example.png",
height = 2000,
width = 2000,
@MarcinKosinski
MarcinKosinski / RTCGA_info_gain.R
Last active April 11, 2016 19:46
information gain for RTCGA BRCA
library(devtools)
devtools::install_github('RTCGA/RTCGA')
library(RTCGA)
installTCGA('RTCGA.clinical')
installTCGA('RTCGA.rnaseq')
library(RTCGA.rnaseq)
library(RTCGA.clinical)
library(dplyr)
survivalTCGA(BRCA.clinical) %>%