【R言語】統計解析フリーソフトＲ 第４章【GNU R】
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2011/01/09(日) 17:31:49統計計算で重宝するデータ型や、複数要素を処理する演算や関数、
解析結果を表示するグラフィックなど、多彩な機能を提供します。
●関連サイト
The R Project
http://www.r-project.org/
RjpWiki
http://www.okada.jp.org/RWiki/
リンク集
http://www.okada.jp.org/RWiki/?%A5%EA%A5%F3%A5%AF%BD%B8
0696１３２人目の素数さん
2012/12/25(火) 21:10:27.10[,1] [,2] [,3]
[1,] 1 4 7
[2,] 1 4 8
[3,] 1 5 7
[4,] 1 5 8
[5,] 1 6 7
[6,] 1 6 8
[7,] 2 4 7
[8,] 2 4 8
[9,] 2 5 7
[10,] 2 5 8
[11,] 2 6 7
[12,] 2 6 8
0697695
2012/12/25(火) 21:25:57.68それなら、expand.grid()を再帰的に使えば良い。
入れ子にするのは面倒なので二項演算子を定義すると、
> "%+%" <- function(x, y) apply(expand.grid(x, y), 1, paste, collapse = "")
> a[[1]] %+% a[[2]] %+% a[[3]]
[1] "147" "247" "157" "257" "167" "267" "148" "248" "158" "258" "168" "268"
あとはRecall()を使って、これを一般化すれば良い。
> f696 <- function(x) { f <- function(i){ if(i > 0) return(Recall(i - 1) %+% x[[i]])
+ else return("")}
+ return(f(length(x)))
+ }
> f696(a)
[1] "147" "247" "157" "257" "167" "267" "148" "248" "158" "258" "168" "268"
次に質問するときは、最初から>>696のように書いた方が良いよ。
または、
> forを3重にすればできるのですが、
この具体的なコードを書くとか。抽象的に説明すると、解釈の余地に幅が生まれる。
0698１３２人目の素数さん
2012/12/25(火) 21:29:26.760699695
2012/12/25(火) 21:31:46.70n桁の整数は連結させるのではなく、
桁をばらばらにして行列にしたかったのか。
それなら話はもっと簡単だ。
pasteの処理をやめれば良い
0700１３２人目の素数さん
2012/12/25(火) 22:31:04.88(function() as.matrix((expand.grid(a)->x)[do.call(order, x),]))()
0701１３２人目の素数さん
2012/12/26(水) 09:37:32.22あぁ、しまった。expand.grid()は3次元次元でも使えたのか。
>>700さんの別解。rev()を2重に使う。
as.matrix(rev(expand.grid(rev(a))))
0702令嬢
2012/12/26(水) 21:17:40.81, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/
. | ＼ ∠イ ,ｲイ| ,`-' |
| l＾,人| ` `-' ゝ | このスレには馬と鹿と豚さんしかいないのね。
| ` -'＼ ｰ' 人
| /(l ＿＿／ ヽ、
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0703１３２人目の素数さん
2012/12/26(水) 23:44:03.23revの二重使いは柔軟な思考ですね。参考になりました。
0704狢 ◆yEy4lYsULH68
2012/12/27(木) 09:03:34.27>増田哲也こそ笑い者。
>俺が逮捕されて懲戒免職させる日本こそ沈めって、一発逆転をねらっている愚民そのもの。
>
0705１３２人目の素数さん
2012/12/29(土) 23:41:13.98SASでQ01-Q50などとやるようなことをしたいのですが、方法がわかりません。ご存知の方、お教えください。
50項目で主成分分析をしようなどという場合を考えると、気が遠くなります。
0706１３２人目の素数さん
2012/12/30(日) 07:20:47.800707１３２人目の素数さん
2012/12/30(日) 08:39:23.19outcome ~ .
って指定するとすべての変数が入る
0708１３２人目の素数さん
2012/12/30(日) 11:33:43.75できるけど、"Q01-Q50"のような記法ではできない。
>>707が助言するようにまるごと入れることもできるし、
paste()とas.formula()をつかって長いformulaを作成することもできる。
formulaのヘルプにある例
> xnam <- paste0("x", 1:25)
> (fmla <- as.formula(paste("y ~ ", paste(xnam, collapse= "+"))))
y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9 + x10 + x11 +
x12 + x13 + x14 + x15 + x16 + x17 + x18 + x19 + x20 + x21 +
x22 + x23 + x24 + x25
0709１３２人目の素数さん
2012/12/30(日) 12:25:38.54ありがとうございます。
これで、エクセルのオートフィルで変数リスト作成してRのコンソールに貼り付ける作業から解放されます。
0710１３２人目の素数さん
2012/12/30(日) 14:53:50.700711１３２人目の素数さん
2012/12/30(日) 14:55:00.08その中で、irisデータを用いてコマンドの学習をしているのですがエラーが出てしまい先に進めなくなってしまいました。
以下にコマンドを載せますものが私の用いたものです。
library(rpart)
library(adabag)
data <- iris
ndata <- nrow(data)
#乱数指定
set.seed(101)
#学習データ(data.learn)とテストデータ(data.test)に分ける
ridx <- sample(ndata, ndata * 0.5)
data.learn <- data[ridx,]
data.test <- data[-ridx,]
#3-fold crossvalidation、弱識別器の数を10とし、学習データに対しboostingを行ったものをdata.adaCvに代入
data.adaCv <- boosting.cv(Species ~ .,data = data.learn,v=3, mfinal = 10)
#テストデータに学習データを照らし合わせる。
resultPredict <- predict(data.adaCv, newdata = data.test, type="class")
というコマンドですが「predict」を行った際に
以下にエラー UseMethod("predict") :
'predict' をクラス "character" のオブジェクトに適用できるようなメソッドがありません
というエラーが出てしまい、ネットや参考文献を見てもわからず投稿させていただきました。
もしこのエラーの解決方法がお分かりのかたがいらっしゃいましたらよろしくお願いします。
0712695
2012/12/30(日) 20:21:19.02>以下にエラー UseMethod("predict") :
>'predict' をクラス "character" のオブジェクトに適用できるようなメソッドがありません
分かるも何も、書いてあるままだよ。
> class(data.adaCv)
[1] "list"
boosting.cv()は、predictメソッドをもつクラスを吐かず、単純にlistクラスを吐く。
それだけの話。
0713１３２人目の素数さん
2012/12/30(日) 20:22:36.64名前欄にゴミが入ったorz
0714１３２人目の素数さん
2012/12/30(日) 23:08:32.14実行してみたところおっしゃたように"list"と出力されました。
ということはboosting.cvで得られた変数(data.adaCv)を用いての
predictはできないということでしょうか？
また、もしできるのであればそのコマンドと使用例を教えていただきたいです。
宜しくお願いします。
0715１３２人目の素数さん
2012/12/31(月) 15:27:18.84, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/
. | ＼ ∠イ ,ｲイ| ,`-' |
| l＾,人| ` `-' ゝ | このスレには馬と鹿と豚さんしかいないのね。
| ` -'＼ ｰ' 人
| /(l ＿＿／ ヽ、
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0716１３２人目の素数さん
2013/01/02(水) 11:10:09.30プログラミング言語の解説というかRをアプリケーションとしてみた解説って感じで
言語的なことはほとんど勉強になりませんでしたね。
やっぱ俺はS言語の本読むべきでしたね。
0717１３２人目の素数さん
2013/01/04(金) 09:40:38.25ttp://www.amazon.co.jp/dp/4873115795
ttp://www.amazon.co.jp/dp/4621061313
0718１３２人目の素数さん
2013/01/04(金) 10:07:16.37Rjpwikiの方に書いたからもういいと思うのだけど、念のため。
boosting.cvの中身は見た？
bootsting()とpredict.boosting()を中で使っているよね。
> boosting.cv
function (formula, data, v = 10, boos = TRUE, mfinal = 100, coeflearn = "Breiman",
control)
{
vardep <- data[, as.character(formula[[2]])]
n <- length(vardep)
if (v > n)
stop(" v should be in [2, n]")
if (v < 2)
stop(" v should be in [2, n]")
predclass <- rep("O", n)
for (i in 1:v) {
test <- v * (0:floor(n/v)) + i
test <- test[test < n + 1]
fit <- boosting(formula, data[-test, ], boos, mfinal,
coeflearn, control = control)
fit.predict <- predict.boosting(fit, data[test, ])
predclass[test] <- fit.predict$class
cat("i: ", c(i, date()), "\n")
}
tabla <- table(predclass, vardep, dnn = c("Predicted Class",
"Observed Class"))
error <- 1 - sum(predclass == vardep)/n
output <- list(class = predclass, confusion = tabla, error = error)
}
<environment: namespace:adabag>
このfit.predictが欲しいなら、自分でbooting.cv()をコピーして編集すればよい。
そうじゃないなら、何をしたいのか、もっとはっきり伝えること。
0719１３２人目の素数さん
2013/01/05(土) 03:55:18.11Rノウハウ本の購入層はやりたい解析があってそこまでたどり着ければいい
って人が大多数なんじゃないかな、んでそれにあわせた内容になってると
＋αを求める人は多くは無いと思う
0720１３２人目の素数さん
2013/01/05(土) 12:25:03.16質問です。
x<-c(1,2,3,4,5,6)
y<-c(1,4,9,7,6,8)
というデータがあったとして、
sp<-smooth.spline(x,y)
predict(sp,x=2.5)
とすると、xが2.5のときのyの値は下記のように予測できます。
$x
[1] 2.5
$y
[1] 6.923206
逆にyからxを求める方法がお分かりの方、いらっしゃいましたらお教えください。
0721１３２人目の素数さん
2013/01/05(土) 14:02:30.39stats:::predict.smooth.spline.fitにざっと目を通した限りは、
できなさそうだけど。
GAMFIT (http://lib.stat.cmu.edu/general/gamfit)の中までは読んでいない。
sp$fitの中身を使って数学的に解くか、xを細かい間隔でpredict()に与えて、
与えられたyに最も近いものを選ぶしかなさそうだが、
ちゃんとした識者に答えてもらった方がよいかも。
後者なら、
predict.y <- function(y){
x0 <- c(1, 2, 3, 4, 5, 6)
y0 <- c(1, 4, 9, 7, 6, 8)
sp <- smooth.spline(x0, y0)
i <- seq(from = min(x0), to = max(x0), by = 0.0001)
Y <- predict(sp, i)
j <- which.min((Y$y - y)^2)
return(list(y = y, x = i[j], fitted = Y$y[j]))
}
こんな感じでよいのでは。
> predict.y(2.345)
$y
[1] 2.345
$x
[1] 1.4487
$fitted
[1] 2.345004
0722１３２人目の素数さん
2013/01/05(土) 14:47:31.66, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/ 私は只の数ヲタなんかとは付き合わないわ。
. | ＼ ∠イ ,ｲイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。
| l＾,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。
| ` -'＼ ｰ' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて
| /(l ＿＿／ ヽ、 良い論文の出版を遅らせるお馬鹿な人。
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、 あなたの論文が Ann of Math に accept される確率は?
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼ それとも最近は Inv. Math. の方が上かしら？
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0723１３２人目の素数さん
2013/01/05(土) 16:46:53.71ありがとうございました。参考になりました。
0724１３２人目の素数さん
2013/01/05(土) 20:24:09.52, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/
. | ＼ ∠イ ,ｲイ| ,`-' |
| l＾,人| ` `-' ゝ | このスレには馬と鹿と豚さんしかいないのね。
| ` -'＼ ｰ' 人
| /(l ＿＿／ ヽ、
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0725１３２人目の素数さん
2013/01/07(月) 16:01:51.28いまさらだけどビギニングス ガイド イン Rのミススペルな。
0726１３２人目の素数さん
2013/01/07(月) 18:48:16.16ttp://www.amazon.com/dp/0387938362
これのことか？
でも「ビギニングス ガイド イン R」じゃないぞ。
0727１３２人目の素数さん
2013/01/09(水) 09:49:00.54正確には
ビギニングス ガイド イン R
でした。
0728１３２人目の素数さん
2013/01/11(金) 23:37:12.84コード 日付1 日付2 日付3 データ1 データ2 データ3
0001 201201 201103 201003 52 33 34
0002 201203 NA NA 22 NA 3
0003 201103 201001 200901 32 2 73
0004 201003 200905 200804 98 87 88
0005 NA NA NA NA NA NA
0006 201105 201102 201009 NA 32 55
日付XとデータXが対応
日付1>日付2>日付3となっている
この時、日付が201101よりも大きい日付に対応するNAでないデータを取得したい
即ち望まれる出力は以下の通り
コード 日付 データ
0001 201201 52
0002 201203 22
0003 201103 32
0004 NA NA
0005 NA NA
0006 201102 32
このように取得したい場合どのような操作をすればよいでしょうか?
実際には数千のコードがあるのでループで一行ずつ回すのはあまり現実的ではありません。
よろしくお願いします
0729１３２人目の素数さん
2013/01/12(土) 03:32:35.10宿題じゃないよね？
> a <- scan(text="0001 201201 201103 201003 52 33 34
+ 0002 201203 NA NA 22 NA 3
+ 0003 201103 201001 200901 32 2 73
+ 0004 201003 200905 200804 98 87 88
+ 0005 NA NA NA NA NA NA
+ 0006 201105 201102 201009 NA 32 55")
Read 42 items
> b <- data.frame(t(matrix(a, ncol = 6)))
> names(b) <- c("コード", paste0("日付", 1:3), paste0("データ", 1:3))
で、ヒントは次の通り。
> m <- ((b[, 2:4] > 201101) & (b[, 2:4] & b[, 5:7])) * b[, 2:4]
> m[(m == 0)] <- NA
> i <- apply(m, 1, which.max)
> unlist(sapply(1:length(i), function(j){ifelse(length(i[j][[1]]) == 0, NA, (b[j, 5:7])[i[j][[1]]])}))
[1] 52 22 32 NA NA 32
> unlist(sapply(1:length(i), function(j){ifelse(length(i[j][[1]]) == 0, NA, (b[j, 2:4])[i[j][[1]]])}))
[1] 201201 201203 201103 NA NA 201102
0730１３２人目の素数さん
2013/01/13(日) 13:37:10.99確認したいのですが、↓だと何故か２になってしまいます。
どこがダメなんだしょう？
x = rnorm(500)
hist(x, prob=T)
sum(hist(x, prob=T)$density)
0731１３２人目の素数さん
2013/01/13(日) 13:52:12.00乱数の分散を変えてみると、理解出来るのでは？
> set.seed(1234); sum(hist(rnorm(500, sd = 2), prob = TRUE)$density)
[1] 1
> set.seed(1234); sum(hist(rnorm(500, sd = 5), prob = TRUE)$density)
[1] 0.2
> set.seed(1234); sum(hist(rnorm(500, sd = .1), prob = TRUE)$density)
[1] 20
0732１３２人目の素数さん
2013/01/13(日) 15:13:55.16わからないよ (´・ω・｀)
sum(hist(x, prob=T)$density)
と
sum(hist(x, prob=T)$count/500)
は同じになると思ったんだけど、違うってことですかね？？
0733１３２人目の素数さん
2013/01/14(月) 13:24:07.24730 の最後の式では、確率密度の積分 ∫f(x) dx = 1 の dx の部分が
抜けてしまっています。
hist を用いる場合、dx に相当する値は
diff(hist(x, prob=T)$breaks) で求めることができます。
x <- rnorm(500)
h <- hist(x, prob=T)
sum(h$density * diff(h$breaks))
0734１３２人目の素数さん
2013/01/14(月) 22:54:04.42ありがとうございます！
0735１３２人目の素数さん
2013/01/16(水) 12:15:17.92言語を英語にすると、立ち上げた際のメッセージなどはきちんと英語になるのですが、
cor.test()などの関数を実行すると出力が日本語メッセージでかえってきます
日本語のwinOSでRの返答を完全に英語にする方法をどなたかご存知ないでしょうか？
ちなみに↓サイトの起動アイコンカスタマイズの方法で英語化しています
ttp://minato.sip21c.org/swtips/R.html
0736１３２人目の素数さん
2013/01/16(水) 12:41:27.37妙な挙動ですね。WindowsだとLANG=Cがどこかで強制的に書き換えられるのでしょうか。
Sys.getenv()で、メッセージが日本語になる前後でLANGやLC_ALLを確認してみてはどうでしょうか。
LANGを確認するには、Sys.getenv("LANG")とすればOKです。
なお、メッセージの言語を指定する環境変数はLC_MESSAGESです。
環境変数によるローケルの管理を知っておいても損はないと思います。
例えば、ttp://www.itmedia.co.jp/enterprise/articles/0605/16/news027.html
の「環境変数の階層」とかに目を通すとか。
0737736
2013/01/16(水) 12:54:26.20ttp://www.itmedia.co.jp/enterprise/articles/0605/16/news027.html
を読んで思い出したが、昔、起動アイコンに、LANG=Cじゃなくて、language=enをつけていた。
「R language=en」でググるとCRANのR for Windows FAQがヒット。
ttp://cran.r-project.org/bin/windows/base/rw-FAQ.html
... if you want to ensure that menus and messages are in (American) English, LANGUAGE=en.
とか、
You can ensure that R uses English messages by appending LANGUAGE=en to the shortcut
you use to start R, or setting it in the Rconsole file.
って書いてあるよ。FAQに。
0738１３２人目の素数さん
2013/01/16(水) 23:32:51.99遅くなりましたがありがとうございました。
無事動作しました。
0739728
2013/01/16(水) 23:41:01.80再び質問です。
横型ではなく縦型の場合はいかがでしょうか?
以下のようなデータフレームを仮定します
コード 日付 データ
0001 20120101 332
0001 20120111 542
0002 20120201 5554
0002 20120301 52334
0002 20120511 15234
0003 20121101 35234
0003 20121201 NA
0004 20120501 75234
0005 20120601 NA
各コードに対してデータがNAでない日付が最大の値を取得したいのです。
NAしかない場合にはその値を取得します。
上記の場合以下の値が期待されます。
コード 日付 データ
0001 20120111 542
0002 20120511 15234
0003 20121101 35234
0004 20120501 75234
0005 20120601 NA
ヒントだけでも結構ですのでよろしくお願いいたします。
0740１３２人目の素数さん
2013/01/17(木) 00:37:07.40ローカルに保存したCSVファイルではなく、直接ダウンロード、読込みをするファイルの文字化けを解消したいのですが、
エンコードを変える記述では上手くいきません。
ご教示いただけないでしょうか？
例)
url.data <- "http://maxis.muam.jp/e/fund/download/";
eMAXIS <- read.table(url.data, sep=",", stringsAsFactor = FALSE)
0741１３２人目の素数さん
2013/01/17(木) 14:18:27.23fileEncoding = "SJIS"
オプションとかではダメってこと？
0742１３２人目の素数さん
2013/01/17(木) 14:25:14.37plyrにあるddply()でIDで分割してwhich.max()で日付が最大のところのindexをとって、その場所の数字を取ればいいかな。
NAの処理を条件分岐を使わないでやるなら、ID日付でソートしたうえで、ddply()の中でzooのna.locf()を使えばいいか？
結構重くなってしまうかも。
0743１３２人目の素数さん
2013/01/17(木) 19:30:55.74, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/
. | ＼ ∠イ ,ｲイ| ,`-' |
| l＾,人| ` `-' ゝ | このスレは馬と鹿と豚さんばかりね。
| ` -'＼ ｰ' 人
| /(l ＿＿／ ヽ、
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0744１３２人目の素数さん
2013/01/17(木) 22:40:29.47ありがとうございます。
ddplyで無事いけました。
20万行程度のデータフレームですが、数秒で完了するので問題ありませんでした。
0745１３２人目の素数さん
2013/01/18(金) 01:06:24.60すみません。勘違いしてました。
オプションでいけました。
ご指摘ありがとうございます。
0746１３２人目の素数さん
2013/01/26(土) 23:13:16.48> install.packages("RMySQL", type="source")
Installing package(s) into ‘D:/program_files/R-2.15.2/library’
(as ‘lib’ is unspecified)
URL 'http://cran.rstudio.com/src/contrib/RMySQL_0.9-3.tar.gz' を試しています
Content type 'application/x-gzip' length 165363 bytes (161 Kb)
開かれた URL
downloaded 161 Kb
* installing *source* package 'RMySQL' ...
** パッケージ 'RMySQL' の解凍およびMD5サムの検証に成功しました
checking for $MYSQL_HOME... D:program_filesMySQL
MYSQL_HOME is set but does not point to a directory
ERROR: configuration failed for package 'RMySQL'
* removing 'D:/program_files/R-2.15.2/library/RMySQL'
Warning in install.packages :
命令 'D:/PROGRA~1/R-215~1.2/bin/x64/R CMD INSTALL -l "D:/program_files/R-2.15.2/library" C:\Users\<<User>>\AppData\Local\Temp\RtmpmAoBu3/downloaded_packages/RMySQL_0.9-3.tar.gz' の実行は状態 1 を持ちました
Warning in install.packages :
installation of package ‘RMySQL’ had non-zero exit status
The downloaded source packages are in
‘C:\Users\<<User>>\AppData\Local\Temp\RtmpmAoBu3\downloaded_packages’
が出て、MYSQL_HOMEの値がダメなんだと思いますが、mysql自体にはパスが通るので合ってるはずです。どこが悪いんでしょうか？
0747１３２人目の素数さん
2013/01/27(日) 08:48:14.43これってソースしかないの？
0748746
2013/01/27(日) 23:38:33.32やっぱりMYSQL_HOMEの値がおかしかったようです。
0749１３２人目の素数さん
2013/02/01(金) 10:37:40.43Rjpwikiの時間の数値変換を、使っているデータの時間フォーマットにあわせてアレンジしてるんですが
date2num <- function(x, ...) {
.day <- as.numeric(substr(x,10,11))
.hour <- as.numeric(substr(x,13,14))
.minute <- as.numeric(substr(x,16,17))
.second <- as.numeric(substr(x,19,20))
.year <- as.numeric(substr(x,2,5))
.mon <- as.numeric(substr(x,7,8))
a <- .day - 1
b <- .year - 1970
c <- 1
if(.mon==2) c <- 31
else if(.mon==3) c <- 59
else if(.mon==4) c <- 90
else if(.mon==5) c <- 120
else if(.mon==6) c <- 151
else if(.mon==7) c <- 181
else if(.mon==8) c <- 212
else if(.mon==9) c <- 243
else if(.mon==10) c <- 273
else if(.mon==11) c <- 304
else if(.mon==12) c <- 334
else if(.mon==1) c <-0
else m <-1
0750１３２人目の素数さん
2013/02/01(金) 10:40:10.94cat(file) # 画面に出力
cat("\n") # 改行コード
file<-sprintf("%03d",a)
cat(file) # 画面に出力
cat("\n") # 改行コード
file<-sprintf("%03d",c)
cat(file) # 画面に出力
cat("\n") # 改行コード
c<-c+a
file<-sprintf("%03d",c)
cat(file) # 画面に出力
cat("\n") # 改行コード
return(b*31556926+c*86400+a*86400+.hour*3600+.minute*60+.second)
}
これに" 2012/11/30 16:50:00"のような時刻データを
date,sum
2012/12/05 20:30:00,17012
2012/12/05 20:30:00,16512
2012/12/05 19:53:00,16512
2012/12/05 19:53:00,16012
以上のようなデータフレームで入力しているのですが
月を調べているif文が1度しか動きません
if文が問題なのだと思っているのですがどう直せばよいものでしょうか
0751１３２人目の素数さん
2013/02/01(金) 11:19:52.70よく分からない。date2numの引数はxだけなのに、なんで、dateとsumのデータフレームを入れるの？
それに、月を7文字目と8文字目ってしているけど、6文字目と7文字目じゃないの？
> date2num("2012/12/05 20:30:00")
以下にエラー if (.mon == 2) c <- 31 else if (.mon == 3) c <- 59 else if (.mon == (from #11) :
TRUE/FALSE が必要なところが欠損値です
[snip]
再現できるように質問しようよ。
もっと言えば、
> as.numeric(as.POSIXlt("2012/12/05 20:30:00"))
[1] 1354707000
とかを使うのは嫌いなの？
0752１３２人目の素数さん
2013/02/01(金) 12:21:01.54自分の力が及ばないせいでデータファイルを作るとき時刻データの前にどうしても消えない半角スペースがあって
それで6番目と7番目ではなく7番目8番目になっていますorz
それと
> as.numeric(as.POSIXlt("2012/12/05 20:30:00"))
[1] 1354707000
これをみて目からうろこでした、自分のやってたことの馬鹿さ加減にあきれました
御助言ありがとうございました
0753１３２人目の素数さん
2013/02/01(金) 13:26:36.91>時刻データの前にどうしても消えない半角スペース
関数を自作するレベルなのになぜw
例えば、
> (a <- c(" 2013-02-01 13:21:17", "2013-01-01 09:21:17", " 2013-03-01 11:21:17"))
[1] " 2013-02-01 13:21:17" "2013-01-01 09:21:17" " 2013-03-01 11:21:17"
こんな感じで、時刻前に空白があったりなかったりするなら、
> library(stringr)
> str_trim(a)
[1] "2013-02-01 13:21:17" "2013-01-01 09:21:17" "2013-03-01 11:21:17"
これでよいのでは？
オプションなしでは、前後の空白を取り除くから、前だけならそのように指定する。
0754１３２人目の素数さん
2013/02/01(金) 13:49:46.73なにからなにまでありがとうございます
0755１３２人目の素数さん
2013/02/01(金) 14:11:14.62>Rの方でtrimするって手もありましたね
別にRの外でも、どこででも空白を削除すれば良いと思う。
例えば下記のようなcsvがあったら、
$ cat tmp.csv
"A","B","C"
"BE DAC","Jan",0.35
"AC BDE","Feb",0.34
" EBACD","Mar",0.66
" DABEC","Apr",0.44
"BEA CD","May",0.84
"CDA EB","Jun",0.3
"EDCAB ","Jul",0.59
"CBEA D","Aug",0.53
"ECDA B","Sep",0.44
"ABDCE ","Oct",0.33
$ sed -e 's/¥" ¥([^ ,]*¥)/¥"¥1/g' tmp.csv
"A","B","C"
"BE DAC","Jan",0.35
"AC BDE","Feb",0.34
"EBACD","Mar",0.66
"DABEC","Apr",0.44
"BEA CD","May",0.84
"CDA EB","Jun",0.3
"EDCAB ","Jul",0.59
"CBEA D","Aug",0.53
"ECDA B","Sep",0.44
"ABDCE ","Oct",0.33
とすると3,4行目にあった文字列前の空白を消せる。
0756１３２人目の素数さん
2013/02/02(土) 10:13:07.27http://cse.naro.affrc.go.jp/takezawa/r-tips/r/43.html
上記サイトの複数の列に対してソートを行う場合は関数を使用して
同じように
a b c
11 1 2 11
12 1 1 12
9 2 4 9
10 2 3 10
7 3 6 7
8 3 5 8
5 5 8 5
6 5 7 6
4 7 9 4
1 9 12 1
2 9 11 2
3 9 10 3
というデータフレームをaを昇順に整列しaで同じものがあった場合はcの小さいものを上にするようにしたいのですが
In Ops.factor(mmm, each) : < 因子に対しては無意味です というエラーが出てうまくいきません
どこが問題でしょうか？
0757１３２人目の素数さん
2013/02/02(土) 10:18:28.45a b c
1 2 11
1 1 12
2 4 9
2 3 10
3 6 7
3 5 8
5 8 5
5 7 6
7 9 4
9 12 1
9 11 2
9 10 3
です
実際にはデータフレーム名test1で
sortlist <- order(test1$a, pmax(test1$a,test1$c))
test1[sortlist,]
以上に実行しました
0758756,757
2013/02/02(土) 10:45:26.080759１３２人目の素数さん
2013/02/02(土) 11:01:54.12複数カラムをキーにしたソートなんて考えたことがなかった。
いったいどんなときに必要になるの？
# 表計算ソフトのように「見せる」目的なら理解できるけど。
0760１３２人目の素数さん
2013/02/02(土) 11:09:52.570761１３２人目の素数さん
2013/02/02(土) 12:03:35.44a,b
3,5
2,3
2,2
2,1
1,3
1,1
a列は必ず降順になっているときに
b列が降順になっていないとき
この例では5行目1,3ですが
これを削除したいのですが 関数化して自動で判別して削除したいです
どのようにしたら良いでしょうか
0762１３２人目の素数さん
2013/02/02(土) 12:16:56.00doByパッケージのorderBy()が使いやすい。
0763１３２人目の素数さん
2013/02/02(土) 12:36:52.26pmax()というのを初めて知った。
bの列を一行下げて頭にInfをいれたものをcとしてデータフレームにいれる、bとcに対してpmax()して、結果をdとしてデータフレームにいれる。dがbと一致しているところを落とす。
でどうだろう？
0764761
2013/02/02(土) 13:18:11.28御助言ありがとうございます
しかし当方最近Rを使いだしたもので何をどうすればいいのか分かりません・・・
あつかましいのですが具体的にどんな感じに書けばいいでしょうか
0765１３２人目の素数さん
2013/02/02(土) 14:23:56.55dat = read.csv(text = "a,b
3,5
2,3
2,2
2,1
1,3
1,1")
dat
dat$c = c(Inf , head(dat$b, -1))
dat
with(dat, pmax(b, c))
dat$d = with(dat, pmax(b, c)) == dat$c
dat
dat[dat$d,]
を想定。関数化されてないけど。もっと大きい例で動くはチェックしていない。
0766１３２人目の素数さん
2013/02/02(土) 15:20:57.47> dat = read.csv(text = "a,b
+ 3,5
+ 2,3
+ 2,2
+ 2,1
+ 1,3
+ 1,1")
> dat[c(TRUE, diff(dat$b) < 0), ]
a b
1 3 5
2 2 3
3 2 2
4 2 1
6 1 1
0767１３２人目の素数さん
2013/02/02(土) 15:35:34.61bが例えば、10, 8, 6, 9, 8, 5とかになっていると、
>>766の方法を1回だけ適用したのでは通用しない。
> d2
a b
1 3 10
2 2 8
3 2 6
4 2 9
5 1 8
6 1 5
> d2[c(TRUE, diff(d2$b) < 0), ]
a b
1 3 10
2 2 8
3 2 6
5 1 8
6 1 5
そのため、十分な回数を繰り返し適用することになる。
> sum(diff(d2$b) > 0)
[1] 1
これが0になるのが、終了条件。
0768１３２人目の素数さん
2013/02/02(土) 19:57:32.38を使えばいいか？
0769１３２人目の素数さん
2013/02/03(日) 12:03:27.26dat = read.csv(text = "a,b
3,5
2,3
2,2
2,1
1,9
1,8
1,1")
dat
dat$cummin <- cummin(dat$b)
dat$correct <- with(dat, b == cummin)
dat
dat[dat$correct == TRUE,]
0770１３２人目の素数さん
2013/02/03(日) 13:34:45.070771１３２人目の素数さん
2013/02/03(日) 22:46:23.58膨大なデータを統計処理して活用するために
RやJavaのhadoopが使える人材が求められてるって書いてあったな
0772１３２人目の素数さん
2013/02/05(火) 08:32:45.91Baseパッケージで何か良い方法ありませんでしょうか？
今は↓のように書いていますが、ablineで線を引くとはみ出してしまう等手間がかかるので、
簡単にかける方法をご存知の方おりましたらお願いします
plot(x=0,y=0,xlim=c(0,6),ylim=c(0,3),lwd=2,ann=F,cex.axis=1,las=1,type="n",axes=F)
axis(side=1,seq(0,6),lab=F,pos=0,tck=0); axis(side=2,seq(0,3),lab=F,tck=0,pos=0)
axis(side=3,at=0:6,lab=F,tck=0,pos=3); axis(side=4,at=0:3,tck=0,lab=F,pos=6)
axis(side=1,at=0,tck=0,pos=0)
>>735-37
レス忘れていました、申し訳ありません
Rconsoleが2つ存在し妙な読み込みがされていたのが原因でした
アドバイスありがとうございました
0773１３２人目の素数さん
2013/02/05(火) 09:04:30.55> plot.new()
> box()
とすると、左下に原点が来るけど。
> points(0, 0)
0774１３２人目の素数さん
2013/02/05(火) 09:17:22.01R的なアプローチとして、複雑な手続きを関数化して簡単にする方法がある。
f772 <- function(){
plot(x=0,y=0,xlim=c(0,6),ylim=c(0,3),lwd=2,ann=F,cex.axis=1,las=1,type="n",axes=F)
axis(side=1,seq(0,6),lab=F,pos=0,tck=0); axis(side=2,seq(0,3),lab=F,tck=0,pos=0)
axis(side=3,at=0:6,lab=F,tck=0,pos=3); axis(side=4,at=0:3,tck=0,lab=F,pos=6)
axis(side=1,at=0,tck=0,pos=0)
}
とスクリプトの最初に書いておけば、
> f772()
以降は、これだけでよい。Baseパッケージのコマンドを探すよりも簡単。
0775１３２人目の素数さん
2013/02/06(水) 03:22:55.78plot(1,1, xaxs = "i", yaxs = "i", xlim = c(0,1), ylim = c(0,1), tck = 0)
?parでbaseの細かいオプションの説明がみられる。
これも役立つかも
www.stat.auckland.ac.nz/~paul/RG2e/chapter3.html
0776１３２人目の素数さん
2013/02/06(水) 06:43:29.05, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/ 私は只の数ヲタなんかとは付き合わないわ。
. | ＼ ∠イ ,ｲイ| ,`-' | 頭が良くて数学が出来てかっこいい人。それが必要条件よ。
| l＾,人| ` `-' ゝ | さらに Ann.of Math に論文書けば十分条件にもなるわよ。
| ` -'＼ ｰ' 人 一番嫌いなのは論文数を増やすためにくだらない論文を書いて
| /(l ＿＿／ ヽ、 良い論文の出版を遅らせるお馬鹿な人。
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、 あなたの論文が Ann of Math に accept される確率は?
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼ それとも最近は Inv. Math. の方が上かしら？
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0777１３２人目の素数さん
2013/02/06(水) 17:31:58.49レスありがとうございます
非常に参考になりました
0778１３２人目の素数さん
2013/02/15(金) 09:19:22.40軸(目盛の位置)がずれてしまうのですが、
揃える方法ってあるのでしょうか？
barplot(1:10, xlim=c(0,10), ylim=c(0,10))
par(new=T)
plot(1:10, col=2, xlim=c(0,10), ylim=c(0,10))
0779１３２人目の素数さん
2013/02/15(金) 10:52:07.52よく分からないが、
> a <- barplot(1:10); points(a, 1:10, col = 2)
こうしたいの？
0780１３２人目の素数さん
2013/02/15(金) 11:18:55.04縦軸の0の位置が、barplot と plot で違うんで揃えたいんです。
0781１３２人目の素数さん
2013/02/15(金) 11:34:08.72いや、そもそも揃っているけど。
> a <- barplot(1:10); points(a, rep(0, 10), col = 2)
0782781
2013/02/15(金) 11:39:38.94あっ、分かった。plot()を無理矢理使っているんだ。
addオプションをサポートしているクラス以外では、plot()を重ねるのは禁止。
想定外の使用形態。
描画を重ねるとき、普通はpoints()やlines()やaxis()を使う。
0783１３２人目の素数さん
2013/02/15(金) 11:44:27.170784781
2013/02/15(金) 11:53:12.37例えば、棒グラフと折れ線グラフを重ねたかったら、次のようにする。
x <- 1:10
n <- length(x)
y <- jitter(rep(10, n), 10)
op <- par(mai = c(1.02, 0.82, 0.82, 0.82))
a <- barplot(x, ylim = c(0, 13))
points(a, y, type = "l", lty = 2)
axis(4, label = (0:6) * 2 / 100, at = (0:6) * 2, las = 2)
par(op)
0785１３２人目の素数さん
2013/03/09(土) 12:17:39.65plyr : 中規模以上のデータで遅すぎ（idataframe使ってもそこまで変わらず）
aggregate : 普通
sqldf : まとも、でもRっぽくない。
他にお勧めがあったら教えて
0786１３２人目の素数さん
2013/03/13(水) 01:03:30.88rawデータをcharacterにしたいです。
chromeの localstorageのデータを整理しようとしています。
ドメインごとのsqlite形式？ファイルがあり、ItemTableというテーブル名で、
key,valueの項目でデータが入っています。
valueのデータ形式がblob型です。
RSQLiteで、
ret <- dbGetGuery(con, "select value ItemTable")
でデータフレームでデータが返って来ました。
class(ret[1,1][[1]])
>> raw
それで、これから先、
ret[1,1][[1]] を characterにしたのですが、
どうすればいいのでしょうか？
0787１３２人目の素数さん
2013/03/13(水) 12:10:10.92内部でSQLiteを使っているソフトが最近多いなと思っていたが、
Chromeもそうなんだ。
再現して検証できない質問なので、適当な回答になるけど、
as.character(ret[1,1][[1]])でどのようなエラーになりますか？
0788786
2013/03/13(水) 12:57:52.15ありがとうございます。
chrome側のconsoleで、 localstorage.hoge = "abc" として、
sqlteデータを入れた時に、
R側でそのblob型のデータ(abcというデータ, 他のsqlite viewerで確認済み) を読むと、
> ret[3,2][[1]]
[1] 32 00 30 00 31 00 33 00 2d 00 30 00 33 00 2d 00 31 00 32 00 54 00 31 00 34
[26] 00 3a 00 33 00 35 00 3a 00 35 00 30 00 2e 00 31 00 34 00 33 00 5a 00
となり、
as.character(ret[3,2][[1]])
[1] "32" "00" "30" "00" "31" "00" "33" "00" "2d" "00" "30" "00" "33" "00" "2d"
[16] "00" "31" "00" "32" "00" "54" "00" "31" "00" "34" "00" "3a" "00" "33" "00"
[31] "35" "00" "3a" "00" "35" "00" "30" "00" "2e" "00" "31" "00" "34" "00" "33"
[46] "00" "5a" "00"
となります。
こういうのを16進数のbyteデータというのですよね？
これをcharacterにencode?してやりたいです。
0789１３２人目の素数さん
2013/03/13(水) 13:54:01.89元データのエンコードが分からないと厳しいと思います。日本語でしょうか。
英語(ASCII)なら、
rawToChar(ret[3,2][[1]])とするとうまく行くかも知れません。
0790789
2013/03/13(水) 15:25:46.98ちゃんと読んでいませんでした。申し訳ない。
これは「abc」という文字列なんですね。
ちなみに"abc"の8進数表記は次の通り
$ echo abc | od -to1c
0000000 141 142 143 012
a b c \n
10進数表記
$ echo abc | od -td1c
0000000 97 98 99 10
a b c \n
16進数表記
$ echo abc | od -tx1c
0000000 61 62 63 0a
a b c \n
>>788 を見る限り、abcじゃなさそうに思える。
16進数と仮定すると、
2 NUL 0 NUL 1 NUL 3 NUL - NUL 3 NUL -
NUL 1 NUL 2 NUL T NUL 1 NUL 4 NUL : NUL 3 NUL
5 NUL : NUL 5 NUL 3 NUL . NUL 1 NUL 4 NUL 3
NUL Z NUL
となるが、NULを飛ばすと、2013-3-12T14:35:53.143Zで、日付時間っぽいな。
これ以上は分からないので、識者の方にパス。
0791１３２人目の素数さん
2013/03/13(水) 15:28:55.56ありがとうございます。
かなり前進しました。
さきほど、データの提示は、間違え(一つずれたものでした)
"abc"というデータを入れたときは、
ret[4,2][[1]]
[1] 61 00 62 00 63 00
rawToChar(ret[4,2][[1]],multiple=T)
[1] "a" "" "b" "" "c" ""
なぜか、空文字が入りますがOKです。
また、日本語の場合ですが、
私のwindows環境の文字コードは、cp932だと思います。
>Sys.getlocale()
[1] "LC_COLLATE=Japanese_Japan.932;LC_CTYPE=Japanese_Japan.932;......."
それで、データとして "あい" と localstorageに入れて、
RSQLiteでとりだしたところ、
ret[5,2][[1]]
[1] 42 30 44 30
となり、
これは、文字コード表を見ると、
unicode (あい) => UTF16形式 3042 3044 と、
２バイトの間隔で前後入れ替えしたものに見えます。
mapply(function(a,b)paste(b,a,sep=""), x[seq(1,3,2)], x[seq(2,4,2)])
[1] "3042" "3044"
これを、"あい" という文字で、
できれば、UTF8形式でファイルにおとしたいのですが、、
0792１３２人目の素数さん
2013/03/13(水) 15:31:21.55返事が入れ違いになりました。
"abc"ではなく、日付でした。
間違えてました。
0793１３２人目の素数さん
2013/03/13(水) 21:35:55.82下記のようなCSVファイル(,区切り)が存在します。ID、NAME、MONTHは各情報の最初の行にしかデータとして入っておりません。
---------- データ部 ----------
ID, NAME, MONTH, PRICE, NUMBER
AAA1, Apple, 1, 100, 10
, , , 98, 12
, , , 102, 8
AAA2, Orange, 2, 120, 3
, , , 130, 4
, , , 112, 8
------------------------------
これをデータフレームとして読み込みます。
ID NAME MONTH PRICE NUMBER
AAA1 Apple 1 100 10
"" "" "" 98 12
"" "" "" 102 8
AAA2 Orange 2 120 3
"" "" "" 130 4
"" "" "" 112 8
.......
となります。
ここでID NAME MONTHを上の行から補間したいのですが、どのようにすればよいでしょうか?
ID NAME MONTH PRICE NUMBER
AAA1 Apple 1 100 10
AAA1 Apple 1 98 12
AAA1 Apple 1 102 8
AAA2 Orange 2 120 3
AAA2 Orange 2 130 4
........
としたいです。補間さえできれば形式はデータフレームでなくても結構です。よろしくお願いいたします。
0794786
2013/03/13(水) 22:48:40.11http://stackoverflow.com/questions/2776135/last-observation-carried-forward-in-a-data-frame
を参考にするといいよ。
""をNAに変えておかないといけないので、
df[ df == ""] <- NA
としておいて、
library(zoo)
na.locf(df)
まだ、rawデータからcharacterに変えられてない、、、
0795１３２人目の素数さん
2013/03/14(木) 09:59:33.96単純にコピーするだけでは駄目なの？
> sample <- data.frame(matrix(c("AAA1","Apple",1,100,10,
+ "","","",98,12,
+ "","","",102,8,
+ "AAA2","Orange",2,120,3,
+ "","","",130,4,
+ "","","",112,8), ncol = 5, byrow = TRUE))
> names(sample) <- c("ID", "NAME", "MONTH", "PRICE", "NUMBER")
> sample
ID NAME MONTH PRICE NUMBER
1 AAA1 Apple 1 100 10
2 98 12
3 102 8
4 AAA2 Orange 2 120 3
5 130 4
6 112 8
> for(i in 1:(nrow(sample) - 1)){
+ if(sample[i+1, 1] == ""){
+ sample[i+1, 1:3] <- sample[i, 1:3]
+ }
+ }
> sample
ID NAME MONTH PRICE NUMBER
1 AAA1 Apple 1 100 10
2 AAA1 Apple 1 98 12
3 AAA1 Apple 1 102 8
4 AAA2 Orange 2 120 3
5 AAA2 Orange 2 130 4
6 AAA2 Orange 2 112 8
0796793
2013/03/14(木) 22:00:08.24ありがとうございます。
無事うまくいけました。
質問に答えることができればよかったのですが、
残念ながら私にはわかりません。申し訳ありません。
>>795
実際のデータは100万行ほどありますので、
ループを回すと試してはいませんが時間がかかるかと。。。。
■ このスレッドは過去ログ倉庫に格納されています