ybenjo

library('ggplot2')

data <- read.table("/tmp/ret.csv", header = T, sep = ",")
data$topic_id <- as.factor(data$topic_id)
qplot(timestamp, gamma, data = data, color = topic_id)

ybenjo

# -*- coding: utf-8 -*-
require 'nokogiri'
require 'open-uri'
require 'logger'

$log  = Logger.new('./crawl.log')

def get(num)
  ret = [ ]
  url = "http://shasetsu.ps.land.to/index.cgi/event/#{num}/"

ybenjo

ybenjo

ybenjo

# -*- coding: utf-8 -*-
require "kconv"
io = IO.popen("mecab", "r+")

io.puts "吾輩は猫である"
io.each{|e|
  puts e
  break if e.chomp == "EOS"
}

ybenjo

> # 決定木を使いたい
> library(mvpart)
> 
> # 本来なら直接
> # rpart(iris$Species~., iris)
> # とやりたいけど iris の予測すべきラベル名がわからず
> # データフレームの最後にあることだけはわかる，みたいな状況
> 
> # ここで最終行の列名取得
> label <- colnames(iris)[ncol(iris)]

ybenjo

釘宮理恵,200007,0.31963575947887185
釘宮理恵,200010,0.31578446348916184
釘宮理恵,200012,0.31564458559250186
釘宮理恵,200101,0.31531536347356165
釘宮理恵,200102,0.3152865223926949
釘宮理恵,200103,0.31479491969312356
釘宮理恵,200104,0.5039254644570497
釘宮理恵,200105,1.0558031621952395
釘宮理恵,200107,1.005958003523609
釘宮理恵,200110,0.9323008071023162

ybenjo

釘宮理恵 さんは 2001年4月頃にブレイクした…?
  0.31479491969312356 => 0.5039254644570497, sd: 0.0014502546594198278
釘宮理恵 さんは 2001年5月頃にブレイクした…?
  0.5039254644570497 => 1.0558031621952395, sd: 0.07085797611505189
釘宮理恵 さんは 2004年10月頃にブレイクした…?
  3.352655091760327 => 5.818529948330466, sd: 0.1639852185544332
釘宮理恵 さんは 2006年4月頃にブレイクした…?
  6.181845671016907 => 8.363148049010016, sd: 0.33924422910288793
子安武人 さんは 1994年4月頃にブレイクした…?
  1.238813633722919 => 2.0370318925579873, sd: 0.14508261990005888

ybenjo

(defn -main [& [args]]
  (println args))

ybenjo

# -*- coding: utf-8 -*-
# Practice of Viterbi Algorithm
# Data from
# http://www.chokkan.org/lectures/2011nlp/03nlp.pdf

class HMM
  def initialize
    @all_x = [ ]
    @all_y = [ ]
    @n_x_y = Hash.new{0.0}