Rスクリプト：SSM職業分類を基にした職業威信スコア

やっぱりweb上にはないようなので、自作しました。
とは言え、ほんとうは、寺沢の自作と言える部分はごくわずかです。
立教大・村瀬洋一先生のサイトで公開されているSPSS用シンタクスを流用させて頂きました。
ここにあるのは95年版を元にしたものなので、ご注意を。詳しくは以下の村瀬先生のページをご覧下さい。
http://www2.rikkyo.ac.jp/web/murase/ssm/prscore.html

あと、800番台の職業に対応していない点もご注意下さい（念のため）

## 職業威信スコア
## 立教大学・村瀬先生作成SPSSシンタクス【95年版】をR用に改変）

  library(car)   # recode() 関数使用のため

# 例：「Dataset」というデータセットに含まれている
# 「syokugyou」という職業の変数をもとに
# 「job.isin」という連続変数を新たに作る


Dataset $ job.isin <- recode(Dataset$syokugyou,  "

  c(501,502)=72.0;
  c(503)=66.3;
  c(504,505)=72.0;
  c(506,507)=66.3;
  c(508,509)=90.1;
  c(510)=65.7;
  c(511,512,513,514,515,516)=59.7;
  c(517)=86.9;
  c(518,519)=70.8;
  c(520)=58.3;
  c(521,512,523)=63.6;
  c(524)=84.3;
  c(525)=63.6;
  c(526)=58.3;
  c(527)=60.3;
  c(528)=66.6;
  c(529)=52.2;
  c(530)=66.6;
  c(531,532)=64.6;
  c(533)=66.6;
  c(534,535)=69.0;
  c(536)=65.7;
  c(537)=52.9;
  c(538)=59.7;
  c(539)=58.3;
  c(540,541)=70.8;
  c(542)=70.0;
  c(543,544)=63.6;
  c(545)=67.2;
  c(546)=74.9;
  c(547)=67.2;
  c(548)=78.1;
  c(549)=67.2;
  c(550)=59.7;
  c(551,552)=56.9;
  c(553)=56.1;
  c(554)=52.2;
  c(555,556)=48.2;
  c(557,558)=52.2;
  c(559)=52.9;
  c(560)=47.8;
  c(561,562)=42.4;
  c(563)=47.8;
  c(564,565)=48.2;
  c(566)=51.3;
  c(567)=52.9;
  c(568)=51.3;
  c(569,570,571)=42.4;
  c(572,573)=47.2;
  c(574)=44.3;
  c(575)=47.2;
  c(576,577)=42.4;
  c(578)=38.1;
  c(579)=49.7;
  c(580)=42.0;
  c(581)=51.6;
  c(582,583)=38.1;
  c(584)=70.0;
  c(585,586)=38.1;
  c(587)=48.9;
  c(588)=38.1;
  c(589)=51.3;
  c(590)=39.9;
  c(591)=64.6;
  c(592)=38.1;
  c(593,594,595)=57.9;
  c(596,597)=39.9;
  c(598)=57.9;
  c(599,600,601,602,603)=45.6;
  c(604,605)=46.5;
  c(606)=51.3;
  c(607)=48.9;
  c(608,609,610)=82.5;
  c(611,612,613)=47.8;
  c(614)=39.0;
  c(615,616)=51.3;
  c(617,618)=46.2;
  c(619)=51.3;
  c(620,621,622)=36.7;
  c(623,624)=52.1;
  c(625,626)=44.0;
  c(627,628,629,630,631,632,633,634,635,636,637,638,639,640,641,642)=48.9;
  c(643)=42.2;
  c(644)=44.6;
  c(645,646)=42.2;
  c(647)=44.6;
  c(648,649,650)=42.0;
  c(651)=52.1;
  c(652,653)=42.0;
  c(654)=52.1;
  c(655)=53.1;
  c(656)=52.1;
  c(657)=42.0;
  c(658)=44.0;
  c(659,660)=48.9;
  c(661)=44.0;
  c(662,663,664,665,666)=52.1;
  c(667,668,669)=44.0;
  c(670,671)=52.1;
  c(672)=42.0;
  c(673,674,675)=46.8;
  c(676,677)=50.4;
  c(678,679)=53.1;
  c(680)=46.1;
  c(681)=53.1;
  c(682,683)=39.0;
  c(684)=47.9;
  c(685,688)=39.0;
  c(701)=52.9;
  c(702)=53.1;
  c(703)=63.6;
  c(704)=42.0;
  else=NA           ")