source('S.R'); Start()

散布図

gPoint(E=D,x=1,y=2,g=0,lm=’‘,lx=’‘,ly=’‘,lg=’‘,r=F,s=1,w=F,f=12,
Lx=F,Ly=F,z=F,a=0,p=’vh’,d=‘both’)
E:,x:x軸,y:y軸,g:グループ軸,lm,lx,ly,lg:ラベル,r:行名,s:ポイントサイズ,
w: 白黒,f:フォントサイズ,Lx:対数軸(x),Ly:対数軸(y),z:両軸ゼロ開始,a:行名角度,
p:配置(‘vh’:平面,‘v’:縦図,‘h’:横図,d:行名の位置調整(‘both’,‘x’,‘y’)

Dt(iris) #データ例:アヤメの萼片と花弁 (df) nr=150
## Class: data.frame / Rows: 150 / Columns: 5
##  
##    Sepal.Length Sepal.Width Petal.Length Petal.Width Species
## 1  5.1          3.5         1.4          0.2         setosa 
## 2  4.9          3.0         1.4          0.2         setosa 
## 3  4.7          3.2         1.3          0.2         setosa 
## 4  4.6          3.1         1.5          0.2         setosa 
## 5  5.0          3.6         1.4          0.2         setosa 
## 6  5.4          3.9         1.7          0.4         setosa 
## 7  4.6          3.4         1.4          0.3         setosa 
## 8  5.0          3.4         1.5          0.2         setosa 
## 9  4.4          2.9         1.4          0.2         setosa 
## 10 4.9          3.1         1.5          0.1         setosa
gPoint(iris,'Sepal.L','Sepal.W') #x,y:アヤメの萼片の長さ・幅

gPoint(iris,1,2)                 #x,y:列番号で指定(同)

gPoint(iris,1,2,5)               #g:グループ軸=5

gPoint(iris,1,2,'Sp')            #g:グループ軸='Sp'(同)

gPoint(iris,1,2,5,z=T)           #z:両軸ゼロ開始

gPoint(iris,1,2,5,s=2)           #s:ポイントサイズ=2

gPoint(iris,1,2,5,f=15)          #f:フォントサイズ=15

gPoint(iris,1,2,5,w=T)           #w:白黒

gPoint(iris,1,2,5,lm='Iris', lx='Sepal length',ly='Sepal width') #ラベル

Dt(diamonds)          #ダイアモンド
## Class: tbl_df, tbl, data.frame / Rows: 53940 / Columns: 10
##  
##    carat cut       color clarity depth table price x    y    z   
## 1  0.23  Ideal     E     SI2     61.5  55    326   3.95 3.98 2.43
## 2  0.21  Premium   E     SI1     59.8  61    326   3.89 3.84 2.31
## 3  0.23  Good      E     VS1     56.9  65    327   4.05 4.07 2.31
## 4  0.29  Premium   I     VS2     62.4  58    334   4.20 4.23 2.63
## 5  0.31  Good      J     SI2     63.3  58    335   4.34 4.35 2.75
## 6  0.24  Very Good J     VVS2    62.8  57    336   3.94 3.96 2.48
## 7  0.24  Very Good I     VVS1    62.3  57    336   3.95 3.98 2.47
## 8  0.26  Very Good H     SI1     61.9  55    337   4.07 4.11 2.53
## 9  0.22  Fair      E     VS2     65.1  61    337   3.87 3.78 2.49
## 10 0.23  Very Good H     VS1     59.4  61    338   4.00 4.05 2.39
gPoint(diamonds,x='ca',y='pr',g='co') #7 sec.

D=decathlon2; Dt()    #データ例:十種競技(df) nr=27
## Class: data.frame / Rows: 27 / Columns: 13
##  
##           X100m Long.jump Shot.put High.jump X400m X110m.hurdle Discus
## SEBRLE    11.04 7.58      14.83    2.07      49.81 14.69        43.75 
## CLAY      10.76 7.40      14.26    1.86      49.37 14.05        50.72 
## BERNARD   11.02 7.23      14.25    1.92      48.93 14.99        40.87 
## YURKOV    11.34 7.09      15.19    2.10      50.42 15.31        46.26 
## ZSIVOCZKY 11.13 7.30      13.48    2.01      48.62 14.17        45.67 
## McMULLEN  10.83 7.31      13.76    2.13      49.91 14.38        44.41 
## MARTINEAU 11.64 6.81      14.57    1.95      50.14 14.93        47.60 
## HERNU     11.37 7.56      14.41    1.86      51.10 15.06        44.99 
## BARRAS    11.33 6.97      14.09    1.95      49.48 14.48        42.10 
## NOOL      11.33 7.27      12.68    1.98      49.20 15.29        37.92 
##           Pole.vault Javeline X1500m Rank Points Competition
## SEBRLE    5.02       63.19    291.7   1   8217   Decastar   
## CLAY      4.92       60.15    301.5   2   8122   Decastar   
## BERNARD   5.32       62.77    280.1   4   8067   Decastar   
## YURKOV    4.72       63.44    276.4   5   8036   Decastar   
## ZSIVOCZKY 4.42       55.37    268.0   7   8004   Decastar   
## McMULLEN  4.42       56.37    285.1   8   7995   Decastar   
## MARTINEAU 4.92       52.33    262.1   9   7802   Decastar   
## HERNU     4.82       57.19    285.1  10   7733   Decastar   
## BARRAS    4.72       55.40    282.0  11   7708   Decastar   
## NOOL      4.62       57.44    266.6  12   7651   Decastar
gPoint()          #デフォルト(g=0:グループなし,r=F:行名なし)

gPoint(g='C')     #g=C:グループ: [C]ompetition

gPoint(r=T)       #r=T:行名あり

gPoint(g='C',r=T) #g=C:グループ: [C]ompetition, r=T:行名あり

gPoint(p='v')           #p:配置='v':縦図

gPoint(p='v',r=T)       #p:配置='v':縦図, r=T:行名あり

gPoint(p='v',g='C')     #p:配置='v':縦図, g=C:グループ: [C]ompetition

gPoint(p='v',g='C',r=T) #p:配置='v':縦図, g=C:グループ: [C]ompetition, r=T:行名あり

gPoint(w=T)                    #両軸 (v1=1,v2=2, w=T:白黒)

gPoint(g='C',lg='Competition') #両軸 (x=1,y=2, lg:凡例)

gPoint(g='C',lg='Competition',s=2)    #s:ポイントサイズ=2

gPoint(g='C',lg='Competition',s=2.5)  #s:ポイントサイズ=2.5

gPoint(g='C',lg='Competition',s=3)    #s:ポイントサイズ=3

D=USArrests; Dt()      #データ例:米国の暴力犯罪率(df) nr=50
## Class: data.frame / Rows: 50 / Columns: 4
##  
##             Murder Assault UrbanPop Rape
## Alabama     13.2   236     58       21.2
## Alaska      10.0   263     48       44.5
## Arizona      8.1   294     80       31.0
## Arkansas     8.8   190     50       19.5
## California   9.0   276     91       40.6
## Colorado     7.9   204     78       38.7
## Connecticut  3.3   110     77       11.1
## Delaware     5.9   238     72       15.8
## Florida     15.4   335     80       31.9
## Georgia     17.4   211     60       25.8
gPoint()                    #デフォルト(x=1,y=2,r=F:行名なし)

gPoint(x='M',y='A',r=T)     #x軸g=[M]urder, y軸=[A]ssault, r=T:行名あり

gPoint(x='M',y='R',r=T)     #x軸g=[M]urder, y軸=[R]ape, r=T:行名あり

gPoint(x='M',y='U',r=T)     #x軸g=[M]urder, y軸=[U]rbanPop, r=T:行名あり

ボックスグラフ

gBox(df=D,y=1,g=NULL,p=NULL,lm=’‘,lx=’‘,ly=’‘,lg=’’,y1=0,y2=0,a=0,f=12,l=F,z=F,w=F)
df,y:y軸,g:グループ列,p:ポイント表示,lx,ly,lg:ラベル,y1,y2:y軸範囲,w:白黒
lx:ラベルx,ly:ラベルy,lg:ラベルg,y1:y軸下限,y2:y軸上限,
a:横軸名角度,f:フォントサイズ,l:対数軸,z:軸ゼロ開始w: 白黒

gBox(iris,'Sepal.Length') #データ: iris, 列:'Sepal.Length'

gBox(iris,1) #同

gBox(iris,1,5) #同+g:グループ列=5

gBox(iris,1,5,p=T) #同+p:ポイント表示+T

gBox(iris,1,5,p=T,w=T) #同+w:白黒+T

gBox(iris,1,5,p=T,lx='がくの長さ',ly='単位: cm') #同+lx:ラベル(x)

gBox(iris,1,5,p=T,lm='Iris',lx='Sepal length',ly='Unit : cm',lg='Iris species',y1=0,y2=10) #同+範囲(y1,y2)

Dt(diamonds) #データ: diamonds
## Class: tbl_df, tbl, data.frame / Rows: 53940 / Columns: 10
##  
##    carat cut       color clarity depth table price x    y    z   
## 1  0.23  Ideal     E     SI2     61.5  55    326   3.95 3.98 2.43
## 2  0.21  Premium   E     SI1     59.8  61    326   3.89 3.84 2.31
## 3  0.23  Good      E     VS1     56.9  65    327   4.05 4.07 2.31
## 4  0.29  Premium   I     VS2     62.4  58    334   4.20 4.23 2.63
## 5  0.31  Good      J     SI2     63.3  58    335   4.34 4.35 2.75
## 6  0.24  Very Good J     VVS2    62.8  57    336   3.94 3.96 2.48
## 7  0.24  Very Good I     VVS1    62.3  57    336   3.95 3.98 2.47
## 8  0.26  Very Good H     SI1     61.9  55    337   4.07 4.11 2.53
## 9  0.22  Fair      E     VS2     65.1  61    337   3.87 3.78 2.49
## 10 0.23  Very Good H     VS1     59.4  61    338   4.00 4.05 2.39
gBox(diamonds,1,2,p=T) #7 sec.

グラフ保存

gSave(p=G,f=‘x.jpg’,w=30,h=20,d=600)
p:グラフオブジェクト,f:ファイル名,w:幅(cm),h:高さ(cm),d:dpi

G=gPoint(iris,'Sepal.L','Sepal.W') #x,y:アヤメの萼片の長さ・幅
G

gSave(f='@1.jpg',w=20,h=15,d=300) #Gを保存, f:ファイル名,w:幅(cm),h:高さ(cm)
gSave(G,'@1.jpg',20,15,300) #同

—–

参考 Referencia

トップページ Portada

(東京大学 上田博人 Hiroto Ueda, Universidad de Tokio, 2022)