[1] "libdir: /xchip/tcga/gdac_prod/applications/process_mgmt/firehose_task_registry/tcga-gdac/ClinicalAnalysisAllGenes_new/broadinstitute.org/cancer.genome.analysis/00333/74/"
[1] "op: MUTATION_RATE"
[1] "dfn: /xchip/cga/gdac-prod/tcga-gdac/jobResults/MutSigRun/SKCM-TM/22985898/iteration1/SKCM-TM.patients.counts_and_rates.txt"
[1] "cfn: /xchip/cga/gdac-prod/tcga-gdac/jobResults/Append_Data/SKCM-TM/22507034/SKCM-TM.merged_data.txt"
[1] "gv: ALL"
[1] "gfn: "
[1] "sfn: "
[1] "fv: ALL"
[1] "ofn: "
[1] "dx: "
[1] "cfn"
[1] "/xchip/cga/gdac-prod/tcga-gdac/jobResults/Append_Data/SKCM-TM/22507034/SKCM-TM.merged_data.txt"
[1] "ok3"

nSamples in clinical file=368, in expression file=290, common to both=290
Number of genes in original expression dataset=2
[1] "ALL"
[1] "data2feature, selection=ALL"
 [1] "YEARS_TO_BIRTH"                      
 [2] "VITAL_STATUS"                        
 [3] "DAYS_TO_DEATH"                       
 [4] "DAYS_TO_LAST_FOLLOWUP"               
 [5] "DAYS_TO_SUBMITTED_SPECIMEN_DX"       
 [6] "PATHOLOGIC_STAGE"                    
 [7] "PATHOLOGY_T_STAGE"                   
 [8] "PATHOLOGY_N_STAGE"                   
 [9] "PATHOLOGY_M_STAGE"                   
[10] "MELANOMA_ULCERATION"                 
[11] "MELANOMA_PRIMARY_KNOWN"              
[12] "BRESLOW_THICKNESS"                   
[13] "GENDER"                              
[14] "DATE_OF_INITIAL_PATHOLOGIC_DIAGNOSIS"
[15] "RADIATION_THERAPY"                   
[16] "RACE"                                
[17] "ETHNICITY"                           

Input Data has 17 rows and 290 columns.

[1] "Last Follow UP"
Variable 1:'YEARS_TO_BIRTH':	nDistinctValues=66,	numeric=TRUE,	binary=FALSE,	exclude=FALSE.
Variable 2:'VITAL_STATUS':	nDistinctValues=2,	numeric=TRUE,	binary=TRUE,	exclude=FALSE.
Variable 3:'DAYS_TO_DEATH':	nDistinctValues=158,	numeric=TRUE,	binary=FALSE,	exclude=TRUE.
[1] "exclude grep('DAYS_?TO', vnms) to deal with survival parameters seperately"
Variable 4:'DAYS_TO_LAST_FOLLOWUP':	nDistinctValues=122,	numeric=TRUE,	binary=FALSE,	exclude=TRUE.
[1] "exclude grep('FOLLOWUP', vnms) to deal with survival parameters seperately"
Variable 5:'DAYS_TO_SUBMITTED_SPECIMEN_DX':	nDistinctValues=231,	numeric=TRUE,	binary=FALSE,	exclude=TRUE.
[1] "exclude grep('DAYS_?TO', vnms) to deal with survival parameters seperately"
Variable 6:'PATHOLOGIC_STAGE':	nDistinctValues=14,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 7:'PATHOLOGY_T_STAGE':	nDistinctValues=15,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 8:'PATHOLOGY_N_STAGE':	nDistinctValues=10,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 9:'PATHOLOGY_M_STAGE':	nDistinctValues=5,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 10:'MELANOMA_ULCERATION':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 11:'MELANOMA_PRIMARY_KNOWN':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 12:'BRESLOW_THICKNESS':	nDistinctValues=99,	numeric=TRUE,	binary=FALSE,	exclude=FALSE.
Variable 13:'GENDER':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 14:'DATE_OF_INITIAL_PATHOLOGIC_DIAGNOSIS':	nDistinctValues=31,	numeric=TRUE,	binary=FALSE,	exclude=TRUE.
[1] "exclude grep('DATE', vnms) to deal with survival parameters seperately"
Variable 15:'RADIATION_THERAPY':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 16:'RACE':	nDistinctValues=3,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 17:'ETHNICITY':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
[1] "## **** detect survival parameters (defined in index such as ind_OS, ind_MFS, ind_RFS, ind_RFS, ind_BCR and ind_d2ssd) *** ##"
[1] "detected survival parameters using [ind_OS, overall_survival]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "detected survival parameters using [ind_OS, curated_overall_survival]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "detected survival parameters using [ind_TCGAOS]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "detected survivial parameters using [ind_MFS]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "detected survival parameters using [ind_RFS]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "detected survival parameters using [ind_BCR]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "detected survival parameters using [ind_Progression]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "detected survival parameters using [index_additional_survival_time]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "detected survival parameters using condition: [is.null(surv.mat)&&(selection=='SURV')]"
[1] "survival parameters accumulated so far"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "************ conversion from categorical data to rank data ********** "
[1] "PATHOLOGY_T_STAGE is converted to numeric rank data using modified categoies"
[1] "PATHOLOGY_N_STAGE is converted to numeric rank data using modified categoies"
[1] "PATHOLOGY_M_STAGE is converted to rank data using modified categoies"
[1] "****** SUMMARY ***** "
Output Data has 290 columns, 2 survival variables, and 12 non-survival variables.
[1] "* survival variables: "
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[2] "VITAL_STATUS"                              
[3] "DAYS_TO_DEATH_OR_LAST_FUP"                 
[4] "VITAL_STATUS"                              
[1] "* non-survival variables: "
 [1] "YEARS_TO_BIRTH"         "PATHOLOGIC_STAGE"       "PATHOLOGY_T_STAGE"     
 [4] "PATHOLOGY_N_STAGE"      "PATHOLOGY_M_STAGE"      "MELANOMA_ULCERATION"   
 [7] "MELANOMA_PRIMARY_KNOWN" "BRESLOW_THICKNESS"      "GENDER"                
[10] "RADIATION_THERAPY"      "RACE"                   "ETHNICITY"             
[1] "changed to 13 non-survival variables adding another age variable for linear regression analysis on mutaion rate."
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[1] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[1] "S"                                         
[2] "TIME_FROM_SPECIMEN_DX_TO_DEATH_OR_LAST_FUP"
[3] "Month"                                     
[1] "check if there is any case_to_report in survival time data or not"
[1] "DAYS_TO_DEATH_OR_LAST_FUP"
[1] "DAYS_TO_DEATH_OR_LAST_FUP"
[1] "D"                         "DAYS_TO_DEATH_OR_LAST_FUP"
[3] "Month"                    
[1] "check if there is any case_to_report in survival time data or not"
[1] "alarming case(s) exist!"
[1] "[  1  ] case_to_report(s) is(are) excluded in survival analysis"
AGE, nv=66, binary=FALSE, numeric=TRUE
AGE_mutation.rate, nv=66, binary=FALSE, numeric=TRUE
$MUTATIONRATE_NONSYNONYMOUS

Call:
lm(formula = as.numeric(i) ~ vv)

Residuals:
       Min         1Q     Median         3Q        Max 
-2.085e-05 -1.186e-05 -5.630e-06  1.460e-06  4.213e-04 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)  
(Intercept) 2.500e-06  6.810e-06   0.367   0.7138  
vv          2.371e-07  1.174e-07   2.021   0.0443 *
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 3.073e-05 on 283 degrees of freedom
  (5 observations deleted due to missingness)
Multiple R-squared: 0.01422,	Adjusted R-squared: 0.01074 
F-statistic: 4.083 on 1 and 283 DF,  p-value: 0.04426 


$MUTATIONRATE_SILENT

Call:
lm(formula = as.numeric(i) ~ vv)

Residuals:
       Min         1Q     Median         3Q        Max 
-1.112e-05 -6.142e-06 -2.915e-06  7.640e-07  2.412e-04 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)  
(Intercept) 7.988e-07  3.791e-06   0.211   0.8333  
vv          1.280e-07  6.534e-08   1.959   0.0511 .
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 1.711e-05 on 283 degrees of freedom
  (5 observations deleted due to missingness)
Multiple R-squared: 0.01338,	Adjusted R-squared: 0.009897 
F-statistic: 3.839 on 1 and 283 DF,  p-value: 0.05106 


PATHOLOGIC_STAGE, nv=14, binary=FALSE, numeric=FALSE
PATHOLOGY_T_STAGE, nv=5, binary=FALSE, numeric=TRUE
[1] "PATHOLOGY_T_STAGE"
[1] "num to class table"
vv1
T0 T1 T2 T3 T4 
23 30 62 60 55 
PATHOLOGY_N_STAGE, nv=4, binary=FALSE, numeric=TRUE
[1] "PATHOLOGY_N_STAGE"
[1] "num to class table"
vv1
 N0  N1  N2  N3 
143  53  34  35 
PATHOLOGY_M_STAGE, nv=2, binary=FALSE, numeric=TRUE
MELANOMA_ULCERATION, nv=2, binary=FALSE, numeric=FALSE
MELANOMA_PRIMARY_KNOWN, nv=2, binary=FALSE, numeric=FALSE
BRESLOW_THICKNESS, nv=99, binary=FALSE, numeric=TRUE
GENDER, nv=2, binary=FALSE, numeric=FALSE
RADIATION_THERAPY, nv=2, binary=FALSE, numeric=FALSE
RACE, nv=3, binary=FALSE, numeric=FALSE
ETHNICITY, nv=2, binary=FALSE, numeric=FALSE
[1] "saved param, results, example.expr in analysis.result.Rdata "
