約 2,105,543 件
https://w.atwiki.jp/thecockrockshockpop/pages/1617.html
http //www.reflexmusic.jp/betawolf/ Beta Wolf EP Beta Wolf EP 2010年8月24日 ( HP ) 1. Undertow 2. Los Angeles 3. Say It To Me Now 4. Dig 5. Girl On Fire 6. Undertow [ acoustic ]
https://w.atwiki.jp/ce00582/pages/1758.html
Private Sub Command1_Click() Dim byear As Single Dim age As Single Dim car As Single Dim mis(15 To 64, 1 To 49) As Single Dim mos(15 To 64, 1 To 49) As Single Dim wis(15 To 64, 1 To 49) As Single Dim wos(15 To 64, 1 To 49) As Single Dim m2by(1900 To 2100, 15 To 69) As Single Dim f2by(1900 To 2100, 15 To 69) As Single Dim alpha(15 To 69) As Single Dim beta(15 To 69) As Single Dim gamma(15 To 69) As Single Dim theta(15 To 64) As Single Dim mw(15 To 69) As Single Dim fw(15 To 69) As Single Dim mde(1985 To 2100, 0 To 99) As Single Dim avgmw(1985 To 2100) As Single Dim z1 As Single Dim z2 As Single Dim z3 As Single Dim avg As Single Dim g(5 To 99) As Single Dim slide(5 To 99) As Single Dim year As Single Dim v(15 To 69) As Single Open "c /eli/data/男子脱退力.txt " For Input As #2 Do Until EOF(2) Input #2, a1, a2, a3, a4 age = a1 gamma(age) = a2 alpha(age) = a3 beta(age) = a4 Loop Close #2 Open "c /eli/data/再加入率.txt " For Input As #3 Do Until EOF(3) Input #3, a1, a2, a3, a4 age = a1 theta(age) = a2 Loop Close #3 Open "c /eli/gdata/変形厚生年金被保険者.txt " For Input As #5 Do Until EOF(5) Input #5, a1, a2, a3, a4 byear = a1 age = a2 m2by(byear, age) = a3 f2by(byear, age) = a4 Loop Close #5 age = 15 Open "c /eli/data/標準報酬.txt " For Input As #12 Do Until EOF(12) Input #12, a1, a2 mw(age) = a1 fw(age) = a2 age = age + 1 Loop Close #12 Open "c /eli/gdata/比例スライド.txt " For Input As #42 Do Until EOF(42) Input #42, a1, a2 year = a1 slide(year) = a2 Loop Close #42 g1 = 0.025 For byear = 1990 To 2034 age = 15 mis(age, 1) = m2by(byear, age) For car = 1 To 49 year = byear + age - 2000 wis(age, car) = mw(age) wos(age, car) = mw(age) Next For age = 16 To 64 z1 = (1 - gamma(age - 1)) * m2by(byear, age - 1) z2 = m2by(byear, age) - z1 If z2 0 Then z2 = 0 z3 = (1 - theta(age)) * z2 z4 = theta(age) * z2 mis(age, 1) = z3 z5 = 0 For car = 1 To 49 z5 = z5 + mos(age - 1, car) Next If z5 = 0 Then z5 = 1 v(age - 1) = z4 / z5 If v(age - 1) 1 Then v(age - 1) = 1 For car = 2 To 49 mis(age, car) = (1 - gamma(age - 1)) * mis(age - 1, car - 1) + v(age - 1) * mos(age - 1, car - 1) Next For car = 1 To 49 mos(age, car) = (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) Next year = byear + age - 2000 wis(age, 1) = mw(age) For car = 2 To 49 z6 = (1 - gamma(age - 1)) * mis(age - 1, car - 1) + v(age - 1) * mos(age - 1, car - 1) z7 = wis(age - 1, car - 1) * (1 - gamma(age - 1)) * mis(age - 1, car - 1) + wos(age - 1, car - 1) * v(age - 1) * mos(age - 1, car - 1) If z6 = 0 Then z7 = 0 If z6 = 0 Then z6 = 1 z8 = z7 / z6 year = byear + age - 2000 wis(age, car) = (mw(age) + (1 + g1 - slide(year)) * (car - 1) * z8 / (1 + g1)) / car Next For car = 2 To 49 z9 = (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) z10 = wis(age - 1, car) * (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + wos(age - 1, car) * (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) If z9 = 0 Then z10 = 0 If z9 = 0 Then z9 = 1 year = byear + age - 2000 wos(age, car) = (1 + g1 - slide(year)) * z10 / (z9 * (1 + g1)) Next Next z1 = 0 z2 = 0 For car = 25 To 49 z1 = z1 + mis(64, car) + mos(64, car) z2 = z2 + wis(64, car) * mis(64, car) + wos(64, car) * mos(64, car) Next year = byear + 65 - 2000 avgmw(byear) = (1 + g1 - slide(year)) * z2 / ((1 + g1) * z1) Next Open "c /eli/gdata/男子賃金1.txt " For Output As #4 For byear = 1990 To 2100 Write #4, byear, avgmw(byear) Next Close #4 End Sub
https://w.atwiki.jp/ce00582/pages/1903.html
Private Sub Command1_Click() Dim byear As Single Dim age As Single Dim car As Single Dim mis(15 To 64, 1 To 49) As Single Dim mos(15 To 64, 1 To 49) As Single Dim wis(15 To 64, 1 To 49) As Single Dim wos(15 To 64, 1 To 49) As Single Dim m2by(1900 To 2100, 15 To 69) As Single Dim f2by(1900 To 2100, 15 To 69) As Single Dim alpha(15 To 69) As Single Dim beta(15 To 69) As Single Dim gamma(15 To 69) As Single Dim theta(15 To 64) As Single Dim mw(15 To 69) As Single Dim fw(15 To 69) As Single Dim mde(1985 To 2100, 0 To 99) As Single Dim avgmw(1985 To 2100) As Single Dim z1 As Single Dim z2 As Single Dim z3 As Single Dim v(15 To 69) As Single Dim avg As Single Dim g(5 To 99) As Single Dim year As Single Dim slide(5 To 99) As Single Open "c /dig/gdata/比例スライド.txt " For Input As #42 Do Until EOF(42) Input #42, a1, a2 year = a1 slide(year) = a2 Loop Close #42 Open "c /dig/data/女子脱退力.txt " For Input As #2 Do Until EOF(2) Input #2, a1, a2, a3, a4 age = a1 gamma(age) = a2 alpha(age) = a3 beta(age) = a4 Loop Close #2 Open "c /dig/data/再加入率.txt " For Input As #3 Do Until EOF(3) Input #3, a1, a2, a3, a4 age = a1 theta(age) = a3 Loop Close #3 Open "c /dig/gdata/変形厚生年金被保険者.txt " For Input As #5 Do Until EOF(5) Input #5, a1, a2, a3, a4 byear = a1 age = a2 m2by(byear, age) = a3 f2by(byear, age) = a4 Loop Close #5 age = 15 Open "c /dig/data/標準報酬.txt " For Input As #12 Do Until EOF(12) Input #12, a1, a2 mw(age) = a1 fw(age) = a2 age = age + 1 Loop Close #12 g1 = 0.025 For byear = 1990 To 2034 age = 15 mis(age, 1) = f2by(byear, age) For car = 1 To 49 year = byear + age - 2000 wis(age, car) = fw(age) wos(age, car) = fw(age) Next For age = 16 To 64 z1 = (1 - gamma(age - 1)) * f2by(byear, age - 1) z2 = f2by(byear, age) - z1 If z2 0 Then z2 = 0 z3 = (1 - theta(age)) * z2 z4 = theta(age) * z2 mis(age, 1) = z3 z5 = 0 For car = 1 To 49 z5 = z5 + mos(age - 1, car) Next If z5 = 0 Then z5 = 1 v(age - 1) = z4 / z5 If v(age - 1) 1 Then v(age - 1) = 1 For car = 2 To 49 mis(age, car) = (1 - gamma(age - 1)) * mis(age - 1, car - 1) + v(age - 1) * mos(age - 1, car - 1) Next For car = 1 To 49 mos(age, car) = (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) Next year = byear + age - 2000 wis(age, 1) = fw(age) For car = 2 To 49 z6 = (1 - gamma(age - 1)) * mis(age - 1, car - 1) + v(age - 1) * mos(age - 1, car - 1) z7 = wis(age - 1, car - 1) * (1 - gamma(age - 1)) * mis(age - 1, car - 1) + wos(age - 1, car - 1) * v(age - 1) * mos(age - 1, car - 1) If z6 = 0 Then z7 = 0 If z6 = 0 Then z6 = 1 z8 = z7 / z6 year = byear + age - 2000 wis(age, car) = (fw(age) + (1 + g1 - slide(year)) * (car - 1) * z8 / (1 + g1)) / car Next For car = 2 To 49 z9 = (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) z10 = wis(age - 1, car) * (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + wos(age - 1, car) * (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) If z9 = 0 Then z10 = 0 If z9 = 0 Then z9 = 1 year = byear + age - 2000 wos(age, car) = (1 + g1 - slide(year)) * z10 / (z9 * (1 + g1)) Next Next z1 = 0 z2 = 0 For car = 25 To 49 z1 = z1 + mis(64, car) + mos(64, car) z2 = z2 + wis(64, car) * mis(64, car) + wos(64, car) * mos(64, car) Next year = byear + 65 - 2000 avgmw(byear) = (1 + g1 - slide(year)) * z2 / ((1 + g1) * z1) Next Open "c /dig/gdata/女子賃金1.txt " For Output As #4 For byear = 1990 To 2100 Write #4, byear, avgmw(byear) Next Close #4 End Sub
https://w.atwiki.jp/javadsge/pages/301.html
package h2; public class pro { int[] y=new int[200]; int[] p=new int[200]; double[][][] v=new double[200][21][21]; int[][][] goton=new int[200][21][21]; int[][][] gotom=new int[200][21][21]; int s; double phi; int[] opn=new int[200]; int[] opm=new int[200]; public static void main(String[] args) { pro test=new pro(); } pro(){ phi=0.5; y[0]=100; p[0]=100; for(s=1;s 102;s++){ y[s]=100; p[s]=100; } double theta,lam,beta; theta=0.5; lam=0.5; beta=0.5; int n,m; int gap,yy; int gg; gg=0; int t=1; for(n=1;n 21;n++){ for(m=1;m 21;m++){ v[1][m][n]=0; }} int nx,mx,nxs,mxs; int g; for(t=2;t 101;t++){ for(n=1;n 21;n++){ for(m=1;m 21;m++){ gap=n-10; g=m-10; yy=y[t]-gap; double maxv=9999; nxs=0; mxs=0; for(nx=1;nx 21;nx++){ for(mx=1;mx 21;mx++){ int gap2=nx-10; int yy2=y[t-1]-gap2; double e1=gap-phi*gap2; double e2=yy-yy2-g; int g2=mx-10; double e3=g-theta*gg-(1-theta)*g2; double e4=p[t]-(1-lam)*p[t+1]-lam*p[t-1]-beta*gap; double v1=e1*e1+e2*2+e3*e3+e4*e4+v[t-1][mx][nx]; if(v1 maxv)nxs=nx; if(v1 maxv)mxs=mx; if(v1 maxv)maxv=v1; }} goton[t][m][n]=nxs; gotom[t][m][n]=mxs; v[t][m][n]=maxv; }}} double maxv=9999; mx=0; nx=0; for(n=1;n 21;n++){ for(m=1;m 21;m++){ if(v[100][m][n] maxv)mx=m; if(v[100][m][n] maxv)nx=n; if(v[100][m][n] maxv)maxv=v[100][m][n]; }} opn[100]=nx; opm[100]=mx; for(s=1;s 100;s++){ t=101-s; m=opm[t]; n=opn[t]; opm[t-1]=gotom[t][m][n]; opn[t-1]=goton[t][m][n]; } for(s=1;s 101;s++){ System.out.println(y[s]+","+(opn[s]-10)); } } }
https://w.atwiki.jp/lotroquest/pages/204.html
RETURN THE CASK/ 感謝の一樽 ブレガドール(Bregedur)との会話 「ありがとう[名前]。そして私の思いを汲んでくれた友ブレシルウェン(BRETHILWEN)にも感謝しよう。私は ブレシルウェン(BRETHILWEN)が汚らわしいゴブリンどもに悩まされていることを知っていたのだよ。そして あなたはこの素晴らしきワインと実りのために勇敢に戦ってくれたのだな。それは我々にとって意味ある大切なことなのだ。」 ●GIFT FOR A FRIEND /友への贈り物につづく
https://w.atwiki.jp/hmiku/pages/1896.html
【検索用 ONIYAGURAtheMetalおにやくら 登録タグ 2008年 O VOCALOID シンP 曲 曲英 神威がくぽ】 + 目次 目次 曲紹介 代表的なPV紹介 歌詞 コメント 作詞:シンP 作曲:シンP 編曲:シンP 唄:神威がくぽ 曲紹介 曲名:『ONIYAGURA the Metal -鬼櫓-』(ONIYAGURA the Metal -おにやぐら-) がくぽの公式絵を見た瞬間「これはやるしかない」と思いました。(作者コメ転載) 代表的なPV紹介 歌詞 山颪荒ぶる 雨月の夜 数多の高楼を仰ぎ 杯を独り掲げ踊るは 幽世と顕世の神楽 沖津風遠く呼び 灯す棕櫚の篝火よ 誰そ彼に 魍魎率いて 唐草を踏みしめ 葛茂る頂望む 前に聳えるは 万丈の山 後に支うは 千尋の谷 白河に臨みて催すは 終わることなき曲水の宴 車輪のごとき夢幻の回廊を輪廻す 笹の舟よ どこまでも狂々と ここは黄泉道繋ぐ 鬼櫓 我は踊らん 鬼神楽 花吹雪舞い散る 荘園に 祀るは八百万の御霊 蜻蛉を独り 集め遊ぶは 星屑と望月の蹴鞠 岩清水かき寄せて 掬い映す姿には 何者も 囁くことなし 天竺を拝賀す 屍者の嘆きは闇に溶けて わっしょい わっしょい わっしょい わっしょい わっしょい わっしょい わっしょい わっしょい そいや そいや そいや そいや そいや そいや そいや そいや らっせー らっせー らっせー らっせー らっせー らっせー らっせー らっせー そーりゃ そーりゃ そーりゃ そーりゃ そーりゃ そーりゃ そーりゃ そーりゃ 黒曜の剣で 穿つとも 始まることなき 月次祭 神威の雙六を たどる霊魂に 六文銭の賽よ いつまでも空頃と ここは追儺きたる 鬼神楽 我は鬼の子 鬼櫓 コメント 格好良すぎて惚れた -- 名無しさん (2009-05-01 20 30 37) がくっぽくどはいいよ! -- 名無しさん (2009-09-07 18 17 26) シンPの曲はうろたんとジャスティスに埋もれてる。名作なのに世に出てないの多いんだよね・・・ -- 名無しさん (2009-09-17 20 30 53) 神降臨…マジヤバス -- 名無しさん (2010-08-08 13 56 52) シンP神!! -- たんたん (2010-10-30 13 57 25) 醒とか破とか也とか格好いい! -- 逆上 瞳 (2014-10-24 22 08 02) 名前 コメント
https://w.atwiki.jp/swtfu2/pages/111.html
THE BATTLE OF ENDOR A long time ago in a galaxy far, far away.... 遠い昔、はるか彼方の銀河系で・・・ STAR WARS スター・ウォーズ THE FORCE UNLEASHED II THE BATTLE OF ENDOR フォース・アンリーシュド2 エンドアの戦い The Rebellion is on the brink of destruction. Mounting a desperate attack on the second DEATH STAR, the Rebels hope to restore freedom to the Galaxy. 反乱同盟軍は滅亡の瀬戸際にあった。第2デス・スターへの絶望的な攻撃を仕掛けることで、反乱軍は銀河系に自由を取り戻そうとしていた。 While the battle rages in space, a small band of Rebels land on the Forest Moon of ENDOR to assault the shield generator protecting the battle station. 宇宙空間で激しい戦いが繰り広げられている間に、反乱軍の小隊はバトル・ステーションを保護しているシールド・ジェネレーターを攻撃する為、森林の月エンドアに着陸した。 Little do the Rebels know, Darth Vader has dispatched his personal assassin, a dark clone of Starkiller, to finish the Rebellion once and for all.... 反乱軍のほとんどの者は知らなかったが、これを最後に反乱同盟軍の息の根を止める為、ダース・ヴェイダーは彼の個人的な暗殺者、スターキラーの邪悪なクローンを派遣していた.... Darth Vader The Rebels must be stopped before they can destroy the shield generator. Leave no survivors. 反逆者どもがシールド・ジェネレーターを破壊する前に、彼らを止めるのだ。生存者を残すな。 Dark Apprentice Understood, my Master. 了解しました、マスター。 Darth Vader There will be no one to stop us this time. 今や我々を止められる者は誰もいない。 Dark Apprentice Yes, my Master. はい、マスター。 Rebel Soldier Everyone take positions, we have incoming! 全員配置につくんだ、来たぞ! Rebel Soldier Don’t let him through! 絶対に奴を通すな! Rebel Soldier Activate those turrets. 砲塔を起動せよ。 Rebel Soldier He’s over here! 彼はこちらにいます! Han Solo Chewie, we got company. チューイー、お客さんが来たようだ! Han Solo Stay down, Chewie. I’ve got you covered. そこにいろ、チューイー。俺がお前を援護する。 Han Solo Can’t you figure out how to control that thing!? It’s getting a little hot down here! 何とかしてそれをコントロール出来ないのか!?少し暑くなってきてるようだぞ! Han Solo Well think of something! えっと、何か考えろ! Han Solo We have to buy the princess more time! 俺達はプリンセスの為にもっと時間を稼いでやる必要がある! Han Solo CHEWIE!!! チューイー!!! Rebel Soldier Princess, you go on! We’ll cover you. You have to set those charges! プリンセス、進んで下さい!我々が援護します。あなたは爆薬を設置する必要があります! Princess Leia Understood... and thank you. 了解しました…ありがとう。 Rebel Soldier Warn the others! 他の者に警告しろ! Rebel Soldier Don’t let him through! 絶対に奴を通すな! Dark Apprentice I don’t know why you even bother... 貴様が何故悩んでいるのか分からないが… Princess Leia There is much you fail to see, Sith. あなたでは見られないものが沢山あります、シス。 Princess Leia After my brother fell on Hoth, I fulfilled the destiny that he could not. 私の兄がホスで倒れた後、私は彼が出来なかった運命を成し遂げました。 Princess Leia Now I am a Jedi. 今や私がジェダイです。 Dark Apprentice You show promise, but you have much to learn. 貴様は見込みがあるが、まだ学ぶべき事も沢山あるぞ。 Dark Apprentice It s over, you rebel scum. これで終わりだ、反逆者のクズめ。 The emperor A great day, Lord Vader. Today marks the death of the Rebellion, and the birth of a new era of peace for the Empire. There is but one last detail to take care of. 記念すべき日であるな、ヴェイダー卿。この日、反乱軍の死と帝国の新しき平和の時代の誕生を祝うのだ。だが、最後に一つだけ片付けねばならぬことがあるぞ。 The emperor Do you take me for a fool? Your feeble schemes to overthrow me using that pathetic clone of your former apprentice were pitifully obvious from their inception. 余を愚か者だと思うておるのか?かつての弟子のみすぼらしいクローンを使って余を倒すというお前の貧相な計画など、初めから分かりきっておったわ。 The emperor Now Skywalker… you will die. スカイウォーカーよ…お前は死ぬのだ。 The emperor Ha ha ha! ハー、ハッ、ハッ! The emperor Captain Sarcli... サークリィ大尉… Captain Sarcli My Liege. はい陛下。 The emperor Take your troops and destroy the clone. 隊を送りクローンを抹殺せよ。
https://w.atwiki.jp/tryo/pages/18.html
1362%5B1%5D.gif 重回帰モデル回帰係数の推定量(最尤推定量など)を求めて,最小値であることを示すときに困りました>< 対数尤度関数からについて2次偏導関数をもとめ, は 次正方行列であり,逆行列をもつ ⇔ ,各固有値の積 なので,固有値はすべて 0でない。 さらに,各固有値の和 なので,・・・ あぁぁぁ!!やっぱり固有値はすべて正とは言えない><できなーーい↓↓ 名前 コメント TeX 微分符号は⇒\delta ではなく,⇒\partial である は\theta \hatではなく,\hat {\theta }である 「(ページ数)/(総ページ数)ページ」は,\usepackage{lastpage} とした上で\rfoot{\thepage{} ~/ \pageref{LastPage}{} ~ページ}である 場合分け:\begin{eqnarray} %改行 f(x)=\left\{ \begin{array}{ll} %改行 2x (0 x 1) \\ %改行 0 (その他) \\ %改行 \end{array} \right. %改行 \end{eqnarray} は \infty である は \xrightarrow[f_t]{ } である。←うまく表示できないが矢印の下にを書かせてます>< 統計 R 統計知識 (覚えておこう~!!^口^) 27_4200_92.gif 標準誤差 標準偏差(><) 標準誤差【standard error, 】: 標準偏差をとしたとき,を標準誤差という。はデータの個数である。 論文では,のように使われることが多い。 推定値の標準偏差のことである。 標準偏差【standard deviation, 】: √分散 である データに対して,を,標本標準偏差という。 分散の単位は平均値の単位の2乗になるので,データのばらつきを見るには標準偏差のほうが便利である。 信頼空間【confidence interval】 確率95%や99%で母集団のパラメータを推定するとき,その区間を信頼空間という。 そのときの確率を信頼係数(信頼水準?!)という。 信頼空間の両端をそれぞれ信頼限界という。
https://w.atwiki.jp/ce00582/pages/1902.html
Private Sub Command1_Click() Dim byear As Single Dim age As Single Dim car As Single Dim mis(15 To 64, 1 To 49) As Single Dim mos(15 To 64, 1 To 49) As Single Dim wis(15 To 64, 1 To 49) As Single Dim wos(15 To 64, 1 To 49) As Single Dim m2by(1900 To 2100, 15 To 69) As Single Dim f2by(1900 To 2100, 15 To 69) As Single Dim alpha(15 To 69) As Single Dim beta(15 To 69) As Single Dim gamma(15 To 69) As Single Dim theta(15 To 64) As Single Dim mw(15 To 69) As Single Dim fw(15 To 69) As Single Dim mde(1985 To 2100, 0 To 99) As Single Dim avgmw(1985 To 2100) As Single Dim z1 As Single Dim z2 As Single Dim z3 As Single Dim avg As Single Dim g(5 To 99) As Single Dim slide(5 To 99) As Single Dim year As Single Dim v(15 To 69) As Single Open "c /dig/data/男子脱退力.txt " For Input As #2 Do Until EOF(2) Input #2, a1, a2, a3, a4 age = a1 gamma(age) = a2 alpha(age) = a3 beta(age) = a4 Loop Close #2 Open "c /dig/data/再加入率.txt " For Input As #3 Do Until EOF(3) Input #3, a1, a2, a3, a4 age = a1 theta(age) = a2 Loop Close #3 Open "c /dig/gdata/変形厚生年金被保険者.txt " For Input As #5 Do Until EOF(5) Input #5, a1, a2, a3, a4 byear = a1 age = a2 m2by(byear, age) = a3 f2by(byear, age) = a4 Loop Close #5 age = 15 Open "c /dig/data/標準報酬.txt " For Input As #12 Do Until EOF(12) Input #12, a1, a2 mw(age) = a1 fw(age) = a2 age = age + 1 Loop Close #12 Open "c /dig/gdata/比例スライド.txt " For Input As #42 Do Until EOF(42) Input #42, a1, a2 year = a1 slide(year) = a2 Loop Close #42 g1 = 0.025 For byear = 1990 To 2034 age = 15 mis(age, 1) = m2by(byear, age) For car = 1 To 49 year = byear + age - 2000 wis(age, car) = mw(age) wos(age, car) = mw(age) Next For age = 16 To 64 z1 = (1 - gamma(age - 1)) * m2by(byear, age - 1) z2 = m2by(byear, age) - z1 If z2 0 Then z2 = 0 z3 = (1 - theta(age)) * z2 z4 = theta(age) * z2 mis(age, 1) = z3 z5 = 0 For car = 1 To 49 z5 = z5 + mos(age - 1, car) Next If z5 = 0 Then z5 = 1 v(age - 1) = z4 / z5 If v(age - 1) 1 Then v(age - 1) = 1 For car = 2 To 49 mis(age, car) = (1 - gamma(age - 1)) * mis(age - 1, car - 1) + v(age - 1) * mos(age - 1, car - 1) Next For car = 1 To 49 mos(age, car) = (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) Next year = byear + age - 2000 wis(age, 1) = mw(age) For car = 2 To 49 z6 = (1 - gamma(age - 1)) * mis(age - 1, car - 1) + v(age - 1) * mos(age - 1, car - 1) z7 = wis(age - 1, car - 1) * (1 - gamma(age - 1)) * mis(age - 1, car - 1) + wos(age - 1, car - 1) * v(age - 1) * mos(age - 1, car - 1) If z6 = 0 Then z7 = 0 If z6 = 0 Then z6 = 1 z8 = z7 / z6 year = byear + age - 2000 wis(age, car) = (mw(age) + (1 + g1 - slide(year)) * (car - 1) * z8 / (1 + g1)) / car Next For car = 2 To 49 z9 = (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) z10 = wis(age - 1, car) * (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + wos(age - 1, car) * (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) If z9 = 0 Then z10 = 0 If z9 = 0 Then z9 = 1 year = byear + age - 2000 wos(age, car) = (1 + g1 - slide(year)) * z10 / (z9 * (1 + g1)) Next Next z1 = 0 z2 = 0 For car = 25 To 49 z1 = z1 + mis(64, car) + mos(64, car) z2 = z2 + wis(64, car) * mis(64, car) + wos(64, car) * mos(64, car) Next year = byear + 65 - 2000 avgmw(byear) = (1 + g1 - slide(year)) * z2 / ((1 + g1) * z1) Next Open "c /dig/gdata/男子賃金1.txt " For Output As #4 For byear = 1990 To 2100 Write #4, byear, avgmw(byear) Next Close #4 End Sub
https://w.atwiki.jp/ce00582/pages/1912.html
Private Sub Command1_Click() Dim byear As Single Dim age As Single Dim car As Single Dim mis(15 To 64, 1 To 49) As Single Dim mos(15 To 64, 1 To 49) As Single Dim wis(15 To 64, 1 To 49) As Single Dim wos(15 To 64, 1 To 49) As Single Dim m2by(1900 To 2100, 15 To 69) As Single Dim f2by(1900 To 2100, 15 To 69) As Single Dim alpha(15 To 69) As Single Dim beta(15 To 69) As Single Dim gamma(15 To 69) As Single Dim theta(15 To 64) As Single Dim mw(15 To 69) As Single Dim fw(15 To 69) As Single Dim mde(1985 To 2100, 0 To 99) As Single Dim avgtmw(1985 To 2100) As Single Dim z1 As Single Dim z2 As Single Dim z3 As Single Dim v(15 To 69) As Single Dim year As Single Dim slide(5 To 99) As Single Open "c /dig/data/女子脱退力.txt " For Input As #2 Do Until EOF(2) Input #2, a1, a2, a3, a4 age = a1 gamma(age) = a2 alpha(age) = a3 beta(age) = a4 Loop Close #2 Open "c /dig/data/再加入率.txt " For Input As #3 Do Until EOF(3) Input #3, a1, a2, a3, a4 age = a1 theta(age) = a3 Loop Close #3 Open "c /dig/gdata/変形厚生年金被保険者.txt " For Input As #5 Do Until EOF(5) Input #5, a1, a2, a3, a4 byear = a1 age = a2 m2by(byear, age) = a3 f2by(byear, age) = a4 Loop Close #5 age = 15 Open "c /dig/data/標準報酬.txt " For Input As #12 Do Until EOF(12) Input #12, a1, a2 mw(age) = a1 fw(age) = a2 age = age + 1 Loop Close #12 Open "c /dig/gdata/比例スライド.txt " For Input As #42 Do Until EOF(42) Input #42, a1, a2 year = a1 slide(year) = a2 Loop Close #42 g1 = 0.025 For byear = 1990 To 2034 age = 15 mis(age, 1) = f2by(byear, age) For car = 1 To 49 wis(age, car) = fw(age) wos(age, car) = fw(age) Next For age = 16 To 64 z1 = (1 - gamma(age - 1)) * f2by(byear, age - 1) z2 = f2by(byear, age) - z1 If z2 0 Then z2 = 0 z3 = (1 - theta(age)) * z2 z4 = theta(age) * z2 mis(age, 1) = z3 z5 = 0 For car = 1 To 49 z5 = z5 + mos(age - 1, car) Next If z5 = 0 Then z5 = 1 v(age - 1) = z4 / z5 If v(age - 1) 1 Then v(age - 1) = 1 For car = 2 To 49 mis(age, car) = (1 - gamma(age - 1)) * mis(age - 1, car - 1) + v(age - 1) * mos(age - 1, car - 1) Next For car = 1 To 49 mos(age, car) = (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) Next wis(age, 1) = fw(age) For car = 2 To 49 z6 = (1 - gamma(age - 1)) * mis(age - 1, car - 1) + v(age - 1) * mos(age - 1, car - 1) z7 = wis(age - 1, car - 1) * (1 - gamma(age - 1)) * mis(age - 1, car - 1) + wos(age - 1, car - 1) * v(age - 1) * mos(age - 1, car - 1) If z6 = 0 Then z7 = 0 If z6 = 0 Then z6 = 1 z8 = z7 / z6 year = byear + age - 2000 wis(age, car) = (fw(age) + (1 + g1 - slide(year)) * (car - 1) * z8 / (1 + g1)) / car Next For car = 2 To 49 z9 = (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) z10 = wis(age - 1, car) * (gamma(age - 1) - alpha(age - 1) - beta(age)) * mis(age - 1, car) + wos(age - 1, car) * (1 - alpha(age - 1) - v(age - 1)) * mos(age - 1, car) If z9 = 0 Then z10 = 0 If z9 = 0 Then z9 = 1 year = byear + age - 2000 wos(age, car) = (1 + g1 - slide(year)) * z10 / ((1 + g1) * z9) Next Next z1 = 0 z2 = 0 For car = 1 To 24 z1 = z1 + mis(64, car) + mos(64, car) z2 = z2 + wis(64, car) * mis(64, car) + wos(64, car) * mos(64, car) Next year = byear + 65 - 2000 avgtmw(byear) = (1 + g1 - slide(year)) * z2 / (z1 * (1 + g1)) Debug.Print byear, avgtmw(byear) Next Open "c /dig/gdata/女子通算賃金1.txt " For Output As #4 For byear = 1990 To 2100 Write #4, byear, avgtmw(byear) Next Close #4 End Sub