| 使用OS | PowerPoint | Word |
|---|---|---|
| Windows | PowerPoint 2016 | Word 2016 |
| Mac | PowerPoint for Mac ver. 15.27 | Word for Mac ver. 15.27 |
import bpy
import numpy as np
# 一旦すべてのオブジェクトを消す
for item in bpy.context.scene.objects:
if item.type == 'MESH':
bpy.context.scene.objects.unlink(item)
for item in bpy.data.objects:
if item.type == 'MESH':
bpy.data.objects.remove(item)
for item in bpy.data.meshes:
bpy.data.meshes.remove(item)
for item in bpy.data.materials:
bpy.data.materials.remove(item)
theta = 0
# mesh arrays
verts = []
faces = []
# mesh variables
numX = 200
numY = 200
# wave variables ここは好きに変更してください
amp = 15 # 振幅
lam = 3 # 波長
scale = .15 # 数値が小さければ細かいグリッドになる
#fill verts array
for i in range (-numX, numX):
for j in range(-numY, numY):
x = scale * i
y = scale * j
R1 = amp * np.cos(2 * np.arctan(y/(x + 0.00001)) - theta + lam * np.sqrt(x**2 + y**2))
# 2は中心から出る波形の数
R2 = 60 + np.sqrt(x**2 + y**2)
z = R1/R2
vert = (x,y,z)
verts.append(vert)
#fill faces array
count = 0
for i in range (0, 2*numY * (2*numX-1) ):
if count < 2*numY-1:
A = i # 1番目の点
B = i+1 # 2番目の点
C = (i+2*numY)+1 # 3番目の点
D = (i+2*numY) # 4番目の点
face = (A,B,C,D)
faces.append(face)
count = count + 1
else:
count = 0
#create mesh and object
mesh = bpy.data.meshes.new("wave")
object = bpy.data.objects.new("wave",mesh)
#set mesh location
object.location = bpy.context.scene.cursor_location
bpy.context.scene.objects.link(object)
#create mesh from python data
mesh.from_pydata(verts,[],faces)
mesh.update(calc_edges=True)
add_theme_support( 'custom-header', apply_filters( 'twentysixteen_custom_header_args', array(
'default-text-color' => $default_text_color,
'width' => 1200,
'height' => 280,
'flex-height' => true,
'wp-head-callback' => 'twentysixteen_header_style',
) ) );
names(iris)
[1] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width" "Species"
xnam<-c("Sepal.Length","Sepal.Width","Petal.Length")plots<-list()
aes(x=Species, y=Sepal.Length)
aes_string(x="Species", y="Sepal.Length")
j <- 1
for (i in 1:length(xnam)) {
plots[[j]]<-ggplot(iris, aes_string(x="Species",y=xnam[i])) + geom_boxplot()
j <- j+1
}
args <- c(plots, list(ncol=3, nrow=1)) do.call(grid.arrange, args)grid.arrangeはgridExtraパッケージの関数です。問題なく実行されると以下の様なグラフが出力されるはずです。
install.packages("reshape")
library(reshape)
iris2 <- iris[,-4] #Petal.Widthを削除（上の例に合わせてます） iris2.melt <- melt.data.frame(iris2,id.vars="Species",variable_name="Region")実行してやると、ほらこの通り。魔法みたい。
> head(iris2.melt) Species Region value 1 setosa Sepal.Length 5.1 2 setosa Sepal.Length 4.9 3 setosa Sepal.Length 4.7 4 setosa Sepal.Length 4.6 5 setosa Sepal.Length 5.0 6 setosa Sepal.Length 5.4列名のvalueはデフォルトで決まっているようなので後で変えたいときは以下のようなのを。
names(iris2.melt) <- c("Species","Region","Length")
で、プロット。
ggplot(iris2.melt,aes(Species,Length)) + geom_boxplot() + facet_wrap(~Region, scales="free_y")結果はこんな感じ。
 |
| （クリックすると拡大します） |
 |
| （クリックすると拡大します） |
d <- read.csv("path to file", header = T or F) ワークディレクトリはデフォルトではホームフォルダに設定されてる。変更したければRstudioの環境設定から変える。＜注意：csvファイルのカンマ区切りはスペースを含まないようにする。スペースも文字としてカウントされてしまうため。＞| Species | Season | Value |
|---|---|---|
| A | Spring | 12 |
| A | Summer | 36 |
| B | Spring | 41 |
| B | Summer | 19 |
d[1,] （行）
d[,1] （列）
d[,-c(2,3,4)] #（2,3,4列目を削除）
rm(a,b,c)
d.spring <- subset(d, Season=="Spring")
#ベクトル作成
> v<-c(2,3,4,5,6)
#普通に表示
> v
[1] 2 3 4 5 6
#2で割って余りが0になる＝TRUE
> v%%2==0
[1] TRUE FALSE TRUE FALSE TRUE
#TRUEになるデータだけピックアップ
> v[v%%2 ==0]
[1] 2 4 6
#TRUEのデータが何番目かを表示
> which(v%%2==0)
[1] 1 3 5
#1,3,5番目のデータをピックアップ、二つ上のコードと結果は同じだけど。
> v[which(v%%2==0)]
[1] 2 4 6
plot(d$Season, d$Value, main="title of figure")
ちなみに、タイトルが長すぎて改行を入れたい場合は、 main="title of figure\nsubtitle"
のように改行したいところに\nを入れれば、"title of figure" と "subtitle"の間に改行が入る。 plot(data, xaxt="n")
とxaxtパラメータでx軸の名前を入れないようにして一旦描画。次に、 axis(side=1, at=c(1:n), labels=LabelName)
としてx軸に上書き。sideパラメータは1：下，2：左，3：上，4：右。atはx軸の項目数。その数と同じだけのLabelNameを先に作っておいて、x軸の項目名を置き換えるわけです。このLabelNameは文字列（character）でもベクターでもOKなようです。ここで、項目名が長いために90度回転させたいならlasというパラメータを入れる。こちらの有り難いページによれば以下とのこと。可能な値は 0 (標準・軸に平行), 1 (x軸，y軸とも水平), 2 (x軸，y軸とも軸に垂直)，3 (x軸，y軸とも垂直（半時計回りで90度回転)） の４つ
axis(side=1, at=c(1:n), labels=LabelName, las=3)
hist(d$Value)
install.packages("ggplot2")
読み込み library("ggplot2")
qplot(data=hoge.rule1, Season, Value, main="title of figure") 分類（Species）に基づいて色分けしたいときは qplot(data=hoge.rule1, Season, Value, main="title of figure", color = Species)
+ facet_wrap( ~ Species, scales = "free_y") レイヤーの発想で点を重ねることも可能。以下の例ではsizeが3の点の上に、sizeが1.5の薄いグレーの点を重ねることで赤い境界線だけに見えるプロットになる。 + geom_point(size=3, color="red") + geom_point(size=1.5, color="grey90")
p <- ggplot(d, aes ( Season, Value ) )
p + geom_boxplot ( aes ( fill = Species )
sink(“results.txt”)
と書くと、これ以降の出力をresults.txtに書き出してくれる。止めるときは、
sink()
「（リアルな）自分とは別の存在として、しかしあたかも他人にとっては自分であるかのように振舞う存在があるのが嫌」というもの。
> install.packages('outliers')
> library(outliers)
> grubbs.test(df)
Grubbs test for one outlier
data: df
G = 2.5060, U = 0.5194, p-value = 0.03119
alternative hypothesis: highest value 138.08 is an outlier
type=20 → 大小どちらか片側から最大2つまで外れ値を表示
type=11 → 大小両方から最大2つまでの外れ値を表示
type=10 → 大小どちらか片側から外れ値を1つ表示
>SG(df)
[[1]]
スミルノフ・グラブス検定
data: min(df) = 72.99
t = 1.275, df = 13, p-value = 1
[[2]]
スミルノフ・グラブス検定
data: max(df) = 138.08
t = 2.506, df = 13, p-value = 0.03119
attr(,"class")
[1] "SG"
ここで有意水準を5％とすると、p < 0.05 の場合はその最大値または最小値が外れ値であると判断できます。リンク先の注意書きにもあるように外れ値は1個ずつ除外すべきですので、最大値または最小値を除外して再度SG関数で計算してみる必要があります。
SG(x)$xqplot(carat, price, data = diamonds) + facet_wrap (clarity~cut)
qplot(carat, price, data = diamonds) + facet_wrap (clarity~cut, scales="free")
qplot(carat, price, data = diamonds) + facet_wrap (clarity~cut) + coord_cartesian(xlim = c(0, 1), ylim = c(250,1000))
「調べてみると、どうやら "--resize" オプションが有効なのは可変サイズ (Dynamically expanding storage) の VDI または VHD 形式で、今回失敗した固定サイズ (Fixed-size storage) ではダメらしい。また (上限を小さくする) 縮小もできない。 」固定サイズにしている場合は、上記ブログに方法が載っています。感謝です。
VirtualBox: 仮想ディスクのサイズを拡張する (VDI) - Natz's Digital 漂流記
VboxManage clonehd <path to your original drive> <path of the copy>VBoxManage modifyhd <path to your vdi> --resize <new size in Megabyte>
tell application "Google Chrome" set myURL to get URL of active tab of first window set myTitle to get title of active tab of first window end tell tell application "Safari" to add reading list item myURL with title myTitle
-- my gui_check()
set temp to the clipboard as text
tell application "Firefox"
activate
set myTitle to the name of front window
tell application "System Events"
keystroke "l" using {command down}
keystroke "c" using {command down}
delay 0.5
set myURL to the clipboard as text
end tell
end tell
tell application "Safari" to add reading list item myURL with title myTitle
delay 0.5
set the clipboard to temp as text
(* to gui_check()
tell application "System Events"
if UI elements enabled is false then
tell application "System Preferences"
activate
set current pane to pane "com.apple.preference.universalaccess"
set msg to "GUI scripting is not available." & return & "Do you put the check in" & return & "\"Enable access for assistive devices\" ?"
--display dialog msg buttons {"OK"} default button "OK" with icon note --giving up after 10
display dialog msg buttons {"cansel", "check and try"} with icon note
end tell
set UI elements enabled to true
delay 0.5
tell application "System Preferences" to quit
delay 0.5
end if
end tell
end gui_check *)
tell application "Safari" to quit
set temp to the clipboard as text
tell application "Sleipnir" to activate
tell application "System Events"
keystroke "l" using {command down, option down}
keystroke "c" using {command down}
key code 53
delay 0.5
set myURL to the clipboard as text
end tell
tell application "Safari" to add reading list item myURL with title "from Sleipnir"
delay 0.5
set the clipboard to temp as text
tell application "Safari" to quit