2019-11-14

ＰＣＡ（主成分分析を勉強してみた）

１．データーの準備

生徒NO <- seq(1, 10, 1)
kubun1 <- c(1,1,1,2,2,2,2,3,3,3)
kubun2 <- c("a","a","a","b","b","b","b","c","c","c")
国語 <- c(96,91,89,78,67,63,53,38,31,39)
英語 <- c(96,91,89,78,67,53,63,36,24,32)
数学 <- c(67,68,69,58,61,55,49,82,96,93)

d01 <- data.frame(生徒NO,kubun1,kubun2,国語,英語,数学)
d01
sapply(d01, class)

(出力)

> d01
   生徒NO kubun1 kubun2 国語 英語 数学
1       1      1      a   96   96   67
2       2      1      a   91   91   68
3       3      1      a   89   89   69
4       4      2      b   78   78   58
5       5      2      b   67   67   61
6       6      2      b   63   53   55
7       7      2      b   53   63   49
8       8      3      c   38   36   82
9       9      3      c   31   24   96
10     10      3      c   39   32   93
> sapply(d01, class)
   生徒NO    kubun1    kubun2      国語      英語      数学 
"numeric" "numeric"  "factor" "numeric" "numeric" "numeric"

２．pcaを実行してみる

> pca <- prcomp(d01[,c(4,5,6)], scale = T)
> print(pca)
Standard deviations (1, .., p=3):
[1] 1.5658770 0.7287046 0.1304565

Rotation (n x k) = (3 x 3):
            PC1       PC2         PC3
国語 -0.6071171 0.4078050  0.68198526
英語 -0.6210608 0.2918497 -0.72739756
数学  0.4956735 0.8651698 -0.07608478
> biplot(pca)

f:id:buriburipipi2:20191114055912j:plain
３．おしゃれな表示方法

library("FactoMineR")
library("factoextra")
fviz_pca_biplot(pca, repel = T)

f:id:buriburipipi2:20191114062107j:plain

fviz_pca_ind(pca,   habillage=d01$kubun2,
     addEllipses=TRUE, ellipse.level=0.95, palette = "Dark2")

f:id:buriburipipi2:20191114062525j:plain

library(ggplot2)
library(ggfortify)
library("ggrepel")
autoplot(pca,data = d01, loadings = T, loadings.colour =  "black", loadings.label = T,
         loadings.label.colour = "black") +   geom_label_repel(label=rownames(d01))

f:id:buriburipipi2:20191114063052j:plain

2017-12-06

library(dplyr)

library(beeswarm)

x <- rnorm(100, 100, 100)
y <- rnorm(100, 300, 100)
z <- rnorm(100, 500, 100)

#boxplot(x, y, z)
b01 <- beeswarm(list(x, y, z),
col = c("orange", "blue", "magenta"),
method = "square", pch = 20)
# add = TRUE)

library(reshape2)
head(reshape2::french_fries)
molten = reshape2::melt(reshape2::french_fries,
id.vars=c("time", "treatment", "subject", "rep"),
variable.name="flavor", na.rm=TRUE)
head(molten)
tail(molten)

library(dplyr)

molten2 <- reshape2::french_fries %>%
tidyr::gather(key = time , value = flavor,potato, buttery ,grassy ,rancid ,painty)

#head(reshape2::french_fries)
#head (molten)
head (molten2)

2017-08-31

ちょっと毒があったかしらん

某ＱＡサイトで、人事担当らしきひとが、人事考課をまじめにやる人と、適当にみんな普通みたいにつける人とが混在してどう甘辛調整したもんだかというご質問がありそれへの回答

・考課者の研修を定期的におこなって基準を合わせるしかないと思いますけど。
・そうですね。変に調整しちゃうより、一度考課を差し戻す方が無難じゃないでしょうか。
・調整って言っても、不真面目な考課というなら、数式では何ともならない気もしますが。
・まあ、どんな考課をしてもどのグループにいるかで処遇は概ね決まっているみたいな、組織に不公平感が蔓延した結果でそんな風になっているなら、その辺を変える問題意識も必要でしょうね。
・一般論でいうと、人事のその手の問題は組織を滅ぼしかねませんね。どこかで役員が聞きかじったかっこいい仕組みを作ってみたものの魂が入っていないのでどうしようもないでしょう。経営陣と相談されて、適正運用のコンサルを入れたらどうでしょうかね。

2017-07-13

20179713のＬＩＮＫ

韓国の政権交代アジア経済研究所
 【アンケート速報】小学生のお父さん・お母さん3,000人に聞きました　子どもの夏休みの自由研究「親にとって負担」約７割日本能率協会
 ISA540の改訂に係るIAASB公開草案の仮訳日本公認会計士協会
 米国のインフラ整備・投資は進むのか国際協力銀行
 高度自動運転における権限委譲方法の基礎的検討（第2報）～運転以外の作業種類による比較日本自動車研究所
 ショートコメント vol.77 地方圏の人口減少の主因は人口の流出ではない～人口流出が加速のイメージも、実際は増えておらず：関西景気レポートりそな総合研究所
 りゅうぎん調査7月号～経営トップに聞く：生活協同組合コープおきなわ / 沖縄県内における2017年プロ野球春季キャンプの経済効果... りゅうぎん総合研究所
 事例に学ぶ効果的なデータヘルスの実践厚生労働省
 2018年卒マイナビ大学生就職内定率調査＜6月末＞～2017年6月末時点の内々定率は73.3％と、前月（5月末）よりも20pt高い割合になった。採用サポネット
 【キーワード】『街角景気』からわかる、景気に対する実感は？三井住友アセットマネジメント
 微小粒子状物質（PM2.5）中のn-アルカンに関する文献調査日本自動車研究所
 クレーム（苦情）は宝物か～不満感情のしくみを理解する、クレーム対応の5ステップ... JAGAT
日本企業での就業を見据えた留学生支援～OJTにおける外国人社員のコンフリクトへの対処方略に着目して日本学生支援機構
 事業継続計画（BCP）に対する埼玉県企業の意識調査（2017年5月）～BCP策定意向企業は半数以下にとどまる。BCP策定による効果、業務の効率化や取引先の信頼向上も帝国データバンク
 深化するコーポレートガバナンス改革～コーポレートガバナンスを巡る動向を紹介し、取組みを深化させていくためのポイントを解説... KPMG
【事例紹介】香川県における留学生の就職支援～「留学生と地域企業をつなぐ就職関連支援」を中心とした各種事業の紹介：多文化共生社会で生きる日本学生支援機構
 トランプ政権のパリ協定離脱を整理する～その真意と影響は?：誤解だらけのエネルギー・環境問題国際環境経済研究所
 大学生は働くことに何を求めるのか：その特徴と変化採用サポネット
 データを読む：2017年上半期（1-6月）「太陽光関連事業者」の倒産状況～原因別　「販売不振」が半数東京商工リサーチ
 金融緩和頼みの経済成長は持続不可能静岡県立大学グローバル地域センター
 不十分なボーナス還元の理由～構造変化が下押し圧力：Economic Trends 第一生命経済研究所
 「職場の蘇生」～先手を打つ管理者こそが、活力ある職場を創り出す～【第３回管理者が『管理』すること】産業能率大学総合研究所
 ウェブマガジン『留学交流』～特集：多文化共生社会で生きる日本学生支援機構
 進化を続けるリバースモーゲージとヴィアジェ～超高齢社会に向けた英米仏のチャレンジ：ニッセイ基礎研所報ニッセイ基礎研究所
 物価指標のインパクト～３回利上げ＋バランスシート縮小でも金利低下：Market Flash 第一生命経済研究所
 内外経済の中長期展望 2017-2030年度三菱総合研究所
 2017年度　高校生の保護者のためのキャリアガイダンス～保護者から自立支援者へ高校生になったわが子とどう向き合うか？リクルート進学総研
 有機ハイドライド由来不純物が燃料電池性能に及ぼす影響日本自動車研究所
 平成28年度食品産業における取引慣行の実態調査報告書食品産業センター
 デザインを巡る現状と論点：産業競争力とデザインを考える研究会（第1回）配布資料経済産業省
 バーゼル銀行監督委員会による市中協議文書「簡素で、透明性が高く、比較可能な短期証券化商品の自己資本規制上の取扱い」の公表について日本銀行
 リサーチTODAY：日本経済の中長期展望、日本は2020年代に衰退しないみずほ総合研究所
 競合国輸出環境調査（中国・緑茶）報告書（2017年3月）ジェトロ
 教えてＱ＆Ａ：米国企業の業績～米国主要企業の4-6月期決算の見通しは？三井住友アセットマネジメント
 地方経済天気図 2017年7月分～持ち直し全国地方銀行協会
 グローバル企業の温暖化対策第4回：フィリップス～残り5年で「再エネ100%」へアジア拠点向けの調達が鍵みずほ情報総研
 遅すぎる出口戦略の副作用～BISはリーマン・ショック以来の金融危機を警告：欧州経済大和総研
 アジア通貨　備え強化～通貨危機から20年... 国際通貨研究所
 けいざい早わかり 2017年度第3号：大枠合意に達した日ＥＵ・ＥＰＡ三菱UFJリサーチ＆コンサルティング
 迅速かつ大胆な改革が急がれる成長戦略～未来投資戦略2017のポイントとその課題：経済・社会構造分析レポート大和総研
 トランプ政権発足から半年を経て：「トランプは変わるか？」：DCインサイトレポート三菱東京UFJ銀行
 宮田レポート（マンスリー）2017年7月号～原油相場、金相場、長期金利、米国株、ドル円... 三菱UFJモルガン・スタンレー証券
 高速衝突まで対応した軽乗用車前面のエネルギ吸収特性日本自動車研究所
 経済は最悪期を脱出、高まるトルコ国債の投資魅力：Market eyes 大和投資信託
 【書評】福祉国家再編に向けたヒント　田中拓道著『福祉政治史～格差に抗するデモクラシー』東京財団
 日EU・EPA交渉大枠合意の意義～日本の通商戦略の再起動：みずほインサイト政策みずほ総合研究所
 BCBS及びIOSCOによる市中協議文書「簡素で、透明性が高く、比較可能な短期証券化商品を特定する要件」の公表について日本銀行
 企業物価指数（2017年6月）～国内企業物価指数は、前月比0.0％（前年比＋2.1％）。日本銀行
 サービス品質の日米比較～日本のサービス品質は、宅配便や地下鉄、コンビニなど今回対象とした28のサービス分野のほとんどで米国を上回る。日本生産性本部
 嶋中雄二の月例景気報告 No.87～「第3の超景気」下の日本経済～ゴールデン・サイクルは２度来る三菱UFJモルガン・スタンレー証券
 婚活実態調査2017 ～2016年に結婚した人の11.3％が婚活サービスで結婚。婚活サービスを使った結婚が広がりつつある... リクルートマーケティングパートナーズ
 中国科学10大進展について中国総合研究センター
 消費税増税を確実に行うための政策的枠組みが必要大和総研
 長時間労働に関する実態調査 20～30代正社員の労働時間についての実態と意識～男性の長時間労働が目立ち、過労死ライン越えも約1割リクルートマネジメントソリューションズ
 世界経済・金融市場の見通し（2017年7月号）アライアンス・バーンスタイン
 多文化共生2.0の時代へ～総務省の取り組みを中心に：多文化共生社会で生きる日本学生支援機構
 ジャーナルに注目した主要国の論文発表の特徴～オープンアクセス、出版国、使用言語の分析科学技術政策研究所
 大手企業も組織の枠を超え、連携を模索～横つながりでイノベーション加速：ダイバーシティ2.0 vol.2 経済産業省
 ロボット産業市場の現状について～2025年の世界の産業用ロボット市場は2.4兆円に中部圏社会経済研究所
 2017年6月の転職求人倍率を発表。6月は1.87倍（前年同月差＋0.07）。【業界別】コンサルティング業界(5.62倍)、インターネット業界(4.96倍）が高いリクルートキャリア
 日本が支援した缶詰はターゲットに届いているか？（ビエンチャン・ラオス）アジア経済研究所
 ID管理システム導入における現状把握チェックリスト（第１版）日本ネットワークセキュリティ協会

2017-07-13

link test

2017-04-02

ラジオ録音データの再生環境

再生環境
　以前はウォークマンの安いのを使っていたが4年弱使ったらイヤフォンジャックが接触不良になってきて、使うのをあきらめた。
Amazon CAPTCHA
実家で小型エンジンものや電動工具等はノイズキャンセラーで不思議なくらい静かになるので、ウォークマンは便利なのだが高いからね（以前は逆輸入の安いのがチラホラあったのだが）。

今は、AGPtEK A02(ブラック)にている。
www.amazon.co.jp
8GBあるから、記憶容量は十分。1週間分はらくらく入る。
正直、音はそれなり。音楽用としてはどうなんだろうと思うが、なれるからいいのかな。
こいつが致命的なのは、ファイルをコピーした順番で再生してしまうということで、つまり、エクスプローラで複数選択してコピー＆ペーストでは、再生順がぐじゃぐじゃになる。zcopy（by　c.mosさん、、vzで開発してた時代が懐かしい）を使うなどして、行儀よく複写する必要がある。
Software - c.mos
飛ばしで再生が簡単にできるように、rrtで10分刻みで録音しています。音がつなぎ目で重なったりするけど、この辺は我慢ですね。
http://www.geocities.jp/buriburipipi/index.html

イヤフォンは、さすがに低価格プレーヤー付属のはそれなりなので、wp150というのにしている
https://www.amazon.co.jp/dp/B01B2PSIBE
もっと安いイヤフォンでも良いのだが、ケーブルが細いのはどうしても断線しやすい。この辺は使い方というか性格というかなんだろうけど。

2017-03-25

RRTの一括録音用予約ファイル生成プログラム作成

会社でdelphiは定期的に使ってはいるものの、プログラムって、既存資産のCut & Paste　というのが多いので、（基本はファイル処理だから、パターン化される）、自宅で作ってみると、相当手間取ってしまう。
参考書はpdf化してしまっているのだが、検索ができるようにOCR処理をいれたものの使いにくいなあ。
pdfは順読みなら耐えられるけど、ランダムアクセス＆しおりで飛ばすというは、まどろっこしい。
google先生になっちゃうなあ。
　スケジュールファイルを柔軟に汎用的にプログラム生成するって、作ってみると大変かもしれない。今は、1日2回電源断（ＰＣを電源タイマー操作、PCのシャットダウンはrrt側でさすがにやらせる）。電源onの開始タイミング操作もおこなっているので、予約ファイルに書いてあっても、実際には録音していない時間帯もある。
　今は汎用性がほとんどない、自分専用のプロトタイプ状態ですね。

f:id:buriburipipi2:20170325080329j:plain

unit Unit1;

interface

uses
  Winapi.Windows, Winapi.Messages, System.SysUtils, System.Variants, System.Classes, Vcl.Graphics,
  Vcl.Controls, Vcl.Forms, Vcl.Dialogs, Vcl.ComCtrls, Vcl.StdCtrls,dateUtils;

type
  TForm1 = class(TForm)
    DateTimePicker1: TDateTimePicker;
    DateTimePicker2: TDateTimePicker;
    DateTimePicker3: TDateTimePicker;
    Label1: TLabel;
    Label2: TLabel;
    DateTimePicker4: TDateTimePicker;
    Label3: TLabel;
    DateTimePicker5: TDateTimePicker;
    DateTimePicker6: TDateTimePicker;
    DateTimePicker7: TDateTimePicker;
    Edit1: TEdit;
    Button1: TButton;
    procedure Button1Click(Sender: TObject);
    procedure DateTimePicker1Change(Sender: TObject);
  private
    { Private 宣言 }
  public
    { Public 宣言 }
  end;

var
  Form1: TForm1;

  str_date :Tdate;
  end_date :Tdate;
  wk_datetime:Tdatetime;
var
  Hour_c:  Word;
  Min_c:  Word;
  Sec_c: Word;
  MSec_c: Word;
  Hour_w:  Word;
  Min_w:  Word;
  Sec_w: Word;
  MSec_w: Word;

  Hour_s: array[1..4] of Word;
  Min_s: array[1..4] of Word;
  Sec_s: array[1..4] of Word;
  MSec_s: array[1..4] of Word;
  Hour_e: array[1..4] of Word;
  Min_e: array[1..4] of Word;
  Sec_e: array[1..4] of Word;
  MSec_e: array[1..4] of Word;
  Year_s:Word;
  Month_s:Word;
  day_s:Word;
  Year_e:Word;
  Month_e:Word;
  day_e:Word;

//  Hour02s,Min02s,Sec02s,MSec02s:Word;
//  Hour01e,Min01e,Sec01e,MSec01e:Word;
//  Hour02e,Min02e,Sec02e,MSec02e:Word;
  add_min:Word;
implementation

{$R *.dfm}

procedure TForm1.Button1Click(Sender: TObject);
var
  F:textfile;
  buf:array[0..5000] of char;
  wStr,OffStr:string;
  atEndFlag,PrintFlag,PrintFlag_next:boolean;
  base_point:integer;

begin
//   str_date:=Date;
//   datetimePicker1.Date:=str_date;
//   DateTimePicker1.Show;
   AssignFile(F,'d:\test\text.csv');
   System.SetTextBuf(F, Buf);
   ReWrite(F);

   end_date:=datetimePicker4.Date;
   if str_date > end_date then
   begin
      end_date:=str_date +1;
      datetimePicker4.Date:=end_date;
      DateTimePicker4.Show;

   end;
   DecodeDate(Str_date,Year_s,Month_s,day_s);
   DecodeDate(end_date,Year_e,Month_e,day_e);

   DecodeTime(DateTimePicker2.time,Hour_s[1],Min_s[1],Sec_s[1],MSec_s[1]);
   DecodeTime(DateTimePicker5.time,Hour_e[1],Min_e[1],Sec_e[1],MSec_e[1]);

   DecodeTime(DateTimePicker6.time,Hour_s[2],Min_s[2],Sec_s[2],MSec_s[2]);
   DecodeTime(DateTimePicker7.time,Hour_e[2],Min_e[2],Sec_e[2],MSec_e[2]);

   DecodeTime(DateTimePicker3.time,Hour_c,Min_c,Sec_c,MSec_c);
   add_min:= (Hour_c * 60) + Min_c;
//   wk_datetime:= str_date + DateTimePicker2.time;
   wk_dateTime := EncodeDateTime(Year_s, Month_s, Day_s, Hour_s[1], Min_s[1], 0, 0);

//   # 日時,分数,フォルダ,タイトル,電源断,
//   2014/10/01 06:00:00,45,c:\rokuon,福祉政策の課題('14)_01_,0,
   wStr:='# 日時,分数,フォルダ,タイトル,電源断,';
   Writeln(F,wStr);
   atEndFlag:=False;
   base_point:=1;
   PrintFlag:=True;
   PrintFlag_next:=True;

   while (not atEndFlag) do
   begin
     DecodeTime(wk_DateTime,Hour_w,Min_w,Sec_w,MSec_w);
     OffStr:='0';
     if printFlag then
     begin
       if (Hour_e[base_point] = Hour_w) and ( (Min_e[base_point] = Min_w)) then
       begin
         PrintFlag_next:=false;
         inc(base_point);
         OffStr:='1';
       end;
     end
     else
       if (Hour_s[base_point] = Hour_w) and ( (Min_s[base_point] = Min_w)) then
       begin
         PrintFlag_next:=True;
         PrintFlag:=True;
       end;
     wStr:= FormatDateTime('yyyy/MM/DD hh:mm:ss',wk_dateTime) + ',' +
           format('%d',[add_min]) + ',' +
           'c:\rokuon' +',' +
           'tbs' + ',' +
           OffStr + ',';
     if PrintFlag then
        Writeln(F,wStr);
     PrintFlag:=PrintFlag_next;
     wk_dateTime := IncMinute(Wk_dateTime, add_min);
     if Base_Point > 2 then
        Base_Point:=1;
     if Wk_DateTime > end_date then
        atEndFlag:=True;
   end;


   Flush(F);
   CloseFile(F);

end;

procedure TForm1.DateTimePicker1Change(Sender: TObject);
begin
     str_date := datetimePicker1.Date;

end;

end.