doraneko94

import 'package:latlong2/latlong.dart';
import 'dart:math';

const f = 298.257222;
const r = 6378137.0;
const m0 = 0.9999;
const n1 = 1 / (2 * f - 1);
const alpha = [
  1.0,
  0.0008377318250120707,

doraneko94

doraneko94

// This code was inspired by https://github.com/rp-rs/rp-hal/blob/main/rp2040-hal/examples/watchdog.rs

#![no_std]
#![no_main]

use panic_halt as _;
use rp2040_hal as hal;
use hal::pac;

use embedded_hal::digital::v2::OutputPin;

doraneko94

from matplotlib import pyplot as plt
from matplotlib.patches import Ellipse
import numpy as np

fig, ax = plt.subplots(1)
ax.set_aspect('equal')

m = np.array([8, 12])
S = np.array([[16, np.sqrt(78)], [np.sqrt(78), 9]])
r = np.random.multivariate_normal(m, S, size=100) # random points to plot

doraneko94

import numpy as np
import pandas as pd
from scipy import stats

class TwoWayAnovaRM:
    
    def __init__(self, data: pd.DataFrame, subject, samples: list):
        
        self.data = data
        self.subject = subject

doraneko94

doraneko94

import numpy as np

class PolarDiff:

    def __init__(self, n, r2=1, x=None, clip=None, epsilon=1e-7):

        if x is not None:
            if n != len(x):
                raise ValueError
            self.theta = np.zeros(n-1)

doraneko94

# ２種類の気体分子が円形の容器内で動く様子のシミュレーションを行い、
# その様子をアニメーションとして保存する。
# 詳細は https://ushitora.net/archives/2197 を参照

import numpy as np
from matplotlib import pyplot as plt
from matplotlib import animation as animation

""" データの準備 """
# 乱数で２種類の気体の極座標生成

doraneko94

import numpy as np
from matplotlib import pyplot as plt
from scipy.stats import norm

def bic_raw(x, k):
    n = len(x)
    return -2 * np.sum(np.log(norm.pdf(x))) + k * np.log(n)

def bic_rss(x, k):
    n = len(x)

doraneko94

import numpy as np

def tau_num(X: np.ndarray, Y: np.ndarray):
    XpY = X + Y
    XmY = X - Y
    A = (XpY == 2).sum()
    S = (XpY == 0).sum()
    dX = (XmY == 1).sum()
    dY = (XmY == -1).sum()
    return A * S - dX * dY