散点图线性拟合

everglow

为什么我的线性拟合y轴不对应呢？求问？

massage

看着像是两张图叠在一起了，用了overlay吗？

everglow

massage 发表于 2024-4-17 11:09
看着像是两张图叠在一起了，用了overlay吗？

是的，因为我要把线性拟合的线overlay在散点图上，不知道是拟合出了问题吗？为什么Y轴和散点图的Y轴不对应，数字也不对

zglnice

   resr@gsLineThicknessF= 3
   resr@gsLineDashPattern = 16
   rc_pos    = regline(xx_pos, yy_pos)          ; slope            
   df_pos    = rc_pos@nptxy-2
   prob_pos  = (1 - betainc(df_pos/(df_pos+rc_pos@tval^2), df_pos/2.0, 0.5) )

   xx1a = (/0,0.9/)
   yy1a = (/0*rc_pos+rc_pos@yintercept,0.9*rc_pos+rc_pos@yintercept/)

   rc_neg    = regline(xx_neg, yy_neg)          ; slope            
   df_neg    = rc_neg@nptxy-2
   prob_neg  = (1 - betainc(df_neg/(df_neg+rc_neg@tval^2), df_neg/2.0, 0.5) )


     xx2a = (/-0.9,0/)
     yy2a = (/-0.9*rc_neg+rc_neg@yintercept,0*rc_neg+rc_neg@yintercept/)

     if (prob_pos.ge.0.95) then
       resr@gsLineDashPattern = 0
     else
       resr@gsLineDashPattern = 16
     end if

    aadd_line_pos   = gsn_add_polyline(wks, plot1(1), xx1a, yy1a, resr)

     if (prob_neg.ge.0.95) then
       resr@gsLineDashPattern = 0
     else
       resr@gsLineDashPattern = 16
     end if
    aadd_line_neg   = gsn_add_polyline(wks, plot1(1), xx2a, yy2a, resr)

everglow

zglnice 发表于 2024-4-25 09:34
resr@gsLineThicknessF= 3
   resr@gsLineDashPattern = 16
   rc_pos    = regline(xx_pos, yy_pos) ...

谢谢你的回复，我去试试！

everglow

zglnice 发表于 2024-4-25 09:34
resr@gsLineThicknessF= 3
   resr@gsLineDashPattern = 16
   rc_pos    = regline(xx_pos, yy_pos) ...

你好，请问这是对两条线性拟合进行显著性检验吗

everglow

everglow 发表于 2024-6-1 15:04
你好，请问这是对两条线性拟合进行显著性检验吗

应该是的

everglow

zglnice 发表于 2024-4-25 09:34
resr@gsLineThicknessF= 3
   resr@gsLineDashPattern = 16
   rc_pos    = regline(xx_pos, yy_pos) ...

请问这个是过95%的显著性检验吗？那xx2a为什么（/-0.9,0/）这不太理解