edb_data_residual_analysis.go 23 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698
  1. package models
  2. import (
  3. "encoding/json"
  4. "eta/eta_index_lib/global"
  5. "eta/eta_index_lib/utils"
  6. "fmt"
  7. "github.com/shopspring/decimal"
  8. "gorm.io/gorm"
  9. "math"
  10. "sort"
  11. "strconv"
  12. "time"
  13. )
  14. type edbDataResidualAnalysis struct {
  15. EdbDataId int `gorm:"column:edb_data_id;type:int(11);primaryKey;not null;"`
  16. EdbInfoId int `gorm:"column:edb_info_id;type:int(11);comment:指标id;default:NULL;"` // 指标id
  17. EdbCode string `gorm:"column:edb_code;type:varchar(190);comment:指标编码;default:NULL;"` // 指标编码
  18. DataTime string `gorm:"column:data_time;type:date;comment:数据日期;default:NULL;"` // 数据日期
  19. Value float64 `gorm:"column:value;type:double;comment:数据值;default:NULL;"` // 数据值
  20. CreateTime time.Time `gorm:"column:create_time;type:datetime;comment:创建时间;default:NULL;" ` // 创建时间
  21. ModifyTime time.Time `gorm:"column:modify_time;type:datetime;comment:修改时间;default:NULL;"` // 修改时间
  22. DataTimeStamp int64 `gorm:"column:data_timestamp;type:bigint(20);comment:数据日期时间戳;default:0;"` // 数据日期时间戳
  23. }
  24. // AfterFind 在该模型上设置钩子函数,把日期转成正确的string,所以查询函数只能用Find函数,First或者Scan是不会触发该函数的来获取数据
  25. func (m *edbDataResidualAnalysis) AfterFind(db *gorm.DB) (err error) {
  26. m.DataTime = utils.GormDateStrToDateStr(m.DataTime)
  27. return
  28. }
  29. // ConvertTimeStr
  30. // @Description: 转成需要输出的格式
  31. // @receiver m
  32. func (m *edbDataResidualAnalysis) ConvertTimeStr() {
  33. m.DataTime = utils.GormDateStrToDateStr(m.DataTime)
  34. return
  35. }
  36. // AddResidualAnalysisData 新增指标数据
  37. func AddResidualAnalysisData(dataList []edbDataResidualAnalysis) (err error) {
  38. err = global.DEFAULT_DB.CreateInBatches(dataList, utils.MultiAddNum).Error
  39. if err != nil {
  40. return err
  41. }
  42. return nil
  43. }
  44. // RefreshAllCalculateResidualAnalysis 刷新残差分析
  45. func RefreshAllCalculateResidualAnalysis(edbInfoId, source, subSource, formulaInt, moveType int, fromEdbInfo *EdbInfo, edbCode, startDate, endDate, moveFrequency string) (err error) {
  46. to := global.DEFAULT_DB
  47. /*to, err := o.Begin()
  48. if err != nil {
  49. return
  50. }
  51. defer func() {
  52. if err != nil {
  53. fmt.Println("RefreshAllCalculateResidualAnalysis,Err:" + err.Error())
  54. _ = to.Rollback()
  55. } else {
  56. _ = to.Commit()
  57. }
  58. }()*/
  59. //只清空残差分析数据自身指标的数据
  60. configMapping, err := GetConfigMappingListByConditionNotBase(edbInfoId)
  61. //if err != nil {
  62. // return err
  63. //}
  64. //
  65. //var edbInfoIdList []int64
  66. //for _, v := range configMapping {
  67. // edbInfoIdList = append(edbInfoIdList, v.EdbInfoId)
  68. //}
  69. //清空原有数据
  70. //sql := ` DELETE FROM edb_data_residual_analysis WHERE edb_info_id in (` + utils.GetOrmInReplace(len(edbInfoIdList)) + `) `
  71. //var params []interface{}
  72. //for _, i := range edbInfoIdList {
  73. // params = append(params, i)
  74. //}
  75. var params []interface{}
  76. sql := ` DELETE FROM edb_data_residual_analysis WHERE edb_info_id = ?`
  77. params = append(params, edbInfoId)
  78. //var params []interface{}
  79. //for _, i := range edbInfoIdList {
  80. // params = append(params, i)
  81. //}
  82. err = to.Exec(sql, params...).Error
  83. if err != nil {
  84. return
  85. }
  86. //计算数据
  87. err = refreshAllCalculateResidualAnalysis(edbInfoId, source, subSource, formulaInt, moveType, fromEdbInfo, edbCode, startDate, endDate, moveFrequency, configMapping)
  88. return
  89. }
  90. // refreshAllCalculateResidualAnalysis 刷新所有残差分析
  91. func refreshAllCalculateResidualAnalysis(edbInfoId, source, subSource, formulaInt, moveType int, fromEdbInfo *EdbInfo, edbCode, startDate, endDate, moveFrequency string, configMapping []CalculateResidualAnalysisConfigMapping) (err error) {
  92. fmt.Println("refreshAllCalculateResidualAnalysis startDate:", startDate)
  93. configMap := make(map[int64]CalculateResidualAnalysisConfigMapping, len(configMapping))
  94. for _, v := range configMapping {
  95. configMap[v.EdbInfoId] = v
  96. }
  97. calculateMappingList, err := GetCalculateMappingListByEdbInfoId(edbInfoId)
  98. if err != nil {
  99. return err
  100. }
  101. //这里顺序和之前的不一致导致的计算错误
  102. var edbInfoIdA, edbInfoIdB int
  103. if configMap[int64(calculateMappingList[0].FromEdbInfoId)].IndexType == 3 {
  104. edbInfoIdA = calculateMappingList[0].FromEdbInfoId
  105. edbInfoIdB = calculateMappingList[1].FromEdbInfoId
  106. } else {
  107. edbInfoIdA = calculateMappingList[1].FromEdbInfoId
  108. edbInfoIdB = calculateMappingList[0].FromEdbInfoId
  109. }
  110. mappingList, err := GetEdbInfoListByIds([]int{edbInfoIdA, edbInfoIdB})
  111. if err != nil {
  112. return err
  113. }
  114. var edbInfoMappingA, edbInfoMappingB *EdbInfoList
  115. for _, v := range mappingList {
  116. if v.Unit == "无" {
  117. v.Unit = ""
  118. }
  119. if v.EdbInfoId == edbInfoIdA {
  120. edbInfoMappingA = v
  121. }
  122. if v.EdbInfoId == edbInfoIdB {
  123. edbInfoMappingB = v
  124. }
  125. }
  126. if edbInfoMappingA == nil {
  127. return fmt.Errorf("指标A不存在")
  128. }
  129. if edbInfoMappingB == nil {
  130. return fmt.Errorf("指标B不存在")
  131. }
  132. // 从配置中取出时间范围
  133. analysisConfig, err := GetResidualAnalysisConfigById(edbInfoId)
  134. if err != nil {
  135. return err
  136. }
  137. configString := analysisConfig.Config
  138. var config ResidualAnalysisConfig
  139. err = json.Unmarshal([]byte(configString), &config)
  140. // 时间处理
  141. switch config.DateType {
  142. case -1:
  143. startDate = config.StartDate
  144. endDate = config.EndDate
  145. case 0:
  146. startDate = config.StartDate
  147. endDate = ""
  148. default:
  149. startDate = utils.GetPreYearTime(config.DateType)
  150. endDate = ""
  151. }
  152. // 原始图表信息
  153. originalEdbList := make([]EdbInfoList, 0)
  154. originalEdbList, fullADataList, fullBDataList, err := fillOriginalChart(config, fromEdbInfo, mappingList, startDate, endDate, edbInfoMappingA, edbInfoMappingB, originalEdbList)
  155. if err != nil {
  156. return err
  157. }
  158. indexADataMap := map[string]*EdbData{}
  159. for _, indexData := range edbInfoMappingA.DataList {
  160. indexADataMap[indexData.DataTime] = indexData
  161. }
  162. // 映射图表信息
  163. mappingEdbList, _, _, _, err := fillMappingChartInfo(config, fromEdbInfo, edbInfoMappingA, edbInfoMappingB, originalEdbList, indexADataMap, startDate, endDate, fullADataList, fullBDataList)
  164. if err != nil {
  165. return err
  166. }
  167. // 残差图表信息
  168. residualEdbList, _, err := fillResidualChartInfo(config, fromEdbInfo, edbInfoMappingA, edbInfoMappingB, mappingEdbList)
  169. // 映射指标 与 残差指标 同步刷新
  170. //for _, mapping := range configMapping {
  171. mapping := configMap[int64(edbInfoId)]
  172. var edbDataResidualAnalysisList []edbDataResidualAnalysis
  173. if mapping.IndexType == 1 {
  174. edbInfo, err := GetEdbInfoById(int(mapping.EdbInfoId))
  175. if err != nil {
  176. return err
  177. }
  178. for _, edbData := range mappingEdbList[0].DataList {
  179. value, _ := strconv.ParseFloat(edbData.Value, 64)
  180. edbDataResidualAnalysisList = append(edbDataResidualAnalysisList, edbDataResidualAnalysis{
  181. EdbInfoId: int(mapping.EdbInfoId),
  182. EdbCode: edbInfo.EdbCode,
  183. DataTime: edbData.DataTime,
  184. Value: value,
  185. CreateTime: time.Now(),
  186. ModifyTime: time.Now(),
  187. DataTimeStamp: edbData.DataTimestamp,
  188. })
  189. }
  190. err = AddResidualAnalysisData(edbDataResidualAnalysisList)
  191. if err != nil {
  192. return err
  193. }
  194. } else if mapping.IndexType == 2 {
  195. edbInfo, err := GetEdbInfoById(int(mapping.EdbInfoId))
  196. if err != nil {
  197. return err
  198. }
  199. for _, edbData := range residualEdbList[1].DataList {
  200. value, _ := strconv.ParseFloat(edbData.Value, 64)
  201. edbDataResidualAnalysisList = append(edbDataResidualAnalysisList, edbDataResidualAnalysis{
  202. EdbInfoId: int(mapping.EdbInfoId),
  203. EdbCode: edbInfo.EdbCode,
  204. DataTime: edbData.DataTime,
  205. Value: value,
  206. CreateTime: time.Now(),
  207. ModifyTime: time.Now(),
  208. DataTimeStamp: edbData.DataTimestamp,
  209. })
  210. }
  211. err = AddResidualAnalysisData(edbDataResidualAnalysisList)
  212. if err != nil {
  213. return err
  214. }
  215. }
  216. return
  217. }
  218. func fillResidualChartInfo(config ResidualAnalysisConfig, req *EdbInfo, edbInfoMappingA *EdbInfoList, edbInfoMappingB *EdbInfoList, mappingEdbList []EdbInfoList) ([]EdbInfoList, float64, error) {
  219. // 计算公式 映射残差 = 因变量指标 - 映射指标
  220. var edbInfoA, edbInfoB EdbInfoList
  221. if mappingEdbList[0].EdbInfoId == edbInfoMappingA.EdbInfoId {
  222. edbInfoA = mappingEdbList[0]
  223. edbInfoB = mappingEdbList[1]
  224. } else {
  225. edbInfoA = mappingEdbList[1]
  226. edbInfoB = mappingEdbList[0]
  227. }
  228. dataAList := edbInfoA.DataList
  229. edbData := make([]*EdbDataList, len(dataAList))
  230. for i, data := range dataAList {
  231. //f, _ := strconv.ParseFloat(data.Value, 64)
  232. f, _ := decimal.NewFromString(data.Value)
  233. edbData[i] = &EdbDataList{
  234. Value: f.InexactFloat64(),
  235. DataTimestamp: data.DataTimestamp,
  236. DataTime: data.DataTime,
  237. EdbInfoId: data.EdbInfoId,
  238. EdbDataId: data.EdbDataId,
  239. }
  240. }
  241. dataBList := edbInfoB.DataList
  242. // 映射指标开始时间
  243. var startTime string
  244. if len(dataBList) > 0 {
  245. startTime = dataBList[0].DataTime
  246. }
  247. var indexDataBMap = make(map[string]*EdbData)
  248. for _, data := range dataBList {
  249. indexDataBMap[data.DataTime] = data
  250. }
  251. // 求R2
  252. var valueB, sumValueA, averageValueA, residualQuadraticSum, totalQuadraticSum, R2 float64
  253. for _, indexData := range edbData {
  254. // 因变量的值总和
  255. sumValueA += indexData.Value
  256. }
  257. // 因变量平均值
  258. averageValueA = sumValueA / float64(len(edbData))
  259. var indexMax, indexMin float64
  260. var edbDataResp []*EdbDataList
  261. if len(edbData) > 0 {
  262. indexMax = edbData[0].Value
  263. indexMin = edbData[0].Value
  264. for _, indexData := range edbData {
  265. if dataB, ok := indexDataBMap[indexData.DataTime]; ok {
  266. f, _ := strconv.ParseFloat(dataB.Value, 64)
  267. valueB = f
  268. } else {
  269. continue
  270. }
  271. // 总因变量平方和
  272. totalQuadraticSum += math.Pow(indexData.Value-averageValueA, 2)
  273. // 补全残差值
  274. indexData.Value = math.Round((indexData.Value-valueB)*10000) / 10000
  275. // 残差平方和
  276. residualQuadraticSum += math.Pow(indexData.Value, 2)
  277. if indexData.Value > indexMax {
  278. indexMax = indexData.Value
  279. }
  280. if indexData.Value < indexMin {
  281. indexMin = indexData.Value
  282. }
  283. // 获取映射指标之后的数据
  284. if startTime != "" && utils.CompareDate(startTime, indexData.DataTime) {
  285. edbDataResp = append(edbDataResp, indexData)
  286. }
  287. }
  288. }
  289. // 计算R2 公式:R2=1-SSE/SST R2越大,越符合线性 R2 = 1 - 残差平方和/总平方和
  290. R2 = 1 - residualQuadraticSum/totalQuadraticSum
  291. mappingEdb := make([]EdbInfoList, len(mappingEdbList))
  292. copy(mappingEdb, mappingEdbList)
  293. for i, mapping := range mappingEdb {
  294. if mapping.EdbInfoId != edbInfoMappingA.EdbInfoId {
  295. toEdbData := convertEdbDataListToEdbData(edbDataResp)
  296. mappingEdb[i].DataList = toEdbData
  297. mappingEdb[i].EdbName = edbInfoMappingA.EdbName + "映射残差/" + edbInfoMappingB.EdbName
  298. if config.IndexType == 2 {
  299. if config.LeadValue > 0 {
  300. mappingEdb[i].EdbName = edbInfoMappingA.EdbName + "映射残差/" + edbInfoMappingB.EdbName + "(领先" + strconv.Itoa(config.LeadValue) + config.LeadFrequency + ")"
  301. }
  302. }
  303. }
  304. }
  305. //reverse(mappingEdb)
  306. return mappingEdb, R2, nil
  307. }
  308. // 双指针法
  309. func reverse(slice []EdbInfoList) {
  310. left, right := 0, len(slice)-1
  311. for left < right {
  312. // 交换元素
  313. slice[left], slice[right] = slice[right], slice[left]
  314. left++
  315. right--
  316. }
  317. }
  318. func fillMappingChartInfo(config ResidualAnalysisConfig, req *EdbInfo, edbInfoMappingA *EdbInfoList, edbInfoMappingB *EdbInfoList, originalEdbList []EdbInfoList, indexADataMap map[string]*EdbData, startDate string, endDate string, fullADataList []*EdbDataList, fullBDataList []*EdbDataList) ([]EdbInfoList, float64, float64, float64, error) {
  319. // 计算公式:Y=aX+b,Y为映射后的指标,X为自变量指标
  320. // 正序:a=(L2-L1)/(R2-R1) b=L2-R2*a
  321. // 逆序:a=(L2-L1)/(R1-R2) b=L2-R1*a
  322. // L2:左轴下限 R2:右轴上限 L1:左轴上限 R1:右轴下限
  323. var a, b, r float64
  324. // 映射残差 计算a,b
  325. if config.ResidualType == 1 {
  326. if config.IsOrder {
  327. a = (config.LeftIndexMax - config.LeftIndexMin) / (config.RightIndexMin - config.RightIndexMax)
  328. b = config.LeftIndexMax - config.RightIndexMin*a
  329. } else {
  330. a = (config.LeftIndexMax - config.LeftIndexMin) / (config.RightIndexMax - config.RightIndexMin)
  331. b = config.LeftIndexMax - config.RightIndexMax*a
  332. }
  333. }
  334. dataList := edbInfoMappingB.DataList
  335. // 指标B数据补充
  336. // 新建一个切片来保存补充的数据
  337. var replenishDataList []*EdbData
  338. for index := 0; index < len(dataList)-1; index++ {
  339. // 获取当前数据和下一个数据
  340. beforeIndexData := dataList[index]
  341. afterIndexData := dataList[index+1]
  342. // 从最早时间开始,补充时间为自然日
  343. for utils.IsMoreThanOneDay(beforeIndexData.DataTime, afterIndexData.DataTime) {
  344. // 创建补充数据
  345. nextDay := utils.GetNextDay(beforeIndexData.DataTime)
  346. toTime := utils.StringToTime(nextDay)
  347. replenishIndexData := EdbData{
  348. DataTime: nextDay, // 计算下一个自然日
  349. DataTimestamp: toTime.UnixMilli(),
  350. Value: beforeIndexData.Value, // 可以选择使用前一天的值,或者其他逻辑来计算值
  351. }
  352. // 将补充数据加入补充数据列表
  353. replenishDataList = append(replenishDataList, &replenishIndexData)
  354. // 更新 beforeIndexData 为新创建的补充数据
  355. beforeIndexData = &replenishIndexData
  356. }
  357. }
  358. // 将补充数据插入原始数据列表
  359. dataList = append(dataList, replenishDataList...)
  360. // 排序
  361. sort.Sort(ByDataTime(dataList))
  362. // 拟合残差 计算a,b
  363. var coordinateList []utils.Coordinate
  364. var replenishADataList []*EdbDataList
  365. var replenishBDataList []*EdbDataList
  366. if config.ResidualType == 2 {
  367. //
  368. // 因变量指标也转换为日度
  369. for index := 0; index < len(fullADataList)-1; index++ {
  370. // 获取当前数据和下一个数据
  371. beforeIndexData := fullADataList[index]
  372. afterIndexData := fullADataList[index+1]
  373. replenishADataList = append(replenishADataList, beforeIndexData)
  374. // 从最早时间开始,补充时间为自然日
  375. if utils.IsMoreThanOneDay(beforeIndexData.DataTime, afterIndexData.DataTime) {
  376. for {
  377. // 创建补充数据
  378. nextDay := utils.GetNextDay(beforeIndexData.DataTime)
  379. toTime := utils.StringToTime(nextDay)
  380. replenishIndexData := EdbDataList{
  381. DataTime: nextDay, // 计算下一个自然日
  382. DataTimestamp: toTime.UnixMilli(),
  383. Value: beforeIndexData.Value, // 可以选择使用前一天的值,或者其他逻辑来计算值
  384. }
  385. // 将补充数据加入补充数据列表
  386. replenishADataList = append(replenishADataList, &replenishIndexData)
  387. // 更新 beforeIndexData 为新创建的补充数据
  388. beforeIndexData = &replenishIndexData
  389. // 检查是否还需要继续补充数据
  390. if !utils.IsMoreThanOneDay(beforeIndexData.DataTime, afterIndexData.DataTime) {
  391. break
  392. }
  393. }
  394. }
  395. }
  396. replenishADataList = append(replenishADataList, fullADataList[len(fullADataList)-1])
  397. // 自变量指标也转换为日度
  398. for index := 0; index < len(fullBDataList)-1; index++ {
  399. // 获取当前数据和下一个数据
  400. beforeIndexData := fullBDataList[index]
  401. afterIndexData := fullBDataList[index+1]
  402. replenishBDataList = append(replenishBDataList, beforeIndexData)
  403. // 从最早时间开始,补充时间为自然日
  404. if utils.IsMoreThanOneDay(beforeIndexData.DataTime, afterIndexData.DataTime) {
  405. for {
  406. // 创建补充数据
  407. nextDay := utils.GetNextDay(beforeIndexData.DataTime)
  408. toTime := utils.StringToTime(nextDay)
  409. replenishIndexData := EdbDataList{
  410. DataTime: nextDay, // 计算下一个自然日
  411. DataTimestamp: toTime.UnixMilli(),
  412. Value: beforeIndexData.Value, // 可以选择使用前一天的值,或者其他逻辑来计算值
  413. }
  414. // 将补充数据加入补充数据列表
  415. replenishBDataList = append(replenishBDataList, &replenishIndexData)
  416. // 更新 beforeIndexData 为新创建的补充数据
  417. beforeIndexData = &replenishIndexData
  418. // 检查是否还需要继续补充数据
  419. if !utils.IsMoreThanOneDay(beforeIndexData.DataTime, afterIndexData.DataTime) {
  420. break
  421. }
  422. }
  423. }
  424. }
  425. replenishBDataList = append(replenishBDataList, fullBDataList[len(fullBDataList)-1])
  426. // replenishADataList --> map
  427. replenishADataMap := make(map[string]*EdbDataList)
  428. for _, indexData := range replenishADataList {
  429. if (utils.StringToTime(indexData.DataTime).After(utils.StringToTime(startDate)) || utils.StringToTime(indexData.DataTime).Equal(utils.StringToTime(startDate))) && (endDate == "" || utils.StringToTime(indexData.DataTime).Before(utils.StringToTime(endDate))) {
  430. replenishADataMap[indexData.DataTime] = indexData
  431. }
  432. }
  433. for _, indexData := range replenishBDataList {
  434. if _, ok := replenishADataMap[indexData.DataTime]; ok {
  435. coordinate := utils.Coordinate{
  436. X: indexData.Value,
  437. Y: replenishADataMap[indexData.DataTime].Value,
  438. }
  439. coordinateList = append(coordinateList, coordinate)
  440. }
  441. }
  442. a, b = utils.GetLinearResult(coordinateList)
  443. r = utils.ComputeCorrelation(coordinateList)
  444. }
  445. // 填充映射指标值 使得时间长度一致
  446. dataList = FillDataBList(dataList, edbInfoMappingA)
  447. // 根据指标A的时间key,在B的映射指标中筛选出对应的值
  448. var dataBList []*EdbDataList
  449. var indexMax, indexMin string
  450. if len(dataList) > 0 {
  451. indexMax = dataList[0].Value
  452. indexMin = dataList[0].Value
  453. for _, indexData := range dataList {
  454. if _, ok := indexADataMap[indexData.DataTime]; ok {
  455. indexDataCopy := *indexData
  456. // 计算指标B映射值
  457. f, _ := strconv.ParseFloat(indexData.Value, 64)
  458. indexDataCopy.Value = fmt.Sprintf("%f", math.Round((a*f+b)*10000)/10000)
  459. // 比较最大值
  460. if indexData.Value > indexMax {
  461. indexMax = indexData.Value
  462. }
  463. // 比较最小值
  464. if indexData.Value < indexMin {
  465. indexMin = indexData.Value
  466. }
  467. // 将副本添加到 dataBList
  468. copyValue, _ := strconv.ParseFloat(indexDataCopy.Value, 64)
  469. dataBList = append(dataBList, &EdbDataList{
  470. DataTime: indexDataCopy.DataTime,
  471. DataTimestamp: indexDataCopy.DataTimestamp,
  472. EdbDataId: indexDataCopy.EdbDataId,
  473. EdbInfoId: indexDataCopy.EdbInfoId,
  474. Value: copyValue,
  475. })
  476. }
  477. }
  478. }
  479. mappingEdbList := make([]EdbInfoList, len(originalEdbList))
  480. copy(mappingEdbList, originalEdbList)
  481. for i, mapping := range mappingEdbList {
  482. if mapping.EdbInfoId != edbInfoMappingA.EdbInfoId {
  483. mappingEdbList[i].EdbInfoId = 0
  484. mappingEdbList[i].EdbCode = ""
  485. mappingEdbList[i].EdbName = edbInfoMappingB.EdbName + "映射" + edbInfoMappingA.EdbName
  486. if config.IndexType == 2 {
  487. if config.LeadValue > 0 {
  488. mappingEdbList[i].EdbName = edbInfoMappingB.EdbName + "映射" + edbInfoMappingA.EdbName + "(领先" + strconv.Itoa(config.LeadValue) + config.LeadFrequency + ")"
  489. }
  490. }
  491. edbData := convertEdbDataListToEdbData(dataBList)
  492. mappingEdbList[i].DataList = edbData
  493. }
  494. }
  495. return mappingEdbList, a, b, r, nil
  496. }
  497. // FillDataBList 填充B的数据 使得与A的时间保持一致
  498. func FillDataBList(dataList []*EdbData, edbInfoMappingA *EdbInfoList) []*EdbData {
  499. dataAList := edbInfoMappingA.DataList
  500. if len(dataList) > 0 && len(dataAList) > 0 {
  501. for utils.StringToTime(dataList[len(dataList)-1].DataTime).Before(utils.StringToTime(dataAList[len(dataAList)-1].DataTime)) {
  502. // 使用A的时间填充时间差
  503. timeDiff := utils.GetNextDayN(dataList[len(dataList)-1].DataTime, 1)
  504. // 创建新的数据点并填充 前值填充
  505. newDataPoint := &EdbData{
  506. DataTime: timeDiff,
  507. Value: dataList[len(dataList)-1].Value,
  508. DataTimestamp: utils.StringToTime(timeDiff).UnixMilli(),
  509. }
  510. // 将新数据点添加到dataList末尾
  511. dataList = append(dataList, newDataPoint)
  512. }
  513. }
  514. return dataList
  515. }
  516. func fillOriginalChart(config ResidualAnalysisConfig, req *EdbInfo, mappingList []*EdbInfoList, startDate string, endDate string, edbInfoMappingA *EdbInfoList, edbInfoMappingB *EdbInfoList, originalEdbList []EdbInfoList) ([]EdbInfoList, []*EdbDataList, []*EdbDataList, error) {
  517. var fullADataList, fullBDataList []*EdbDataList
  518. for _, v := range mappingList {
  519. var edbInfoMapping EdbInfoList
  520. edbInfoMapping.EdbName = v.EdbName
  521. // 获取图表中的指标数据
  522. dataList, err := GetEdbDataList(v.Source, v.SubSource, v.EdbInfoId, startDate, endDate)
  523. if err != nil {
  524. return nil, nil, nil, fmt.Errorf("获取指标数据失败,Err:%s", err.Error())
  525. }
  526. data := convertEdbDataListToEdbData(dataList)
  527. // 重新获取指标数据 产品要求需要和计算指标-拟合残差逻辑保持一致
  528. fullDataList, err := GetEdbDataList(v.Source, v.SubSource, v.EdbInfoId, "", "")
  529. if err != nil {
  530. return nil, nil, nil, fmt.Errorf("获取指标数据失败,Err:%s", err.Error())
  531. }
  532. if v.EdbInfoId == edbInfoMappingB.EdbInfoId {
  533. // 领先指标 dataList进行数据处理
  534. if config.IndexType == 2 {
  535. if config.LeadValue < 0 {
  536. return nil, nil, nil, fmt.Errorf("领先值不能小于0")
  537. } else if config.LeadValue > 0 {
  538. edbInfoMapping.EdbName = v.EdbName + "(领先" + strconv.Itoa(config.LeadValue) + config.LeadFrequency + ")"
  539. for _, indexData := range dataList {
  540. switch config.LeadFrequency {
  541. case "天":
  542. indexData.DataTime = utils.GetNextDayN(indexData.DataTime, config.LeadValue)
  543. case "周":
  544. indexData.DataTime = utils.GetNextDayN(indexData.DataTime, config.LeadValue*7)
  545. case "月":
  546. indexData.DataTime = utils.TimeToString(utils.AddDate(utils.StringToTime(indexData.DataTime), 0, config.LeadValue), utils.YearMonthDay)
  547. case "季":
  548. indexData.DataTime = utils.TimeToString(utils.AddDate(utils.StringToTime(indexData.DataTime), 0, config.LeadValue*3), utils.YearMonthDay)
  549. case "年":
  550. indexData.DataTime = utils.TimeToString(utils.AddDate(utils.StringToTime(indexData.DataTime), config.LeadValue, 0), utils.YearMonthDay)
  551. }
  552. indexData.DataTimestamp = utils.StringToTime(indexData.DataTime).UnixMilli()
  553. }
  554. }
  555. }
  556. edbInfoMappingB.DataList = data
  557. fullBDataList = fullDataList
  558. } else {
  559. edbInfoMappingA.DataList = data
  560. fullADataList = fullDataList
  561. }
  562. edbInfoMapping.EdbInfoId = v.EdbInfoId
  563. edbInfoMapping.EdbCode = v.EdbCode
  564. edbInfoMapping.Unit = v.Unit
  565. edbInfoMapping.Frequency = v.Frequency
  566. edbInfoMapping.Source = v.Source
  567. edbInfoMapping.SourceName = v.SourceName
  568. edbInfoMapping.LatestDate = v.LatestDate
  569. edbInfoMapping.LatestValue = v.LatestValue
  570. edbInfoMapping.DataList = data
  571. originalEdbList = append(originalEdbList, edbInfoMapping)
  572. }
  573. return originalEdbList, fullADataList, fullBDataList, nil
  574. }
  575. func convertEdbDataListToEdbData(edbDataLists []*EdbDataList) []*EdbData {
  576. var edbDataList []*EdbData
  577. for _, edbData := range edbDataLists {
  578. data := &EdbData{
  579. DataTime: edbData.DataTime,
  580. DataTimestamp: edbData.DataTimestamp,
  581. EdbDataId: edbData.EdbDataId,
  582. EdbInfoId: edbData.EdbInfoId,
  583. Value: fmt.Sprintf("%f", edbData.Value),
  584. }
  585. edbDataList = append(edbDataList, data)
  586. }
  587. return edbDataList
  588. }
  589. type ByDataTime []*EdbData
  590. func (a ByDataTime) Len() int {
  591. return len(a)
  592. }
  593. func (a ByDataTime) Swap(i, j int) {
  594. a[i], a[j] = a[j], a[i]
  595. }
  596. func (a ByDataTime) Less(i, j int) bool {
  597. t1 := utils.StringToTime(a[i].DataTime)
  598. t2 := utils.StringToTime(a[j].DataTime)
  599. return t1.Before(t2)
  600. }