【R言語】統計解析フリーソフトR 第5章【GNU R】at MATH
【R言語】統計解析フリーソフトR 第5章【GNU R】 - 暇つぶし2ch760:¥ ◆2VB8wsVUoo
16/09/25 12:16:10.32 VBoukDxw.net


>349 名前:132人目の素数さん :2016/09/22(木) 01:38:42.56 ID:CSzeBKOI
> 337 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 338 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 339 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 340 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 341 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 342 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 343 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 344 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 345 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 346 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>

761:132人目の素数さん
16/09/25 12:29:48.98 b94OmM0w.net
r2norm <- function(n, mu, sigma, rho) {
tmp <- rnorm(n)
x <- mu+sigma*tmp
y <- rho*x + sqrt(1-rho^2)*rnorm(n)
return(data.frame(x=x,y=y))
}
mydata <- r2norm(100, 0, 1, 0.5)
cor(mydata$x,mydata$y)

x1 <- numeric(length=10000)
for (i in 1:10000) {mydata <- r2norm(10, 0, 1, 0)
s<- cor(mydata$x,mydata$y)
x1[i] <- s}
x2 <- numeric(length=10000)
for (i in 1:10000) {mydata <- r2norm(10, 0, 1, .25)
s2<- cor(mydata$x,mydata$y)
x2[i] <- s2}
x3 <- numeric(length=10000)
for (i in 1:10000) {mydata <- r2norm(10, 0, 1, .5)
s3<- cor(mydata$x,mydata$y)
x3[i] <- s3}

plot(density(x1),xlim=c(-1.5,1.5),ylim=c(0,1.6))
lines(density(x2),col=2)
lines(density(x3),col=3)

762:132人目の素数さん
16/09/25 12:31:31.25 b94OmM0w.net
>>727は、>>724に触発されて作ってみました。
複数のサイトのコピペで作っただけですが、なんだか、スッキリしました。

>>724さん、改めてありがとうございました。

763:¥ ◆2VB8wsVUoo
16/09/25 12:38:12.82 VBoukDxw.net


>349 名前:132人目の素数さん :2016/09/22(木) 01:38:42.56 ID:CSzeBKOI
> 337 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 338 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 339 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 340 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 341 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 342 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 343 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 344 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 345 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 346 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>

764:132人目の素数さん
16/09/25 12:39:28.00 b94OmM0w.net
図々しいのですが、加えて質問をさせてください。
ランダム生成からのplotではなく、
確率密度関数によって、きれいにplotするプログラミンも教えてもらえませんか。

765:¥ ◆2VB8wsVUoo
16/09/25 12:43:06.99 VBoukDxw.net


>349 名前:132人目の素数さん :2016/09/22(木) 01:38:42.56 ID:CSzeBKOI
> 337 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 338 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 339 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 340 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 341 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 342 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 343 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 344 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 345 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 346 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>

766:132人目の素数さん
16/09/25 12:57:25.55 R70tUm2r.net
>>730
x<-seq(from=-1.5,to=1.5,length.out = 1000)
d1<-dPearson(x = x,N = 10,rho = 0)
d2<-dPearson(x = x,N = 10,rho = 0.25)
d3<-dPearson(x = x,N = 10,rho = 0.5)
plot(x,d1,type="l",col=1,ylim=c(0,1.6))
lines(x,d2,col=2)
lines(x,d3,col=3)

767:¥ ◆2VB8wsVUoo
16/09/25 13:14:40.61 VBoukDxw.net


>349 名前:132人目の素数さん :2016/09/22(木) 01:38:42.56 ID:CSzeBKOI
> 337 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 338 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 339 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 340 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 341 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 342 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 343 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 344 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 345 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 346 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>

768:132人目の素数さん
16/09/25 13:19:00.04 b94OmM0w.net
>>732
まことに、ありがとうございました。
教えていただいたプログラミングの意味を勉強します。

769:¥ ◆2VB8wsVUoo
16/09/25 13:21:05.27 VBoukDxw.net


>349 名前:132人目の素数さん :2016/09/22(木) 01:38:42.56 ID:CSzeBKOI
> 337 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 338 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 339 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 340 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 341 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 342 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 343 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 344 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 345 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 346 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>

770:132人目の素数さん
16/09/25 15:32:03.95 BL+yjTYo.net
みんな頭いいんだなぁ
勉強中だけどさっぱり分からん

771:¥ ◆2VB8wsVUoo
16/09/25 16:35:30.09 VBoukDxw.net


>349 名前:132人目の素数さん :2016/09/22(木) 01:38:42.56 ID:CSzeBKOI
> 337 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 338 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 339 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 340 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 341 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 342 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 343 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 344 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 345 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 346 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>

772:¥ ◆2VB8wsVUoo
16/09/25 17:13:10.81 VBoukDxw.net


773:¥ ◆2VB8wsVUoo
16/09/25 17:13:27.67 VBoukDxw.net


774:¥ ◆2VB8wsVUoo
16/09/25 17:13:44.32 VBoukDxw.net


775:¥ ◆2VB8wsVUoo
16/09/25 17:13:59.96 VBoukDxw.net


776:¥ ◆2VB8wsVUoo
16/09/25 17:14:16.36 VBoukDxw.net


777:¥ ◆2VB8wsVUoo
16/09/25 17:14:31.55 VBoukDxw.net


778:¥ ◆2VB8wsVUoo
16/09/25 17:14:46.46 VBoukDxw.net


779:¥ ◆2VB8wsVUoo
16/09/25 17:15:03.49 VBoukDxw.net


780:¥ ◆2VB8wsVUoo
16/09/25 17:15:20.02 VBoukDxw.net


781:¥ ◆2VB8wsVUoo
16/09/25 17:15:37.70 VBoukDxw.net


782:¥ ◆2VB8wsVUoo
16/09/25 18:17:03.92 VBoukDxw.net
馬鹿板遊びはもうヤメレ。頭の悪い奴が跋扈したらアカンやろ。東京都庁
みたいにナルぞ。



783:132人目の素数さん
16/10/01 14:58:06.90 uxqnkz1u.net
year month income expense diff
2016 07 50 15 35
2016 08 30 70 -40

このようなdata.frameがあります。
その月毎のincomeとexpenseをこの順に左からならべてbarplotにしたいです。
2016年7月分はincomeを(2016-07,-15)と(2016-07, 35)の間の縦棒、expenseは(2016-07, -15)と(2016-07,0)の間の縦棒。この2本はくっつけて表示させたいです。
次の2016年8月分はincome(2016-08, 0)と(2016-08, 30)の間の縦棒、expenseは(2016-08, -40)と(2016-08,30)の間の縦棒を描画したいです。この2つは7月分とは少し離して描画できないでしょうか。
よろしくお願いします。

784:¥ ◆2VB8wsVUoo
16/10/01 16:14:10.42 hX8kEwRO.net
馬鹿板遊びはもうヤメレ。頭の悪い奴が跋扈したらアカンやろ。東京都庁
みたいにナルぞ。



785:¥ ◆2VB8wsVUoo
16/10/01 17:42:19.82 hX8kEwRO.net


786:¥ ◆2VB8wsVUoo
16/10/01 17:42:37.77 hX8kEwRO.net


787:¥ ◆2VB8wsVUoo
16/10/01 17:42:55.21 hX8kEwRO.net


788:¥ ◆2VB8wsVUoo
16/10/01 17:43:11.49 hX8kEwRO.net


789:¥ ◆2VB8wsVUoo
16/10/01 17:43:28.22 hX8kEwRO.net


790:¥ ◆2VB8wsVUoo
16/10/01 17:43:43.75 hX8kEwRO.net


791:¥ ◆2VB8wsVUoo
16/10/01 17:43:58.43 hX8kEwRO.net


792:¥ ◆2VB8wsVUoo
16/10/01 17:44:13.48 hX8kEwRO.net


793:¥ ◆2VB8wsVUoo
16/10/01 17:44:29.89 hX8kEwRO.net


794:¥ ◆2VB8wsVUoo
16/10/01 17:44:48.85 hX8kEwRO.net
馬鹿板遊びはもうヤメレ。頭の悪い奴が跋扈したらアカンやろ。東京都庁
みたいにナルぞ。



795:¥ ◆2VB8wsVUoo
16/10/02 05:35:39.87 gY266DhX.net


796:¥ ◆2VB8wsVUoo
16/10/02 10:15:40.80 gY266DhX.net


797:132人目の素数さん
16/10/02 21:17:03.85 jWdphjbd.net
>>749
とりあえず、ヒントだけ。
> d <- data.frame(i = 1:5, e = 1:5)
> barplot(t(d[, c("i", "e")]), beside = TRUE, names.arg = paste(7:11, "月"))

798:132人目の素数さん
16/10/04 00:03:38.42 Eu7Cy3cF.net
> 757

ヒントありがとうございます。

year<- c(2016,2016)
month<-c(07,08)
income<-c(50,30)
expense<-c(15,70)
diff<-c(35,-40)
(dtt<-data.frame(YEAR=year, MONTH=month, INCOME=income, EXPENSE=expense, DIFF=diff))
dttm<-as.matrix(dtt)
barplot(t(dttm[,c(3,4)]), beside=TRUE)

教えてもらったここまではできました。

799:¥ ◆2VB8wsVUoo
16/10/04 00:15:40.99 Z


800:aAz6bOT.net



801:¥ ◆2VB8wsVUoo
16/10/04 00:21:26.14 ZaAz6bOT.net


802:¥ ◆2VB8wsVUoo
16/10/04 00:21:43.34 ZaAz6bOT.net


803:¥ ◆2VB8wsVUoo
16/10/04 00:22:00.28 ZaAz6bOT.net


804:¥ ◆2VB8wsVUoo
16/10/04 00:22:17.04 ZaAz6bOT.net


805:¥ ◆2VB8wsVUoo
16/10/04 00:22:31.56 ZaAz6bOT.net


806:¥ ◆2VB8wsVUoo
16/10/04 00:22:47.99 ZaAz6bOT.net


807:¥ ◆2VB8wsVUoo
16/10/04 00:23:08.19 ZaAz6bOT.net


808:¥ ◆2VB8wsVUoo
16/10/04 00:23:26.84 ZaAz6bOT.net


809:¥ ◆2VB8wsVUoo
16/10/04 00:23:45.56 ZaAz6bOT.net


810:¥ ◆2VB8wsVUoo
16/10/04 00:24:05.13 ZaAz6bOT.net


811:132人目の素数さん
16/10/04 23:17:06.05 5lO37yPD.net
rect()とかでbarを描画していくしかないでしょうか。
barplotのままでbarをshiftさせていくのは難しいか。

812:¥ ◆2VB8wsVUoo
16/10/04 23:21:18.58 ZaAz6bOT.net
若い奴を相手にして、『人前ではメッキで遣り過ごせ!』と教える芳雄が、
何と「研究者としての基本的態度」を申し述べ、毎晩の様に連呼して、そ
の重要性を訴えてた。所が芳雄は一向に「何をどうスル事なのか」を説明
せず、子供の私は日々頭を悩ませた。普通に考えれば、ソレは:
★★★「研究対象に向かう時にどういう風に頭を使って考えて創造的になるか」★★★
であろう。そして試行錯誤を繰り返しながら失敗を重ね、その度毎に芳雄
に罵倒され、そして「釜ヶ崎へ行け!」と恐喝された。

でも芳雄を良く観察し、そして:
★★★『芳雄と一緒になって母親を無根拠に罵倒したら「芳雄は大喜び」した!』★★★
ので、やっと理解した。芳雄が言う「研究者としての基本的態度」とは:
1.その場の自分の損得を考える。
2.何も考えずに「芳雄に同調」スル。
3.発言の中身を空虚にして、表現にだけ注意して敬語を使う。
という様な事であり、要は『顔色を窺って、その場を繕え』という事だ。

コレを理解した時はホンマに嬉しかった。芳雄のオツムが『サル並である』
という定理を証明した瞬間だった。



813:¥ ◆2VB8wsVUoo
16/10/04 23:32:08.30 ZaAz6bOT.net


814:¥ ◆2VB8wsVUoo
16/10/04 23:32:24.07 ZaAz6bOT.net


815:¥ ◆2VB8wsVUoo
16/10/04 23:32:40.35 ZaAz6bOT.net


816:¥ ◆2VB8wsVUoo
16/10/04 23:40:51.35 ZaAz6bOT.net


817:¥ ◆2VB8wsVUoo
16/10/04 23:41:06.10 ZaAz6bOT.net


818:¥ ◆2VB8wsVUoo
16/10/04 23:41:21.72 ZaAz6bOT.net


819:¥ ◆2VB8wsVUoo
16/10/04 23:41:38.00 ZaAz6bOT.net


820:¥ ◆2VB8wsVUoo
16/10/04 23:41:53.13 ZaAz6bOT.net


821:¥ ◆2VB8wsVUoo
16/10/04 23:42:11.56 ZaAz6bOT.net


822:¥ ◆2VB8wsVUoo
16/10/04 23:42:28.26 ZaAz6bOT.net


823:¥ ◆2VB8wsVUoo
16/10/04 23:49:23.12 ZaAz6bOT.net


824:132人目の素数さん
16/10/05 00:56:33.07 PcC8Rgne.net
久しぶりに来たらなんだこれ
なんかRに恨みでもあんのか

825:132人目の素数さん
16/10/05 01:21:59.62 0MSpG4tB.net
構うな。NGワードへ登録

826:132人目の素数さん
16/10/08 00:01:46.18 CjOpE3tl.net
>763さん
stack overflowで質問してみました。
Jpegのグラフも載せたのでご覧いただけたら幸いです。
URLリンク(stackoverflow.com)

827:¥ ◆2VB8wsVUoo
16/10/08 00:04:21.68 tA6BJFvJ.net
馬鹿板遊びはもうヤメレ。頭の悪い奴が跋扈したらアカンやろ。東京都庁
みたいにナルぞ。



828:¥ ◆2VB8wsVUoo
16/10/08 00:08:15.98 tA6BJFvJ.net
ワシが炎上さしたろか。



829:¥ ◆2VB8wsVUoo
16/10/08 06:24:51.60 tA6BJFvJ.net


830:¥ ◆2VB8wsVUoo
16/10/08 06:25:08.33 tA6BJFvJ.net


831:¥ ◆2VB8wsVUoo
16/10/08 06:25:25.88 tA6BJFvJ.net


832:¥ ◆2VB8wsVUoo
16/10/08 06:25:42.80 tA6BJFvJ.net


833:¥ ◆2VB8wsVUoo
16/10/08 06:25:58.59 tA6BJFvJ.net


834:¥ ◆2VB8wsVUoo
16/10/08 06:26:14.08 tA6BJFvJ.net


835:¥ ◆2VB8wsVUoo
16/10/08 06:26:31.11 tA6BJFvJ.net


836:¥ ◆2VB8wsVUoo
16/10/08 06:26:47.22 tA6BJFvJ.net


837:¥ ◆2VB8wsVUoo
16/10/08 06:27:05.57 tA6BJFvJ.net


838:¥ ◆2VB8wsVUoo
16/10/08 06:27:24.78 tA6BJFvJ.net


839:¥ ◆2VB8wsVUoo
16/10/08 22:36:33.94 tA6BJFvJ.net


840:132人目の素数さん
16/10/09 09:55:08.20 T1ly88NG.net
$

841:¥ ◆2VB8wsVUoo
16/10/09 10:18:25.03 t75Uv93o.net


842:132人目の素数さん
16/10/09 18:03:54.15 iuUISe+C.net
>>791
あなたが何をやりたいのがこちらに伝わったら、
多分回答できると思うけど、やっぱり分からない。
丁寧に説明しようと努力しているの


843:はひしひしと伝わるのだが。 とりあえず、水色とクリーム色の棒グラフを再現すれば次のようようになる。 1. 描画する数値を計算する m1 <- matrix(c(30,0,30,30), 2) m2 <- matrix(c(-15,-15,0,-40), 2) 2. 描く barplot(m1, beside = TRUE, ylim = c(-40, 40)) barplot(m2, beside = TRUE, add = TRUE)



844:132人目の素数さん
16/10/09 18:06:32.85 iuUISe+C.net
>>807
abline(h=0)
を忘れていた。追加。

845:132人目の素数さん
16/10/10 10:47:31.22 nDVpiw1f.net
基準点をあわせるように自分でベクトルを平行移動させないと多分ムリ
何をやりたいのかは俺もさっぱり分からないけどw

846:743
16/10/13 18:40:49.69 BRFThXT1.net
x

month income_expense
1 2016-07 50
2 2016-07 -15
3 2016-08 30
4 2016-08 -75

x$id <- seq_along(x$income_expense)
x$type <- ifelse(x$income_expense > 0,

847:743
16/10/13 18:42:39.50 BRFThXT1.net
"in",
"out")
x[x$month %in% c("Starting Cash", "End Cash"),
"type"] <- "net"
month_type <- paste(x$month, x$type, sep="-")
x.paste<-cbind(x, month_type)

for (i in 1:(length(x.paste)/2)) {
x.paste$start5[2*i-1] <- ifelse(abs(x.paste$income_expense[2*i-1]) > abs(x.paste$income_expense[2*i]) , x.paste$income_expense[2*i],
0)
x.paste$end5[2*i-1] <- ifelse(abs(x.paste$income_expense[2*i-1]) > abs(x.paste$income_expense[2*i]) , x.paste$income_expense[2*i-1]+x.paste$income_expense[2*i],
x.paste$income_expense[2*i-1])
x.paste$start5[2*i] <- ifelse(abs(x.paste$income_expense[2*i-1]) > abs(x.paste$income_expense[2*i]) ,x.paste$income_expense[2*i] ,
x.paste$income_expense[2*i-1]+x.paste$income_expense[2*i])
x.paste$end5[2*i] <- ifelse(abs(x.paste$income_expense[2*i-1]) > abs(x.paste$income_expense[2*i]) , 0,
x.paste$income_expense[2*i-1])
}

png("Waterfall_bs01.png", width = 500, height = 300, pointsize = 12, bg = "white")
library(ggplot2)
ggplot(x.paste, aes(month_type, fill = type)) + geom_rect(aes(x = month_type,
xmin = id - 0.45, xmax = id + 0.45, ymin = end5,
ymax = start5))

graphics.off()

>> 801さん, 803さん
返答ありがとうございました。
基準線を合わせるように自分でベクトル平行移動させました。
geom_rectで描くことができました。
ありがとうございました。

848:¥ ◆2VB8wsVUoo
16/10/13 19:19:35.31 vumtC5Zn.net


>349 名前:132人目の素数さん :2016/09/22(木) 01:38:42.56 ID:CSzeBKOI
> 337 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 338 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 339 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 340 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 341 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 342 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 343 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 344 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 345 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>
> 346 名前:あぼ~ん[NGName:¥ </b>◆2VB8wsVUoo <b>] 投稿日:あぼ~ん
>

849:¥ ◆2VB8wsVUoo
16/10/13 19:24:07.69 vumtC5Zn.net


850:¥ ◆2VB8wsVUoo
16/10/13 19:24:26.79 vumtC5Zn.net


851:¥ ◆2VB8wsVUoo
16/10/13 19:24:43.92 vumtC5Zn.net


852:¥ ◆2VB8wsVUoo
16/10/13 19:24:59.90 vumtC5Zn.net


853:¥ ◆2VB8wsVUoo
16/10/13 19:25:17.55 vumtC5Zn.net


854:¥ ◆2VB8wsVUoo
16/10/13 19:25:35.26 vumtC5Zn.net


855:¥ ◆2VB8wsVUoo
16/10/13 19:25:52.29 vumtC5Zn.net


856:¥ ◆2VB8wsVUoo
16/10/13 19:26:09.67 vumtC5Zn.net


857:¥ ◆2VB8wsVUoo
16/10/13 19:26:27.29 vumtC5Zn.net


858:¥ ◆2VB8wsVUoo
16/10/13 19:26:44.51 vumtC5Zn.net


859:132人目の素数さん
16/10/13 22:21:29.74 nyVxQQ1F.net


860:¥ ◆2VB8wsVUoo
16/10/14 03:34:44.81 WA0SUenS.net


861:132人目の素数さん
16/10/19 20:13:01.33 eayAaE2q.net
R初心者の者です。
> c <- 10

> c(c = c)
 c
10
となる理由がよく分かりません。
2行目のうち、最初のcが関数で、次のcがただの文字、最後のcが1行目のc(=10)ということですか?
よろしくお願いします。

862:¥ ◆2VB8wsVUoo
16/10/19 20:13:48.01 2Wa6Uw62.net


>191 名前:132人目の素数さん :2016/10/12(水) 14:26:06.21 ID:5WGW5M4m
> 「筑波大学准教授の強制わいせつ痴漢行為」に思う | 月光院璋子の日記
> URLリンク(plaza.rakuten.co.jp)
>
> 芳雄「哲也の話にするな、忘れろ!奴は増田家の恥、もはや養子でさえない」
>

863:¥ ◆2VB8wsVUoo
16/10/19 20:22:29.83 2Wa6Uw62.net


864:¥ ◆2VB8wsVUoo
16/10/19 20:22:47.59 2Wa6Uw62.net


865:¥ ◆2VB8wsVUoo
16/10/19 20:23:03.68 2Wa6Uw62.net


866:¥ ◆2VB8wsVUoo
16/10/19 20:23:18.13 2Wa6Uw62.net


867:¥ ◆2VB8wsVUoo
16/10/19 20:23:32.30 2Wa6Uw62.net


868:¥ ◆2VB8wsVUoo
16/10/19 20:23:47.54 2Wa6Uw62.net


869:¥ ◆2VB8wsVUoo
16/10/19 20:24:03.30 2Wa6Uw62.net


870:¥ ◆2VB8wsVUoo
16/10/19 20:24:18.86 2Wa6Uw62.net


871:¥ ◆2VB8wsVUoo
16/10/19 20:24:34.19 2Wa6Uw62.net


872:¥ ◆2VB8wsVUoo
16/10/19 20:24:49.00 2Wa6Uw62.net


873:¥ ◆2VB8wsVUoo
16/10/20 23:12:45.01 4i85UFaq.net


874:¥ ◆2VB8wsVUoo
16/10/20 23:25:04.98 4i85UFaq.net


875:¥ ◆2VB8wsVUoo
16/10/20 23:29:49.21 4i85UFaq.net


876:¥ ◆2VB8wsVUoo
16/10/21 19:17:48.82 bLSYGql/.net


877:132人目の素数さん
16/10/21 22:15:37.66 eqGGKWzQ.net
>>825
合ってるよ。

878:¥ ◆2VB8wsVUoo
16/10/21 22:18:35.89 bLSYGql/.net


879:132人目の素数さん
16/10/21 22:53:24.77 kk6Osvdu.net
>>841
ご回答ありがとうございます!

880:¥ ◆2VB8wsVUoo
16/10/22 03:14:05.61 dMde0TVU.net


>679 名前:132人目の素数さん :2016/10/20(木) 08:48:57.82 ID:1+lfflhP
> 阪大ごときで研究者目指したらアカンやろw
>
>680 名前:¥ ◆2VB8wsVUoo :2016/10/20(木) 08:53:56.11 ID:4i85UFaq
> ホウ、なるほどナ。
>
> ¥
>
>681 名前:¥ ◆2VB8wsVUoo :2016/10/20(木) 09:21:36.53 ID:4i85UFaq
> そやし東大と京大以外の大学院は全部閉鎖せんとアカンわ。無駄やさかいナ。
>
> ¥
>
>682 名前:¥ ◆2VB8wsVUoo :2016/10/20(木) 09:28:07.86 ID:4i85UFaq
> ほんで宮廷以外の数学科かて全部閉鎖せんとアカンわ。馬鹿板人みたいな
> 低能ゾンビばっかし居っても税金の無駄にナルだけで役に立たへんしナ。
>
> ¥
>

881:¥ ◆2VB8wsVUoo
16/10/22 03:18:56.51 dMde0TVU.net


882:¥ ◆2VB8wsVUoo
16/10/22 03:19:15.84 dMde0TVU.net


883:¥ ◆2VB8wsVUoo
16/10/22 03:19:35.54 dMde0TVU.net


884:¥ ◆2VB8wsVUoo
16/10/22 03:19:53.64 dMde0TVU.net


885:¥ ◆2VB8wsVUoo
16/10/22 03:20:12.39 dMde0TVU.net


886:¥ ◆2VB8wsVUoo
16/10/22 03:20:31.04 dMde0TVU.net


887:¥ ◆2VB8wsVUoo
16/10/22 03:20:50.12 dMde0TVU.net


888:¥ ◆2VB8wsVUoo
16/10/22 03:21:19.32 dMde0TVU.net


889:¥ ◆2VB8wsVUoo
16/10/22 03:21:37.86 dMde0TVU.net


890:¥ ◆2VB8wsVUoo
16/10/22 03:21:54.27 dMde0TVU.net


891:132人目の素数さん
16/10/22 17:06:34.82 XZKF8OyD.net
>>825
これ何かの本で見た覚えがあるんだけど思い出せない

892:132人目の素数さん
16/10/22 21:37:00.71 M7jVWcyA.net
>>855
R言語徹底解説です

893:132人目の素数さん
16/10/22 23:20:02.84 XZKF8OyD.net
>>856
ああ。今確認した。エクササイズ6.2.5の1か。
c(c=c)の最初のcがcombine関数の関数(省略)名cで、
二番目のc(c関数の引数の左辺値)がベクトルの構成要素の名前を定義するcで、
三番目のc(c関数の引数の右辺値)が前の行でcという名前で定義されているオブジェクト(の値)。

894:¥ ◆2VB8wsVUoo
16/10/23 02:03:11.21 3ZIBPF5r.net


>269 名前:132人目の素数さん :2016/09/24(土) 11:16:21.64 ID:ERaem3b8
> 生活保護受給者の三割がメタボだそうだ
> マス哲もそうなのか?
>
>270 名前:¥ ◆2VB8wsVUoo :2016/09/24(土) 12:04:07.63 ID:gc7JIRwz
> 粗食してますんで、太ってはいますがまあ大丈夫そうですわ。今は精神科
> 系の強い薬物が一切必要が無いので、在任中よりも猛烈に健康になったと
> 思いますね。ストレスが何も無いっちゅうんはホンマに有り難いですわ。
>
> ¥
>
>271 名前:132人目の素数さん :2016/10/22(土) 11:25:03.06 ID:EFTkFzKF
> どぴゅ
>
>506 名前:132人目の素数さん :2016/06/11(土) 12:35:37.39 ID:Ise/AxZk
>  
> 哲也はコンヌの黒歴史
>
>507 名前:132人目の素数さん :2016/10/22(土) 11:24:37.17 ID:EFTkFzKF
> どぴゅ
>

895:¥ ◆2VB8wsVUoo
16/10/23 02:15:59.47 3ZIBPF5r.net


896:¥ ◆2VB8wsVUoo
16/10/23 02:16:16.15 3ZIBPF5r.net


897:¥ ◆2VB8wsVUoo
16/10/23 02:16:33.13 3ZIBPF5r.net


898:¥ ◆2VB8wsVUoo
16/10/23 02:16:51.23 3ZIBPF5r.net


899:¥ ◆2VB8wsVUoo
16/10/23 02:17:08.24 3ZIBPF5r.net


900:¥ ◆2VB8wsVUoo
16/10/23 02:17:24.20 3ZIBPF5r.net


901:¥ ◆2VB8wsVUoo
16/10/23 02:17:40.94 3ZIBPF5r.net


902:¥ ◆2VB8wsVUoo
16/10/23 02:17:58.13 3ZIBPF5r.net


903:¥ ◆2VB8wsVUoo
16/10/23 02:18:16.48 3ZIBPF5r.net


904:¥ ◆2VB8wsVUoo
16/10/23 02:18:34.98 3ZIBPF5r.net


905:132人目の素数さん
16/11/13 17:45:07.40 AJbUUUy0V
R-3.3.1   MagaiHin
R-3.2.5   Last

906:132人目の素数さん
16/11/14 05:39:42.16 LliUCPaqU
無料誹謗中傷対策のブッダワークス
URLリンク(www.buddha-works.net)

907:132人目の素数さん
17/03/18 22:10:26.69 9lKmiZM8.net
盛り上げていこうぜ!

908:132人目の素数さん
17/03/19 08


909::25:18.19 ID:BUbMna52.net



910:132人目の素数さん
17/03/30 22:10:28.44 K8rGc65l.net
Rってどれくらい使えたら役に立つんだ?ベイズ階層位はできんとだめ?

911:132人目の素数さん
17/03/31 16:49:02.43 LyGz7dL7.net
目的によるんじゃない? 考えてみ

912:132人目の素数さん
17/04/23 08:56:55.19 qEAq5gtK.net
R 3.4.0
URLリンク(www.r-project.org)

913:132人目の素数さん
17/04/24 23:49:44.87 38+Bilc8.net
>>873
Stanを使えばベイズ階層も苦にならない。
Samplingの待ち時間が苦痛かも。

914:132人目の素数さん
17/04/25 18:50:16.56 M0sS+FAa.net
>>875
おかげでアップデートに気づいたわ

ちなみにyou stupid darknessってどういう意味やろか?

915:132人目の素数さん
17/04/25 19:39:43.28 24RIix2u.net
URLリンク(kurtisscaletta.com)

916:132人目の素数さん
17/04/25 19:46:55.22 24RIix2u.net
1. なにそれ?
2. 暗いと不平を言うよりも、すすんで灯りをつけましょうって聞いたんだ。
3. それはそうだけど、世の常としてそう思わない人もいるだろうね。
4. 暗いぞバカヤロー!

って感じかな?

917:132人目の素数さん
17/04/26 11:52:55.96 tqLIcCJs.net
取引先の、統計をやってるという人と話をした。
名刺には「シニアスタティスティシャン」と書いてある。

話しをしてみると、統計学のダメな入門書に書いてあるような
大ウソを堂々と主張する。

だが相手は取引先であるから、「馬鹿かてめーは!」とか言えない。

言えないからここに書くのだが、
名刺に「シニアスタティスティシャン」とか書くなよ馬鹿たれが!

918:132人目の素数さん
17/04/26 11:59:53.14 tqLIcCJs.net
>>876
Stanって知らんかった。
便利そうだな使ってみる。
ありがとう!

919:869
17/04/26 12:31:23.28 Vq01HD2H.net
>>878
>>879
ありがとう
ピーナッツが読めるぐらいに英語が堪能ならなぁ

920:ま
17/04/28 19:19:21.81 10ctknwL.net
a+b=C

921:ま
17/04/28 20:57:46.81 10ctknwL.net
>>883今の現代数学は間違ってできてるため、すうじでフっ素をあらわしてる。
そのため、人をすいとる。養分の話しか、してない。数学が人を、おしつぶして、時をすいとってしまった
数式をかえよう
X(f)が、あらたな、人を救う。数式。だということになるよね

922:ま
17/04/28 21:05:59.12 10ctknwL.net
4代元素が、すべて、支払いをし、責任をおうべきだと、思う。この世は、人間をつくり、だまし、最後には、すう式で、すいとってしまい、4代元素は次第には、元素になり、死にいたらしめる、などフリーメーソンなどつくった、元素からの安易な暗号の挑戦状なんだと思う

923:ま
17/04/29 02:13:36.09 PXwKZAGU.net
>>883 a+b=abでいい。c+c=ccでいいと思うんだけど…
同じものを+んだから、同じものになればいいわけで。a×90乗=aで、いいいわけじゃない。なんか、変化して元素に、なるなら、回りのものも、全部元素にならないと、いけないのに、自分たちは、固くのこってて、これを、やりたいなら、人が死んでから、やるべき

924:ま
17/04/29 02:39:27.48 PXwKZAGU.net
>>886固くのこされた人の気持ちを考えるべき。だE=MC2にしたって
E=C2で、やるべきなのに、添字にMを、そえて、そのエネルギーの座標に、人間をさしかえてるだけの、話しでね、
時間一次元と空間3次元の4成分だけ、もたされて、それを、吸う、すう、しき、とか、うけとるだけやないか。元素って。最初からE=C2で、この世を回すべき。
E=Mにかきかえるべきだ。それが、正しいと思うし。正当的。アルケーが歩け

元素はまちがった空間をつくっている
E=C2、または、E=M、で、場をつくるべき
相対性理論はまちがっている
死後の世界観を座標に4つ、ひきつける空間をもたせて、座標に吸わせることを、表しており、正当性のかけらもない
死後なら死後で座標をCがやればいい。わけで、それを有限にしており、できなかった、のだ
だから、元素が悪いのだすべて火、土、水、空に原因がある

925:132人目の素数さん
17/04/29 02:58:57.49 OVQ3EEWI.net
馬鹿が現れた

926:132人目の素数さん
17/05/01 17:04:16.76 xPREFfVb.net
-δC2こうで

927:132人目の素数さん
17/05/01 17:09:14.52 xPREFfVb.net
>>883я

928:132人目の素数さん
17/05/01 17:12:13.05 xPREFfVb.net
タダシさ

929:132人目の素数さん
17/05/01 17:19:22.91 xPREFfVb.net
タダシさ=こたえ

930:¥
17/05/08 20:06:42.83 OR+quqWp.net


931:¥
17/05/08 20:14:24.76 OR+quqWp.net


932:¥
17/05/08 20:14:49.35 OR+quqWp.net


933:¥
17/05/08 20:15:09.69 OR+quqWp.net


934:¥
17/05/08 20:15:31.60 OR+quqWp.net


935:¥
17/05/08 20:15:52.66 OR+quqWp.net


936:¥
17/05/08 20:16:14.65 OR+quqWp.net


937:¥
17/05/08 20:16:41.12 OR+quqWp.net


938:¥
17/05/08 20:17:02.39 OR+quqWp.net


939:¥
17/05/08 20:17:24.76 OR+quqWp.net


940:132人目の素数さん
17/05/18 12:12:04.69 b7yGK+hkS
ウィンドゥズにインストールされたRではR検定に合格できません。
R専用のスパコンでなければ、すべての解析ができないからです。
Fortran2008の関数が導入してあるコンパイラーがあれば、検定用Rのパッケージが構築可能だとか
GCCはバージョン4.12.5で導入予定。完成は未定です。


R Developersより

941:132人目の素数さん
17/05/25 19:49:56.38 zLanYWay.net
数値の20170501を2017-05-01ってな感じで日付に変換したいんですが、何か良い方法はないでしょうか。
力技は考えたんですが、100万行くらい変換しないといけないので、力技以外で変換したいんです。

942:132人目の素数さん
17/05/25 20:11:21.15 t+im+FUk.net
as.Date("20170501",format = "%Y%m%d")

943:¥
17/05/26 09:14:55.04 iCP5fMHR.net


944:¥
17/05/26 09:15:13.28 iCP5fMHR.net


945:¥
17/05/26 09:15:33.73 iCP5fMHR.net


946:¥
17/05/26 09:15:55.04 iCP5fMHR.net


947:¥
17/05/26 09:16:19.41 iCP5fMHR.net


948:¥
17/05/26 09:16:43.19 iCP5fMHR.net


949:¥
17/05/26 09:17:05.17 iCP5fMHR.net


950:¥
17/05/26 09:17:26.87 iCP5fMHR.net


951:¥
17/05/26 09:17:48.95 iCP5fMHR.net


952:¥
17/05/26 09:18:09.48 iCP5fMHR.net


953:132人目の素数さん
17/06/16 15:45:46.35 r+tPpOy0.net
R version 3.1.1ユーザです。

> a = 5
> is.double(a)
[1] TRUE
> is.integer(a)
[1] FALSE

のようにTRUEとFALSEが逆になってしまうのはなぜなんでしょう。
整数を代入するだけでは整数だと予測してくれないんでしょうか?

954:¥
17/06/16 16:12:39.43 TDPjrtUP.net


>112 名前:¥ ◇2VB8wsVUoo 2017/06/15(木) 08:39:06.70 ID:7Rm0/VD6
> ★★★数学徒は情熱的な霊感により主観的に暮らし、唯ひたすら自己の世界に沈潜すべき。★★★
>
> 佐藤幹夫を見よ、ラマヌジャンを見よ
>
> 本物の¥
>

955:132人目の素数さん
17/06/16 17:21:40.28 nDRgYZZV.net
>>916
その代入だと実数として扱われる

a <- 5L

とすれば整数として扱われるようになる

956:¥
17/06/16 17:44:03.74 TDPjrtUP.net


>112 名前:¥ ◇2VB8wsVUoo 2017/06/15(木) 08:39:06.70 ID:7Rm0/VD6
> ★★★数学徒は情熱的な霊感により主観的に暮らし、唯ひたすら自己の世界に沈潜すべき。★★★
>
> 佐藤幹夫を見よ、ラマヌジャンを見よ
>
> 本物の¥
>

957:¥
17/06/16 19:31:56.18 TDPjrtUP.net


958:¥
17/06/16 19:32:16.87 TDPjrtUP.net


959:¥
17/06/16 19:32:37.07 TDPjrtUP.net


960:¥
17/06/16 19:32:57.87 TDPjrtUP.net


961:¥
17/06/16 19:33:17.21 TDPjrtUP.net


962:¥
17/06/16 19:33:35.76 TDPjrtUP.net


963:¥
17/06/16 19:33:55.49 TDPjrtUP.net


964:¥
17/06/16 19:34:15.63 TDPjrtUP.net


965:¥
17/06/16 19:34:36.83 TDPjrtUP.net


966:¥
17/06/16 19:34:57.09 TDPjrtUP.net


967:907
17/06/17 10:07:37.47 oboZrWLI.net
>>918
ユーザーが明示的に指定してやらないとすべて実数になっちゃうんですね?
お答えくださってどうも有難うございました。

968:¥
17/06/17 10:26:55.03 x2f6D4gs.net
★★★忖度と処世術に汚染された日本人:権威主義的な支配と損したくない人達★★★
  ~~~芳雄氏が言う『研究者としての基本的態度』とは一体何だろうか~~~

佐藤幹夫:自分自身の素朴な疑問に真剣に耳を傾ける。⇒不滅の金字塔を打ち立てる。
糞父芳雄:人間関係を駆使し他人を操り根回しを行う。⇒ハリボテお教授として君臨。

隠蔽の財務省、嘘吐きの文科省、そして問答無用に屈服させる官邸。コレでも先進国?

(佐藤師がしてたのは本物の研究だ。だが)芳雄氏がしてたのはケケケ、ケンキュウ。
外見を繕って偉そう見せさえすれば何でもヨロシ。ほんで教授になりさえすれば研究の
中身なんて何でもヨロシ。そもそも論文なんてモンは、外国の権威ある雑誌に掲載され
さえすれば、その中身のギロンなんて何でもヨロシ。そやし適当に書いてしまえ~~~
中身がダメだと知ってて、ソレでもSTAP論文を外国に投稿して受理される。発覚したら
適当に言い逃れる醜い態度。オツムのダメな大学院生に「虚偽の良品ラベル」を貼って
世間に出荷するハリボテ大学は詐欺行為そのもの。世間に媚びを売って客商売に徹し、
『売れさえすれば学生の脳の質なんて何でもヨロシ』と居直る大学。そしてブランド名
だけを見て仕入れる世間。●●は一流大学やさかい、きっと優秀なエリートやろwww

中身を何も説明しないで、問答無用に上から押し付ける。ソレをイチャモンで騒いで、
そして邪魔して潰そうとする周囲の下々。大学教員も国会議事堂も、そして馬鹿板人の
遣ってる事も皆同じだ。日本人はバカ民族であり、今は外国にもちゃんとバレてるので
海外からも軽蔑されるだけであり、そのうちにどの国からも信用されなくなるだろう。

近視眼的で打算的な人生観を息子に押し付ける父親と、大脳に栄養が足りてない連中が
跋扈する永田町や霞が関に支配される国に住む不幸、一体どうしてくれるというのか。

☆☆☆数学徒が馬鹿板をしたらダメ。さもないと国家議事堂みたいになります。☆☆☆



969:¥
17/06/17 12:35:00.60 x2f6D4gs.net


970:¥
17/06/17 12:35:20.71 x2f6D4gs.net


971:¥
17/06/17 12:35:41.80 x2f6D4gs.net


972:¥
17/06/17 12:36:00.84 x2f6D4gs.net


973:¥
17/06/17 12:36:21.78 x2f6D4gs.net


974:¥
17/06/17 12:36:42.46 x2f6D4gs.net


975:¥
17/06/17 12:37:03.93 x2f6D4gs.net


976:¥
17/06/17 12:37:23.71 x2f6D4gs.net


977:¥
17/06/17 12:37:59.24 x2f6D4gs.net


978:¥
17/06/17 12:38:23.00 x2f6D4gs.net


979:132人目の素数さん
17/07/23 16:07:33.29 AqYrJa+q.net
下のy式をもっとスマートに書く方法を教えてください
つまり、yの式は実際には+1, +2 …, +100までしなくてはいけないのです。
x <- 1:10
y <- c(sum((x+1)^2), sum((x+2)^2), sum((x+3)^2), …, sum((x+100)^2))

よろしくお願いします。

980:132人目の素数さん
17/07/23 16:33:37.21 FF49HrjI.net
z<-c()
for(i in 1:100){
z<-c(z, sum(x+i)^2))
}
でもできるけどダメ?

981:¥
17/07/23 16:35:16.97 r8UcZ6ry.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



982:132人目の素数さん
17/07/23 17:11:12.91 AqYrJa+q.net
>>943
ありがとうございます。一つ解決しました。
しかし、続いて教わった方法で、平均の最小2乗法による最小値を図示しようと思ったのですが
seq(0, 10, .01)を入れると、x軸がなぜか100倍されてしまします。
なぜなのか分かりませんか?おそらく、このforを使った解決策が
その後の応用(最小値を図示する)に適切でないのでしょうが
それがなんでなのか・・・(私はforについて、勉強し直します)

x <- c(3,5,6,8,2,9)
mean(x)

z<-c()
for(i in seq(from=0, to=10, by=.01)){
z<-c(z,sum(x-i)^2)
}
plot(z, type ="l") # x=5.5が最小のはずなのに、550が最小になっている

983:¥
17/07/23 17:45:58.37 r8UcZ6ry.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



984:132人目の素数さん
17/07/23 23:11:56.44 k8lcunzY.net
横軸の数値用ベクトルwにサンプリング数値を保存しておいて、
グラフ化時のx軸の値として指定してみてはどうか?

x <- c(3,5,6,8,2,9)
mean(x)

z<-c()
w<-c()
for(i in seq(from=0, to=10, by=.01)){
z<-c(z,sum(x-i)^2)
w<-append(w, i)
}
plot(w, z, type ="l")

985:132人目の素数さん
17/07/23 23:14:46.35 k8lcunzY.net
w<-append(w, i)

w<-c(w, i)
の方が簡潔で上との整合性も保てるか。

986:132人目の素数さん
17/07/24 00:56:57.68 7iMASKV5.net
>>942
z <- outer(x, 1:100, function(x, y) (x + y)^2)
apply(z, 2, sum)

987:132人目の素数さん
17/07/24 01:07:34.75 7iMASKV5.net
>>945
y軸の値しか指定しないから、x軸の値は1から始まる整数になる。
xv <- seq(0, 10, by=0.01)
z <- outer(x, xv, function(x, y) (x - y)^2)
yv <- apply(z, 2, sum)
plot(xv, yv, type="l")

988:132人目の素数さん
17/07/24 05:19:21.68 0E8vffzB.net
>>947-950
たくさんの返信ありがとうございます。
仕事から帰ったら、試してみます。

989:¥
17/07/24 09:36:16.88 /2fQM0DG.net


990:¥
17/07/24 09:36:35.46 /2fQM0DG.net


991:¥
17/07/24 09:36:55.11 /2fQM0DG.net


992:¥
17/07/24 09:37:12.58 /2fQM0DG.net


993:¥
17/07/24 09:37:31.25 /2fQM0DG.net


994:¥
17/07/24 09:37:48.24 /2fQM0DG.net


995:¥
17/07/24 09:38:05.80 /2fQM0DG.net


996:¥
17/07/24 09:38:23.56 /2fQM0DG.net


997:¥
17/07/24 09:38:48.61 /2fQM0DG.net


998:¥
17/07/24 09:39:05.54 /2fQM0DG.net


999:描
17/07/24 09:53:25.77 3peLAFj6.net
$

1000:¥
17/07/24 11:30:57.34 /2fQM0DG.net


1001:132人目の素数さん
17/07/24 17:15:25.78 RJfWHfZU.net
この熱いあぼーん達はなんなんだろ…

1002:¥
17/07/24 17:33:19.17 /2fQM0DG.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



1003:132人目の素数さん
17/07/24 21:41:07.75 J3GDpKdM.net
>>950
添削。
> system.time(for(i in 1:10){yv <- apply(z, 2, sum)})
ユーザ システム 経過
0.032 0.005 0.037
> system.time(for(i in 1:10){yv <- colSums(z)})
ユーザ システム 経過
0.002 0.000 0.001
ということで、apply()ではなくcolSums()を使うべき

1004:132人目の素数さん
17/07/25 00:40:01.16 w3drJ+gw.net
>>966
初心者には、速度を追及するよりは、応用の効くやりかたを学んでもらったほうが良いのでは?

1005:132人目の素数さん
17/07/25 00:41:44.53 BWjSvXL+.net
数学板はなんで荒らし投稿ばかりなの?
ロボットのしわざかな?
常連さんのPC�


1006:ェ乗っ取られているの?



1007:132人目の素数さん
17/07/25 06:50:18.42 7svw9fzv.net
>>947, 939
ありがとうございます。x軸の値の指定が間違っていたのですね。
y軸の値をx軸に指定していても図の形が同じなため気がつきませんでした。
しかし、y軸の値をx軸にあてはめても図の形が同じなのはなぜだろうか?
というか、y軸を指定しないでx軸だけ指定したのに、エラーが出ず
それなりの図(plot)が書けているのはなぜだ?
世の中分からんことばかりです。

>>949, 941, 957
outer()とapply()も実は使ったことなくて、いろいろいじってみて勉強します。
この分野は計算速度とか測って極めていくのか。
colSums()とか初めて見ました。

>>967
働き出してから、統計を勉強し直している初心者です。
ご迷惑おかけします。
基本的には論文(心理学・医療系)を理解するための勉強です。
ふと、このスレで、とんちんかんな質問をしていたら私かもしれません。

1008:¥
17/07/25 08:34:59.00 1OMr9h78.net


1009:¥
17/07/25 08:35:16.52 1OMr9h78.net


1010:¥
17/07/25 08:35:33.91 1OMr9h78.net


1011:¥
17/07/25 08:35:50.65 1OMr9h78.net


1012:¥
17/07/25 08:36:07.97 1OMr9h78.net


1013:¥
17/07/25 08:36:25.20 1OMr9h78.net


1014:¥
17/07/25 08:36:41.70 1OMr9h78.net


1015:¥
17/07/25 08:36:58.72 1OMr9h78.net


1016:¥
17/07/25 08:37:19.20 1OMr9h78.net


1017:¥
17/07/25 08:37:35.98 1OMr9h78.net


1018:132人目の素数さん
17/07/25 21:48:57.22 PpGdU5wz.net
>>969
x軸の値を指定せずにplotすると、横軸はデータのインデックス(順番)が表示される。
インデックスに対応する、本当の値をx軸の値としてplotしても、当然同じ図形になる。

1019:¥
17/07/25 22:08:54.95 1OMr9h78.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



1020:132人目の素数さん
17/07/25 22:22:42.60 ylT3F2qf.net
次スレは>>990あたり?

1021:¥
17/07/25 22:24:57.98 1OMr9h78.net
〒〒〒馬鹿板は悪い習慣であり、この行為は脳を悪くする。そやし足を洗いなさい。〒〒〒



1022:¥
17/07/26 05:47:15.01 e+L8tFA6.net


1023:¥
17/07/26 05:47:35.53 e+L8tFA6.net


1024:¥
17/07/26 05:47:52.90 e+L8tFA6.net


1025:¥
17/07/26 05:48:08.88 e+L8tFA6.net


1026:¥
17/07/26 05:48:25.31 e+L8tFA6.net


1027:¥
17/07/26 05:48:42.11 e+L8tFA6.net


1028:¥
17/07/26 05:48:59.88 e+L8tFA6.net


1029:¥
17/07/26 05:49:26.67 e+L8tFA6.net


1030:¥
17/07/26 05:49:42.86 e+L8tFA6.net


1031:¥
17/07/26 05:50:00.74 e+L8tFA6.net


1032:132人目の素数さん
17/08/02 19:58:43.15 CTkpktyp.net
自由度10のt分布に従う乱数を50個作りたい時は
rt(50,10)
でokでしょうか

1033:¥
17/08/02 20:01:22.80 6XC/wDpJ.net


1034:¥
17/08/02 22:21:30.33 6XC/wDpJ.net


1035:¥
17/08/02 22:21:48.18 6XC/wDpJ.net


1036:¥
17/08/02 22:22:05.33 6XC/wDpJ.net


1037:¥
17/08/02 22:22:21.58 6XC/wDpJ.net


1038:¥
17/08/02 22:22:39.25 6XC/wDpJ.net


1039:¥
17/08/02 22:22:57.70 6XC/wDpJ.net


1040:¥
17/08/02 22:23:16.14 6XC/wDpJ.net


1041:¥
17/08/02 22:23:36.18 6XC/wDpJ.net


1042:¥
17/08/02 22:23:53.86 6XC/wDpJ.net


1043:¥
17/08/02 22:24:11.53 6XC/wDpJ.net


1044:¥
17/08/02 22:24:30.81 6XC/wDpJ.net


1045:¥
17/08/02 22:24:47.55 6XC/wDpJ.net


1046:¥
17/08/02 22:25:05.15 6XC/wDpJ.net


1047:猫
17/08/02 22:25:35.27 6XC/wDpJ.net


1048:1001
Over 1000 Thread.net
このスレッドは1000を超えました。
もう書けないので、新しいスレッドを立ててくださいです。。。
life time: 1406日 9時間 18分 13秒

1049:過去ログ ★
[過去ログ]
■ このスレッドは過去ログ倉庫に格納されています


最新レス表示
レスジャンプ
類似スレ一覧
スレッドの検索
話題のニュース
おまかせリスト
オプション
しおりを挟む
スレッドに書込
スレッドの一覧
暇つぶし2ch