CamDavidsonPilon

# -*- coding: utf-8 -*-
from __future__ import annotations
import signal
from pioreactor.automations import LEDAutomationContrib
from pioreactor.automations import events
from pioreactor.actions.led_intensity import LED_Channel

__plugin_summary__ = "An LED automation for ON/OFF cycles"
__plugin_version__ = "0.0.2"
__plugin_name__ = "Light Dark Cycle LED automation"

CamDavidsonPilon

#!/bin/bash

sudo apt update
sudo apt install -y git

git clone https://github.com/Pioreactor/pioreactor.git --depth 1
cd pioreactor

make install-leader

CamDavidsonPilon

#!/bin/bash

sudo apt update
sudo apt install -y git

git clone https://github.com/Pioreactor/pioreactor.git --depth 1
cd pioreactor

make install-worker-from-args newHostname=$1

CamDavidsonPilon

#!/bin/bash

sudo apt update
sudo apt install -y git

git clone https://github.com/Pioreactor/pioreactor.git --depth 1
cd pioreactor

make install-leader-as-worker

CamDavidsonPilon

CamDavidsonPilon

[
    {
        "id": "328cb8d0.ad53d8",
        "type": "tab",
        "label": "Flow 1",
        "disabled": false,
        "info": ""
    },
    {
        "id": "e064c4c9.ed0b18",

CamDavidsonPilon

import time
import board
import adafruit_dht
from  paho.mqtt import publish

# Initial the dht device, with data pin connected to:
dhtDevice = adafruit_dht.DHT22(board.D21)


while True:

CamDavidsonPilon

CamDavidsonPilon

"""
Our model is

 log(cell/ml) = alpha * (secchi stick depth) + intercept + Noise

However, our cell/mL comes from a noisy hemocytometer, and our
depth measured with the secchi stick is also noisy. Our goal is
still to infer alpha & intercept as well as possible, so we can
dump the hemocytometer.

CamDavidsonPilon

# test previous algorithm

actuals = pd.read_csv("https://gist.githubusercontent.com/csaid/a57c4ebaa1c7b0671cdc9692638ea4c4/raw/ad1709938834d7bc88b62ff0763733502eb6a329/shower_problem_tau_samples.csv")


DELTA = 0.1

def survival_function(t, lambda_=50., rho=1.5):
    # Assume simple Weibull model
    return np.exp(-(t/lambda_) ** rho)