[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/UCS-TP/22539817/iteration1/UCS-TP.patients.counts_and_rates.txt"
[1] "cfn: /xchip/cga/gdac-prod/tcga-gdac/jobResults/Append_Data/UCS-TP/22507158/UCS-TP.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/UCS-TP/22507158/UCS-TP.merged_data.txt"
[1] "ok3"

nSamples in clinical file=57, in expression file=57, common to both=57
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] "TUMOR_TISSUE_SITE"                   
 [6] "PATHOLOGIC_STAGE"                    
 [7] "PATHOLOGY_T_STAGE"                   
 [8] "PATHOLOGY_N_STAGE"                   
 [9] "PATHOLOGY_M_STAGE"                   
[10] "GENDER"                              
[11] "DATE_OF_INITIAL_PATHOLOGIC_DIAGNOSIS"
[12] "RADIATION_THERAPY"                   
[13] "HISTOLOGICAL_TYPE"                   
[14] "RACE"                                
[15] "ETHNICITY"                           

Input Data has 15 rows and 57 columns.

[1] "Last Follow UP"
Variable 1:'YEARS_TO_BIRTH':	nDistinctValues=30,	numeric=TRUE,	binary=FALSE,	exclude=FALSE.
Variable 2:'VITAL_STATUS':	nDistinctValues=2,	numeric=TRUE,	binary=TRUE,	exclude=FALSE.
Variable 3:'DAYS_TO_DEATH':	nDistinctValues=34,	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=22,	numeric=TRUE,	binary=FALSE,	exclude=TRUE.
[1] "exclude grep('FOLLOWUP', vnms) to deal with survival parameters seperately"
Variable 5:'TUMOR_TISSUE_SITE':	nDistinctValues=1,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "TUMOR_TISSUE_SITE is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 6:'PATHOLOGIC_STAGE':	nDistinctValues=0,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "PATHOLOGIC_STAGE is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 7:'PATHOLOGY_T_STAGE':	nDistinctValues=0,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "PATHOLOGY_T_STAGE is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 8:'PATHOLOGY_N_STAGE':	nDistinctValues=0,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "PATHOLOGY_N_STAGE is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 9:'PATHOLOGY_M_STAGE':	nDistinctValues=0,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "PATHOLOGY_M_STAGE is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 10:'GENDER':	nDistinctValues=1,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "GENDER is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 11:'DATE_OF_INITIAL_PATHOLOGIC_DIAGNOSIS':	nDistinctValues=11,	numeric=TRUE,	binary=FALSE,	exclude=TRUE.
[1] "exclude grep('DATE', vnms) to deal with survival parameters seperately"
Variable 12:'RADIATION_THERAPY':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 13:'HISTOLOGICAL_TYPE':	nDistinctValues=3,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 14:'RACE':	nDistinctValues=3,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 15:'ETHNICITY':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "ETHNICITY excluded in the analysis because there is no case of (both >= 3) in the table below"
              HISPANIC OR LATINO NOT HISPANIC OR LATINO
freq.values   "1"                "43"                  
freq.contrast "43"               "1"                   
both >= 3     "FALSE"            "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] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "detected survival parameters using [ind_OS, curated_overall_survival]"
[1] "survival parameters accumulated so far"
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "detected survival parameters using [ind_TCGAOS]"
[1] "survival parameters accumulated so far"
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "detected survivial parameters using [ind_MFS]"
[1] "survival parameters accumulated so far"
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "detected survival parameters using [ind_RFS]"
[1] "survival parameters accumulated so far"
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "detected survival parameters using [ind_BCR]"
[1] "survival parameters accumulated so far"
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "detected survival parameters using [ind_Progression]"
[1] "survival parameters accumulated so far"
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "detected survival parameters using [index_additional_survival_time]"
[1] "survival parameters accumulated so far"
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "detected survival parameters using condition: [is.null(surv.mat)&&(selection=='SURV')]"
[1] "survival parameters accumulated so far"
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "************ conversion from categorical data to rank data ********** "
[1] "****** SUMMARY ***** "
Output Data has 57 columns, 1 survival variables, and 4 non-survival variables.
[1] "* survival variables: "
[1] "DAYS_TO_DEATH_OR_LAST_FUP" "VITAL_STATUS"             
[1] "* non-survival variables: "
[1] "YEARS_TO_BIRTH"    "RADIATION_THERAPY" "HISTOLOGICAL_TYPE"
[4] "RACE"             
[1] "changed to 5 non-survival variables adding another age variable for linear regression analysis on mutaion rate."
[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=30, binary=FALSE, numeric=TRUE
AGE_mutation.rate, nv=30, binary=FALSE, numeric=TRUE
$MUTATIONRATE_NONSYNONYMOUS

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

Residuals:
       Min         1Q     Median         3Q        Max 
-7.874e-06 -4.010e-06 -1.525e-06 -2.090e-07  9.900e-05 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)
(Intercept) -1.611e-05  1.414e-05  -1.139    0.260
vv           2.841e-07  2.011e-07   1.413    0.163

Residual standard error: 1.394e-05 on 55 degrees of freedom
Multiple R-squared: 0.03501,	Adjusted R-squared: 0.01747 
F-statistic: 1.996 on 1 and 55 DF,  p-value: 0.1634 


$MUTATIONRATE_SILENT

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

Residuals:
       Min         1Q     Median         3Q        Max 
-1.894e-06 -7.762e-07 -4.213e-07 -8.300e-08  1.999e-05 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)
(Intercept) -3.322e-06  2.893e-06  -1.148    0.256
vv           6.024e-08  4.115e-08   1.464    0.149

Residual standard error: 2.852e-06 on 55 degrees of freedom
Multiple R-squared: 0.03751,	Adjusted R-squared: 0.02001 
F-statistic: 2.143 on 1 and 55 DF,  p-value: 0.1489 


RADIATION_THERAPY, nv=2, binary=FALSE, numeric=FALSE
HISTOLOGICAL_TYPE, nv=3, binary=FALSE, numeric=FALSE
RACE, nv=3, binary=FALSE, numeric=FALSE
[1] "saved param, results, example.expr in analysis.result.Rdata "
