Hatena::ブログ(Diary)

Triad sou. このページをアンテナに追加 RSSフィード Twitter

2013-09-05 SAS, Getting data set information from the sashelp library

sashelp.vcolumn と sashelp.vtable によるデータセット情報の抽出

sashelp.vcolumn や sashelp.vtable などには、SAS にロードされた全データセットの情報が格納されている。

そのため、フォーマットだったり、オブザベーション数だったり、変数の数だったり、他にも非常に細かい情報を取り出すことができる。

大胆だなぁ。

proc format;
  value myfmt 1 = "Yes" 2 = "No";
run;

data test;
  x = 1;
  y = 2;
  z = 1;
  output;
  x = 2;
  y = 2;
  z = 2;
  output;
  format x best32. y best8. z myfmt.;
data test2;
  x = 1;
  y = 2;
  z = 1;
  output;
  format x best32. y best8. z myfmt.;
run;

proc print data = sashelp.vcolumn;
  where libname = "WORK" & memtype = "DATA";
  var libname memname format;
run;

proc print data = sashelp.vtable;
  where libname = "WORK" & memtype = "DATA";
  var libname memname nobs nvar;
run;
 Obs    libname    memname    format

6468     WORK       TEST      BEST32.
6469     WORK       TEST      BEST8.
6470     WORK       TEST      MYFMT.
6471     WORK       TEST2     BEST32.
6472     WORK       TEST2     BEST8.
6473     WORK       TEST2     MYFMT.
Obs    libname    memname    nobs    nvar

806     WORK       TEST        2       3
807     WORK       TEST2       1       3

2013-01-24 RcmdrPlugin.KMggplot2_0.2-0 is on CRAN now

RcmdrPlugin.KMggplot2_0.2-0 is on CRAN now

I posted a new version of the "RcmdrPlugin.KMggplot2" package, which is an Rcmdr plug-in for a "ggplot2" GUI front-end.

This package assists you to make "ggplot2" graphics.


NEWS

Changes in version 0.2-0 (2013-01-23)
  • Added themes: classic, minimal (ggplot2-0.9.3).
  • Added themes: tufte, economist, solarized, stata, excel, igray, few, wsj (ggthemes).
  • Kaplan-Meier plot: added a option to show a p-value of the log-rank test (Thanks to Sharma Vivek).
  • Kaplan-Meier plot: added a option to draw a line at median survival (Thanks to Sharma Vivek).
  • Scatter plot matrix: fixed bugs and reimplemented, because the ggplot2::plotmatrix() is deprecated.
  • Scatter plot matrix: added stratified plots.
  • Modified start-up for the problem that Commander window open twice when a plug-in is loaded via library() (see Rcmdr 1.9-3 NEWS). In RcmdrPlugin.KMggplot2 0.2-0, you can directly load by library(RcmdrPlugin.KMggplot2).

themes
  • theme_classic (ggplot2)

f:id:triadsou:20130124110741p:image

  • theme_minimal (ggplot2)

f:id:triadsou:20130124110748p:image

  • theme_tufte (ggthemes)

f:id:triadsou:20130124110751p:image

  • theme_economist (ggthemes)

f:id:triadsou:20130124110742p:image

  • theme_solarized (ggthemes)

f:id:triadsou:20130124110749p:image

  • theme_stata (ggthemes)

f:id:triadsou:20130124110750p:image

  • theme_excel (ggthemes)

f:id:triadsou:20130124110743p:image

  • theme_igray (ggthemes)

f:id:triadsou:20130124110745p:image

  • theme_few (ggthemes)

f:id:triadsou:20130124110744p:image

  • theme_wsj (ggthemes)

f:id:triadsou:20130124110753p:image


Kaplan-Meier plot: log-rank p-value and lines at median survival

f:id:triadsou:20130124110746p:image


Scatter plot matrix: stratified plots

f:id:triadsou:20130124110747p:image


Related Posts

2012-12-22 Firefox のタブが接続中のままになるエラー

Firefox のタブが接続中のままになるエラー

エラー: ReferenceError: oldSetTabTitle is not defined
ソースファイル: chrome://tabmixplus/content/utils.js
行: 347

こんなエラーがでてタブの名前が全部接続中に・・・

Tab Mix Plus をオフにするとなおるけどどうしたもんか。



アドオンのソースを調べると IE Tab 2 と干渉していたらしい事が分かった。

IE Tab 2 の方はあまり使わないので切っておくことにしたら無事解決した。

IE Tab 2 修正版 (4.12.22.2) にアップデートしてもなおるみたいです。


エラーが起きた環境

2012-11-22 Time-dependent (varying) covariate / effect

Time-dependent (varying) covariate / Time-dependent (varying) effect

比例ハザードモデルとその拡張のちょっとしたメモ。

  • 色々調べると、time-dependent と time-varying は特に区別していない文献ばかりだった
  • 当然 covariate と effect では意味が違うので注意
    • Proportional hazard
    • Time-dependent covariate
    • Time-dependent effect
    • Time-dependent covariate + effect

2012-09-26 Recoding a factor variable in a data.frame by the levels() funct

Recoding a factor variable in a data.frame by the levels() function

When you want to change the coding of a factor variable in R, you can use the levels() function.

Recoding z as 1 = A, 2-6 = B, 7 = C

d <- data.frame(z = gl(7, 1, 14))
d$zc <- factor(d$z)
levels(d$zc) <- list(A = 1, B = 2:6, C = 7)
d
Result
> levels(d$zc) <- list(A = 1, B = 2:6, C = 7)
> d
   z zc
1  1  A
2  2  B
3  3  B
4  4  B
5  5  B
6  6  B
7  7  C
8  1  A
9  2  B
10 3  B
11 4  B
12 5  B
13 6  B
14 7  C

Recoding z as 1 = "1", 2-6 = "2-6", 7 = "7"

levels(d$zc) <- list("1" = 1, "2-6" = 2:6, "7" = 7)
d
Result
> levels(d$zc) <- list("1" = 1, "2-6" = 2:6, "7" = 7)
> d
   z  zc
1  1   1
2  2 2-6
3  3 2-6
4  4 2-6
5  5 2-6
6  6 2-6
7  7   7
8  1   1
9  2 2-6
10 3 2-6
11 4 2-6
12 5 2-6
13 6 2-6
14 7   7