-
Notifications
You must be signed in to change notification settings - Fork 1
/
阿里报告 - II.Rmd
743 lines (598 loc) · 27.6 KB
/
阿里报告 - II.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
---
title: "阿里彩票"
subtitle: "分析阿里数据 (第II部)"
author: "雷欧 <img src='文艺坊图库/QQ图片20200720121653.png' height='14'>"
date: "`r lubridate::today('Asia/Tokyo')`"
output:
html_document:
mathjax: https://cdn.jsdelivr.net/npm/mathjax@3/es5/tex-mml-chtml.js
number_sections: yes
toc: yes
toc_depth: 4
toc_float:
collapsed: yes
smooth_scroll: yes
code_folding: hide
css: CSSBackgrounds.css
---
# 主题
<img src='文艺坊图库/商场如战场.png' width='560'>
[[**主题曲**]{style="color:blue"}](https://github.com/englianhu/report/discussions)
🚄东方快车,🚄一带一路。
<audio controls loop autoplay src="文艺坊歌曲库/東方快車合唱團 Oriental Express - 紅紅青春敲呀敲 《黑松沙士》廣告主題曲.mp3" controls></audio>
<br>
# 设定
## SCSS 设置
<style>
pre {
overflow-x: auto;
}
pre code {
word-wrap: normal;
white-space: pre;
}
.table-hover > tbody > tr:hover {
background-color: #8D918D;
}
</style>
```{r load-sass, class.source='bg-success', class.output='bg-primary'}
# install.packages('remotes', dependencies = TRUE, INSTALL_opts = '--no-lock')
library('BBmisc', 'rmsfuns')
#remotes::install_github("rstudio/sass")
lib('sass')
## https://support.rstudio.com/hc/en-us/articles/200532197
## https://community.rstudio.com/t/r-does-not-display-korean-chinese/30889/3?u=englianhu
#Sys.setlocale("LC_CTYPE", "en_US.UTF-8")
#Sys.setlocale("LC_CTYPE", "zh_CN.UTF-8")
#Sys.setlocale(category = "LC_CTYPE", "Chinese (Simplified)_China.936")
#Sys.setlocale(locale = "Chinese")
#Sys.setlocale(locale = "Japanese")
#Sys.setlocale(locale = "English")
# rmarkdown::render('/home/englianhu/Documents/owner/ryo-cn.Rmd', encoding = 'UTF-8')
#Sys.setlocale("LC_CTYPE", "UTF-8")
#Sys.setlocale(locale = "UTF-8")
#Sys.setlocale(category = "LC_ALL", locale = "chs")
#Sys.setlocale(category = "LC_ALL", locale = "UTF-8")
#Sys.setlocale(category = "LC_ALL", locale = "Chinese")
#Sys.setlocale(category = "LC_ALL", locale = "zh_CN.UTF-8")
Sys.setlocale("LC_ALL", "en_US.UTF-8")
```
```{scss set-scss, class.source='bg-success', class.output='bg-primary'}
/* https://stackoverflow.com/a/66029010/3806250 */
h1 { color: #002C54; }
h2 { color: #2F496E; }
h3 { color: #375E97; }
h4 { color: #556DAC; }
h5 { color: #92AAC7; }
/* ----------------------------------------------------------------- */
/* https://gist.github.com/himynameisdave/c7a7ed14500d29e58149#file-broken-gradient-animation-less */
.hover01 {
/* color: #FFD64D; */
background: linear-gradient(155deg, #EDAE01 0%, #FFEB94 100%);
transition: all 0.45s;
&:hover{
background: linear-gradient(155deg, #EDAE01 20%, #FFEB94 80%);
}
}
.hover02 {
color: #FFD64D;
background: linear-gradient(155deg, #002C54 0%, #4CB5F5 100%);
transition: all 0.45s;
&:hover{
background: linear-gradient(155deg, #002C54 20%, #4CB5F5 80%);
}
}
.hover03 {
color: #FFD64D;
background: linear-gradient(155deg, #A10115 0%, #FF3C5C 100%);
transition: all 0.45s;
&:hover{
background: linear-gradient(155deg, #A10115 20%, #FF3C5C 80%);
}
}
```
```{r gb-opts, class.source='hover01', class.output='hover02'}
## Set the timezone but not change the datetime
Sys.setenv(TZ = 'Asia/Shanghai')
## Setting to omit all warnings
## https://stackoverflow.com/a/36846793/3806250
## Set width
## options(knitr.table.format = 'html') will set all kableExtra tables to be 'html', otherwise need to set the parameter on every single table.
options(warn = -1, width = 999, knitr.table.format = 'html')#, digits.secs = 6)
## https://stackoverflow.com/questions/39417003/long-vectors-not-supported-yet-abnor-in-rmd-but-not-in-r-script
## https://yihui.org/knitr/options
knitr::opts_chunk$set(
class.source = 'hover01', class.output = 'hover02', class.error = 'hover03',
message = FALSE, warning = FALSE, error = TRUE,
autodep = TRUE, aniopts = 'loop',
progress = TRUE, verbose = TRUE,
cache = FALSE, cache.lazy = FALSE, result = 'asis')
```
<br><br>
## 设置
```{r libs, eval=FALSE}
## Setup Options, Loading Required Libraries and Preparing Environment
## Loading the packages and setting adjustment
source('函数/libs.R')
```
读取程序包
```{r load-pkgs}
## Setup Options, Loading Required Libraries and Preparing Environment
## 3210448065@qq.com
## leiou123
## 2849108450@qq.com
## leiou123
## https://rstudio.cloud/project/1198888
## Loading the package 'BBmisc'
if(suppressMessages(!require('BBmisc'))){
install.packages('BBmisc', dependencies = TRUE, INSTALL_opts = '--no-lock')
suppressMessages(library('BBmisc'))
}
if (suppressMessages(!require('rmsfuns'))) {
install.packages('rmsfuns', dependencies = TRUE, INSTALL_opts = '--no-lock')
suppressMessages(library('rmsfuns'))
}
if(!require('REmap')) devtools::install_github('lchiffon/REmap')
## Loading multiple packages at once
#pkgs <- c('readr', 'plyr', 'dplyr', 'magrittr', 'tidyverse', 'devtools', 'zoo', 'lubridate', 'stringr', 'rvest', 'markdown', 'googleVis', 'knitr', 'rmarkdown', 'htmltools', 'knitr', 'kableExtra', 'formattable', 'echarts4r', 'radarchart', 'MASS', 'htmlwidgets', 'maps', 'REmap', 'ggmap', 'vembedr')
pkgs <- c('plyr', 'dplyr', 'magrittr', 'stringr', 'knitr', 'kableExtra', 'lubridate', 'broom', 'purrr', 'readxl', 'tibble', 'DT')
suppressAll(lib(pkgs))
load_pkg(pkgs)
rm(pkgs)
## Set the googleVis options first to change the behaviour of plot.gvis, so that
## only the chart component of the HTML file is written into the output file.
op <- options(gvis.plot.tag = 'chart')
## <audio src='music/bigmoney.mp3' autoplay controls loop></audio>
```
```{r set-prefer}
conflict_prefer('filter', 'dplyr')
conflict_prefer('select', 'dplyr')
conflict_prefer('mutate', 'dplyr')
conflict_prefer('rename', 'dplyr')
```
<br><br>
# 介绍
今天尝试分析阿里**样本阿里数据**,分析`总添加`、`总注册`、`总首冲`、`总投注`、`总流水`、`总充值`、`总盈利`、`净投注`、`总体留存率`、`人均充值`能产生多少`负盈利`。需要预测公司整体`三个月`、`半年`、`一年`报表数据最佳,不过只有两个月数据也可以试试。
# 数据
读取样本数据。
```{r warning=FALSE}
## 读取数据
smp <- suppressAll(read_excel('.xlsx') %>% tbl_df)
names(smp) <- smp[1,]; smp %<>% .[-1,]
# 总盈利, 总添加, 总注册, 总首冲, 总投注, 总流水, 总充值
# names(smp) <- c('Date', 'PL', 'Add', 'RG', 'FD', 'TB', 'TO', 'TD')
smp %<>% mutate_if(is.character, as.numeric)
smp$日期 %<>% as.numeric %>% as.Date()
smp %>%
kable(caption = '阿里数据') %>%
kable_styling(bootstrap_options = c('striped', 'hover', 'condensed', 'responsive')) %>%
scroll_box(width = '100%', height = '400px')
```
# 线性模型
统计建模。
```{r warning=FALSE}
## 模型比较
m1 <- lm(总盈利 ~ 总添加 + 总注册 + 总首冲 + 总投注 + 总流水 + 总充值, data = smp)
m2 <- lm(总盈利 ~ 总添加 + 总注册 + 总首冲 + 总投注 + 总流水, data = smp)
m3 <- lm(总盈利 ~ 总添加 + 总注册 + 总首冲 + 总投注, data = smp)
m4 <- lm(总盈利 ~ 总添加 + 总注册 + 总首冲, data = smp)
m <- list(m1 = m1, m2 = m2, m3 = m3, m4 = m4)
rm(m1, m2, m3, m4)
```
# 移动模型
## 测试1个月移动数据
```{r warning=FALSE, message=FALSE}
## 转换数据类型
nm <- c('Date', 'PL', 'Add', 'RG', 'FD', 'TB', 'TO', 'TD', 'NB', 'NBR', 'AD')
names(smp) <- nm
smp %<>% na.omit
datt <- llply(smp[,-1], function(x) {
x %<>% xts(order.by = smp$Date)
x
})
## 阿里数据 1个月
dbt <- smp$Date[month(1):length(smp$Date)]
datt <- llply(dbt, function(x) {
dbtt <- smp$Date[smp$Date <= x]
dbtt <- dbtt[(length(dbtt) - month(1)):length(dbtt)]
z <- llply(datt, function(y) {
y <- y[index(y) %in% dbtt]
xx <- auto.arima(y, D=1)
fst <- forecast(xx, h=1) %>%
tbl_df %>%
select('Point Forecast') %>%
rename(`FP` = `Point Forecast`)
names(y) %<>% str_extract_all('[A-Z]{1,}')
y <- y[nrow(y)]
zz <- data.frame(Date = index(y), y, fst) %>% tbl_df
zz
}) %>% bind_cols
}) %>% bind_rows
datt %<>% select(-Date1, -Date2, -Date3, -Date4, -Date5,
-Date6, -Date7, -Date8, -Date9) %>%
mutate(FP = c(NA, FP[-length(FP)]), FP1 = c(NA, FP1[-length(FP1)]), FP2 = c(NA, FP2[-length(FP2)]), FP3 = c(NA, FP3[-length(FP3)]), FP5 = c(NA, FP5[-length(FP5)]), FP6 = c(NA, FP6[-length(FP6)]), FP7 = c(NA, FP7[-length(FP7)]), FP8 = c(NA, FP8[-length(FP8)]), FP9 = c(NA, FP9[-length(FP9)]))
datt$FP4 = c(NA, datt$FP4[-length(datt$FP4)])
names(datt) <- c('日期', '总添加', '预测总添加',
'总注册', '预测总注册', '总首冲', '预测总首冲',
'总投注', '预测总投注', '总流水', '预测总流水',
'总充值', '预测总充值', '总盈利', '预测总盈利',
'净投注', '预测净投注', '总体留存率', '预测总体留存率',
'人均充值', '预测人均充值')
```
### 绘图
```{r warning=FALSE}
## 绘图
## 预测总添加 - 总添加
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总添加, colour = 总添加)) +
geom_line(aes(y = 预测总添加, colour = 预测总添加))
```
上图显示预测总添加。
```{r warning=FALSE}
## 预测总注册 - 总注册
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总注册, colour = 总注册)) +
geom_line(aes(y = 预测总注册, colour = 预测总注册))
```
上图显示预测总注册。
```{r warning=FALSE}
## 预测总投注 - 总投注
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总投注, colour = 总投注)) +
geom_line(aes(y = 预测总投注, colour = 预测总投注))
```
上图显示预测总投注。
```{r warning=FALSE}
## 预测总流水 - 总流水
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总流水, colour = 总流水)) +
geom_line(aes(y = 预测总流水, colour = 预测总流水))
```
上图显示预测总流水。
```{r warning=FALSE}
## 预测总充值 - 总充值
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总充值, colour = 总充值)) +
geom_line(aes(y = 预测总充值, colour = 预测总充值))
```
上图显示预测总充值。
```{r warning=FALSE}
## 预测总盈利 - 总盈利
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总盈利, colour = 总盈利)) +
geom_line(aes(y = 预测总盈利, colour = 预测总盈利))
```
上图显示预测总盈利。
```{r warning=FALSE}
## 预测净投注 - 净投注
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 净投注, colour = 净投注)) +
geom_line(aes(y = 预测净投注, colour = 预测净投注))
```
上图显示预测总投注。
```{r warning=FALSE}
## 预测总体留存率 - 总体留存率
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总体留存率, colour = 总体留存率)) +
geom_line(aes(y = 预测总体留存率, colour = 预测总体留存率))
```
上图显示预测总体留存率。
```{r warning=FALSE}
## 预测人均充值 - 人均充值
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 人均充值, colour = 人均充值)) +
geom_line(aes(y = 预测人均充值, colour = 预测人均充值))
```
上图显示预测总人均充值。
### 比较预测值精准度
```{r warning=FALSE}
datt %>% na.omit %>%
kable(caption = '预测样本数据') %>%
kable_styling(bootstrap_options = c('striped', 'hover', 'condensed', 'responsive')) %>%
scroll_box(width = '100%', height = '400px')
```
```{r warning=FALSE}
## MSE 比较预测值精准度
mse <- datt %>% na.omit %>%
mutate(
均方误差.总添加 = mean((预测总添加 - 总添加)^2),
均方误差.总注册 = mean((预测总注册 - 总注册)^2),
均方误差.总投注 = mean((预测总投注 - 总投注)^2),
均方误差.总流水 = mean((预测总流水 - 总流水)^2),
均方误差.总充值 = mean((预测总充值 - 总充值)^2),
均方误差.总盈利 = mean((预测总盈利 - 总盈利)^2),
均方误差.净投注 = mean((预测净投注 - 净投注)^2),
均方误差.总体留存率 = mean((预测总体留存率 - 总体留存率)^2),
均方误差.人均充值 = mean((预测人均充值 - 人均充值)^2))
mse1 <- mse %>%
select(均方误差.总添加, 均方误差.总注册, 均方误差.总投注,
均方误差.总流水, 均方误差.总充值, 均方误差.总盈利,
均方误差.净投注, 均方误差.总体留存率, 均方误差.人均充值) %>%
unique
suppressAll(mse1 %>%
kable(booktabs = T) %>%
kable_styling())
```
## 测试40天移动数据
```{r warning=FALSE, message=FALSE}
## 转换数据类型
nm <- c('Date', 'PL', 'Add', 'RG', 'FD', 'TB', 'TO', 'TD', 'NB', 'NBR', 'AD')
names(smp) <- nm
smp %<>% na.omit
datt <- llply(smp[,-1], function(x) {
x %<>% xts(order.by = smp$Date)
x
})
## 阿里数据 40天
dbt <- smp$Date[40:length(smp$Date)]
datt <- llply(dbt, function(x) {
dbtt <- smp$Date[smp$Date <= x]
dbtt <- dbtt[(length(dbtt) - 40):length(dbtt)]
z <- llply(datt, function(y) {
y <- y[index(y) %in% dbtt]
xx <- auto.arima(y, D=1)
fst <- forecast(xx, h=1) %>%
tbl_df %>% select('Point Forecast') %>%
rename(`FP` = `Point Forecast`)
names(y) %<>% str_extract_all('[A-Z]{1,}')
y <- y[nrow(y)]
zz <- data.frame(Date = index(y), y, fst) %>% tbl_df
zz
}) %>% bind_cols
}) %>% bind_rows
datt %<>% select(-Date1, -Date2, -Date3, -Date4, -Date5,
-Date6, -Date7, -Date8, -Date9) %>%
mutate(FP = c(NA, FP[-length(FP)]), FP1 = c(NA, FP1[-length(FP1)]), FP2 = c(NA, FP2[-length(FP2)]), FP3 = c(NA, FP3[-length(FP3)]), FP5 = c(NA, FP5[-length(FP5)]), FP6 = c(NA, FP6[-length(FP6)]), FP7 = c(NA, FP7[-length(FP7)]), FP8 = c(NA, FP8[-length(FP8)]), FP9 = c(NA, FP9[-length(FP9)]))
datt$FP4 = c(NA, datt$FP4[-length(datt$FP4)])
names(datt) <- c('日期', '总添加', '预测总添加',
'总注册', '预测总注册', '总首冲', '预测总首冲',
'总投注', '预测总投注', '总流水', '预测总流水',
'总充值', '预测总充值', '总盈利', '预测总盈利',
'净投注', '预测净投注', '总体留存率', '预测总体留存率',
'人均充值', '预测人均充值')
```
### 绘图
```{r warning=FALSE}
## 绘图
## 预测总添加 - 总添加
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总添加, colour = 总添加)) +
geom_line(aes(y = 预测总添加, colour = 预测总添加))
```
上图显示预测总添加。
```{r warning=FALSE}
## 预测总注册 - 总注册
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总注册, colour = 总注册)) +
geom_line(aes(y = 预测总注册, colour = 预测总注册))
```
上图显示预测总注册。
```{r warning=FALSE}
## 预测总投注 - 总投注
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总投注, colour = 总投注)) +
geom_line(aes(y = 预测总投注, colour = 预测总投注))
```
上图显示预测总投注。
```{r warning=FALSE}
## 预测总流水 - 总流水
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总流水, colour = 总流水)) +
geom_line(aes(y = 预测总流水, colour = 预测总流水))
```
上图显示预测总流水。
```{r warning=FALSE}
## 预测总充值 - 总充值
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总充值, colour = 总充值)) +
geom_line(aes(y = 预测总充值, colour = 预测总充值))
```
上图显示预测总充值。
```{r warning=FALSE}
## 预测总盈利 - 总盈利
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总盈利, colour = 总盈利)) +
geom_line(aes(y = 预测总盈利, colour = 预测总盈利))
```
上图显示预测总盈利。
```{r warning=FALSE}
## 预测净投注 - 净投注
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 净投注, colour = 净投注)) +
geom_line(aes(y = 预测净投注, colour = 预测净投注))
```
上图显示预测总投注。
```{r warning=FALSE}
## 预测总体留存率 - 总体留存率
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 总体留存率, colour = 总体留存率)) +
geom_line(aes(y = 预测总体留存率, colour = 预测总体留存率))
```
上图显示预测总体留存率。
```{r warning=FALSE}
## 预测人均充值 - 人均充值
datt %>% ggplot(aes(日期)) +
geom_line(aes(y = 人均充值, colour = 人均充值)) +
geom_line(aes(y = 预测人均充值, colour = 预测人均充值))
```
上图显示预测总人均充值。
### 比较预测值精准度
```{r warning=FALSE}
datt %>% na.omit %>%
kable(caption = '预测样本数据') %>%
kable_styling(bootstrap_options = c('striped', 'hover', 'condensed', 'responsive')) %>%
scroll_box(width = '100%', height = '400px')
```
```{r warning=FALSE}
## MSE 比较预测值精准度
mse <- datt %>% na.omit %>%
mutate(均方误差.总添加 = mean((预测总添加 - 总添加)^2),
均方误差.总注册 = mean((预测总注册 - 总注册)^2),
均方误差.总投注 = mean((预测总投注 - 总投注)^2),
均方误差.总流水 = mean((预测总流水 - 总流水)^2),
均方误差.总充值 = mean((预测总充值 - 总充值)^2),
均方误差.总盈利 = mean((预测总盈利 - 总盈利)^2),
均方误差.净投注 = mean((预测净投注 - 净投注)^2),
均方误差.总体留存率 = mean((预测总体留存率 - 总体留存率)^2),
均方误差.人均充值 = mean((预测人均充值 - 人均充值)^2))
mse2 <- mse %>%
select(均方误差.总添加, 均方误差.总注册, 均方误差.总投注,
均方误差.总流水, 均方误差.总充值, 均方误差.总盈利,
均方误差.净投注, 均方误差.总体留存率, 均方误差.人均充值) %>%
unique
suppressWarnings(mse2 %>%
kable(booktabs = T) %>%
kable_styling())
```
# 结论
## 筛选最优线性模型
筛选最优线性模型。
```{r warning=FALSE}
aic <- suppressAll(llply(1:length(m), function(i) {
y <- data.frame(paste0('m', i), t(data.frame(extractAIC(m[[i]]))))
names(y) <- (c('模型', '自由度', 'aic'))
y
}) %>% bind_rows)
## 筛选最佳模型
aic %>%
kable(booktabs = T) %>%
kable_styling() %>%
row_spec(aic$aic %>% which.min, bold = T, color = "white", background = 'goldenrod')
```
```{r warning=FALSE}
suppressWarnings(m[[aic %>% .$aic %>% which.min]] %>%
tidy %>%
kable(booktabs = T) %>%
kable_styling())
```
阐述最佳模型。
```
- `截距(Intercept)`=是在没有任何变量的默认下的负盈利,负盈利为`r coef(m[[aic %>% .$aic %>% which.min]])[1]`。
- `总添加`=每添加一个`总添加`,即可产生`r coef(m[[aic %>% .$aic %>% which.min]])[2]`负盈利。
- `总注册`=每添加一个`总注册`,即可产生`r coef(m[[aic %>% .$aic %>% which.min]])[3]`负盈利。
- `总首冲`=每添加一个`总首冲`,即可产生`r coef(m[[aic %>% .$aic %>% which.min]])[4]`负盈利。
- `总投注`=每添加一个`总投注`,即可产生`r coef(m[[aic %>% .$aic %>% which.min]])[5]`负盈利。
- `总流水`=每添加一个`总流水`,即可产生`r coef(m[[aic %>% .$aic %>% which.min]])[6]`负盈利。
- `总充值`=每添加一个`总充值`,即可产生`r coef(m[[aic %>% .$aic %>% which.min]])[7]`负盈利。
- `总添加`+`总注册`,即可产生`r coef(m[[aic %>% .$aic %>% which.min]])[2]`+`r coef(m[[aic %>% .$aic %>% which.min]])[3]`=`r coef(m[[aic %>% .$aic %>% which.min]])[2] + coef(m[[aic %>% .$aic %>% which.min]])[3]`负盈利,以此类推。
```
## 筛选最优移动模型
筛选最优移动模型。
```{r warning=FALSE}
mse3 <- data.frame('数据量大小' = c('14天', '21天'), rbind(mse1, mse2))
mse3 <- data.frame(mse3, '总体均方误差' = rowMeans(mse3[-1]))
mse3 %>%
kable(booktabs = T) %>%
kable_styling()
```
根据以上均方误差比较,`...`数据会比`...`数据来得精准。此报表只是个样本数据。
## 总结论
- 分析好公司报表后,
- 还需要分析`彩种`、`赔率`、`流水`、`充值`、`负盈利`数据。
- 还需要分析`进粉成本`、`租金`等等其他成本,`负盈利`数据。
```
中博手游
https://zbcp.yayaread.cn/#/home
已有系统优势:
产品多、优惠多...
- 主打产品置放在显眼的位置,还有走势图方便玩家参考。
- 热门彩种列在左手边,方便玩家们选择。
- 开奖通告、资讯中心、最新中奖榜都在主页,方便玩家一览即可。
- 区分高频与低频彩种等等产品...
额外系统提升(人工智能化):
- 如果会员可以在登录时,根据自己的等级,自动弹出个等级级别,还需要多少流水或充值量才能晋级,刺激玩家充值来玩晋级,刺激消费。
- 如果玩家或代理盈亏到一定的(后台可以自己设置),就会弹出开发代理的流水返点,然后在图表旁边有教程视频或图,方便玩家操作,刺激消费。
- 依据每个玩家平时玩什么彩种,然后自动推荐相关彩种,刺激消费。
- 存取款旁也有教程视频或者图。
- 打开网站5秒未点击任何按钮,就会自动转到注册页面。
- 自动出现默认邀请码,那么就不会空有流量却没有注册人数。(设置个小组专门负责这些会员,额外提升公司业绩)
- 注册页面需要填写微信、QQ和邮箱,那么有优惠活动时可以通知会员。(回访提升公司业绩)
- 注册页面和首次要充值的页面,左右两边会出现首充优惠、晋级奖励和每日嘉奖优惠。
- 已注册尚未首充,一登录账户,就会弹出充值页面,左右两边会出现首充优惠、晋级奖励和每日嘉奖优惠。
- 已充值过的会员,一到充值页面,左右两边会出现转卡优惠(根据周一至周五,周六周日的优惠而自动分辨)、晋级奖励和每日嘉奖优惠。
- 系统自动检测每个会员的盈亏是否达到指定的百分比(比方说赢了60%或亏了40%,可以设置每周或每月),然后会弹出信息提示。可以设置总充值低于多少的话,就不会弹出提示。(自动根据各别会员的盈亏,自动产生不同的盈利图,和宏图,刺激消费)
- 如果亏超过40%的会员,一登录账户就会弹出要不要接受亏损金。(如果亏超过40%的会员,也会自动收到邮件、微信、QQ通知,要不要来网站领取亏损金)。
- 如果赢超过60%,系统自动弹出提示,那么方便我们进一步处理。(包括查看玩什么彩种和几倍流水,卡款等等措施)
- 如果赢超过60%的会员,只要金额都不大(可以设置多少金额)的话,就会收到微信、QQ、邮件建议激活邀请码当代理,然后内容引导怎么操作和利润(难保他会开发到大会员)。
- 风险管理,以足彩为例子,就是一开奖后,赔率就会自动调低,然后目前的最大投注额是100万(可以随着投注金额的大小自动变动最大投注额,或者累计投注额到了60%就会自动调动),可以根据累计投注额,设置赔率调动。
- 开发下级:添加个上载图片视频的功能,比方说单带多少的赚多少、介绍平台、首充优惠、转卡优惠等等图片视频。
- 体验金,智能弹出开发代理视窗,可以直接激活邀请码,并且可以输入所有朋友的QQ或微信号,然后系统提供开发下级教程、优惠等等图和视频,直接群发。
- 吃喝嫖赌都是传统行业,可以尝试和彩票杂志、烟酒、网贷或银行贷款、情色合作,添加个额外服务(凡是顾客买杂志、烟酒或嫖妓,那些商家可以成为我们代理,然后直接通过手机提供应用下载或二维码和邀请码开发下级,并且推销优惠活动和教导如何注册充值。亏钱的会员可以自行通过网贷刺激消费,只是法律上的需要考量)
- 存取款和网址都能智能转换,比方说有备用5到8个
- 第一:使用模拟号投注赢钱录制视频,账号打马赛克;第二:视频赢钱欢呼,然后在视频内有个输入金额和马上投注,方便消费,刺激消费。
- 第一:高流量不注册的话,可以和其它商家合作;第二:注册时填写职业、收入、兴趣、买车买房等等数据(可以添加个黄页,通过supper cookies收集浏览者的爱好的数据),然后根据那些数据自动弹出广告。赚取广告费。
- API
个人觉得比我们目前的平台好很多...
公司目前的商业模式已经可以赚到钱,不过可能公司寿命只有五六年,如果走品牌路线的话,可能可以学习威廉希尔、立博那些网站,靠品牌经营五六十年,不过成本与盈利就是的考量了...
```
<br><br>
# 附录
## 文书明细
以下乃此文书的文件信息。
- 文集建立日:2022-05-31
- 文集最新更新日:`r today('Asia/Shanghai')`
- `r R.version.string`
- [**rmarkdown**](https://github.com/rstudio/rmarkdown) 程序包版本:`r packageVersion('rmarkdown')`
- 文集版本:0.2.1
- 文集作者:[®γσ, ξηg Lιαη Ημ](https://rpubs.com/englianhu/ryo-cn)
- 猫舍:[源代码](https://github.com/englianhu/report)
- 追加附属信息
```{r info}
#suppressMessages(require('formattable', quietly = TRUE))
#suppressMessages(require('knitr', quietly = TRUE))
#suppressMessages(require('kableExtra', quietly = TRUE))
#suppressMessages(require('magittr', quietly = TRUE))
#suppressMessages(require('devtools', quietly = TRUE))
sys1 <- session_info()$platform |>
unlist() |>
{\(.) data.frame(row.names = 1:length(.),
Category = names(.), session_info = .)}()
sys2 <- data.frame(Sys.info()) |>
{\(.) data.frame(Category = row.names(.), Sys.info = .[,1])}()
#remarks, dim(sys1), dim(sys2)
if (nrow(sys1) == 11 && nrow(sys2) == 8) {
sys2 <- sys2 |>
{\(.) rbind(., data.frame(
Category = c('rmarkdown', 'rsconnect', 'Current time'),
Sys.info = c(as.character(getwd()),
as.character(packageVersion('rsconnect')),
paste(as.character(lubridate::now('Asia/Shanghai')), 'CST 🗺'))))}()
} else if (nrow(sys1) == 10 && nrow(sys2) == 8) {
sys1 <- rbind(sys1, data.frame(Category = '', session_info = ''))
sys2 <- sys2 |>
{\(.) rbind(., data.frame(
Category = c('rmarkdown', 'rsconnect', 'Current time'),
Sys.info = c(as.character(getwd()),
as.character(packageVersion('rsconnect')),
paste(as.character(lubridate::now('Asia/Shanghai')), 'CST 🗺'))))}()
}
sys <- cbind(sys1, sys2) |>
{\(.)
kbl(., caption = 'Additional session information:')}() |>
{\(.)
kable_styling(., bootstrap_options = c('striped', 'hover', 'condensed', 'responsive'))}() |>
{\(.)
row_spec(., 0, background = 'DimGrey', color = 'yellow')}() |>
{\(.)
column_spec(., 1, background = 'CornflowerBlue', color = 'red')}() |>
{\(.)
column_spec(., 2, background = 'grey', color = 'black')}() |>
{\(.)
column_spec(., 3, background = 'CornflowerBlue', color = 'blue')}() |>
{\(.)
column_spec(., 4, background = 'grey', color = 'white')}() |>
{\(.)
row_spec(., 11, bold = TRUE, color = 'yellow', background = '#D7261E')}()
rm(sys1, sys2)
sys
```
## 参考文献
1) []()
2) [GitHub : Statistical Rethinking (2022 Edition)](https://github.com/englianhu/stat_rethinking_2022)
3) [GitHub : `rethinking` package](https://github.com/englianhu/rethinking) McElreath 2020. Statistical Rethinking, 2nd edition
4) [GitHub : Stan](https://github.com/stan-dev)
5) [GitHub : Repository for distributing (some) stan-dev R packages](https://github.com/stan-dev/r-packages)
6) [书栈网](https://www.bookstack.cn/user/englianhu)
---
[<img src="文艺坊图库/Scibrokes.png" height="14"/> Sςιβrοκεrs Trαdιηg®](http://www.scibrokes.com)<br>
<span style='color:RoyalBlue'>**[<img src="文艺坊图库/Scibrokes.png" height="14"/> 世博量化®](http://www.scibrokes.com)企业知识产权及版权所有,盗版必究。**</span>