2018-05-01から1ヶ月間の記事一覧
Rで等高線など3軸データのプロット ## d<-data.frame(y=tmp2$LIFE_TIME,x1=tmp2$TINRYOU,x2=tmp2$TIKUNENSUU) d2<-d%>% dplyr::filter(x2<=25 & x1<=10^5) %>% dplyr::group_by(x1=10000*round(x1/10000),x2=2*round(x2/2)) %>% dplyr:: summarise(y=mean(y…
https://www.listendata.com/2016/08/dplyr-tutorial.htmlsummarise_ifとかsummrise_allとかのがよさげ 基本 x<-rnorm(1000) y<-rnorm(1000) z<-rnorm(1000) g<-sample(1:5,1000,replace=T) s<-data.frame(x,y,z,g) s%>% group_by(g) %>% dplyr::summarise_e…
とりあえず。 https://cran.r-project.org/web/packages/grpreg/vignettes/quick-start.pdf grpreg パッケージの利用 ##2018.05.22 ##Dummy test x<-factor(c("M","M","F")) dummy(x) y<-factor(c("A","B","A")) dummy(y) dummy.data.frame(data.frame(x,y),…
Rでgroup lasso xx b c y x #install.packages("grplasso") library("grplasso")dataans coef(ans)
GMPのインストール
ggplot2で時系列的なデータを2つ並べた図を作る例。図を並べる x1a<-cumsum(rnorm(1000,mean=0,sd=1)) x1b<-cumsum(rnorm(1000,mean=0,sd=1)) m1<-seq(0,2,length.out=length(x1a)) y1<-x1a*m1 y2<-x1b z1<-x1a z2<-x1b/m1 x1<-1:length(y1) set1<-factor(c…