Skip to content

Instantly share code, notes, and snippets.

View ymattu's full-sized avatar

Yuya MATSUMURA ymattu

View GitHub Profile
@ymattu
ymattu / conditional probit.sas
Created January 3, 2016 09:27
20150103ブログ用1
proc import out = modelsample
datafile = "C:\Users\User\Desktop\三田論\mitaron_conditional_probit1113.csv"
dbms=csv REPLACE;
getnames=yes;
datarow=2;
run;
proc print data=modelsample(obs=30);
run:
data model(keep=id week year month day brand uicp sicp qicp kicp picp price cm sns price_rihanari cm_rihanari sns_rihanari choice chosen demo_icp age sex marry income1 income2 );
/*切片*/
@ymattu
ymattu / japanr.stan
Created January 3, 2016 09:29
20150103_ブログ用2
data
{
int<lower=0> N; //サンプルサイズ(ユニークな個人)
int<lower=0> NT; //サンプルサイズ(全体)
int<lower=0, upper=1> y[NT]; //目的変数(離反)
int id[NT]; //ユニークな個人の数
vector[2] x[NT]; //切片と価格
vector[10] w[NT]; //環境変数
matrix[N,4] z; //属性
}
setwd("E:/TokyoR/JapanR2015")
library(rstan)
#stanファイル名
scr <- "mitaronstandraft.stan"
d <- read.csv("sample.csv")
#ユニークな個人の情報
@ymattu
ymattu / R_on_SASIML.sas
Created January 12, 2016 18:11
SASでRを動かすサンプル
submit/R;
library(RMeCab);
result <- RMeCabC("私は統計学を勉強しています。");
result
endsubmit;
@ymattu
ymattu / sas_env.sas
Created January 13, 2016 02:35
sasv9.cfg
/*自分で設定*/
-RLANG
@ymattu
ymattu / sasv9.sas
Created January 13, 2016 02:36
sasv9.cfg_default
-config "C:\Program Files\SASHome\SASFoundation\9.4\nls\ja\sasv9.cfg
@ymattu
ymattu / sas_r_error.sas
Created January 13, 2016 02:37
sas_r_error
ERROR: The RLANG system option must be specified in the SAS configuration file or on the SAS invocation command line
to enable the submission of R language statements.
@ymattu
ymattu / R_on_SASSTAT.sas
Last active January 13, 2016 13:52
stat_r_try
PROC IML;
submit/R;
library(RMeCab);
result <- RMeCabC("私は大学に通っています。");
result
endsubmit;
run;
/*推定用データ*/
proc import out = modelsample
datafile = "C:\Users\User\Desktop\三田論\mitaron_conditional_probit1115_suitei.csv"
dbms=csv REPLACE;
getnames=yes;
datarow=2;
run;
proc print data=modelsample(obs=30);
run:
data model(keep=id week year month day /*brand*/ uicp sicp qicp kicp picp price cm sns price_rihanari cm_rihanari sns_rihanari choice chosen demo_icp age sex marry income1 income2 );
proc iml;
submit/R;
iris;
endsubmit;
run ImportDataSetFromR("work.iris","iris");
quit;
proc print data=iris(obs=30) /*データの確認(変数名など)*/
run;