Ｑ＆Ａ
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
COLOR(blue){SIZE(20){R および RjpWiki に関する質問コーナ...
COLOR(red){SIZE(20){このコーナーはスパム対策のため閉鎖し...
//COLOR(red){注意}：[[Ｑ＆Ａ]] コーナーが少し重くなり過ぎ...
COLOR(red){注意}：新規記事用の入力欄は以下の目次の直後に...
もし新規に質問をされる場合は、効率的にコメントを得るため...
- 広い意味で R に関係する話題に限定!
- 学校の宿題は自分で考えるべきです
- 具体的な状況がわかる様に背景説明をけちらない
- [[CRAN の検索エンジン:http://finzi.psych.upenn.edu/sear...
- 自分には馴染みでも、他の人には馴染みの無い手法・概念も...
- もし問題が解決した(結局解決しない場合を含め)場合はその...
- 後で他の人が参照する際の便宜のため、タイトルには簡潔で...
- 回答者は神通力を持ってはいません。「これこれしたが、う...
~
----
アーカイブ：[[Ｑ＆Ａ(旧1)]]~
COLOR(red){SIZE(20){このコーナーはスパム対策のため閉鎖し...
#contents
----
COLOR(red){注意}：投稿失敗した記事をページ先頭の「新規｜...
//#article
**文字式による数式の関数認識
>[[KOSHI]] (2006-01-18 (水) 22:58:51)~
~
expression()の形式で入力した数式を，関数として認識したい...
例えば
> f <- expression(x*y)
を（x,y）の関数に直したいとき，
> F <- function(x,y){}
> body(F) <- f
とやれば
> F <- (x, y) x * y
とできます．しかし，これはfで使われる変数が（x,y）と既知...
もし，fでの変数が未知の場合，これを関数認識するにはどうし...
> g <- G(expression(x*y*z*u))
などとすると,ちゃんとこの4つの変数（x,y,z,u）を使って~
> g <- function(x,y,z,u) x*y*z*u
となるような関数を作ることは可能でしょうか．~
~
ちなみに~
deriv(expr, name, func=T)
という関数を使うと，勝手に書いた数式~
expr=expression(...)
の，変数nameによる1階偏微分を，その文字による関数として認...
~
どなたかお力をお貸しいただけないでしょうか．~
//
-ソースを見ることについて，どこかのページにも書いてありま...
まず，関数名だけをタイプしてみる
> deriv
function (expr, ...)
UseMethod("deriv")
<environment: namespace:stats>
なんだぁ？これだけか？と思ったら，
> methods(deriv)
[1] deriv.default deriv.formula
と，やってみる。ここに示されるのが，実際に呼び出される候...
単に関数名をコンソールにタイプするのではなくて，以前に出...
> stats:::deriv.default
function (expr, namevec, function.arg = NULL, tag = ".ex...
hessian = FALSE, ...)
.Internal(deriv.default(expr, namevec, function.arg, tag...
<environment: namespace:stats>
もう一つ
> stats:::deriv.formula
function (expr, namevec, function.arg = NULL, tag = ".ex...
hessian = FALSE, ...)
{
if ((le <- length(expr)) > 1)
.Internal(deriv.default(expr[[le]], namevec, fun...
tag, hessian))
else stop("invalid formula in deriv")
}
<environment: namespace:stats>
これらの関数のどれがどのように使われるかは，help で検索の...
ここで示されるのは，R の範囲で書かれるプログラムであって...
そのような部分についてのプログラムの詳細は，（えい面倒だ...
てなことかな（はずしたか。。） -- &new{2006-01-18 (水) 2...
-それとも，あなたは，このようなことがお望みなのかな？ ...
> func <- as.function(a <- function(x, y) x+y)
> func
function(x, y) x+y
> func(2,4)
[1] 6
> func <- as.function(a <- function(x, y) x*y)
> func(3,7)
[1] 21
一応そういうことで。 -- &new{2006-01-18 (水) 23:36:17};
-私の場合，ソースを見ても分からない部類だと思います・・・...
-実際，deriv（）はどんな文字を用いて数式を書いても，それ...
-具体的には，どんなものですか？上の，func だって，x,y を...
-＞実際，deriv（）はどんな文字を用いて数式を書いても，そ...
derive を読んでください(^_^;)~
もしくは，ウイザード召喚！！(^_^;) -- &new{2006-01-18 (...
-たとえば，ユーザが f <- expression(x*y) と入力すること...
> G <- function(expr){
F <- function(x,y){ }
body(F) <- f
return(F)
}
というGを用意しておけば(x,yの数式なら何でも)関数に出来ま...
> f <- expression(u*v)
と書いてしまえば，G(f)は(u,v)の関数に出来ないじゃないです...
-関数定義に出てくる名前をパースして，それぞれが想定してい...
-deriv も何が変数（とみなすべき）かを指示する仕様になって...
> a <- 3; b <- 5
> f <- alist(x=,y=a,x*y+b)
> F <- as.function(f)
> F
function (x, y = a)
x * y + b
> F(1)
[1] 8
-＞ alist を文字列から合成する方法.....一時ファイルに書き...
-なるほど。一時ファイル名をその場で作る関数が確かあったは...
-tempfile 関数ですね。。 -- &new{2006-01-20 (金) 10:56:4...
-皆様、いろいろとご助言有難うございます。自分でも少し進展...
> Func <- function(expr){
df <- deriv(expr, all.vars(expr), func=T) #とり...
F <- function(p,q){ #ここを変数２個にしておく。
df(p,q)[1] #derivの[1]には微分しない関数が返...
}
return(F)
}
こうして関数Fを作れば、ユーザが入力する「変数の個数」さえ...
あとは、自然数nを与えられた時に、自動的にn変数の空の関数...
> body(G(n)) <- df(x1,...xn)[1]
で完成するのではないでしょうか？~
しかし、こういうG(n)ってできますでしょうか？~
forとか使っても関数の変数を増やせなくて・・・あと一歩だと...
よろしくお願いします。 -- [[KOSHI]] &new{2006-01-20 (金) ...
-上に示唆されていた方法を使えば（ファイルはテキスト、文字...
> G <- function(n) {
tmpfile <- tempfile() # 一...
.x <- paste("x", 1:n, ",", sep="") # 引...
.x[n] <- substr(.x[n], 1, nchar(.x[n])-1) # 最...
sink(tmpfile)
cat("function (", .x, ") {}?n") # 関...
sink() # フ...
.f <- source(tmpfile)[[1]] # フ...
unlink(tmpfile) # 一...
return(.f)
}
> f5 <- G(5)
> f5
function ( x1, x2, x3, x4, x5 ) {}
> str(f5)
function (x1, x2, x3, x4, x5)
- attr(*, "source")= chr "function ( x1, x2, x3, x4, x...
- G<-function(n)eval(parse(text=paste("function(",paste("...
-↑ 上手い！ しかし、ファイルに一旦書き出す方法は（少し野...
-どうも，このページに粘着質のスパム投稿者が住み着いたよう...
-Q & A (中級) とでも改名して再出発しましょうか。みたとこ...
-スパム対策のためＱ＆Ａコーナーは閉鎖します。以後は [[Ｑ...
**共分散分析
>[[S.O]] (2006-01-13 (金) 03:10:39)~
~
共分散分析のやり方が、いまいち分かりません。~
yが目的変数で、x1,x2は連続量で、x3,x4は離散です。~
共分散分析は、lmよりaovの方がよいと言われていたので、以下...
aov(y~x1+x2+x3+x4)~
この場合、普通の４元配置分散分析だと思い、x1,x2を共変量に...
お手数ですが、どうぞ宜しくお願いいたします。~
//
-「いまいちわからない」ではこっちも「いまいちわからない」...
**Mac版のR2WinBUGS
>[[まさひと]] (2006-01-11 (水) 20:54:49)~
~
パッケージではMac用のR2WinBUGSが配付されいますが、もとも...
//
-example(bugs) で，なんか動いているけど。よくわかりません...
-wine経由で動かす(system関数使って)物です. 私の手元でも(D...
-Darwineをインストールしてみました. GUIのバージョンの関係...
-なかまさん、ありがとうございます。さっそく試してみます。...
-あー、すいません。よく見るとまだまだ普通には動かんようで...
**Rで＠Risk的なパッケージを自作
>[[Tak]] (2006-01-10 (火) 18:43:58)~
~
リスク解析の勉強をかねて、Rで＠Riskの様なパッケージを自作...
http://www.palisade-europe.com/~
~
既存のRのパッケージをできるだけ使って、＠Riskの様な包括的...
~
他にparametric sensitivity analysis、モンテカルロ法、deci...
//
-Palisade社が何らかの特許を申請していた場合、完全移植には...
-確かにPalisade社の名前をあげる必要はなかったですね。＠Ri...
**plotのwidthの種類
>[[浩太]] (2006-01-05 (木) 13:52:55)~
~
barplotの棒の幅を定めるwidthって何種類ぐらいの幅がありま...
具体的にどんな種類の幅がありますか？~
教えてください~
//
-質問の意味がわかりません。widthは数値ベクトル（リサイク...
> x <- c(1,4,8,5,2)
> barplot(x, width=c(1.2, 3.4, 5.321, 2.345))
これにより以下のような図が得られるんだが。。。
#ref(barplot.png)
それとも，
> barplot(x, width=1)
> barplot(x, width=1000)
のような描画では差がないように見えるからということかな？~
ヘルプを読めばわかるように，xlim で横軸の長さを決めない限...
Specifying a single value will no visible effect unless x...
-ｐarを使って図を重ね合わせる時に差をつけるときも数値でや...
-浩太さんがparで図を重ねることができる人ならば、実際に試...
-四つのbarplotをparを使って重ねあわせました。それぞれの違...
-用いたプログラムをここに載せればいかがですか。長すぎたり...
なお，投稿の失敗は，上のほうの「編集」をクリックすること...
-一般的に言えば，幅を変えると必然的に面積も変わるわけで，...
**barplotでNameが表示されない
>[[鍋]] (2006-01-03 (火) 14:29:13)~
~
> mydata
[1] 29 24 20 18 25 22 23 20 20 23 20 30 21 20 19 26 14 2...
[26] 2 2 4 1 3
> main
[1] "Turnover Number"
> names.arg
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 1...
[26] 26 27 28 29 30 31 32
> xlab
[1] "Event Time Number"
> ylab
[1] "Turnover"
> ylim
[1] 0 32
> barplot(mydata, main=main, names.arg=names.arg, xlab=x...
Error in barplot.default(mydata, main = main, names.arg ...
incorrect number of names
なぜこのようなErrorが表示されるのかわかりません。~
incorrectになるのはどういったときなのでしょうか？~
//
-名前の数が合わないと申しておりますし,実際合ってないよう...
-length(mydata) の結果と，length(names.arg) の結果を比べ...
-国際化された R を使えば，エラーメッセージも日本語で出る...
以下にエラーbarplot.default(mydata, main = main, names.a...
名前の数が誤りです
-- &new{2006-01-03 (火) 15:41:29};
-鍋さんの悪い癖は，コメントに反応しないことだね。 -- &ne...
**R. Tierney 氏の R バイトコンパイラーについて
>[[ma]] (2005-12-29 (木) 18:12:39)~
~
何でも掲示板に出ている「R. Tierney 氏の R バイトコンパイ...
~
”bad function”というエラーが出た場合には具体的にどのよう...
漠然とした質問で申し訳ないのですが、どなたかお教え頂けな...
//
-少々説明を追加させて下さい。cmpfun()でバイトコンパイルは...
-もしあまり長くなければ、エラーが起きる実例を紹介するとコ...
**軸の省略
>[[輝]] (2005-12-23 (金) 16:33:08)~
~
barplot()で，図を描いており，一つだけ大きい項目があり，~
途中のY軸(目盛り）を省きたいのですが，~
どのようにすればよいのですか？~
//
-[[Ｑ＆Ａ (初級者コース)#content_1_29]]にある，「ヒストグ...
-先の質問をしたAkiraです。外れ値がある場合は対数変換とい...
**RExcelの起動時エラー
>[[YK]] (2005-12-23 (金) 12:52:29)~
~
RExcelを利用しようとしたのですが，Excel起動時にエラーが発...
~
エラーメッセージは，以下のとおりです．
「実行時エラー 5 プロシージャの呼び出し、または引数が不...
そのエラーの発生箇所は，☆の付いている箇所です．
Private Sub RemoveRMenu()
On Error Resume Next
☆CommandBars(1).Controls("RExcel").Delete
End Sub
RExcelの実行した環境は以下のとおりです．
Windows XP(SP2)
Office Excel 2003
R 2.2.1
RServ 2.0.0（Rserv1.3.5でも同様にエラーが発生しました）
RExcelに添付されていたExcel 02 - Macro Demoを実行した際に...
~
エラー修正，回避などのヒントをご教授願えないでしょうか。~
よろしくお願いいたします。~
//
-思いつきです。3年ほど前の話ですが、RExcelはエラーが多発...
-ありがとうございます．バージョンによるエラーを視野に入れ...
**MASS lda 判別関数と偏F値について
>[[mt]] (2005-12-21 (水) 11:44:53)~
~
MASS ライブラリ中の lda() 関数を実行した結果から、判別関...
R中ではpredictで判別してくれるのでいいのですが、実際どの...
どなたかヒントなどご教授願えないでしょうか。~
よろしくお願いいたします。~
//
-lda は，正準判別分析を行っていますね。ということで，偏F...
**R2.2.0でのeps出力について。
>[[KJ]] (2005-12-20 (火) 21:37:18)~
~
Windows 版の R-2.2.0 で eps 形式での出力を行ったファイル...
それ以前のバージョンの R であれば問題なく表示できるのです...
原因は何かわかりますでしょうか？
GSview 4.7 2005-03-26
Unknown in Comments section at line 2:
%%DocumentNeededResources: font Times-Roman
Unknown in Comments section at line 3:
%%+ font Times-Bold
Unknown in Comments section at line 4:
%%+ font Times-Italic
中略
End offending input ---
file offset = 5390
gsapi_run_string_continue returns -101
~
//
-うーん, gsの環境のような気がするんですが...もし差し支え...
-ありがとうございます。
x <- seq(-3,3,by=0.01)
plot(x,dnorm(x),type="l")
として、右クリックでeps保存したファイルです。~ -- [[KJ]]
#ref(http://www.sfc.keio.ac.jp/~kj198276/test.eps)
-中は闇Ｒ-2.2.0のようですが,多分GSがCID未対応の設定もしく...
setHook(packageEvent("grDevices", "onLoad"),
function(...) grDevices::ps.options(cidfam...
等として日本語環境を無効にするか, gsが日本語を扱えるよう...
あげてみてください。 -- [[なかま]] &new{2005-12-21 (水) 2...
-あ、闇版を使わないと言う手もあります。(でも、2.3.0になる...
-知識不足で、仰っている意味が捉えきれません…ちょっと調べ...
-やはりGSの設定の問題でした。情報ありがとうございました！...
**barplot( )の凡例位置の変更
>[[兄歯]] (2005-12-17 (土) 14:17:53)~
~
barplot( )の凡例位置がグラフと重なってしまうので，凡例の...
//
-legend()はいかがでしょうか？~
x <- matrix(rnorm(100), ncol=10)
dimnames(x) <- list(1:10, 1:10)
barplot(x, beside=TRUE, legend.text=rownames(x), col=1:1...
barplot(x, beside=TRUE, col=1:10)
legend("topright", legend=rownames(x), pch=15, col=1:10)...
legend("topleft", legend=rownames(x), pch=15, col=1:10)#...
-- [[Akira]] &new{2005-12-17 (土) 15:33:28};
-解決しました。どうもありがとうございます。 -- [[兄歯]] &...
**軸の分割数について
>[[内藤]] (2005-12-16 (金) 16:57:40)~
~
plot( )で，軸（目盛り)の分割を任意に与えたいのですが，~
どのようにすればよいのですか？~
//
-? axis をやってみればよろしいかと。
hist(rnorm(1000),xaxt="n")
axis(1,pos=0,at=c(-1,0.5,2.8),labels=TRUE)
などで，遊んでみる。 -- &new{2005-12-16 (金) 18:00:36};
-どうもありがとうございます。助かりました。 -- [[内藤]] &...
**maptoolsの座標の参照方法
>[[traveller]] (2005-12-06 (火) 12:42:34)~
~
maptoolsの各図形を参照する方法、どなたか、教えていただ...
//
-各図形を参照するとは？もう少し具体的に。何をしたいのかも...
-ポイント、ラインおよびポリゴンの座標、BBOXを参照する方法...
-maptoolsで扱うObjectはMap Obejctと呼ばれ、"Shapes"と"att...
> x <- read.shape(system.file("shapes/sids.shp", package...
> x$Shapes
[[1]]
[[1]]$Pstart
[1] 0
[[1]]$verts
[,1] [,2]
[1,] -81.47276 36.23436
[2,] -81.54084 36.27251
[3,] -81.56198 36.27359
[略]
[26,] -81.45289 36.23959
[27,] -81.47276 36.23436
[[1]]$shp.type
[1] 5
[[1]]$nVerts
[1] 27
[[1]]$nParts
[1] 1
[[1]]$bbox
[1] -81.74107 36.23436 -81.23989 36.58965
[[1]]$shpID
[1] 0
attr(,"nVerts")
[1] 27
attr(,"nParts")
[1] 1
attr(,"shp.type")
[1] 5
attr(,"bbox")
[1] -81.74107 36.23436 -81.23989 36.58965
[以下略]
従って、例えば、1つ目のポリゴンの座標を得る場合には
> x$Shapes[[1]]$verts
でOKです。BBOX座標が必要な場合は、
> x$Shapes[[1]]$bbox
ですね。
-谷村さん、ご親切なご回答ありがとうございます。 -- [[trav...
**ベイズクラスタ
>[[ド級初心者]] (2005-11-07 (月) 12:58:10)~
~
Rでベイズクラスタ分析が可能と風の噂で聞きましたが、どのプ...
//
**コーフェン相関係数の求め方
>[[mori]] (2005-11-02 (水) 14:08:03)~
~
階層的クラスタリングでできた樹形図（デンドログラム）の結...
//
-知らないものだったのでググッてみた。~
Sneath and Sokal『数理分類学』（1994年，内田...
ひどい訳本みたいで，三中先生に「だいたいこの訳本，トンデ...
[[http://cse.niaes.affrc.go.jp/minaka/diary2003-09.html:h...
[[http://www.i-juse.co.jp/statistics/support/sympo/m10/m1...
[[http://www.sci.kagoshima-u.ac.jp/~dllsa/jp/dlls/cluscrt...
-ちゃんとした名前がわかれば，さらにググッてみることができ...
-ありがとうございます。その関数のようです。綴りはcoffenだ...
-"どーして「コーフェン」になるの" -> これは multicolinear...
**ロジスティック回帰での変数指定方法とresidualsなどの値
>[[たかくら]] (2005-11-01 (火) 19:00:18)~
~
glmでロジスティック回帰をしていてほとんど同じ（？）計算...
例えば、独立変数xが１，２，３のそれぞれの場合に、50試行の...
このとき、試行数や成功例数の与え方としてweightsパラメータ...
~
resp<-rep(c(1,0),3)~
x<-(1,1,2,2,3,3)~
weight<-c(5,45,20,30,45,5)~
df1<-data.frame(resp,x,weight)~
summary(glm(resp~x,weights=weight,family=binomial,data=df...
~
x<-(1:3)~
y<-c(5,20,45)~
df2<-data.frame(x,y)~
df2$Ymat<-cbind(df2$y,50-df2$y)~
summary(glm(Ymat~x,family=binomial,data=df2))~
~
このどちらの場合もデータとしては同じだと思いますし、推...
でも、residualsの値やAICなど、モデル全体の適合性に関わる...
//
-なぜ同じと思ったのか不思議ですね。後者は，データは3対で...
-それと，weights は本来の意味とも違うのでは？(resp=1,x=1)...
-あなたのやりたかったことは，むしろこちらのほうなんではな...
> x <- rep(1:3, each=50)
> y <- c(rep(1,5), rep(0,45), rep(1,20), rep(0,30), rep(...
> summary(glm(y ~ x, family=binomial))
Call:
glm(formula = y ~ x, family = binomial)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.0559 -0.4062 -0.4062 0.5075 2.2521
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -4.6727 0.7323 -6.381 1.76e-10 ***
x 2.2191 0.3366 6.594 4.29e-11 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1...
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 207.28 on 149 degrees of freedom
Residual deviance: 133.16 on 148 degrees of freedom
AIC: 137.16
Number of Fisher Scoring iterations: 4
違うか？ -- &new{2005-11-01 (火) 21:12:22};
-どうもありがとうございます。ご指摘の通りです。それをやり...
weightsは重み付けだけで、自由度が変化してないのですね。た...
重み付けだけというのが具体例としてぴんと来ません。わたし...
-いくつか例を当たってみましたが、"Introductory Statistics...
rate <- c(5/50, 20/50, 45/50)
total <- c(50,50,50)
x <- 1:3
summary(glm(rate~x,weights=total,binomial))
とする方が、よりスマートですね。１番目、２番目の方法では...
**有効桁数
>[[初心者K]] (2005-10-28 (金) 12:11:29)~
~
コンピュータでは2進数で数字が表現されているため、~
近似的な計算によって、桁落ちや丸め誤差が生じてしまうと聞...
それはコンピュータでの計算に一般的に言えることなのでしょ...
（やはりRでもアルゴリズムによって、同様の事が生じてしまう...
Rではどのような扱いになっているのか教えてください。~
//
-整数は正確に保持されています（どんな大きな整数でもという...
小数部分を持つ数も，小数部が2進数で表されるもの（0.5, 0.2...
それ以外の数は，"近似値"が保持されますが，実際上は有効桁...
ただ，アルゴリズムに気をつけないととんでもない結果が得ら...
よく知られていることではつい最近までの Excel の分散などを...
また別の例では，二次方程式の解を求める際に，単に解の公式...
よくうっかり落ち込むのは，演算結果が特定の定数値に等しい...
このようなことは R でもその他のソフトでも，あるいは C 言...
-桁落ちは，2 進数であろうと10進数であろうと生じます。10進...
-ありがとうございました。それと関連して追加の質問なのです...
-64bit 版もあるようですし。数値演算プロセッサを使って計算...
-数値演算はど素人(実際そうですのでご心配なく)ですが,だい...
-long doubleは処理系依存です. -- [[なかま]] &new{2005-10-...
-long double というのは定義上は4倍精度で34，5桁でしたっけ...
-そうはいってもアルゴリズムを気をつけなければならないとい...
-「何か他のもの」<- 意味不明で良い質問ではないですね。倍...
-ありがとうございました。頂いたコメントを参考に、またさら...
**ARMAモデルの最尤推定
>[[tanabe]] (2005-10-27 (木) 23:30:27)~
~
いつもお世話になっております。~
~
Rでは時系列のARおよびARIMAモデルに関してパラメータの最尤...
また、多変量ARMAモデルのパラメータ最尤推定に対応した機能...
さらに、欠損値を含むデータに対しても対応する機能があるの...
~
以上、よろしく御願いします。~
//
**時系列の成分分解
>[[kkondo]] (2005-09-29 (木) 13:33:10)~
~
１変量時系列データを、2階型トレンド成分、ARMA(p,q)成分、...
最終的には、AIC (もしくはBIC) が最適になるように、成分分...
//
-自己補足説明なのですが、標準的(?)な方法論を用いた方法を...
**プロットで使える文字
>近藤 (2005-09-28 (水) 23:47:54)~
~
expression() で使える文字の一覧というのは存在するのでしょ...
expression(cdots) は存在するのですが expression(cdot) は...
//
-help("plotmath")に一覧があります.demo(plotmath)で一覧が...
-印刷される場合は[[R-Tips:http://cse.naro.affrc.go.jp/tak...
-ありがとうございました。 中点というのは存在しないみたい...
**R-2.1.1 (Windows版)のコマンドラインからインストール
>[[BR103]] (2005-09-22 (木) 16:28:39)~
~
題名のように、GUIを使わずにインストールはできないのでし...
//
-rw2011.exe /VERYSILENT しかし、これはデフォルトでインス...
**read.dbfの不具合？
>[[谷村]] (2005-08-19 (金) 16:11:14)~
~
r-helpでスルーされたので、こちらに相談します。~
KEYCODE (数値型19桁小数0桁)
422010010
42201002101
42201002102
42201002103
42201002104
422010060
422010071
422010072
42201008001
42201008002
という内容のtest.dbfとtest.csvを用意します。&size(8){実物...
> library(foreign)
> cbind(read.csv("test.csv"),read.dbf("test.dbf"))
KEYCODE KEYCODE
1 422010010 422010010
2 42201002101 NA
3 42201002102 NA
4 42201002103 NA
5 42201002104 NA
6 422010060 422010060
7 422010071 422010071
8 422010072 422010072
9 42201008001 NA
10 42201008002 NA
となります。read.shapeで読み込んだ属性値がおかしいので、...
OpenOffice.org Calcまたはtxt2dbfでtest.dbfを作成して試し...
//
-dbfファイルとは、dBASEの形式でしょうか？ -- &new{2005-0...
-dbfはdBASE形式ファイルの拡張子です。 -- [[谷村]] &new{20...
-write.dbf で書き出したものとは違うのですね。そちらが作っ...
-添付できないので、私が管理するサーバにアップロードしまし...
-セキュリティホール対策のため、画像以外の添付機能は無効に...
> cbind(read.csv("t1.csv",header=F),read.dbf("t2.dbf"))
V1 N1
1 422010010 422010010
2 42201002101 42201002101
3 42201002102 42201002102
4 42201002103 42201002103
5 42201002104 42201002104
6 422010060 422010060
7 422010071 422010071
8 422010072 422010072
9 42201008001 42201008001
10 42201008002 42201008002
と問題ないようです。dbfはNeoOffice/J 1.1 RC Patch8で作成...
-上記からダウンロードした test.dbf を読み込んだところ、最...
-予想通り。4バイト整数として読んでいるようなので，Mac で...
> cbind(read.csv("test.csv"),read.dbf("test.dbf"))
KEYCODE KEYCODE
1 422010010.00000000 422010010
2 42201002101.00000000 2147483647
3 42201002102.00000000 2147483647
4 42201002103.00000000 2147483647
5 42201002104.00000000 2147483647
6 422010060.00000000 422010060
7 422010071.00000000 422010071
8 422010072.00000000 422010072
9 42201008001.00000000 2147483647
10 42201008002.00000000 2147483647
> x <- cbind(read.csv("test.csv"),read.dbf("test.dbf"))
> class(x[2,2])
[1] "integer"
岡田さんと同じ環境なんだけど，ファイルは，谷村さんのをダ...
-ご無沙汰です。> 岡田さん 追試験をありがとうございました...
$ dbf2txt test.dbf
#KEYCODE
422010010
42201002101
42201002102
42201002103
42201002104
422010060
422010071
422010072
42201008001
42201008002
dbfファイルは正常に見えます。oocalcで開いても正常に開きま...
-signed integer の桁あふれみたいですね。2147483647 = 0x7F...
-ちなみにppc64 Rでも 2147483647 でした... foreignのソース...
-正常に読めたdbfは [[ここ:http://epidemiology.md.tsukuba....
-岡田さんも青木先生もありがとうございます。原因が絞り込ま...
-下記のような試行をしてみました。-- [[谷村]] &new{2005-08...
> write.dbf(read.csv("test.csv"),file="test2.dbf")
> cbind(read.csv("test.csv"),read.dbf("test2.dbf"))
KEYCODE KEYCODE1
422010010 4220100102
42201002101 422010021013
42201002102 422010021024
42201002103 422010021035
42201002104 422010021046
422010060 4220100607
422010071 4220100718
422010072 4220100729
42201008001 42201008001
42201008002 42201008002
-岡田さんのと谷村さんのでは，同じ dBase 形式でも中身が違...
$ > hexdump t2.dbf
0000000 0305 0813 0a00 0000 4100 0f00 0000 0000
0000010 0000 0000 0000 0000 0000 0000 0000 0000
0000020 4e31 0000 0000 0000 0000 004e 0000 0000
0000030 0e02 0000 0000 0000 0000 0000 0000 0000
0000040 0d20 3432 3230 3130 3031 302e 3030 0000
0000050 2034 3232 3031 3030 3231 3031 2e30 3020
0000060 3432 3230 3130 3032 3130 322e 3030 2034
0000070 3232 3031 3030 3231 3033 2e30 3020 3432
0000080 3230 3130 3032 3130 342e 3030 2034 3232
0000090 3031 3030 3630 2e30 3000 0020 3432 3230
00000a0 3130 3037 312e 3030 0000 2034 3232 3031
00000b0 3030 3732 2e30 3000 0020 3432 3230 3130
00000c0 3038 3030 312e 3030 2034 3232 3031 3030
00000d0 3830 3032 2e30 3000
00000d7
$ > hexdump test.dbf
0000000 0369 0813 0a00 0000 4100 1400 0000 0000
0000010 0000 0000 0000 0000 0000 0000 0000 0000
0000020 4b45 5943 4f44 4500 0000 004e 0000 0000
0000030 1300 0000 0000 0000 0000 0000 0000 0000
0000040 0d20 2020 2020 2020 2020 2020 3432 3230
0000050 3130 3031 3020 2020 2020 2020 2020 3432
0000060 3230 3130 3032 3130 3120 2020 2020 2020
0000070 2020 3432 3230 3130 3032 3130 3220 2020
0000080 2020 2020 2020 3432 3230 3130 3032 3130
0000090 3320 2020 2020 2020 2020 3432 3230 3130
00000a0 3032 3130 3420 2020 2020 2020 2020 2020
00000b0 3432 3230 3130 3036 3020 2020 2020 2020
00000c0 2020 2020 3432 3230 3130 3037 3120 2020
00000d0 2020 2020 2020 2020 3432 3230 3130 3037
00000e0 3220 2020 2020 2020 2020 3432 3230 3130
00000f0 3038 3030 3120 2020 2020 2020 2020 3432
0000100 3230 3130 3038 3030 3200
0000109
dBase ファイルは比較的おとなしい構造をしている（上記）の...
-単純に scan で読んで，その後選択変換する。。。~
> x <- scan("test.dbf", what="")
> x
[1] "?003i?b?023" "" "A" "42201001...
[5] "42201002101" "42201002102" "42201002103" "42201002...
[9] "422010060" "422010071" "422010072" "42201008...
[13] "42201008002"
> as.numeric(x[-(1:3)])
[1] 422010010.00000000 42201002101.00000000
[3] 42201002102.00000000 42201002103.00000000
[5] 42201002104.00000000 422010060.00000000
[7] 422010071.00000000 422010072.00000000
[9] 42201008001.00000000 42201008002.00000000
などとねぇ。 -- [[青木繁伸]] &new{2005-08-19 (金) 17:26:5...
-数値型19桁小数0桁ではなく数値型19桁小数5桁にしてみると、...
-なるほど、scanを使う方法があるのですね。read.shapeの改造...
-無理矢理 foreign のソース自体を書き換えるという方法もあ...
-うんと,こんなんでいいんじゃないかなぁ -- [[なかま]] &new...
--- foreign.orig/src/dbfopen.c 2004-10-27 20:26:43.0000...
+++ foreign/src/dbfopen.c 2005-08-19 17:42:47.5496...
@@ -970,7 +970,8 @@
|| psDBF->pachFieldType[iField] == 'F' )
/* || psDBF->pachFieldType[iField] == 'D' ) D is...
{
- if( psDBF->panFieldDecimals[iField] > 0 )
+ if( psDBF->panFieldDecimals[iField] > 0 ||
+ psDBF->panFieldSize[iField] > 9 )
return( FTDouble );
else
return( FTInteger );
- かぶった(^^;
*** foreign/src/Rdbfread.c 2005-04-06 04:35:39.000000000...
--- foreign.new/src/Rdbfread.c 2005-08-19 17:39:25.00000...
***************
*** 89,95 ****
SET_VECTOR_ELT(df, nRvar, allocVector(STRSXP,nrec...
break;
case 2:
! SET_VECTOR_ELT(df, nRvar, allocVector(INTSXP,nrec...
break;
case 3:
SET_VECTOR_ELT(df, nRvar, allocVector(REALSXP,nre...
--- 89,95 ----
SET_VECTOR_ELT(df, nRvar, allocVector(STRSXP,nrec...
break;
case 2:
! SET_VECTOR_ELT(df, nRvar, allocVector(REALSXP,nre...
break;
case 3:
SET_VECTOR_ELT(df, nRvar, allocVector(REALSXP,nre...
***************
*** 120,129 ****
case 2:
if( DBFIsAttributeNULL( hDBF, iRecord, i ))
! INTEGER(VECTOR_ELT(df, nRvar))[iRecord] = NA_INT...
else
! INTEGER(VECTOR_ELT(df, nRvar))[iRecord] =
! DBFReadIntegerAttribute( hDBF, iRecord, i );
nRvar++;
break;
--- 120,129 ----
case 2:
if( DBFIsAttributeNULL( hDBF, iRecord, i ))
! REAL(VECTOR_ELT(df, nRvar))[iRecord] = NA_INTEGER;
else
! REAL(VECTOR_ELT(df, nRvar))[iRecord] =
! DBFReadDoubleAttribute( hDBF, iRecord, i );
nRvar++;
break;
-- [[岡田]] &new{2005-08-19 (金) 17:49:05};
-dBaseの仕様だと要するに10進数で,Rの方は整数型と浮動小数...
-う、輻輳。さー、どっちがいいでしょう。(笑) -- [[なかま]]...
-とりあえず、なかまさんの方を試しました。うまくいきました...
-ごった煮パッケージですから,普通にSummaryとして投稿しとけ...
-岡田さんのパッチも試しました。こっちは素直にパッチがあた...
-なかまさんのパッチのほうがスマートですから、そちらをフィ...
-r-helpにサマリを投稿しました。ありがとうございました。 -...
**A I D S モデル
>[[sonkun]] (2005-08-09 (火) 17:29:30)~
~
Rで需要分析のAlmost Ideal Demand System（A I D S モデル）...
どなたか詳しい方いらっしゃいましたらご助言よろしくお願い...
//
-トップページにリンクのある R site search を使いキーワー...
-パッケージmicEconに A I D S 関連の関数があります。 -- &...
**GARCHのパラメータ推定について
>[[松本]] (2005-07-23 (土) 15:19:13)~
~
以前も時系列分析について質問させて頂きました。~
GARCHモデルを使って、時系列分析をしようとしているのですが...
Rのバージョンは、2.1.1です。~
ただ、garchFit(データ、order=c(0,1))のように、ARCH推定は...
私の仕方が悪いのかさえ分かりませんので、教えていただけま...
また、ARMAモデルの次数設定の方法を以前教えて頂いたのです...
用いた導出法は、
~
for (j in 1:5)
for (i in 1:5) {
tmp <- arima(data, order=c(j, 0, i))
cat("order1=", j, "order2=", i, "AIC=", tmp$aic,...
}
//
-実際に操作した内容を，全部！！！コンソールからコピー・ペ...
- [[ヘルプ]]を参照して，美しく投稿しましょう。
-貼り付けたエラーメッセージは，文字化けしていませんね。R...
-rate というオブジェクトは，あなたが定義したもの？わたし...
-そうです。私が使っているのは、Windowsです。~
文字化けについては、ご指摘のように対処してみます。~
ただ、文字化けしていないメッセージを見ても分かりません。t...
気づきました。以前、ARIMAで同様の事をした際は、きちんとで...
どうも扱うデータがおかしいのかもしれないので、もう一度チ...
しかし、なぜGARCH推定はだめで、ARCH推定は可能であったので...
-rateというのは確かに私が定義したものです。具体的には、FR...
-何でかは，私にはわかりません。~
いずれにせよ，あなたには分かるべきですが，rate というもの...
一応，エラーの原因と思われるオブジェクトが示されたのです...
-rateというのは、ｔｘｔ形式のファイルです。 -- [[松本]] &...
-t x t 形式のファイルったって，それがどのような情報を持つ...
極端なことを言えば，ゼロ行列かもしれないし，単位行列かも...
先にもコメントしたように，誰も，あなたと同じ条件で追試で...
rate を作った手順そのものも，R コンソールからコピー・ペー...
要するに，あなた以外の人がその人の環境下で，あなたと同じ...
もっとも，みんなが暇で，あなたの問題につきあってあげよう...
そうでなくても，ここが問題じゃないの？くらいのアドバイス...
-答えが欲しいなら，回答があるかないか分からなくても，なる...
レスポンス返しても，それに対する応答が長い間ないと，もう...
-頻繁に回答のチェックはしています。 -- [[松本]] &new{2005...
-ただ、ここの編集機能の使い方が不慣れで、rateの説明を編集...
rateですが、txtファイルに数値を打ち込んだだけのものです。
全ての数値を打ち込めないので、部分的に挙げます。
126.01
126.24
126.06
126.23
127.38
127.57
127.73
128.14
128.52
129.58
129.95
130.52
131.27
131.18
130.92
132.08
131.54
131.1
というように数値が800個程縦にならんでいます。それ以外には...
一切ありません。このtxtファイルのファイル名は、rate6.txt...
というようにしました。長く時間をかけてすみません。-- [[...
-秘密事項でなければ，そのファイルを添付すればいかがですか...
-どうも
終了行:
COLOR(blue){SIZE(20){R および RjpWiki に関する質問コーナ...
COLOR(red){SIZE(20){このコーナーはスパム対策のため閉鎖し...
//COLOR(red){注意}：[[Ｑ＆Ａ]] コーナーが少し重くなり過ぎ...
COLOR(red){注意}：新規記事用の入力欄は以下の目次の直後に...
もし新規に質問をされる場合は、効率的にコメントを得るため...
- 広い意味で R に関係する話題に限定!
- 学校の宿題は自分で考えるべきです
- 具体的な状況がわかる様に背景説明をけちらない
- [[CRAN の検索エンジン:http://finzi.psych.upenn.edu/sear...
- 自分には馴染みでも、他の人には馴染みの無い手法・概念も...
- もし問題が解決した(結局解決しない場合を含め)場合はその...
- 後で他の人が参照する際の便宜のため、タイトルには簡潔で...
- 回答者は神通力を持ってはいません。「これこれしたが、う...
~
----
アーカイブ：[[Ｑ＆Ａ(旧1)]]~
COLOR(red){SIZE(20){このコーナーはスパム対策のため閉鎖し...
#contents
----
COLOR(red){注意}：投稿失敗した記事をページ先頭の「新規｜...
//#article
**文字式による数式の関数認識
>[[KOSHI]] (2006-01-18 (水) 22:58:51)~
~
expression()の形式で入力した数式を，関数として認識したい...
例えば
> f <- expression(x*y)
を（x,y）の関数に直したいとき，
> F <- function(x,y){}
> body(F) <- f
とやれば
> F <- (x, y) x * y
とできます．しかし，これはfで使われる変数が（x,y）と既知...
もし，fでの変数が未知の場合，これを関数認識するにはどうし...
> g <- G(expression(x*y*z*u))
などとすると,ちゃんとこの4つの変数（x,y,z,u）を使って~
> g <- function(x,y,z,u) x*y*z*u
となるような関数を作ることは可能でしょうか．~
~
ちなみに~
deriv(expr, name, func=T)
という関数を使うと，勝手に書いた数式~
expr=expression(...)
の，変数nameによる1階偏微分を，その文字による関数として認...
~
どなたかお力をお貸しいただけないでしょうか．~
//
-ソースを見ることについて，どこかのページにも書いてありま...
まず，関数名だけをタイプしてみる
> deriv
function (expr, ...)
UseMethod("deriv")
<environment: namespace:stats>
なんだぁ？これだけか？と思ったら，
> methods(deriv)
[1] deriv.default deriv.formula
と，やってみる。ここに示されるのが，実際に呼び出される候...
単に関数名をコンソールにタイプするのではなくて，以前に出...
> stats:::deriv.default
function (expr, namevec, function.arg = NULL, tag = ".ex...
hessian = FALSE, ...)
.Internal(deriv.default(expr, namevec, function.arg, tag...
<environment: namespace:stats>
もう一つ
> stats:::deriv.formula
function (expr, namevec, function.arg = NULL, tag = ".ex...
hessian = FALSE, ...)
{
if ((le <- length(expr)) > 1)
.Internal(deriv.default(expr[[le]], namevec, fun...
tag, hessian))
else stop("invalid formula in deriv")
}
<environment: namespace:stats>
これらの関数のどれがどのように使われるかは，help で検索の...
ここで示されるのは，R の範囲で書かれるプログラムであって...
そのような部分についてのプログラムの詳細は，（えい面倒だ...
てなことかな（はずしたか。。） -- &new{2006-01-18 (水) 2...
-それとも，あなたは，このようなことがお望みなのかな？ ...
> func <- as.function(a <- function(x, y) x+y)
> func
function(x, y) x+y
> func(2,4)
[1] 6
> func <- as.function(a <- function(x, y) x*y)
> func(3,7)
[1] 21
一応そういうことで。 -- &new{2006-01-18 (水) 23:36:17};
-私の場合，ソースを見ても分からない部類だと思います・・・...
-実際，deriv（）はどんな文字を用いて数式を書いても，それ...
-具体的には，どんなものですか？上の，func だって，x,y を...
-＞実際，deriv（）はどんな文字を用いて数式を書いても，そ...
derive を読んでください(^_^;)~
もしくは，ウイザード召喚！！(^_^;) -- &new{2006-01-18 (...
-たとえば，ユーザが f <- expression(x*y) と入力すること...
> G <- function(expr){
F <- function(x,y){ }
body(F) <- f
return(F)
}
というGを用意しておけば(x,yの数式なら何でも)関数に出来ま...
> f <- expression(u*v)
と書いてしまえば，G(f)は(u,v)の関数に出来ないじゃないです...
-関数定義に出てくる名前をパースして，それぞれが想定してい...
-deriv も何が変数（とみなすべき）かを指示する仕様になって...
> a <- 3; b <- 5
> f <- alist(x=,y=a,x*y+b)
> F <- as.function(f)
> F
function (x, y = a)
x * y + b
> F(1)
[1] 8
-＞ alist を文字列から合成する方法.....一時ファイルに書き...
-なるほど。一時ファイル名をその場で作る関数が確かあったは...
-tempfile 関数ですね。。 -- &new{2006-01-20 (金) 10:56:4...
-皆様、いろいろとご助言有難うございます。自分でも少し進展...
> Func <- function(expr){
df <- deriv(expr, all.vars(expr), func=T) #とり...
F <- function(p,q){ #ここを変数２個にしておく。
df(p,q)[1] #derivの[1]には微分しない関数が返...
}
return(F)
}
こうして関数Fを作れば、ユーザが入力する「変数の個数」さえ...
あとは、自然数nを与えられた時に、自動的にn変数の空の関数...
> body(G(n)) <- df(x1,...xn)[1]
で完成するのではないでしょうか？~
しかし、こういうG(n)ってできますでしょうか？~
forとか使っても関数の変数を増やせなくて・・・あと一歩だと...
よろしくお願いします。 -- [[KOSHI]] &new{2006-01-20 (金) ...
-上に示唆されていた方法を使えば（ファイルはテキスト、文字...
> G <- function(n) {
tmpfile <- tempfile() # 一...
.x <- paste("x", 1:n, ",", sep="") # 引...
.x[n] <- substr(.x[n], 1, nchar(.x[n])-1) # 最...
sink(tmpfile)
cat("function (", .x, ") {}?n") # 関...
sink() # フ...
.f <- source(tmpfile)[[1]] # フ...
unlink(tmpfile) # 一...
return(.f)
}
> f5 <- G(5)
> f5
function ( x1, x2, x3, x4, x5 ) {}
> str(f5)
function (x1, x2, x3, x4, x5)
- attr(*, "source")= chr "function ( x1, x2, x3, x4, x...
- G<-function(n)eval(parse(text=paste("function(",paste("...
-↑ 上手い！ しかし、ファイルに一旦書き出す方法は（少し野...
-どうも，このページに粘着質のスパム投稿者が住み着いたよう...
-Q & A (中級) とでも改名して再出発しましょうか。みたとこ...
-スパム対策のためＱ＆Ａコーナーは閉鎖します。以後は [[Ｑ...
**共分散分析
>[[S.O]] (2006-01-13 (金) 03:10:39)~
~
共分散分析のやり方が、いまいち分かりません。~
yが目的変数で、x1,x2は連続量で、x3,x4は離散です。~
共分散分析は、lmよりaovの方がよいと言われていたので、以下...
aov(y~x1+x2+x3+x4)~
この場合、普通の４元配置分散分析だと思い、x1,x2を共変量に...
お手数ですが、どうぞ宜しくお願いいたします。~
//
-「いまいちわからない」ではこっちも「いまいちわからない」...
**Mac版のR2WinBUGS
>[[まさひと]] (2006-01-11 (水) 20:54:49)~
~
パッケージではMac用のR2WinBUGSが配付されいますが、もとも...
//
-example(bugs) で，なんか動いているけど。よくわかりません...
-wine経由で動かす(system関数使って)物です. 私の手元でも(D...
-Darwineをインストールしてみました. GUIのバージョンの関係...
-なかまさん、ありがとうございます。さっそく試してみます。...
-あー、すいません。よく見るとまだまだ普通には動かんようで...
**Rで＠Risk的なパッケージを自作
>[[Tak]] (2006-01-10 (火) 18:43:58)~
~
リスク解析の勉強をかねて、Rで＠Riskの様なパッケージを自作...
http://www.palisade-europe.com/~
~
既存のRのパッケージをできるだけ使って、＠Riskの様な包括的...
~
他にparametric sensitivity analysis、モンテカルロ法、deci...
//
-Palisade社が何らかの特許を申請していた場合、完全移植には...
-確かにPalisade社の名前をあげる必要はなかったですね。＠Ri...
**plotのwidthの種類
>[[浩太]] (2006-01-05 (木) 13:52:55)~
~
barplotの棒の幅を定めるwidthって何種類ぐらいの幅がありま...
具体的にどんな種類の幅がありますか？~
教えてください~
//
-質問の意味がわかりません。widthは数値ベクトル（リサイク...
> x <- c(1,4,8,5,2)
> barplot(x, width=c(1.2, 3.4, 5.321, 2.345))
これにより以下のような図が得られるんだが。。。
#ref(barplot.png)
それとも，
> barplot(x, width=1)
> barplot(x, width=1000)
のような描画では差がないように見えるからということかな？~
ヘルプを読めばわかるように，xlim で横軸の長さを決めない限...
Specifying a single value will no visible effect unless x...
-ｐarを使って図を重ね合わせる時に差をつけるときも数値でや...
-浩太さんがparで図を重ねることができる人ならば、実際に試...
-四つのbarplotをparを使って重ねあわせました。それぞれの違...
-用いたプログラムをここに載せればいかがですか。長すぎたり...
なお，投稿の失敗は，上のほうの「編集」をクリックすること...
-一般的に言えば，幅を変えると必然的に面積も変わるわけで，...
**barplotでNameが表示されない
>[[鍋]] (2006-01-03 (火) 14:29:13)~
~
> mydata
[1] 29 24 20 18 25 22 23 20 20 23 20 30 21 20 19 26 14 2...
[26] 2 2 4 1 3
> main
[1] "Turnover Number"
> names.arg
[1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 1...
[26] 26 27 28 29 30 31 32
> xlab
[1] "Event Time Number"
> ylab
[1] "Turnover"
> ylim
[1] 0 32
> barplot(mydata, main=main, names.arg=names.arg, xlab=x...
Error in barplot.default(mydata, main = main, names.arg ...
incorrect number of names
なぜこのようなErrorが表示されるのかわかりません。~
incorrectになるのはどういったときなのでしょうか？~
//
-名前の数が合わないと申しておりますし,実際合ってないよう...
-length(mydata) の結果と，length(names.arg) の結果を比べ...
-国際化された R を使えば，エラーメッセージも日本語で出る...
以下にエラーbarplot.default(mydata, main = main, names.a...
名前の数が誤りです
-- &new{2006-01-03 (火) 15:41:29};
-鍋さんの悪い癖は，コメントに反応しないことだね。 -- &ne...
**R. Tierney 氏の R バイトコンパイラーについて
>[[ma]] (2005-12-29 (木) 18:12:39)~
~
何でも掲示板に出ている「R. Tierney 氏の R バイトコンパイ...
~
”bad function”というエラーが出た場合には具体的にどのよう...
漠然とした質問で申し訳ないのですが、どなたかお教え頂けな...
//
-少々説明を追加させて下さい。cmpfun()でバイトコンパイルは...
-もしあまり長くなければ、エラーが起きる実例を紹介するとコ...
**軸の省略
>[[輝]] (2005-12-23 (金) 16:33:08)~
~
barplot()で，図を描いており，一つだけ大きい項目があり，~
途中のY軸(目盛り）を省きたいのですが，~
どのようにすればよいのですか？~
//
-[[Ｑ＆Ａ (初級者コース)#content_1_29]]にある，「ヒストグ...
-先の質問をしたAkiraです。外れ値がある場合は対数変換とい...
**RExcelの起動時エラー
>[[YK]] (2005-12-23 (金) 12:52:29)~
~
RExcelを利用しようとしたのですが，Excel起動時にエラーが発...
~
エラーメッセージは，以下のとおりです．
「実行時エラー 5 プロシージャの呼び出し、または引数が不...
そのエラーの発生箇所は，☆の付いている箇所です．
Private Sub RemoveRMenu()
On Error Resume Next
☆CommandBars(1).Controls("RExcel").Delete
End Sub
RExcelの実行した環境は以下のとおりです．
Windows XP(SP2)
Office Excel 2003
R 2.2.1
RServ 2.0.0（Rserv1.3.5でも同様にエラーが発生しました）
RExcelに添付されていたExcel 02 - Macro Demoを実行した際に...
~
エラー修正，回避などのヒントをご教授願えないでしょうか。~
よろしくお願いいたします。~
//
-思いつきです。3年ほど前の話ですが、RExcelはエラーが多発...
-ありがとうございます．バージョンによるエラーを視野に入れ...
**MASS lda 判別関数と偏F値について
>[[mt]] (2005-12-21 (水) 11:44:53)~
~
MASS ライブラリ中の lda() 関数を実行した結果から、判別関...
R中ではpredictで判別してくれるのでいいのですが、実際どの...
どなたかヒントなどご教授願えないでしょうか。~
よろしくお願いいたします。~
//
-lda は，正準判別分析を行っていますね。ということで，偏F...
**R2.2.0でのeps出力について。
>[[KJ]] (2005-12-20 (火) 21:37:18)~
~
Windows 版の R-2.2.0 で eps 形式での出力を行ったファイル...
それ以前のバージョンの R であれば問題なく表示できるのです...
原因は何かわかりますでしょうか？
GSview 4.7 2005-03-26
Unknown in Comments section at line 2:
%%DocumentNeededResources: font Times-Roman
Unknown in Comments section at line 3:
%%+ font Times-Bold
Unknown in Comments section at line 4:
%%+ font Times-Italic
中略
End offending input ---
file offset = 5390
gsapi_run_string_continue returns -101
~
//
-うーん, gsの環境のような気がするんですが...もし差し支え...
-ありがとうございます。
x <- seq(-3,3,by=0.01)
plot(x,dnorm(x),type="l")
として、右クリックでeps保存したファイルです。~ -- [[KJ]]
#ref(http://www.sfc.keio.ac.jp/~kj198276/test.eps)
-中は闇Ｒ-2.2.0のようですが,多分GSがCID未対応の設定もしく...
setHook(packageEvent("grDevices", "onLoad"),
function(...) grDevices::ps.options(cidfam...
等として日本語環境を無効にするか, gsが日本語を扱えるよう...
あげてみてください。 -- [[なかま]] &new{2005-12-21 (水) 2...
-あ、闇版を使わないと言う手もあります。(でも、2.3.0になる...
-知識不足で、仰っている意味が捉えきれません…ちょっと調べ...
-やはりGSの設定の問題でした。情報ありがとうございました！...
**barplot( )の凡例位置の変更
>[[兄歯]] (2005-12-17 (土) 14:17:53)~
~
barplot( )の凡例位置がグラフと重なってしまうので，凡例の...
//
-legend()はいかがでしょうか？~
x <- matrix(rnorm(100), ncol=10)
dimnames(x) <- list(1:10, 1:10)
barplot(x, beside=TRUE, legend.text=rownames(x), col=1:1...
barplot(x, beside=TRUE, col=1:10)
legend("topright", legend=rownames(x), pch=15, col=1:10)...
legend("topleft", legend=rownames(x), pch=15, col=1:10)#...
-- [[Akira]] &new{2005-12-17 (土) 15:33:28};
-解決しました。どうもありがとうございます。 -- [[兄歯]] &...
**軸の分割数について
>[[内藤]] (2005-12-16 (金) 16:57:40)~
~
plot( )で，軸（目盛り)の分割を任意に与えたいのですが，~
どのようにすればよいのですか？~
//
-? axis をやってみればよろしいかと。
hist(rnorm(1000),xaxt="n")
axis(1,pos=0,at=c(-1,0.5,2.8),labels=TRUE)
などで，遊んでみる。 -- &new{2005-12-16 (金) 18:00:36};
-どうもありがとうございます。助かりました。 -- [[内藤]] &...
**maptoolsの座標の参照方法
>[[traveller]] (2005-12-06 (火) 12:42:34)~
~
maptoolsの各図形を参照する方法、どなたか、教えていただ...
//
-各図形を参照するとは？もう少し具体的に。何をしたいのかも...
-ポイント、ラインおよびポリゴンの座標、BBOXを参照する方法...
-maptoolsで扱うObjectはMap Obejctと呼ばれ、"Shapes"と"att...
> x <- read.shape(system.file("shapes/sids.shp", package...
> x$Shapes
[[1]]
[[1]]$Pstart
[1] 0
[[1]]$verts
[,1] [,2]
[1,] -81.47276 36.23436
[2,] -81.54084 36.27251
[3,] -81.56198 36.27359
[略]
[26,] -81.45289 36.23959
[27,] -81.47276 36.23436
[[1]]$shp.type
[1] 5
[[1]]$nVerts
[1] 27
[[1]]$nParts
[1] 1
[[1]]$bbox
[1] -81.74107 36.23436 -81.23989 36.58965
[[1]]$shpID
[1] 0
attr(,"nVerts")
[1] 27
attr(,"nParts")
[1] 1
attr(,"shp.type")
[1] 5
attr(,"bbox")
[1] -81.74107 36.23436 -81.23989 36.58965
[以下略]
従って、例えば、1つ目のポリゴンの座標を得る場合には
> x$Shapes[[1]]$verts
でOKです。BBOX座標が必要な場合は、
> x$Shapes[[1]]$bbox
ですね。
-谷村さん、ご親切なご回答ありがとうございます。 -- [[trav...
**ベイズクラスタ
>[[ド級初心者]] (2005-11-07 (月) 12:58:10)~
~
Rでベイズクラスタ分析が可能と風の噂で聞きましたが、どのプ...
//
**コーフェン相関係数の求め方
>[[mori]] (2005-11-02 (水) 14:08:03)~
~
階層的クラスタリングでできた樹形図（デンドログラム）の結...
//
-知らないものだったのでググッてみた。~
Sneath and Sokal『数理分類学』（1994年，内田...
ひどい訳本みたいで，三中先生に「だいたいこの訳本，トンデ...
[[http://cse.niaes.affrc.go.jp/minaka/diary2003-09.html:h...
[[http://www.i-juse.co.jp/statistics/support/sympo/m10/m1...
[[http://www.sci.kagoshima-u.ac.jp/~dllsa/jp/dlls/cluscrt...
-ちゃんとした名前がわかれば，さらにググッてみることができ...
-ありがとうございます。その関数のようです。綴りはcoffenだ...
-"どーして「コーフェン」になるの" -> これは multicolinear...
**ロジスティック回帰での変数指定方法とresidualsなどの値
>[[たかくら]] (2005-11-01 (火) 19:00:18)~
~
glmでロジスティック回帰をしていてほとんど同じ（？）計算...
例えば、独立変数xが１，２，３のそれぞれの場合に、50試行の...
このとき、試行数や成功例数の与え方としてweightsパラメータ...
~
resp<-rep(c(1,0),3)~
x<-(1,1,2,2,3,3)~
weight<-c(5,45,20,30,45,5)~
df1<-data.frame(resp,x,weight)~
summary(glm(resp~x,weights=weight,family=binomial,data=df...
~
x<-(1:3)~
y<-c(5,20,45)~
df2<-data.frame(x,y)~
df2$Ymat<-cbind(df2$y,50-df2$y)~
summary(glm(Ymat~x,family=binomial,data=df2))~
~
このどちらの場合もデータとしては同じだと思いますし、推...
でも、residualsの値やAICなど、モデル全体の適合性に関わる...
//
-なぜ同じと思ったのか不思議ですね。後者は，データは3対で...
-それと，weights は本来の意味とも違うのでは？(resp=1,x=1)...
-あなたのやりたかったことは，むしろこちらのほうなんではな...
> x <- rep(1:3, each=50)
> y <- c(rep(1,5), rep(0,45), rep(1,20), rep(0,30), rep(...
> summary(glm(y ~ x, family=binomial))
Call:
glm(formula = y ~ x, family = binomial)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.0559 -0.4062 -0.4062 0.5075 2.2521
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -4.6727 0.7323 -6.381 1.76e-10 ***
x 2.2191 0.3366 6.594 4.29e-11 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1...
(Dispersion parameter for binomial family taken to be 1)
Null deviance: 207.28 on 149 degrees of freedom
Residual deviance: 133.16 on 148 degrees of freedom
AIC: 137.16
Number of Fisher Scoring iterations: 4
違うか？ -- &new{2005-11-01 (火) 21:12:22};
-どうもありがとうございます。ご指摘の通りです。それをやり...
weightsは重み付けだけで、自由度が変化してないのですね。た...
重み付けだけというのが具体例としてぴんと来ません。わたし...
-いくつか例を当たってみましたが、"Introductory Statistics...
rate <- c(5/50, 20/50, 45/50)
total <- c(50,50,50)
x <- 1:3
summary(glm(rate~x,weights=total,binomial))
とする方が、よりスマートですね。１番目、２番目の方法では...
**有効桁数
>[[初心者K]] (2005-10-28 (金) 12:11:29)~
~
コンピュータでは2進数で数字が表現されているため、~
近似的な計算によって、桁落ちや丸め誤差が生じてしまうと聞...
それはコンピュータでの計算に一般的に言えることなのでしょ...
（やはりRでもアルゴリズムによって、同様の事が生じてしまう...
Rではどのような扱いになっているのか教えてください。~
//
-整数は正確に保持されています（どんな大きな整数でもという...
小数部分を持つ数も，小数部が2進数で表されるもの（0.5, 0.2...
それ以外の数は，"近似値"が保持されますが，実際上は有効桁...
ただ，アルゴリズムに気をつけないととんでもない結果が得ら...
よく知られていることではつい最近までの Excel の分散などを...
また別の例では，二次方程式の解を求める際に，単に解の公式...
よくうっかり落ち込むのは，演算結果が特定の定数値に等しい...
このようなことは R でもその他のソフトでも，あるいは C 言...
-桁落ちは，2 進数であろうと10進数であろうと生じます。10進...
-ありがとうございました。それと関連して追加の質問なのです...
-64bit 版もあるようですし。数値演算プロセッサを使って計算...
-数値演算はど素人(実際そうですのでご心配なく)ですが,だい...
-long doubleは処理系依存です. -- [[なかま]] &new{2005-10-...
-long double というのは定義上は4倍精度で34，5桁でしたっけ...
-そうはいってもアルゴリズムを気をつけなければならないとい...
-「何か他のもの」<- 意味不明で良い質問ではないですね。倍...
-ありがとうございました。頂いたコメントを参考に、またさら...
**ARMAモデルの最尤推定
>[[tanabe]] (2005-10-27 (木) 23:30:27)~
~
いつもお世話になっております。~
~
Rでは時系列のARおよびARIMAモデルに関してパラメータの最尤...
また、多変量ARMAモデルのパラメータ最尤推定に対応した機能...
さらに、欠損値を含むデータに対しても対応する機能があるの...
~
以上、よろしく御願いします。~
//
**時系列の成分分解
>[[kkondo]] (2005-09-29 (木) 13:33:10)~
~
１変量時系列データを、2階型トレンド成分、ARMA(p,q)成分、...
最終的には、AIC (もしくはBIC) が最適になるように、成分分...
//
-自己補足説明なのですが、標準的(?)な方法論を用いた方法を...
**プロットで使える文字
>近藤 (2005-09-28 (水) 23:47:54)~
~
expression() で使える文字の一覧というのは存在するのでしょ...
expression(cdots) は存在するのですが expression(cdot) は...
//
-help("plotmath")に一覧があります.demo(plotmath)で一覧が...
-印刷される場合は[[R-Tips:http://cse.naro.affrc.go.jp/tak...
-ありがとうございました。 中点というのは存在しないみたい...
**R-2.1.1 (Windows版)のコマンドラインからインストール
>[[BR103]] (2005-09-22 (木) 16:28:39)~
~
題名のように、GUIを使わずにインストールはできないのでし...
//
-rw2011.exe /VERYSILENT しかし、これはデフォルトでインス...
**read.dbfの不具合？
>[[谷村]] (2005-08-19 (金) 16:11:14)~
~
r-helpでスルーされたので、こちらに相談します。~
KEYCODE (数値型19桁小数0桁)
422010010
42201002101
42201002102
42201002103
42201002104
422010060
422010071
422010072
42201008001
42201008002
という内容のtest.dbfとtest.csvを用意します。&size(8){実物...
> library(foreign)
> cbind(read.csv("test.csv"),read.dbf("test.dbf"))
KEYCODE KEYCODE
1 422010010 422010010
2 42201002101 NA
3 42201002102 NA
4 42201002103 NA
5 42201002104 NA
6 422010060 422010060
7 422010071 422010071
8 422010072 422010072
9 42201008001 NA
10 42201008002 NA
となります。read.shapeで読み込んだ属性値がおかしいので、...
OpenOffice.org Calcまたはtxt2dbfでtest.dbfを作成して試し...
//
-dbfファイルとは、dBASEの形式でしょうか？ -- &new{2005-0...
-dbfはdBASE形式ファイルの拡張子です。 -- [[谷村]] &new{20...
-write.dbf で書き出したものとは違うのですね。そちらが作っ...
-添付できないので、私が管理するサーバにアップロードしまし...
-セキュリティホール対策のため、画像以外の添付機能は無効に...
> cbind(read.csv("t1.csv",header=F),read.dbf("t2.dbf"))
V1 N1
1 422010010 422010010
2 42201002101 42201002101
3 42201002102 42201002102
4 42201002103 42201002103
5 42201002104 42201002104
6 422010060 422010060
7 422010071 422010071
8 422010072 422010072
9 42201008001 42201008001
10 42201008002 42201008002
と問題ないようです。dbfはNeoOffice/J 1.1 RC Patch8で作成...
-上記からダウンロードした test.dbf を読み込んだところ、最...
-予想通り。4バイト整数として読んでいるようなので，Mac で...
> cbind(read.csv("test.csv"),read.dbf("test.dbf"))
KEYCODE KEYCODE
1 422010010.00000000 422010010
2 42201002101.00000000 2147483647
3 42201002102.00000000 2147483647
4 42201002103.00000000 2147483647
5 42201002104.00000000 2147483647
6 422010060.00000000 422010060
7 422010071.00000000 422010071
8 422010072.00000000 422010072
9 42201008001.00000000 2147483647
10 42201008002.00000000 2147483647
> x <- cbind(read.csv("test.csv"),read.dbf("test.dbf"))
> class(x[2,2])
[1] "integer"
岡田さんと同じ環境なんだけど，ファイルは，谷村さんのをダ...
-ご無沙汰です。> 岡田さん 追試験をありがとうございました...
$ dbf2txt test.dbf
#KEYCODE
422010010
42201002101
42201002102
42201002103
42201002104
422010060
422010071
422010072
42201008001
42201008002
dbfファイルは正常に見えます。oocalcで開いても正常に開きま...
-signed integer の桁あふれみたいですね。2147483647 = 0x7F...
-ちなみにppc64 Rでも 2147483647 でした... foreignのソース...
-正常に読めたdbfは [[ここ:http://epidemiology.md.tsukuba....
-岡田さんも青木先生もありがとうございます。原因が絞り込ま...
-下記のような試行をしてみました。-- [[谷村]] &new{2005-08...
> write.dbf(read.csv("test.csv"),file="test2.dbf")
> cbind(read.csv("test.csv"),read.dbf("test2.dbf"))
KEYCODE KEYCODE1
422010010 4220100102
42201002101 422010021013
42201002102 422010021024
42201002103 422010021035
42201002104 422010021046
422010060 4220100607
422010071 4220100718
422010072 4220100729
42201008001 42201008001
42201008002 42201008002
-岡田さんのと谷村さんのでは，同じ dBase 形式でも中身が違...
$ > hexdump t2.dbf
0000000 0305 0813 0a00 0000 4100 0f00 0000 0000
0000010 0000 0000 0000 0000 0000 0000 0000 0000
0000020 4e31 0000 0000 0000 0000 004e 0000 0000
0000030 0e02 0000 0000 0000 0000 0000 0000 0000
0000040 0d20 3432 3230 3130 3031 302e 3030 0000
0000050 2034 3232 3031 3030 3231 3031 2e30 3020
0000060 3432 3230 3130 3032 3130 322e 3030 2034
0000070 3232 3031 3030 3231 3033 2e30 3020 3432
0000080 3230 3130 3032 3130 342e 3030 2034 3232
0000090 3031 3030 3630 2e30 3000 0020 3432 3230
00000a0 3130 3037 312e 3030 0000 2034 3232 3031
00000b0 3030 3732 2e30 3000 0020 3432 3230 3130
00000c0 3038 3030 312e 3030 2034 3232 3031 3030
00000d0 3830 3032 2e30 3000
00000d7
$ > hexdump test.dbf
0000000 0369 0813 0a00 0000 4100 1400 0000 0000
0000010 0000 0000 0000 0000 0000 0000 0000 0000
0000020 4b45 5943 4f44 4500 0000 004e 0000 0000
0000030 1300 0000 0000 0000 0000 0000 0000 0000
0000040 0d20 2020 2020 2020 2020 2020 3432 3230
0000050 3130 3031 3020 2020 2020 2020 2020 3432
0000060 3230 3130 3032 3130 3120 2020 2020 2020
0000070 2020 3432 3230 3130 3032 3130 3220 2020
0000080 2020 2020 2020 3432 3230 3130 3032 3130
0000090 3320 2020 2020 2020 2020 3432 3230 3130
00000a0 3032 3130 3420 2020 2020 2020 2020 2020
00000b0 3432 3230 3130 3036 3020 2020 2020 2020
00000c0 2020 2020 3432 3230 3130 3037 3120 2020
00000d0 2020 2020 2020 2020 3432 3230 3130 3037
00000e0 3220 2020 2020 2020 2020 3432 3230 3130
00000f0 3038 3030 3120 2020 2020 2020 2020 3432
0000100 3230 3130 3038 3030 3200
0000109
dBase ファイルは比較的おとなしい構造をしている（上記）の...
-単純に scan で読んで，その後選択変換する。。。~
> x <- scan("test.dbf", what="")
> x
[1] "?003i?b?023" "" "A" "42201001...
[5] "42201002101" "42201002102" "42201002103" "42201002...
[9] "422010060" "422010071" "422010072" "42201008...
[13] "42201008002"
> as.numeric(x[-(1:3)])
[1] 422010010.00000000 42201002101.00000000
[3] 42201002102.00000000 42201002103.00000000
[5] 42201002104.00000000 422010060.00000000
[7] 422010071.00000000 422010072.00000000
[9] 42201008001.00000000 42201008002.00000000
などとねぇ。 -- [[青木繁伸]] &new{2005-08-19 (金) 17:26:5...
-数値型19桁小数0桁ではなく数値型19桁小数5桁にしてみると、...
-なるほど、scanを使う方法があるのですね。read.shapeの改造...
-無理矢理 foreign のソース自体を書き換えるという方法もあ...
-うんと,こんなんでいいんじゃないかなぁ -- [[なかま]] &new...
--- foreign.orig/src/dbfopen.c 2004-10-27 20:26:43.0000...
+++ foreign/src/dbfopen.c 2005-08-19 17:42:47.5496...
@@ -970,7 +970,8 @@
|| psDBF->pachFieldType[iField] == 'F' )
/* || psDBF->pachFieldType[iField] == 'D' ) D is...
{
- if( psDBF->panFieldDecimals[iField] > 0 )
+ if( psDBF->panFieldDecimals[iField] > 0 ||
+ psDBF->panFieldSize[iField] > 9 )
return( FTDouble );
else
return( FTInteger );
- かぶった(^^;
*** foreign/src/Rdbfread.c 2005-04-06 04:35:39.000000000...
--- foreign.new/src/Rdbfread.c 2005-08-19 17:39:25.00000...
***************
*** 89,95 ****
SET_VECTOR_ELT(df, nRvar, allocVector(STRSXP,nrec...
break;
case 2:
! SET_VECTOR_ELT(df, nRvar, allocVector(INTSXP,nrec...
break;
case 3:
SET_VECTOR_ELT(df, nRvar, allocVector(REALSXP,nre...
--- 89,95 ----
SET_VECTOR_ELT(df, nRvar, allocVector(STRSXP,nrec...
break;
case 2:
! SET_VECTOR_ELT(df, nRvar, allocVector(REALSXP,nre...
break;
case 3:
SET_VECTOR_ELT(df, nRvar, allocVector(REALSXP,nre...
***************
*** 120,129 ****
case 2:
if( DBFIsAttributeNULL( hDBF, iRecord, i ))
! INTEGER(VECTOR_ELT(df, nRvar))[iRecord] = NA_INT...
else
! INTEGER(VECTOR_ELT(df, nRvar))[iRecord] =
! DBFReadIntegerAttribute( hDBF, iRecord, i );
nRvar++;
break;
--- 120,129 ----
case 2:
if( DBFIsAttributeNULL( hDBF, iRecord, i ))
! REAL(VECTOR_ELT(df, nRvar))[iRecord] = NA_INTEGER;
else
! REAL(VECTOR_ELT(df, nRvar))[iRecord] =
! DBFReadDoubleAttribute( hDBF, iRecord, i );
nRvar++;
break;
-- [[岡田]] &new{2005-08-19 (金) 17:49:05};
-dBaseの仕様だと要するに10進数で,Rの方は整数型と浮動小数...
-う、輻輳。さー、どっちがいいでしょう。(笑) -- [[なかま]]...
-とりあえず、なかまさんの方を試しました。うまくいきました...
-ごった煮パッケージですから,普通にSummaryとして投稿しとけ...
-岡田さんのパッチも試しました。こっちは素直にパッチがあた...
-なかまさんのパッチのほうがスマートですから、そちらをフィ...
-r-helpにサマリを投稿しました。ありがとうございました。 -...
**A I D S モデル
>[[sonkun]] (2005-08-09 (火) 17:29:30)~
~
Rで需要分析のAlmost Ideal Demand System（A I D S モデル）...
どなたか詳しい方いらっしゃいましたらご助言よろしくお願い...
//
-トップページにリンクのある R site search を使いキーワー...
-パッケージmicEconに A I D S 関連の関数があります。 -- &...
**GARCHのパラメータ推定について
>[[松本]] (2005-07-23 (土) 15:19:13)~
~
以前も時系列分析について質問させて頂きました。~
GARCHモデルを使って、時系列分析をしようとしているのですが...
Rのバージョンは、2.1.1です。~
ただ、garchFit(データ、order=c(0,1))のように、ARCH推定は...
私の仕方が悪いのかさえ分かりませんので、教えていただけま...
また、ARMAモデルの次数設定の方法を以前教えて頂いたのです...
用いた導出法は、
~
for (j in 1:5)
for (i in 1:5) {
tmp <- arima(data, order=c(j, 0, i))
cat("order1=", j, "order2=", i, "AIC=", tmp$aic,...
}
//
-実際に操作した内容を，全部！！！コンソールからコピー・ペ...
- [[ヘルプ]]を参照して，美しく投稿しましょう。
-貼り付けたエラーメッセージは，文字化けしていませんね。R...
-rate というオブジェクトは，あなたが定義したもの？わたし...
-そうです。私が使っているのは、Windowsです。~
文字化けについては、ご指摘のように対処してみます。~
ただ、文字化けしていないメッセージを見ても分かりません。t...
気づきました。以前、ARIMAで同様の事をした際は、きちんとで...
どうも扱うデータがおかしいのかもしれないので、もう一度チ...
しかし、なぜGARCH推定はだめで、ARCH推定は可能であったので...
-rateというのは確かに私が定義したものです。具体的には、FR...
-何でかは，私にはわかりません。~
いずれにせよ，あなたには分かるべきですが，rate というもの...
一応，エラーの原因と思われるオブジェクトが示されたのです...
-rateというのは、ｔｘｔ形式のファイルです。 -- [[松本]] &...
-t x t 形式のファイルったって，それがどのような情報を持つ...
極端なことを言えば，ゼロ行列かもしれないし，単位行列かも...
先にもコメントしたように，誰も，あなたと同じ条件で追試で...
rate を作った手順そのものも，R コンソールからコピー・ペー...
要するに，あなた以外の人がその人の環境下で，あなたと同じ...
もっとも，みんなが暇で，あなたの問題につきあってあげよう...
そうでなくても，ここが問題じゃないの？くらいのアドバイス...
-答えが欲しいなら，回答があるかないか分からなくても，なる...
レスポンス返しても，それに対する応答が長い間ないと，もう...
-頻繁に回答のチェックはしています。 -- [[松本]] &new{2005...
-ただ、ここの編集機能の使い方が不慣れで、rateの説明を編集...
rateですが、txtファイルに数値を打ち込んだだけのものです。
全ての数値を打ち込めないので、部分的に挙げます。
126.01
126.24
126.06
126.23
127.38
127.57
127.73
128.14
128.52
129.58
129.95
130.52
131.27
131.18
130.92
132.08
131.54
131.1
というように数値が800個程縦にならんでいます。それ以外には...
一切ありません。このtxtファイルのファイル名は、rate6.txt...
というようにしました。長く時間をかけてすみません。-- [[...
-秘密事項でなければ，そのファイルを添付すればいかがですか...
-どうも
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
