タグ付けされた質問 「regression」

1つ(または複数)の「従属」変数と「独立」変数の間の関係を分析する手法。

2
ロジスティック回帰における高度に歪んだデータセットの重みの追加
入力変数をバイナリ出力変数に合わせるために、標準バージョンのロジスティック回帰を使用しています。 しかし、私の問題では、負の出力(0)が正の出力(1)をはるかに上回ります。比率は20:1です。そのため、分類子をトレーニングすると、正の出力の可能性を強く示唆する機能でも、対応するパラメーターの値が非常に低い(非常に負の)ようです。これは、パラメーターをその方向に引っ張っている否定的な例が多すぎるために発生しているようです。 だから私は私が肯定的な例に重みを追加できるかどうか疑問に思っています(たとえば1ではなく20を使用して)。これはまったくメリットがありますか?もしそうなら、私はどのように重みを追加するべきですか(以下の方程式で)。 次のようなコスト関数ルックス J=(−1/m)⋅∑i=1my⋅log(h(x⋅θ))+(1−y)(1−log(h(x⋅θ)))J=(−1/m)⋅∑i=1my⋅log⁡(h(x⋅θ))+(1−y)(1−log⁡(h(x⋅θ)))J = (-1 / m) \cdot\sum_{i=1}^{m} y\cdot\log(h(x\cdot\theta)) + (1-y)(1 - \log(h(x\cdot\theta))) このコスト関数(WRTの勾配)次のとおりです。θθ\theta grad=((h(x⋅θ)−y)′⋅X)′grad=((h(x⋅θ)−y)′⋅X)′\mathrm{grad} = ((h(x\cdot\theta) - y)' \cdot X)' ここで、 =テストケースの数、x =特徴行列、y =出力ベクトル、h =シグモイド関数、θ =学習しようとしているパラメーター。mmmxxxyyyhhhθθ\theta 最後に、可能な限り低いを見つけるために勾配降下法を実行します。実装は正しく実行されているようです。JJJ

3
ガウス過程回帰おもちゃ問題
ガウスプロセス回帰の直感を得ようとしていたため、簡単な1Dおもちゃの問題を作成して試してみました。入力としてを取り、応答としてをました。(からの「インスピレーション」)バツ私= { 1 、2 、3 }xi={1,2,3}x_i=\{1,2,3\}、Y = X 2y私= { 1 、4 、9 }yi={1,4,9}y_i=\{1,4,9\}y= x2y=x2y=x^2 回帰には、標準的な二乗指数カーネル関数を使用しました。 k (xp、xq)= σ2fexp( − 12 リットル2| バツp− xq|2)k(xp,xq)=σf2exp⁡(−12l2|xp−xq|2)k(x_p,x_q)=\sigma_f^2 \exp \left( - \frac{1}{2l^2} \left|x_p-x_q\right|^2 \right) 標準偏差ノイズがあると仮定したので、共分散行列は次のようになります。σんσn\sigma_n Kp q= k (xp、xq)+ σ2んδp qKpq=k(xp,xq)+σn2δpqK_{pq} = k(x_p,x_q) + \sigma_n^2 \delta_{pq} ハイパー、データの対数尤度を最大化することによって推定しました。ポイントで予測を行うために、平均と分散をそれぞれ次のように見つけましたX ⋆(σん、l 、σf)(σn,l,σf)(\sigma_n,l,\sigma_f)バツ⋆x⋆x_\star σ 2 X ⋆ = …

1
調整済みR二乗は、固定スコアまたはランダムスコアの母集団をR二乗で推定しようとしますか?
母集団の二乗は、固定スコアまたはランダムスコアを想定して定義できます。ρ2ρ2\rho^2 固定スコア:予測子のサンプルサイズと特定の値は固定されます。したがって、は、予測子の値が一定に保たれているときに、母集団回帰式によって結果で説明される分散の割合です。ρ2fρf2\rho^2_f ランダムスコア:予測子の特定の値は、分布から抽出されます。したがって、は、予測子の値が予測子の母集団分布に対応する母集団の結果で説明される分散の割合を指します。ρ2rρr2\rho^2_r この区別が\ rho ^ 2の推定に大きな違いをもたらすかどうかρ2ρ2\rho^2については、以前尋ねました。また、\ rho ^ 2の公平な見積もりを計算する方法 ρ2ρ2\rho^2についても一般的に質問しました。 サンプルサイズが大きくなると、固定スコアとランダムスコアの区別が重要でなくなることがわかります。ただし、調整済みR2R2R^2が固定スコアまたはランダムスコア\ rho ^ 2を推定するように設計されているかどうかを確認しようとしていますρ2ρ2\rho^2。 ご質問 調整済みR2R2R^2 は、固定スコアまたはランダムスコア\ rho ^ 2を推定するように設計されていρ2ρ2\rho^2ますか? 調整済みr二乗の式が\ rho ^ 2の 1つまたは他の形式にどのように関係するかについての原理的な説明はありρ2ρ2\rho^2ますか? 私の混乱の背景 Yin and Fan(2001、p.206)を読んだとき、彼らはこう書いている: 重回帰モデルの基本的な前提の1つは、独立変数の値が既知の定数であり、実験前に研究者によって固定されることです。従属変数のみがサンプルごとに自由に変化します。その回帰モデルは、固定線形回帰モデルと呼ばれます。 ただし、社会科学および行動科学では、独立変数の値が研究者によって固定されることはほとんどなく、ランダムなエラーの影響も受けます。したがって、アプリケーションの2番目の回帰モデルが提案されており、従属変数と独立変数の両方を変化させることができます(Binder、1959; Park&Dudycha、1974)。そのモデルは、ランダムモデル(または修正モデル)と呼ばれます。ランダムモデルと固定モデルから取得した回帰係数の最尤推定値は正規性の仮定の下では同じですが、それらの分布は大きく異なります。ランダムモデルは非常に複雑であるため、一般的に使用される固定線形回帰モデルの代わりに、それを受け入れる前にさらに調査が必要です。したがって、通常は固定モデルが適用され、仮定が完全に満たされていなくても(Claudy、1978)。仮定に違反した固定回帰モデルをこのように適用すると、「オーバーフィッティング」が発生します。これは、不完全なサンプルデータから導入されたランダムエラーがプロセスで大文字にされる傾向があるためです。その結果、そのようにして得られたサンプルの多重相関係数は、真の母集団多重相関を過大評価する傾向があります(Claudy、1978; Cohen&Cohen、1983; Cummings、1982)。 したがって、上記のステートメントが、調整されたがランダムモデルによって導入されたエラーを補償することを言っているのか、またはこれがランダムモデルの存在を示すペーパーの単なる警告であったのかは不明でしたが、固定モデルに焦点を当てます。R2R2R^2 参考文献 Yin、P.&&Fan、X.(2001)。重回帰における収縮の推定:さまざまな分析方法の比較。Journal of Experimental Education、69(2)、203-224。PDFR2R2R^2

2
既知のブレークポイントを持つ区分的線形回帰における勾配の標準誤差
状況 1つの従属と1つの独立変数データセットがあります。発生する既知の/固定されたブレークポイントを使用して、連続的な区分線形回帰を近似したいと思います。ブレイクポインは不確実性なく知られているので、推定したくありません。次に、の形式の回帰(OLS)を これはX K (1、2、... 、K)Y I = β 0 + β 1 X I + β 2マックス(X I - 1、0 )+ β 3マックス(X I - 2、0 )+ ... + β K + 1つのマックス(Xyyyxxxkkk(a1,a2,…,ak)(a1,a2,…,ak)(a_{1}, a_{2}, \ldots, a_{k})yi=β0+β1xi+β2max(xi−a1,0)+β3max(xi−a2,0)+…+βk+1max(xi−ak,0)+ϵiyi=β0+β1xi+β2max⁡(xi−a1,0)+β3max⁡(xi−a2,0)+…+βk+1max⁡(xi−ak,0)+ϵi y_{i} = \beta_{0} + \beta_{1}x_{i} + \beta_{2}\operatorname{max}(x_{i}-a_{1},0) + \beta_{3}\operatorname{max}(x_{i}-a_{2},0) +\ldots+ \beta_{k+1}\operatorname{max}(x_{i}-a_{k},0) +\epsilon_{i} R …

1
通常の最小二乗法は効率的で公平な推定量をどのような仮定の下で与えますか?
ガウスマルコフの仮定の下で、通常の最小二乗法が効率的で不偏の推定量を与えることは本当ですか? そう: E(ut)= 0E(あなたt)=0E(u_t)=0 すべてのについてttt E(utあなたs)= σ2E(あなたtあなたs)=σ2E(u_tu_s)=\sigma^2 fort = st=st=s E(utあなたs)= 0E(あなたtあなたs)=0E(u_tu_s)=0 fort ≠ st≠st\neq s ここで、は残差です。あなたあなたu

3
Rのlm式の相互作用項を解釈するにはどうすればよいですか?
Rではlm()、次のように関数を呼び出すと、 lm.1 = lm(response ~ var1 + var2 + var1 * var2) summary(lm.1) これは私と応答変数の線形モデル与えvar1、var2およびそれらの間の相互作用を。しかし、相互作用の項はどの程度正確に数値で解釈されますか? ドキュメントでは、これはとの間の「クロス」であるvar1と述べていますが、var2正確には「クロス」が何であるかについての説明はありませんでした。 2つの変数間の相互作用を組み込むためにRが計算している正確な数値を知ることは私にとって役立ちます。
9 r  regression 

2
データの正規化に関連する混乱
私は線形回帰モデルを学ぼうとしています。ただし、データの正規化に関して混乱しています。特徴/予測子をゼロ平均および単位分散に正規化しました。ターゲットについても同じようにする必要がありますか?もしそうなら、なぜですか?

2
ガウス・マルコフの定理:BLUE and OLS
私はwikipediaの Guass-Markovの定理について読んでいます。誰かがこの定理の主要な点を理解するのを手伝ってくれることを願っていました。 :私たちは、マトリクス状に、によって与えられ、線形モデルを想定し と私たちはBLUEを探しています、β。y= Xβ+ ηy=Xβ+η y = X\beta +\eta βˆβ^ \widehat\beta この、私はラベルう "残留"と ε = β - β "エラー"。(つまり、Gauss-Markovページでの使用法の反対です)。η= y− Xβη=y−Xβ\eta = y - X\betaε = βˆ- βε=β^−β\varepsilon = \widehat\beta - \beta OLS(通常の最小二乗)推定量は、のargminとして導出できます。| 残差| | 2 2 = | | η | | 2 2。| | 残差 | |22= …

1
タイプIIIの二乗和
私は、1つのカテゴリー変数(男性と女性)と1つの連続変数を持つ線形回帰モデルを持っています。AAABBB Rのコントラストコードをで設定しましたoptions(contrasts=c("contr.sum","contr.poly"))。そして今、私は、タイプIIIの二乗和と、を使用したそれらの相互作用(A:B)があります。AAABBBdrop1(model, .~., test="F") 私が行き詰まっているのは、二乗和の計算方法です。だBBBと思いますsum((predicted y of the full model - predicted y of the reduced model)^2)。縮小モデルはのようになりますy~A+A:B。しかし、を使用するとpredict(y~A+A:B)、Rは完全なモデルの予測値と同じ予測値を返します。したがって、平方和は0になります。 (の二乗和には、の縮小モデルを使用しました。これはと同じです。)AAAy~B+A:By~A:B 以下は、ランダムに生成されたデータのコード例です。 A<-as.factor(rep(c("male","female"), each=5)) set.seed(1) B<-runif(10) set.seed(5) y<-runif(10) model<-lm(y~A+B+A:B) options(contrasts = c("contr.sum","contr.poly")) #type3 sums of squares drop1(model, .~., test="F") #or same result: library(car) Anova(lm(y~A+B+A:B),type="III") #full model predFull<-predict(model) #Calculate sum of squares #SS(A|B,AB) predA<-predict(lm(y~B+A:B)) …

1
ロジスティック回帰モデル変数のp値の意味
だから私はRでロジスティック回帰モデルを使っています。統計はまだ初めてですが、回帰モデルについて少し理解できたように思いますが、それでも気になることがいくつかあります。 リンクされた画像を見ると、私が作成したサンプルモデルのRプリントの概要が表示されています。このモデルは、データセット内の電子メールがrefoundか(バイナリ変数される場合は、予測しようとしているisRefound)とデータセットが密接に関連する2つの変数が含まれているisRefound、すなわち、next24およびnext7daysこれらはまた、バイナリであり、メールが次にクリックされる場合は教えてくれ- 24時間/ログの現在のポイントから次の7日間。 高いp値は、この変数がモデル予測に与える影響がかなりランダムであることを示しているはずですよね?これに基づいて、これらの2つの変数が計算式から外されている場合、モデル予測の精度が10%未満に低下する理由がわかりません。これらの変数の重要度が非常に低い場合、それらをモデルから削除すると大きな影響があるのはなぜですか? よろしくお願いします、リッキーフォックス 編集: 最初に私はnext24だけを削除しましたが、これはcoefがかなり小さいので影響が少ないはずです。予想通り、ほとんど変更されていません-そのための写真をアップロードしないでください。 next7daysを削除すると、モデルに大きな影響がありました:AIC 200kアップ、精度16%まで、再現率73%まで

3
平均パズルへの回帰
Daniel Kahnemanによる「Thinking、Fast and Slow」の「Regression to the Mean」の章では、例が示され、読者は、全体の売上予測と前年の販売数から、個々の店舗の売上を予測するよう求められています。 。たとえば、(この本の例には4つのストアがありますが、ここでは簡単にするために2つ使用しています): Store 2011 2012 1 100 ? 2 500 ? Total 600 660 単純な予測では、店舗1と2の場合はそれぞれ110%と550、10%増加します。しかし、著者はこの素朴なアプローチは間違っていると主張しています。パフォーマンスの低いストアは10%以上増加する可能性が高く、パフォーマンスの高いストアは10%未満増加または減少する可能性があります。したがって、115(15%増加)および535(7%増加)の予測は、単純な予測よりも「より正確」になります。 私が理解していないことは、ストア1の100の売上が必然的にパフォーマンスの低いストアであるとどのように結論できるかです。おそらく、場所の違いが原因で、店舗1と2の真の時系列平均は10と550であり、店舗1は2011年に非常に大きな年となり、店舗2は2011年に悲惨な年となりました。ストア1の減少とストア2の増加を予測するには? 元の例では時系列情報が記載されていなかったのはわかっていますが、「平均への回帰」は断面平均を指しているため、時系列情報は問題ではないという印象を受けています。私は何を誤解していますか?

1
ロジスティック回帰における巨大な係数-それは何を意味し、何をすべきか?
私はロジスティック回帰の間に巨大な係数を取得しkrajULKVます。 > summary(m5) Call: glm(formula = cbind(ml, ad) ~ rok + obdobi + kraj + resid_usili2 + rok:obdobi + rok:kraj + obdobi:kraj + kraj:resid_usili2 + rok:obdobi:kraj, family = "quasibinomial") Deviance Residuals: Min 1Q Median 3Q Max -2.7796 -1.0958 -0.3101 1.0034 2.8370 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -486.72087 …


2
縦断データによるSVM回帰
患者あたり約500の変数があり、各変数には1つの連続値があり、3つの異なる時点(2か月後と1年後)で測定されます。回帰により、新しい患者の治療結果を予測したいと思います。 そのような長期データでSVM回帰を使用することは可能ですか?

4
0と1の間の範囲の結果に対するロジスティック回帰の拡張
結果が厳密に0、1ではなく、含む0から1までのすべての実数の範囲にある回帰問題があります。Y= [ 0 、0.12 、0.31 、。。。、1 ]Y=[0,0.12,0.31,...,1]Y = [ 0, 0.12, 0.31, ..., 1 ] この問題はこのスレッドで既に説明されていますが、私の質問は少し異なります。 ロジスティック回帰が通常使用されるのと同じ理由で、線形回帰を使用できません。線形回帰では、A)非常に大きなIVs値は予測結果を1にスキューし、B)線形回帰の結果は0,1の制限に制限されません。 教科書からこのロジスティックコスト関数を見る方程式が計算するように設計されていることを収集しますとが同じ値0または1を持たない場合にのみ、0より大きいコスト。y xコスト= − yログ(h (x ))− (1 − y)ログ(1 − h (x ))Cost=−ylog⁡(h(x))−(1−y)log⁡(1−h(x))\text{Cost} = -y \log(h(x)) - (1 - y) \log(1-h(x))yyyバツxx コスト関数を変更してすべての仮説エラーを測定することにより、ロジスティック回帰を使用することは可能ですか?

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.