1 2 S-PLUS python MeCab csv S-PLUS csv S- PLUS S-PLUS csv python word1,word2,word3 2 python MeCab MeCab ipadic html 2010/3/26

Size: px

Start display at page:

Download "1 2 S-PLUS python MeCab csv S-PLUS csv S- PLUS S-PLUS csv python word1,word2,word3 2 python MeCab MeCab ipadic html 2010/3/26"

まいえふじつぐ
5 years ago
Views:

1 A 9 B 10 C 10 D 11 E 12 1

2 1 2 S-PLUS python MeCab csv S-PLUS csv S- PLUS S-PLUS csv python word1,word2,word3 2 python MeCab MeCab ipadic html 2010/3/ /10/ URL 3. MeCab % 50% 2

3 3 1: plot lines locator(1) 3

4 : : : 4

5 4.2 (f g)[n] def = m= f [m]g[n + m] (1) n m (lag) S- PLUS acf 2?? 3:??

6 4.3 Lag lag 0 4: 6

7 5 5 5: 5 2:

8 6 [1],S-PLUS, [2] JIN S PAGE, [3] : :, 83%AC%3A%E3%82%A4%E3%82%AB%E5%A8%98%E3%81%8B%E3%82%8F%E3%81%84%E3%81%84 8

9 A zikeiretu <- structure(c(985, 1812, 1624, 2688, 7683, 2592, 13899, , , , 85159, , 68796, 53159, 54207, 39708, 86626, 52512, , ),.Names = c( 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 )) par(xaxt= n ) plot(zikeiretu,type= p,ylab=,xlab=,sub= ) lines(zikeiretu,col=3) par(xaxt= s ) axis(1,at = c(1:20)) from <- locator(1) to <- locator(1) arrows(from$x,from$y,to$x,to$y) text(from, ) from <- locator(1) to <- locator(1) arrows(from$x,from$y,to$x,to$y) text(from, ) from <- locator(1) to <- locator(1) arrows(from$x,from$y,to$x,to$y) text(from, ) from <- locator(1) to <- locator(1) arrows(from$x,from$y,to$x,to$y) text(from, ) 9

10 B par(xaxt= n ) plot(c(1:20),ika$word90,type= p,col=1,xlab=,ylab=,sub= ) par(new=t) plot(c(1\maketitle20),ika$word90,type= l,col=1,ylab=,xlab= ) par(new=t) plot(c(1:20),ika$word29,type= p,col=3,xlab=,ylab=,labels=false,axes=false) par(new=t) plot(c(1\maketitle20),ika$word29,type= l,col=3,xlab=,ylab=, labels=false,axes=false) axis(side=4,col=3) par(xaxt= s ) axis(1,at = c(1:20)) legend(locator(1),c(, ),lty=1,col=c(1,3)) C myccf <- function(a,b,mainname){ aa <- ts(a) bb <- ts(b) x <- ts.intersect(aa,bb) xx <- acf(x, plot=f) dat <- c(rev(xx$acf[-1, 1, 2]), xx$acf[, 2, 1]) d <- data.frame(lag=round(c(rev(xx$lag[-1, 1, 2]), xx$lag[, 2,1]),3),acf=dat) plot(d$lag,d$acf,type="n",ylab="acf",xlab="lag",main=mainname) abline(h=0) for ( i in 1:length(d$acf)){ } } lines(c(d$lag[i],d$lag[i]),c(0,d$acf[i])) return(d) 10

11 D wariai100 <- read.table("100_wariais.txt",sep=",",head=t) sortlist <- order(wariai100$word101,decreasing=false) wariai100 <- wariai100[sortlist,] #rownames(wariai100) <- c(1:nrow(wariai100)) for ( i in 2:101){ for ( j in 2:101){ a <- myccf(wariai100[i],wariai100[j],mainname = paste(names(wariai100[i]),"&",names(wariai100[j]))) if ( max(a$acf[c(10,11,12)]) >= 0.8 && i!= j){ cat(paste(names(wariai100[i]),"&", names(wariai100[j]),"\n")) } if ( min(a$acf[c(10,11,12)]) <= -0.8 && i!= j){ cat(paste(names(wariai100[i])," ", names(wariai100[j]),"\n")) } } } par(mfrow=c(2,2)) myccf(wariai100$word10,wariai100$word58,mainname=" ") myccf(wariai100$word54,wariai100$word73,mainname=" ") myccf(wariai100$word100,wariai100$word68,mainname=" ") myccf(wariai100$word89,wariai100$word29,mainname=" ") 11

12 E wariai100 <- read.table("100_wariais.txt",sep=",",head=t) sortlist <- order(wariai100$word101,decreasing=false) wariai100 <- wariai100[sortlist,] par(mfrow=c(2,2)) ika.pc<-princomp(wariai100[1:100],cor=f) barplot(ika.pc$scores[,1],main=,names=paste(1\maketitle20)) barplot(ika.pc$scores[,2],main=,names=paste(1:20)) sortlist <- order(ika.pc$coef[,1],decreasing=false) dat <- ika.pc$coef[,1] barplot(c(dat[sortlist[100:96]],dat[sortlist[5:1]]),las=2,names=c(,,,,,,,,, ),main= 1 ) sortlist <- order(ika.pc$coef[,2],decreasing=false) dat <- ika.pc$coef[,2] barplot(c(dat[sortlist[100:96]],dat[sortlist[5:1]]),las=2,names=c(,,,,,,,,, ),main= 2 ) 12