kuchaale

kuchaale

kuchaale

kuchaale

kuchaale

kuchaale

from wrf import getvar
from netCDF4 import Dataset
import xarray as xr
import pyproj

# Extract the variables of interest at time index 17
ds = Dataset('../wrfout_d02_2015-07-12_1200.nc')
variables = [getvar(ds, var, 17) for var in ('z', 'dbz', 'pressure', 'ter', 'ua',
                                             'va', 'wa', 'temp', 'rh')]
data = xr.merge(variables)

kuchaale

kuchaale

library(KFAS)
library(rstan)
data(GlobalTemp)

model_dlm1a <- stan_model("../stan/dlm1a.stan")

y <- as.matrix(GlobalTemp)
data <-
    within(list(),
       {

kuchaale

import numba
import numpy as np

@numba.jit(nopython=True, fastmath=True)
def move_regress(X, Y, w, fit_intercept=True):
    """Moving window multi-regressions
    
    Solves the least-squares problems `Y[t-w+1:t+1] = X[t-w+1:t+1] @ B[t] + A[t]`
    for `B` and `A` for all `t` in `[w-1, len(Y)-1]`.
 

kuchaale