2021年4月3日土曜日

tidyverse tidyr dplyr detach search searchpath library library(help= )

 


> search()   
> searchpaths()    # インストール場所の絶対パスを得る

# あるライブラリ中のオブジェクトの一覧を見る library(help=パッケージ名) 
# もちろんそのライブラリが既にインストールされている必要があります。

> library(gstat)       # ライブラリ gstat をロード
> library(help=gstat) # gstat 中のオブジェクト(関数、データセット)の一覧を得る


> library(tidyverse)
─ Attaching packages ─────────────────────────────────────────────────── tidyverse 1.3.0 ─
✓ tibble  3.1.0     ✓ dplyr   1.0.5
✓ tidyr   1.1.3     ✓ stringr 1.4.0
✓ readr   1.4.0     ✓ forcats 0.5.1
✓ purrr   0.3.4     
─ Conflicts ───────────────────────────────────────────────────── tidyverse_conflicts() ─
x stringr::boundary() masks strucchange::boundary()
x dplyr::filter()     masks stats::filter()
x dplyr::first()      masks xts::first()

x dplyr::lag()        masks stats::lag()
x dplyr::last()       masks xts::last()
x dplyr::select()     masks MASS::select()

> data(mtcars)
> mtcars <- as_tibble(mtcars, rownames = "model") %>% mutate(cyl = as.character(cyl))
> g <- ggplot(mtcars, aes(x = mpg, y = wt, color = cyl)) +
+     geom_text(aes(label = model), family = family_sans, fontface = "plain") +
+     labs(x =paste(family_serif, "ボールドを使用"), y = paste(family_serif, "イタリックを使用"),
+          title = paste(family_serif, "ボールドイタリックを使用")) +
+     annotate("text", x = 10, y = 2, label = paste(family_sans, "標準書体を使用"), hjust = 0) +
+     theme(
+         text = element_text(family = family_serif, face = "plain"),
+         title = element_text(face = "bold.italic"),
+         axis.title = element_text(face = "italic"),
+         axis.title.x = element_text(face = "bold")
+     )
> g
Error in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y,  : 
  polygon edge not found
In addition: There were 50 or more warnings (use warnings() to see the first 50)
> g <- ggplot(mtcars, aes(x = mpg, y = wt, color = cyl)) +
+     geom_text(aes(label = model), family = family_sans, fontface = "plain") +
+     labs(x =paste(family_serif, "ボールドを使用"), y = paste(family_serif, "イタリックを使用"),
+          title = paste(family_serif, "ボールドイタリックを使用")) +
+     annotate("text", x = 10, y = 2, label = paste(family_sans, "標準書体を使用"), hjust = 0) +
+     theme(
+         text = element_text(family = family_serif, face = "plain"),
+         title = element_text(face = "bold.italic"),
+         axis.title = element_text(face = "italic"),
+         axis.title.x = element_text(face = "bold")
+     )
> g
Error in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y,  : 
  polygon edge not found
In addition: Warning message:
In grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y,  :
  no font could be found for family "Hiragino Mincho ProN"
>
last(dmdf[,13],10)   # blocked by dplyr
Error in order(order_by)[[n]] : subscript out of bounds
> detach("package:tidyverse", unload=TRUE)
> last(dmdf[,13],10)
Error in order(order_by)[[n]] : subscript out of bounds
> search()
 [1] ".GlobalEnv"           "package:forcats"      "package:stringr"      "package:dplyr"        "package:purrr"        "package:readr"       
 [7] "package:tidyr"        "package:tibble"       "package:RColorBrewer" "package:ggplot2"      "tools:rstudio"        "package:datasets"    
[13] "package:beepr"        "package:forecast"     "package:mondate"      "package:vars"         "package:lmtest"       "package:urca"        
[19] "package:strucchange"  "package:sandwich"     "package:MASS"         "package:utils"        "package:graphics"     "package:grDevices"   
[25] "package:quantmod"     "package:TTR"          "package:xts"          "package:zoo"          "package:stats"        "package:methods"     
[31] "Autoloads"            "package:base"        
> detach("package:dplyr", unload=TRUE)
Warning message:
‘dplyr’ namespace cannot be unloaded:
  namespace ‘dplyr’ is imported by ‘dbplyr’, ‘broom’, ‘tidyr’ so cannot be unloaded 
> detach("package:tidyr", unload=TRUE)
Warning message:
‘tidyr’ namespace cannot be unloaded:
  namespace ‘tidyr’ is imported by ‘broom’ so cannot be unloaded 
> detach("package:broom", unload=TRUE)
Error in detach("package:broom", unload = TRUE) : invalid 'name' argument
> detach("package:tidyr", unload=TRUE)
Error in detach("package:tidyr", unload = TRUE) : invalid 'name' argument
> detach("package:dplyr", unload=TRUE)
Error in detach("package:dplyr", unload = TRUE) : invalid 'name' argument
> last(dmdf[,13],10)
 [1]  7 18  6  7 15 16 20 12 10 23

2021年4月1日木曜日

annotate , xlab , ylab , scale_color_hue, size, legend


df <- data.frame(case_per_capita=as.vector(apply(mdf[,-48],2,sum) / pref_db$x2017),pop_density=pref_db$x2017/pref_db$size,sign=pref_db$x2017,r=pref_db[,2])

p <- ggplot(df, aes(x=pop_density,y=case_per_capita,size=sign,color=r))
p <- p +
xlab("人口密度") + ylab("人口あたり件数")
p <- p + geom_point(alpha=1)
p <- p+annotate("text",label=pref_db[1,3],x=df[1,2], y=df[1,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[11,3],x=df[11,2], y=df[11,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[12,3],x=df[12,2], y=df[12,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[13,3],x=df[13,2], y=df[13,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[14,3],x=df[14,2], y=df[14,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[26,3],x=df[26,2], y=df[26,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[23,3],x=df[23,2], y=df[23,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[27,3],x=df[27,2], y=df[27,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[28,3],x=df[28,2], y=df[28,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[40,3],x=df[40,2], y=df[40,1]+0.1,colour='red',family = "HiraKakuProN-W3")
p <- p+annotate("text",label=pref_db[47,3],x=df[47,2], y=df[47,1]+0.1,colour='red',family = "HiraKakuProN-W3")
# p <- p +scale_color_brewer(palette="Spectral")
# p <- p +scale_color_brewer(palette=rainbow(47))
# p <- p +scale_color_brewer()
p <- p + theme_gray (base_family = "HiraKakuPro-W3")
p <- p + scale_color_hue(name="都道府県",labels=pref_db[,3])
# p <- p + guides(fill = guide_legend(reverse = F,order = 2),label = TRUE)
p <- p + guides(size = guide_legend(title="人口"))
# don't forget to set "color=". otherwise fails to show up.
p <- p + geom_smooth(method = "lm",se=F,color="red",size=1)
# p + scale_colour_manual(values = pref_db[,3])
plot(p)



for the final version, see below.


palette("Alpabet")
df <- data.frame(case_per_capita=as.vector(apply(mdf[,-48],2,sum) / pref_db$x2017),pop_density=pref_db$x2017/pref_db$size,sign=pref_db$x2017,r=pref_db[,2])
# p <- ggplot(df, aes(x=pop_density,y=case_per_capita,size=sign,color=r))
p <- ggplot(df, aes(x=pop_density,y=case_per_capita,size=sign,color=r))
p <- p + xlab("人口密度") + ylab("人口あたり件数")
p <- p + geom_point(alpha=1)
p <- p+annotate("text",label=pref_db[attributes(df[df$case_per_capita > 3 ,])$row.names,3],x=df[attributes(df[df$case_per_capita > 3 ,])$row.names,2], y=df[attributes(df[df$case_per_capita > 3 ,])$row.names,1]+0.1,colour='red',family = "HiraKakuProN-W3")
# p <- p +scale_color_brewer(palette="Spectral")
# p <- p +scale_color_brewer(palette=rainbow(47))
# p <- p +scale_color_brewer()
p <- p + theme_gray (base_family = "HiraKakuPro-W3")
p <- p + scale_color_hue(name="都道府県",labels=pref_db[,3])
# p <- p + guides(fill = guide_legend(reverse = F,order = 2),label = TRUE)
p <- p + guides(size = guide_legend(title="人口"))
# don't forget to set "color=". otherwise fails to show up.
p <- p + geom_smooth(method = "lm",se=F,color="red",size=1)  
plot(p)

df <- data.frame(death_per_capita=as.vector(apply(dmdf[,-48],2,sum) / pref_db$x2017),pop_density=pref_db$x2017/pref_db$size,sign=pref_db$x2017,r=pref_db[,2])
# p <- ggplot(df, aes(x=pop_density,y=case_per_capita,size=sign,color=r))
p <- ggplot(df, aes(x=pop_density,y=death_per_capita,size=sign,color=r))
p <- p + xlab("人口密度") + ylab("人口あたり死者数")
p <- p + geom_point(alpha=1)
p <- p+annotate("text",label=pref_db[attributes(df[df$death_per_capita > 0.07 ,])$row.names,3],x=df[attributes(df[df$death_per_capita > 0.07 ,])$row.names,2], y=df[attributes(df[df$death_per_capita > 0.07 ,])$row.names,1]+0.002,colour='red',family = "HiraKakuProN-W3")

# p <- p +scale_color_brewer(palette="Spectral")
# p <- p +scale_color_brewer(palette=rainbow(47))
# p <- p +scale_color_brewer()
p <- p + theme_gray (base_family = "HiraKakuPro-W3")
p <- p + scale_color_hue(name="都道府県",labels=pref_db[,3])
# p <- p + guides(fill = guide_legend(reverse = F,order = 2),label = TRUE)
p <- p + guides(size = guide_legend(title="人口"))
# don't forget to set "color=". otherwise fails to show up.
p <- p + geom_smooth(method = "lm",se=F,color="red",size=1)  
plot(p)

2021年3月22日月曜日

upgrading to R 4

 # library(magrittr) ??
install.packages("tidyverse")
install.packages("tidyr")
install.packages("beepr")
install.packages("forecast")
install.packages("mondate")
install.packages("vars")
install.packages("quantmod")
install.packages("xts")

upgrade version up

2021年3月14日日曜日

金利 NASDAQ 相関 cor.test() intersect()


choose date both TNX and NDX delta is available.


 d <- intersect(index(diff(TNX[,4])[!is.na(diff(TNX[,4]))]), index(diff(NDX[,4])[!is.na(diff(NDX[,4]))]))

check outputs.


head(as.Date(d))

[1] "2007-01-04" "2007-01-05" "2007-01-09" "2007-01-10" "2007-01-11" "2007-01-12"

Check correlation in 2021.


cor.test(diff(TNX[,4])["2021"][as.Date(d)],diff(NDX[,4])["2021"][as.Date(d)])

Pearson's product-moment correlation
data:  diff(TNX[, 4])["2021"][as.Date(d)] and diff(NDX[, 4])["2021"][as.Date(d)]
t = -2.3358, df = 36, p-value = 0.02519
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 -0.61153962 -0.04874962
sample estimates:
       cor 
-0.3627793 
 警告メッセージ: 
 z + c(-1, 1) * sigma * qnorm((1 + conf.level)/2) で: 
  Recycling array of length 1 in array-vector arithmetic is deprecated.
  Use c() or as.vector() instead.

Check correlation all time.


cor.test(diff(TNX[,4])[as.Date(d)],diff(NDX[,4])[as.Date(d)])

Pearson's product-moment correlation
data:  diff(TNX[, 4])[as.Date(d)] and diff(NDX[, 4])[as.Date(d)]
t = 9.819, df = 2787, p-value < 2.2e-16
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 0.1467388 0.2184915
sample estimates:
      cor 
0.1828587 
 警告メッセージ: 
 z + c(-1, 1) * sigma * qnorm((1 + conf.level)/2) で: 
  Recycling array of length 1 in array-vector arithmetic is deprecated.
  Use c() or as.vector() instead.


> plot.default(diff(TNX[,4])[as.Date(d)],diff(NDX[,4])[as.Date(d)])
> plot.default(diff(TNX[,4])["2021"][as.Date(d)],diff(NDX[,4])["2021"][as.Date(d)])

2021年3月13日土曜日

NASDAQ EPS

download data from here.

modify data like below. use   %s/\s/,/g   in vi.

2020-12-31,1.34
2020-09-30,1.58
2020-06-30,1.45
2020-03-31,1.22
2019-12-31,1.21
       <skip>
2005-06-30,0.13
2005-03-31,0.13

store data into csv file. " ~/nasdaqeps.txt" is used in the sample. change end date accordingly!.

w <- as.xts(read.zoo(read.csv("~/nasdaqeps.txt",header=FALSE)))
ndx_eps <- as.xts(as.vector(na.omit(filter(as.vector(w),rep(1,4)))),seq(as.Date("2005-10-01"),as.Date("2020-10-01"),by='quarters'))

when k2007 is as below.

>  k2007
[1] "2007-01-01::2020-12-31"

summary(lm(apply.quarterly(NDX[,4][k2007],mean) ~ ndx_eps[k2007]+apply.quarterly(PA[k2007],mean)+apply.quarterly(UC[k2007],mean)+apply.quarterly(CS[k2007],mean)))

Call:
lm(formula = apply.quarterly(NDX[, 4][k2007], mean) ~ ndx_eps[k2007] + 
    apply.quarterly(PA[k2007], mean) + apply.quarterly(UC[k2007], 
    mean) + apply.quarterly(CS[k2007], mean))

Residuals:
   Min     1Q Median     3Q    Max 
-787.4 -293.1 -105.6  106.1 2281.4 

Coefficients:
                                   Estimate Std. Error t value Pr(>|t|)    
(Intercept)                      -1.335e+04  8.966e+02 -14.893  < 2e-16 ***
ndx_eps[k2007]                    2.612e+02  7.864e+01   3.321 0.001664 ** 
apply.quarterly(PA[k2007], mean)  7.181e-01  5.575e-02  12.882  < 2e-16 ***
apply.quarterly(UC[k2007], mean) -2.917e+00  9.409e-01  -3.100 0.003144 ** 
apply.quarterly(CS[k2007], mean)  4.479e+01  1.096e+01   4.088 0.000154 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 532.2 on 51 degrees of freedom
Multiple R-squared:  0.9611, Adjusted R-squared:  0.958 
F-statistic: 314.8 on 4 and 51 DF,  p-value: < 2.2e-16


2021年3月10日水曜日

gawk awk 正規表現 エスケープ sed gsed

 正規表現メモ

gawk 正規表現


\w
これは単語を構成する任意のキャラクタ、つまり 文字、数字、それとアンダースコアにマッチする演算子である。 これは [[:alnum:]_] の簡潔な表現とみなして良い。
\W
これは単語を構成する要素にならない任意のキャラクタにマッチする 演算子である。これは [^[:alnum:]_] の簡潔な表現とみなして良い。
\<
これは単語の先頭にある空文字列にマッチする演算子である。 例えば、/\<away/`away'にマッチするが、 `stowaway'にはマッチしない。
\>
これは単語の末尾にある空文字列にマッチする演算子である。 例えば、/stow\>/`stow'にマッチするが、 `stowaway'にはマッチしない。
\y
これは単語の先頭、あるいは末尾の空文字列とマッチする演算子である (つまり語の区切りとマッチするということである)。 例えば、`\yballs?\y'は独立した単語として `ball' にも `balls'にもマッチする。注意! 他のGNU softwareは単語の区切りに\bを使うらしい。
\B
この演算子は単語中の空文字列にマッチする。言い換えると、`\B'は二つ の単語の構成要素文字の間にある空文字列にマッチするということである。例え ば、/\Brat\B/ は`crate'にマッチする。しかし、`dirty rat' にはマッチしない。`B'は簡単にいうと`\y'の反対語である。
\( .... \) () でマッチの塊を扱う
\n  \n は n 番目の () に対応 ([a-z]*) \1*

gsed 正規表現


\< 
        “単語”の先頭。
\> 
        “単語”の末尾。
\b 
        “単語”の先頭、あるいは末尾。言い換えれば、 “語の区切り”とマッチする。
\B 
        “単語”の中にある空文字列にマッチします
\`  
        バッファの先頭。
\'  
        バッファの末尾。
\w 
        “単語”を構成するキャラクタ、つまり文字と数字それにアンダースコア のいずれか。
\W 
        “単語”を構成するキャラクタ以外のキャラクタの いずれか。
\s 
        「空白」一文字とマッチする。localeの影響に注意。
\S
        「空白以外」の一文字とマッチ。

2021年3月9日火曜日

data.frame locale xts format

 
>  last(PI,8)

                PI
2020-06-01 20032.7
2020-07-01 20173.9
2020-08-01 19624.3
2020-09-01 19762.2
2020-10-01 19627.8
2020-11-01 19386.2
2020-12-01 19499.2
2021-01-01 21453.9

上記がある時日付のフォーマットだけを変更する。

Sys.setlocale("LC_TIME","C") 
data.frame(t=(PI  %>% last(.,8) %>%  index() %>% format(., "%Y %b")), d=(PI  %>% last(.,8) %>%  as.vector()))  
Sys.setlocale("LC_TIME", "ja_JP.UTF-8")


YYYY MMM にするのはこちら

Sys.setlocale("LC_TIME","C")

data.frame(t=(PI  %>% last(.,8) %>%  index() %>% format(., "%b %Y")), d=(PI  %>% last(.,8) %>%  as.vector()))

Sys.setlocale("LC_TIME", "ja_JP.UTF-8") 


         t       d
1 2020 Jun 20032.7
2 2020 Jul 20173.9
3 2020 Aug 19624.3
4 2020 Sep 19762.2
5 2020 Oct 19627.8
6 2020 Nov 19386.2
7 2020 Dec 19499.2
8 2021 Jan 21453.9

>   system("locale")

LANG="ja_JP.UTF-8"
LC_COLLATE="ja_JP.UTF-8"
LC_CTYPE="ja_JP.UTF-8"
LC_MESSAGES="ja_JP.UTF-8"
LC_MONETARY="ja_JP.UTF-8"
LC_NUMERIC="ja_JP.UTF-8"
LC_TIME="ja_JP.UTF-8"
LC_ALL=

日付 フォーマット date time format POSIX

 

1.時刻フィールド


%H     時 (00..23)
%I     時 (01..12)
%k     時 (0..23)
%l     時 (1..12)
%M     分 (00..59)
%p     AM あるいは PM のロケール
%r     時刻、12 時間 (hh:mm:ss [AP]M)
%s     1970-01-01 00:00:00 UTC からの秒数 (標準外の拡張)
%S     秒 (00..60)
%T     時刻、24 時間 (hh:mm:ss)
%X     ロケールによる時刻の表現 (%H:%M:%S)
%Z     タイムゾーン (例 EDT)、あるいはタイムゾーンが決定できないならば無し

2.日付フィールド


%a     ロケールの省略形の曜日名 (Sun..Sat)
%A     ロケールの完全表記の曜日名、可変長 (Sunday..Saturday)
%b     ロケールの省略形の月名 (Jan..Dec)
%B     ロケールの完全表記の月名、可変長 (January..December)
%c     ロケールの日付と時刻 (Sat Nov 04 12:02:33 EST 1989)
%d     月内通算日数 (01..31)
%D     日付 (mm/dd/yy)
%h     %b と同じ
%j     年内通算日数 (001..366)
%m     月 (01..12)
%U     日曜日を週の最初の日とした年内通算週 (00..53)
%w     週のうちの曜日 (0..6) (0 が日曜日)
%W     月曜日を週の最初の日とした年内通算週 (00..53)
%x     ロケールの日付表現 (mm/dd/yy)
%y     年の最後の 2 つの数字 (00..99)
%Y     年 (1970...)

Sys.getenv() 環境

 1. 全体

>  Sys.getenv()

__CF_USER_TEXT_ENCODING             0x1F5:0x1:0xE
__CFBundleIdentifier                org.rstudio.RStudio
CLICOLOR_FORCE                      1
COMMAND_MODE                        unix2003
DISPLAY                             /private/tmp/com.apple.launchd.zbStA3g2rP/org.xquartz:0
DYLD_FALLBACK_LIBRARY_PATH          /Library/Frameworks/R.framework/Resources/lib:/Users/honomoto/lib:/usr/local/lib:/usr/lib:::/lib:/Library/Java/JavaVirtualMachines/jdk1.8.0_241.jdk/Contents/Home/jre/lib/server
EDITOR                              vi
GIT_ASKPASS                         rpostback-askpass
HOME                                /Users/honomoto
LANG                                en
LANGUAGE                            en
LC_CTYPE                            en
LC_TYPE                             en

LN_S                                ln -s
LOGNAME                             honomoto
MAKE                                make
MPLENGINE                           tkAgg
PAGER                               /usr/bin/less
PATH                                /usr/bin:/bin:/usr/sbin:/sbin:/usr/local/bin:/Library/TeX/texbin:/opt/X11/bin
R_BROWSER                           /usr/bin/open
R_BZIPCMD                           /usr/bin/bzip2
R_DOC_DIR                           /Library/Frameworks/R.framework/Resources/doc
R_GZIPCMD                           /usr/bin/gzip
R_HOME                              /Library/Frameworks/R.framework/Resources
R_INCLUDE_DIR                       /Library/Frameworks/R.framework/Resources/include
R_LIBS_SITE                         
R_LIBS_USER                         ~/Library/R/4.0/library
R_PAPERSIZE                         a4
R_PDFVIEWER                         /usr/bin/open
R_PLATFORM                          x86_64-apple-darwin17.0
R_PRINTCMD                          lpr
R_QPDF                              /Library/Frameworks/R.framework/Resources/bin/qpdf
R_RD4PDF                            times,inconsolata,hyper
R_SESSION_TMPDIR                    /var/folders/yw/p0wwlz5r8xj4zf001s7n3j780000gn/T//RtmppzWaAC
R_SHARE_DIR                         /Library/Frameworks/R.framework/Resources/share
R_STRIP_SHARED_LIB                  strip -x
R_STRIP_STATIC_LIB                  strip -S
R_SYSTEM_ABI                        macos,gcc,gxx,gfortran,gfortran
R_TEXI2DVICMD                       /usr/local/bin/texi2dvi
R_UNZIPCMD                          /usr/bin/unzip
R_ZIPCMD                            /usr/bin/zip
RMARKDOWN_MATHJAX_PATH              /Applications/RStudio.app/Contents/Resources/resources/mathjax-26
RS_PPM_FD_READ                      9
RS_PPM_FD_WRITE                     10
RS_RPOSTBACK_PATH                   /Applications/RStudio.app/Contents/MacOS/rpostback
RS_SHARED_SECRET                    19408203731238894728102252184
RSTUDIO                             1
RSTUDIO_CONSOLE_COLOR               256
RSTUDIO_CONSOLE_WIDTH               123
RSTUDIO_PANDOC                      /Applications/RStudio.app/Contents/MacOS/pandoc
RSTUDIO_SESSION_PORT                25883
RSTUDIO_USER_IDENTITY               honomoto
RSTUDIO_WINUTILS                    bin/winutils
SED                                 /usr/bin/sed
SHELL                               /bin/bash
SSH_ASKPASS                         rpostback-askpass
SSH_AUTH_SOCK                       /private/tmp/com.apple.launchd.RiShkX8Usr/Listeners
TAR                                 /usr/bin/tar
TERM                                xterm-256color
TMPDIR                              /var/folders/yw/p0wwlz5r8xj4zf001s7n3j780000gn/T/
TZDIR                               macOS
USER                                honomoto
XPC_FLAGS                           0x0
XPC_SERVICE_NAME                    application.org.rstudio.RStudio.262356.262643


2. 日本語

LANG                                ja_JP.UTF-8
LANGUAGE                            jp 
LC_CTYPE                            ja_JP.UTF-8

3. Locale

> Sys.getlocale()
[1] "ja_JP.UTF-8/ja_JP.UTF-8/ja_JP.UTF-8/C/ja_JP.UTF-8/ja_JP.UTF-8"
>   Sys.setlocale("LC_TIME","C")
[1] "C"

>     Sys.getlocale("LC_TIME", "ja_JP.UTF-8/ja_JP.UTF-8/ja_JP.UTF-8/C/ja_JP.UTF-8/ja_JP.UTF-8")

>     system("locale")

LANG="en"
LC_COLLATE="C"
LC_CTYPE="C"
LC_MESSAGES="C"
LC_MONETARY="C"
LC_NUMERIC="C"
LC_TIME="C"
LC_ALL=

> system("locale")

LANG="ja_JP.UTF-8"
LC_COLLATE="ja_JP.UTF-8"
LC_CTYPE="ja_JP.UTF-8"
LC_MESSAGES="ja_JP.UTF-8"
LC_MONETARY="ja_JP.UTF-8"
LC_NUMERIC="ja_JP.UTF-8"
LC_TIME="ja_JP.UTF-8"
LC_ALL=

4. 日本語環境サンプル
system("locale")
LANG="ja_JP.UTF-8"
LC_COLLATE="ja_JP.UTF-8"
LC_CTYPE="ja_JP.UTF-8"
LC_MESSAGES="ja_JP.UTF-8"
LC_MONETARY="ja_JP.UTF-8"
LC_NUMERIC="ja_JP.UTF-8"
LC_TIME="ja_JP.UTF-8"
LC_ALL=

Personal income forecast 個人所得予測 個人所得 personalincome

 2021MAR09



Sys.setlocale("LC_TIME","C")

data.frame(t=(PI  %>% last(.,8) %>%  index() %>% format(., "%b %Y")), d=(PI  %>% last(.,8) %>%  as.vector()))

Sys.setlocale("LC_TIME", "ja_JP.UTF-8") 



1 Jun 2020 20032.7
2 Jul 2020 20173.9
3 Aug 2020 19624.3
4 Sep 2020 19762.2
5 Oct 2020 19627.8
6 Nov 2020 19386.2
7 Dec 2020 19499.2
8 Jan 2021 21453.9

2020JUL31

1 Nov 2020 19435.0
2 Dec 2020 19562.2
3 Jan 2021 21504.5
4 Feb 2021 19955.1
5 Mar 2021 24142.4
6 Apr 2021 20853.2
7 May 2021 20388.2
8 Jun 2021 20414.3

2020SEP15

1 Dec 2020 19562.2
2 Jan 2021 21504.5
3 Feb 2021 19955.1
4 Mar 2021 24142.4
5 Apr 2021 20848.2
6 May 2021 20404.6
7 Jun 2021 20441.8
8 Jul 2021 20667.7

then,

Jun 2020     20032.7
Jul 2020     20173.9
Aug 2020     19624.3
Sep 2020     19762.2
Oct 2020     19628.7
Nov 2020     19386.3
Dec 2020     19491.3
Jan 2021     21462.2
Feb 2021     19945.6
---------------------
Mar 2021     20120 (2021 MAR 21242.06)
Apr 2021     19995
May 2021     20066
Jun 2021     20185
Jul 2021     20254
Aug 2021     20316
Sep 2021     20369



Feb 2021 19945.6
-------------------------
Mar 2021 20840 (2021 MAR 21242.06)
Apr 2021 20150
May 2021 20066
Jun 2021 20185
Jul 2021 20254
Aug 2021 20316
Sep 2021 20369
Oct 2021 20423


May 2021 20804.2 Jun 2021 20659 Jul 2021 20733 Aug 2021 20919 Sep 2021 21229 Oct 2021 21419 Nov 2021 21106 Dec 2021 21261 Jan 2022 21404

2021SEP15

Aug 2021 20740 ±34 Sep 2021 21229 ±39 Oct 2021 21419 ±42 Nov 2021 21106 ±44 Dec 2021 21261 ±46 Jan 2022 21404 ±47 Feb 2022 21469 ±48 Mar 2022 21553 ±49

2021年2月26日金曜日

EPS 2021FEB26

 


eps_year_xts["2020::"]

             [,1]

2020-01-01 116.33
2020-04-01  99.23
2020-07-01  98.22
2020-10-01  96.45
2021-01-01 119.73
2021-04-01 139.51
2021-07-01 146.75
2021-10-01 154.78


$    tac eps.txt | awk '{gsub("\\$","",$NF);print "eps_year_xts[\"2019::\"]["NR"] <- "$NF}'

eps_year_xts["2019::"][1] <- 134.39
eps_year_xts["2019::"][2] <- 135.27
eps_year_xts["2019::"][3] <- 132.90
eps_year_xts["2019::"][4] <- 139.47
eps_year_xts["2019::"][5] <- 116.33
eps_year_xts["2019::"][6] <- 99.23
eps_year_xts["2019::"][7] <- 98.22
eps_year_xts["2019::"][8] <- 91.15
eps_year_xts["2019::"][9] <- 114.36
eps_year_xts["2019::"][10] <- 133.35
eps_year_xts["2019::"][11] <- 141.08
eps_year_xts["2019::"][12] <- 155.56

2021年2月19日金曜日

ゼロサプレス zero supress bash awk split 文字列分割

 bash

$ seq -f %03g 1 10

001
002
003
004
005
006
007
008
009
010


ループ  loop

for i in `seq -f %03g 1 90`; do wget https://blog-imgs-77-origin.fc2.com/u/r/u/urutoraerogazou/matsushima-kaede-166-$i.jpg ; done


for i in `seq 414 453`; do wget  https://blog-imgs-76-origin.fc2.com/d/e/n/densetsuav/1543tachibanarikopin-$i.jpg ; done


二重ループ double loop サンプル

for j in `seq 1 34`
    do for i in `seq 1 12`
      # do wget -r https://javtube.com/javpic/ameri-ichinose/$j/ameri-ichinose-$i.jpg
      do wget -r https://jjgirls.com/japanese/rin-aikawa/$j/rin-aikawa-$i.jpg
  done
done




R


 for(i in seq(1,47,1)){ colnames(mdf)[i] <-  (paste(sprintf("%02d",i),colnames(mdf)[i],sep=""))}


sprintf("%.2f",100.0111)

[1] "100.01"

sprintf("%.3f",100.0111)

[1] "100.011"

sprintf("%.5f",100.0111)

[1] "100.01110"

sprintf("%.6f",100.0111)

[1] "100.011100"



awk


awk '{$3=sprintf("%02d", $3);


find . -print | awk -F\. '{print $2}' | grep k | awk -F\/ '{print $3"-"$2}' | awk '{split($0,a,"-"); print "cp \.\/"a[4]"\/"a[1]"-"a[2]"-"a[3]".jpg  ~\/tmpimage\/kaede\/jjgirls.com\/"a[1]a[2]a[4]"-"a[3]".jpg"}' > kaedecp.txt


split($0,a,"-")で行全体を-で区切って分割する。分割した各要素はa[i]でアクセスできる。

2021年2月16日火曜日

append new entries to xts object.

 

last(tmp.predict,3)

        SP5.Open SP5.High  SP5.Low SP5.Close   SP5.Volume   spline      eps
10 2020 3385.870 3549.850 3233.940  3269.960  89737600000 3475.842 3450.187
11 2020 3296.200 3645.990 3279.740  3621.630 100977880000 3465.431 3396.378
12 2020 3645.870 3760.200 3633.400  3756.070  96056410000 3420.055 3238.727


create the entry for "2021-01-01". first to create single row matrix and convert to the xts.

as.xts(matrix(c(as.vector(apply.monthly(SP5["2021-01"],mean)),3395,3180),nrow=1),as.Date("2021-01-01"))[,-5]

                    [,1]          [,2]            [,3]             [,4]       [,5]                 [,6]    [,7]
2021-01-01 3797.387 3818.136 3768.962 3793.748 5555199474 3395 3180


run append.

append(tmp.predict,as.xts(matrix(c(as.vector(apply.monthly(SP5["2021-01"],mean)),3395,3180),nrow=1),as.Date("2021-01-01"))[,-5]) %>% last(.,6)

        SP5.Open SP5.High  SP5.Low SP5.Close   SP5.Volume   spline      eps
 8 2020 3288.260 3514.770 3284.530  3500.310  84402300000 3270.029 3340.749
 9 2020 3507.440 3588.110 3209.450  3363.000  92084120000 3416.319 3423.604
10 2020 3385.870 3549.850 3233.940  3269.960  89737600000 3475.842 3450.187
11 2020 3296.200 3645.990 3279.740  3621.630 100977880000 3465.431 3396.378
12 2020 3645.870 3760.200 3633.400  3756.070  96056410000 3420.055 3238.727
 1 2021 3797.387 3818.136 3768.962  3793.748   5555199474 3395.000 3180.000


2021年2月5日金曜日

EPS 2021FEB05

 


eps_year_xts["2020::"]

             [,1]

2020-01-01 116.33
2020-04-01  99.23
2020-07-01  96.47
2020-10-01  92.18
2021-01-01 106.60
2021-04-01 125.68
2021-07-01 132.68
2021-10-01 140.24

$ tac eps.txt | awk '{gsub("\\$","",$NF);print "eps_year_xts[\"2019::\"]["NR"] <- "$NF}'

eps_year_xts["2019::"][1] <- 134.39
eps_year_xts["2019::"][2] <- 135.27
eps_year_xts["2019::"][3] <- 132.90
eps_year_xts["2019::"][4] <- 139.47
eps_year_xts["2019::"][5] <- 116.33
eps_year_xts["2019::"][6] <- 99.23
eps_year_xts["2019::"][7] <- 98.22
eps_year_xts["2019::"][8] <- 96.45
eps_year_xts["2019::"][9] <- 119.73
eps_year_xts["2019::"][10] <- 139.51
eps_year_xts["2019::"][11] <- 146.75
eps_year_xts["2019::"][12] <- 154.78