Tuesday, July 23, 2013

Asians Eat Weird Things

Asians Eat Weird Things 

Thursday, July 11, 2013

How to turn on/off and diagnose the ANC (active noise cancellation) system of 2013 Honda Pilot

How to turn off/on the ANC (active noise cancellation) system of 2013 Honda Pilot
  1. Turn the ignition switch on the I (radio) position.
  2. Turn the radio on and check it operates normally, and turn the radio off.
  3. Press and hold the preset buttons #1, #6, and volumn (VOL) push button at the same time.
  4. when you see 'DIAG' appears on the display, you can release the buttons
  5. Press the #1 preset button, you see 'ANC ON' on the display.
  6. Press the #1 preset button again to turn the ANC OFF and hear a humming/booming noise coming from the speakers for about 1 minutes.
  7. You can press the #1 button repeatedly for ON and OFF, but the humming/booming noise will not come up repeatedly except repeating step 1, 3-6 again.

This is a similar approach for Acura

Wednesday, July 03, 2013

Data visulization - a tutorial

Data visulization - a tutorial
Source: Tyler Rinker's blog

This is a nice video about how to visualize your data more effectively.



You can find the whole blog with slides and data here: The Mechanics of Data Visualization.

Monday, July 01, 2013

How to test proportionality assumption in survival analysis using SAS

How to test proportionality assumption in survival analysis using SAS
 
There are several approaches to test proportionality assumption in survival analysis:
  • Graphical Approach is to plot Log–Log survival function by researched major predictor. If the two line is parallel without cross each other, the assumption is considering confirmed.
PROC LIFETEST DATA=ONE; METHOD=KM PLOTS=(S,LLS);
     TIME SURVTIME*EVENT(0);
     STRATA RISK0;
RUN;

PROC PHREG DATA=ONE;
     MODEL SURVTIME*EVENT(0)=RISK1 AGE SEX;
    BASELINE OUT=LLSOUT LOGLOGS=LOGLOGS;
RUN;
PROC GPLOT DATA=LLSOUT;
     PLOT LOGLOGS*SURVTIME=RISK1;
RUN;
  • Or, you can include an interaction term of a predictor and follow up time in the model. If this interaction term is not significant, there is no violation of assumption. For example:
PROC PHREG DATA=ONE;
MODEL SURVTIME*EVENT(0)=RISK1 TIMEX;
        TIMEX=RISK1*(LOG(SURVTIME)-LOG(average followup time));
          * someone also uses simple forms: TIMEX=RISK1*LOG(SURVTIME);
          * or even: TIMEX=RISK1*SURVTIME; (not good enough); 
     PROPORTIONALITY_TEST: TIMEX;
RUN;
  • Possibly, the easiest and powerful approach after SAS 9.2 is to use ASSESS statement, for example:
ODS GRAPHICS ON;
PROC PHREG DATA=ONE;
MODEL SURVTIME*EVENT(0)=RISK1;
     ASSESS PH/resample;
RUN;
ODS GRAPHICS OFF;
 
What’s New of PROC PHREG of SAS 9.2 (pdf)