4/13の進捗 [反省]
チャロ
朝1 奨学金
統計
山木さんのレジメ
中西ブログ
解析で天体と近くの星までの距離を求める
#1シグマを下記加えたグラフ
kyori.ah<-( (data.set$x.H.1-data.set$x2)^2+(data.set$y.H.1-data.set$y2)^2 )^(1/2)
plot(data.set$time.A,kyori.ah,ylab='distance A-H',xlab='time(JD)')
abline(a=mean(kyori.ah),col='red',b=0)
abline(a=mean(kyori.ah)+sd(kyori.ah),col='green',b=0)
abline(a=mean(kyori.ah)-sd(kyori.ah),col='green',b=0)
#1シグマ以上のものを除いた上でΧ2乗を求める
data<- data.frame(data.set, kyori.ah)
###data.1sig.nに星Aと星Hの距離が1シグマのばらつきで収まっているものだけを入れる。
##data.1sig.n <-subset(data, mean(kyori.ah)+sd(kyori.ah)>data$kyori.ah & mean(kyori.ah)-##sd(kyori.ah) data$kyori.ah & mean(kyori.ah)-3*sd(kyori.ah)
朝1 奨学金
統計
山木さんのレジメ
中西ブログ
解析で天体と近くの星までの距離を求める
#1シグマを下記加えたグラフ
kyori.ah<-( (data.set$x.H.1-data.set$x2)^2+(data.set$y.H.1-data.set$y2)^2 )^(1/2)
plot(data.set$time.A,kyori.ah,ylab='distance A-H',xlab='time(JD)')
abline(a=mean(kyori.ah),col='red',b=0)
abline(a=mean(kyori.ah)+sd(kyori.ah),col='green',b=0)
abline(a=mean(kyori.ah)-sd(kyori.ah),col='green',b=0)
#1シグマ以上のものを除いた上でΧ2乗を求める
data<- data.frame(data.set, kyori.ah)
###data.1sig.nに星Aと星Hの距離が1シグマのばらつきで収まっているものだけを入れる。
##data.1sig.n <-subset(data, mean(kyori.ah)+sd(kyori.ah)>data$kyori.ah & mean(kyori.ah)-##sd(kyori.ah) data$kyori.ah & mean(kyori.ah)-3*sd(kyori.ah)
2010-09-15 01:09
nice!(0)
コメント(0)
トラックバック(0)
コメント 0