Python ret_tv_difmeanの例

コード例 #1

ファイル: dif.dens.py プロジェクト: nbykutsumi/dtanl

        #---------------
        a2dens_tmp = ma.masked_where( a2regionmask ==0.0, a2dens)
        a2dens_tmp = ma.masked_invalid(a2dens_tmp)
        #---------------
        # mean dens
        #---------------
        dlv[reg, era, iclass].append( a2dens_tmp.mean() )
      #-----------------
      # calc regional mean
      #-----------------
      dmv[reg, era, iclass]       = mean(dlv[reg, era, iclass])
      #-----------------
    #-----------------
    # calc t-value
    #-----------------
    dtv[reg, iclass] = func.ret_tv_difmean(dlv[reg, "his",iclass], dlv[reg, "fut",iclass])
    print iclass, dtv[reg, iclass]
    #-----------------
    # calc difference of regional mean
    #-----------------
    dfcm[reg, iclass] = (mean(dlv[reg, "fut", iclass]) - mean(dlv[reg, "his", iclass]))/mean(dlv[reg,"his",iclass]) * 100.0

#****************************************
# write csv
#----------------------------------------
csvdir = csvdir_root + "/%04d.%04d-%04d.%04d"%(diyear["his"], deyear["his"], diyear["fut"], deyear["fut"])

func.mk_dir(csvdir)

#------------
for reg in lreg:

コード例 #2

ファイル: dif.prec.stat.py プロジェクト: nbykutsumi/dtanl

         #
         #---------------
         # mean value
         #---------------
         dlp[reg, era, iclass].append( a2p_tmp.mean() )
         dlsp[reg, era, iclass].append( a2sp_tmp.mean() )
       #-----------------
       # calc regional mean
       #-----------------
       dp_reg[reg, era, iclass]       = mean(dlp[reg, era, iclass])
       dsp_reg[reg, era, iclass]       = mean(dlsp[reg, era, iclass])
       #-----------------
     #-----------------
     # calc t-value
     #-----------------
     dtv_p[reg, iclass] = func.ret_tv_difmean(dlp[reg, "fut",iclass], dlp[reg, "his",iclass])
     dtv_sp[reg, iclass] = func.ret_tv_difmean(dlsp[reg, "fut",iclass], dlsp[reg, "his",iclass])
     print iclass, dtv_p[reg, iclass]
     #-----------------
     # calc difference of regional mean
     #-----------------
     dfc_p[reg, iclass] = (mean(dlp[reg, "fut", iclass]) - mean(dlp[reg, "his", iclass]))/mean(dlp[reg,"his",iclass]) * 100.0
     dfc_sp[reg, iclass] = (mean(dlsp[reg, "fut", iclass]) - mean(dlsp[reg, "his", iclass]))/mean(dlsp[reg,"his",iclass]) * 100.0
 
 #****************************************
 # write csv
 #----------------------------------------
 csvdir = csvdir_root 
 
 func.mk_dir(csvdir)