- PDF Free Download

Size: px

Start display at page:

Download ""

ふじよしこいたばし
5 years ago
Views:

7 :46:18

9 ##### Ruby ##### h = Hash.new(0) while line = gets() if line =~ /^.*?\x01(.*?)\x01/ h[$1] += 1 end end h.sort_by { k,v - v }.each { k,v print v, "\t", k, "\n" } ##### R ##### data = read.delim("...", header=false) plot(1: , data[,1], type="o", pch=16, log="xy", xlab="", ylab="")

10 PersonFinder_Mi (15635 ) PersonFinder_ot PersonFinder_Iw PersonFinder_Fk JorudanLive

11 JorudanLive PersonFinder_Mi PersonFinder_ot PersonFinder_Iw PersonFinder_Fk

12 PersonFinder_Na PersonFinder_Mi PersonFinder_Iw PersonFinder_Fk PersonFinder_Ib PersonFinder_Ao PersonFinder_ot

16 @NHK_PR

18 ./histo.rb ? >histo.out./histo3b.rb histo.out >histo3b.out data = read.csv("histo3b.out", header=false) t = as.posixct(data[,1]) n = data[,2] par(las=1) # las: label style par(mgp=c(2,0.8,0)) # : c(3,1,0) r = range(t) plot(t, n/10000, type="h", xlab="", ylab="", ylim=c(0,360), xaxs="i", yaxs="i", xlim=c(r[1]- 1800,r[2]+1800), xaxt="n", lwd=4, lend=1, col="#0068b7") r = as.posixct(c(" "," ")) br = seq(r[1], r[2], by="days") axis.posixct(1, br, at=br, format="%m/%d") axis(2, labels=" ", at=350, padj=- 1, tick=false) title(" ", line=0.5, cex.main=2)

20 data = read.csv("histo.out", header=false) t = as.posixct(data[,1]) n = data[,2] par(las=1) # las: label style par(mgp=c(2,0.8,0)) # : c(3,1,0) plot(t, n/10000, main="", type="l", xlab="", ylab="", xaxt="n") r = as.posixct(c(" "," ")) br = seq(r[1], r[2], by="days") axis.posixct(1, at=br, format="%m/%d") axis(2, labels=" ", at=10, padj=- 1, tick=false) title(" ", line=0.5, cex.main=2)

31 lat =... # lon =... # library(rgooglemaps) source(" u.ac.jp/~okumura/stat/data/getmap.r") MyMap = GetMap(c(35,140), zoom=5, sensor="false", hl="ja", maptype="roadmap") tmp = PlotOnStaticMap(MyMap, lat=lat, lon=lon, pch=16, col="#ff2800")

32 03/11 15:00 16:00

33 data =... # t = as.posixct(data[,1]) x = data[,2] plot(t, x/ , type="o", pch=16, xaxt="n", xlab="", ylab="") axis.posixct(1, at=round(t,"days"), format="%m/%d") axis(2, labels=" ", at=50, padj=- 1, tick=false)

34 15:00 18:00 07:00

35 t0 = as.posixct(" ") t1 = as.posixct(" ") t2 = as.posixct(" ") t3 = as.posixct(" ") plot(t[t>=t0 & t<=t1], x[t>=t0 & t<=t1]/ , type="o", pch=16, xaxt="n", xlab="", ylab="", ylim=range(x/ )) points(t[t>=t0 & t<=t1], x[t>=t1 & t<=t2]/ , type="o", pch=16, col="#f39800") points(t[t>=t0 & t<=t1], x[t>=t2 & t<=t3]/ , type="o", pch=16, col="#0068b7") axis.posixct(1, t[t >= t0 & t <= t1], format="%h:%m") axis(2, labels=" ", at=50, padj=- 1, tick=false)

40 x =... # par(mgp=c(2,0.8,0)) h = hist(x, right=f, plot=f, breaks=0:58781)$counts plot(0:1000, h[1:1001]/10000, type="h", xlab=" ", ylab="", ylim=range(c(0,h/10000))) axis(2, labels=" ", at=30, padj=- 1, tick=false)

42 plot(0:5000, h[1:5001], type="l", pch=1, log="y", xlab="", ylab="")

54 Cs- 134 (605keV) Cs- 137 (662keV) Cs- 134 (796keV) K- 40 (1461keV)

55 K- 40 I keV 81.5% keV 7.16% Cs keV 85.10% Cs keV 8.338% keV % keV 97.62% keV 85.46% keV 8.688% keV 10.66% Bi keV 45.49% keV 4.894% keV 3.107% keV 14.92% keV 5.834% keV 3.988% keV 15.30% Pb keV 7.251% keV 18.42% keV 35.60%

57 Cs- 134 (605keV) Cs- 137 (662keV) Cs- 134 (796keV)

59 Cs- 137 Cs k

63 ch = 0:1023 sp = scan("spectrum.txt") library(nlme) data = data.frame(ch=ch, sp=sp) fit = gnls(sp ~ (c + d*(ch- 350) + b*(ch- 350)^2 + (e1*( /0.9226)*dnorm((ch *m1)/(0.9226*s)) + e1*( /0.9276)*dnorm((ch *m1)/(0.9276*s)) + e1*(0.9762/0.956)*dnorm((ch- m1)/(0.956*s)) + e2*0.851*dnorm((ch- m2)/s) + e1*(0.8546/1.0967)*dnorm((ch- m3)/(1.0967*s)) + e1*( /1.101)*dnorm((ch *m3)/(1.101*s)) ) / s ) * (ch/350)^u, data=data, subset=260:350, start=list(c=200,d=- 3,b=0.02,e1=10000,e2=10000, m1=313,m2=342,m3=407,s=9,u=- 1.2), weights=varpower(fixed=0.5), control=list(nlstol=1e- 4))

71 data = read.csv( " u.ac.jp/~okumura/stat/data/fukushima csv", colclasses=c("posixct","numeric","numeric","numeric","numeric")) attach(data) # plot(,, type="l", xlab="", ylab="", xaxt="n", col="#f39800") axis(2, "μsv/h", at=0.75, padj=- 1) axis.posixct(1,, format="%y %m ") # plot(,, type="l", xlab="", ylab="", xaxt="n", col="#009944") axis(2, " ", at=35, padj=- 1) axis.posixct(1,, format="%y %m ") # cols = colorramp(c("#ff8000","#0080ff")) t = (as.numeric( ) - as.numeric( [1])) / (as.numeric( [1464]) - as.numeric( [1])) plot(,, type="p", pch=16, xlab=" ", ylab=" μsv/h ", col=rgb(cols(t)/255))

73 > lm( ~ + ) (Intercept) 1.280e e e- 09 # plot(,, type="l", xlab="", ylab="", xaxt="n", col="#f39800") axis(2, "μsv/h", at=0.75, padj=- 1) axis.posixct(1,, format="%y %m ") points(, *( - 20), type="l", col="#0068b7")

78 2F

1 R Windows R 1.1 R The R project web R web Download [CRAN] CRAN Mirrors Japan Download and Install R [Windows 9

1 R Windows R 1.1 R The R project web R web Download [CRAN] CRAN Mirrors Japan Download and Install R [Windows 9 1 R 2007 8 19 1 Windows R 1.1 R The R project web http://www.r-project.org/ R web Download [CRAN] CRAN Mirrors Japan Download and Install R [Windows 95 and later ] [base] 2.5.1 R - 2.5.1 for Windows R