[1] "libdir: /xchip/tcga/gdac_prod/applications/process_mgmt/firehose_task_registry/tcga-gdac/ClinicalAnalysisAllGenes_new/broadinstitute.org/cancer.genome.analysis/00333/69/"
[1] "op: MUTATION_RATE"
[1] "dfn: /xchip/cga/gdac-prod/tcga-gdac/jobResults/MutSigRun/PCPG-TP/13386409/PCPG-TP.patients.counts_and_rates.txt"
[1] "cfn: /xchip/cga/gdac-prod/tcga-gdac/jobResults/Append_Data/PCPG-TP/15087180/PCPG-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/PCPG-TP/15087180/PCPG-TP.merged_data.txt"
[1] "ok3"

nSamples in clinical file=179, in expression file=178, common to both=178
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] "DCC_UPLOAD_DATE"                       
 [6] "GENDER"                                
 [7] "RADIATION_THERAPY"                     
 [8] "RADIATIONS_RADIATION_REGIMENINDICATION"
 [9] "RACE"                                  
[10] "ETHNICITY"                             
[11] "BATCH_NUMBER"                          

Input Data has 11 rows and 178 columns.

[1] "Batch" "11"   
[1] "Last Follow UP"
Variable 1:'YEARS_TO_BIRTH':	nDistinctValues=59,	numeric=TRUE,	binary=FALSE,	exclude=FALSE.
Variable 2:'VITAL_STATUS':	nDistinctValues=2,	numeric=TRUE,	binary=TRUE,	exclude=FALSE.
Variable 3:'DAYS_TO_DEATH':	nDistinctValues=6,	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=162,	numeric=TRUE,	binary=FALSE,	exclude=TRUE.
[1] "exclude grep('FOLLOWUP', vnms) to deal with survival parameters seperately"
Variable 5:'DCC_UPLOAD_DATE':	nDistinctValues=1,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "DCC_UPLOAD_DATE is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 6:'GENDER':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 7:'RADIATION_THERAPY':	nDistinctValues=1,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "RADIATION_THERAPY is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 8:'RADIATIONS_RADIATION_REGIMENINDICATION':	nDistinctValues=1,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "RADIATIONS_RADIATION_REGIMENINDICATION is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
Variable 9:'RACE':	nDistinctValues=4,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 10:'ETHNICITY':	nDistinctValues=2,	numeric=FALSE,	binary=FALSE,	exclude=FALSE.
Variable 11:'BATCH_NUMBER':	nDistinctValues=1,	numeric=FALSE,	binary=FALSE,	exclude=TRUE.
[1] "BATCH_NUMBER is excluded in the analysis because there is no more than two cases of (unique non-NA values)"
[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 178 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" "GENDER"         "RACE"           "ETHNICITY"     
[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=59, binary=FALSE, numeric=TRUE
AGE_mutation.rate, nv=59, binary=FALSE, numeric=TRUE
$MUTATIONRATE_NONSYNONYMOUS

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

Residuals:
       Min         1Q     Median         3Q        Max 
-4.098e-07 -1.233e-07 -2.690e-08  1.071e-07  4.918e-07 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) 2.097e-07  4.257e-08   4.926 1.93e-06 ***
vv          4.476e-09  8.589e-10   5.211 5.22e-07 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 1.725e-07 on 176 degrees of freedom
Multiple R-squared: 0.1337,	Adjusted R-squared: 0.1287 
F-statistic: 27.15 on 1 and 176 DF,  p-value: 5.218e-07 


$MUTATIONRATE_SILENT

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

Residuals:
       Min         1Q     Median         3Q        Max 
-1.466e-07 -5.543e-08 -1.644e-08  5.303e-08  2.871e-07 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)   
(Intercept) 6.272e-08  1.961e-08   3.199  0.00164 **
vv          1.310e-09  3.956e-10   3.312  0.00112 **
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 7.944e-08 on 176 degrees of freedom
Multiple R-squared: 0.05866,	Adjusted R-squared: 0.05331 
F-statistic: 10.97 on 1 and 176 DF,  p-value: 0.001125 


GENDER, nv=2, binary=FALSE, numeric=FALSE
RACE, nv=4, binary=FALSE, numeric=FALSE
ETHNICITY, nv=2, binary=FALSE, numeric=FALSE
[1] "saved param, results, example.expr in analysis.result.Rdata "
