edb_data_residual_analysis.go 22 KB

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