## ----------------------------------------------------------------------------- ## Data Structure ## ----include=FALSE------------------------------------------------------------ ## Data Structure library(MethPed) data(MethPed_sample) MethPed_sample_v1<- MethPed_sample ## ----echo=FALSE--------------------------------------------------------------- head(MethPed_sample_v1) ## ----------------------------------------------------------------------------- ## Data class ## ----echo=FALSE--------------------------------------------------------------- class(MethPed_sample_v1) ## ----include=FALSE------------------------------------------------------------ library(Biobase) ## ----------------------------------------------------------------------------- ## Data Structure ## ----echo=FALSE--------------------------------------------------------------- MethPed_sample_v2<-Biobase::exprs(MethPed_sample) head(MethPed_sample_v2) ## ----------------------------------------------------------------------------- ## Data class ## ----echo=FALSE--------------------------------------------------------------- class(MethPed_sample_v2) ## ----------------------------------------------------------------------------- ## Data Structure ## ----echo=FALSE--------------------------------------------------------------- MethPed_sample_v3 <- data.frame(MethPed_sample_v2) MethPed_sample_v3$TargetID <- rownames(MethPed_sample_v3) rownames(MethPed_sample_v3) <-NULL head(MethPed_sample_v3) ## ----------------------------------------------------------------------------- ## Data class ## ----echo=FALSE--------------------------------------------------------------- class(MethPed_sample_v3) ## ----------------------------------------------------------------------------- # List of 900 probes in predictors library(MethPed) data(MethPed_900probes) head(MethPed_900probes) ## ----Loading_package_and_data------------------------------------------------- # Loading package library(MethPed) # Loading data set data(MethPed_sample) head(MethPed_sample) class(MethPed_sample) ## ----Check_missing_probe_data, echo=TRUE-------------------------------------- ## Check for missing value missingIndex <- checkNA(MethPed_sample) missingIndex ## ----include=FALSE------------------------------------------------------------ set.seed(1000) ## ----Run_MethPed_classifier, echo=TRUE, message=TRUE, warning=FALSE----------- # Run the MethPed classifier myClassification <- MethPed(MethPed_sample) ## ----Output_of_Methped_classifier--------------------------------------------- myClassification ## ----MethPed_classifier_output_description_1-3-------------------------------- # First part myClassification$target_id # Second part myClassification$probes # Third part myClassification$sample ## ----MethPed_classifier_output_description_4---------------------------------- # Fourth part myClassification$probes_missing ## ----MethPed_classifier_output_description_5---------------------------------- # Fifth part myClassification$oob_err.rate ## ----MethPed_classifier_output_description_6---------------------------------- # Sixth part myClassification$predictions ## ----Summary_of_MethPed_classifier_1------------------------------------------ summary(myClassification) ## ----include=FALSE------------------------------------------------------------ set.seed(1000) ## ----include=FALSE------------------------------------------------------------ myClassification_max <- MethPed(MethPed_sample,prob=FALSE) ## ----echo=FALSE--------------------------------------------------------------- summary(myClassification_max) ## ----MethPed_bar_plot_custom_code--------------------------------------------- par(mai = c(1, 1, 1, 2), xpd=TRUE) mat<-t(myClassification$predictions) mycols <- c("green",rainbow(nrow(mat),start=0,end=1)[nrow(mat):1],"red") barplot(mat, col = mycols, beside=FALSE,axisnames=TRUE, ylim=c(0,1),xlab= "Sample",ylab="Probability") legend( ncol(mat)+0.5,1, legend = rownames(mat),fill = mycols,xpd=TRUE, cex = 0.6) ## ----MethPed_bar_plot_generic_function_TRUE, echo=FALSE----------------------- plot(myClassification) ## ----MethPed_bar_plot_generic_function_FALSE, echo=FALSE---------------------- plot(myClassification_max) ## ----Function_to_get_missing_probes------------------------------------------- probeMis(myClassification) ## ----------------------------------------------------------------------------- # Loading dataset data(MethPed_sample) class(MethPed_sample) ## ----include=FALSE------------------------------------------------------------ library(Biobase) ## ----echo=FALSE--------------------------------------------------------------- MethPed_sample_matrix<-Biobase::exprs(MethPed_sample) head(MethPed_sample_matrix) ## ----------------------------------------------------------------------------- class(MethPed_sample_matrix) ## ----------------------------------------------------------------------------- checkNA(MethPed_sample_matrix) ## ----------------------------------------------------------------------------- MethPed_sample_missing<-MethPed_sample_matrix MethPed_sample_missing[c(1,10,200),2]<-NA MethPed_sample_missing[c(4,600,500,1000),1]<-NA head(MethPed_sample_missing,10) ## ----------------------------------------------------------------------------- checkNA(MethPed_sample_missing) ## ----Delete_random_seed------------------------------------------------------- if(exists(".Random.seed")) rm(.Random.seed) ## ----include=FALSE------------------------------------------------------------ library(impute) imputedData <- impute.knn(MethPed_sample_missing) ## ----Combine_probe_name_with_impoted_data, include=FALSE---------------------- imputedData <- imputedData$data ## ----echo=FALSE--------------------------------------------------------------- checkNA(imputedData) ## ----include=FALSE------------------------------------------------------------ set.seed(1000) ## ----include=FALSE------------------------------------------------------------ myClassification <- MethPed(imputedData) ## ----echo=FALSE--------------------------------------------------------------- summary(myClassification) ## ----------------------------------------------------------------------------- citation('MethPed') ## ----------------------------------------------------------------------------- sessionInfo()