スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

Rデー

多重の条件に当てはまる行の抽出
今のところ下記以上にスマートに実行する方法は知らない。

> swiss
Fertility Agriculture Examination Education Catholic Infant.Mortality
Courtelary 80.2 17.0 15 12 9.96 22.2
Delemont 83.1 45.1 6 9 84.84 22.2
Franches-Mnt 92.5 39.7 5 5 93.40 20.2
  ・
  ・
  ・


> x <-swiss[swiss["Fertility"]<=72.0,]
> x
Fertility Agriculture Examination Education Catholic Infant.Mortality
Aigle 64.1 62.0 21 12 8.52 16.5
Aubonne 66.9 67.5 14 7 2.27 19.1
Avenches 68.9 60.7 19 12 4.43 22.7
  ・
  ・
  ・


> y <-x[x["Agriculture"]<=60.0,]
> y
Fertility Agriculture Examination Education Catholic Infant.Mortality
Grandson 71.7 34.0 17 8 3.30 20.0
Lausanne 55.7 19.4 26 28 12.11 20.2
La Vallee 54.3 15.2 31 20 2.15 10.8
  ・
  ・
  ・



逆行列の計算と行列の積

> A <- matrix(nrow=3, ncol=3)
> A

[,1] [,2] [,3]
[1,] NA NA NA
[2,] NA NA NA
[3,] NA NA NA

> A[,1] <- round(runif(3, min = 0, max = 9))
> A[,2] <- round(runif(3, min = 0, max = 9))
> A[,3] <- round(runif(3, min = 0, max = 9))
> A

[,1] [,2] [,3]
[1,] 6 4 7
[2,] 2 1 7
[3,] 7 1 2

> B <- solve(A)
> B
[,1] [,2] [,3]
[1,] -0.04347826 -0.008695652 0.1826087
[2,] 0.39130435 -0.321739130 -0.2434783
[3,] -0.04347826 0.191304348 -0.0173913

> round(A %*% B) # 行列の積

[,1] [,2] [,3]
[1,] 1 0 0
[2,] 0 1 0
[3,] 0 0 1
スポンサーサイト

コメントの投稿


非公開コメント

コメント

Profile

zoa

Author:zoa
Metropolis (1927)

Calendar
05 | 2017/06 | 07
- - - - 1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30 -
Labels
Admin
Previous Posts
Recent Comments
Recent Trackback
Archive
このページのトップへ
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。