# Psilomethoxin, My Case Report (Addendum #2)

##### October 9, 2023

Purely based on subjective experience, I estimate the pharmacodynamics as:

```
pm <- function(x) 7.5*dchisq(4*(x-.5), 6)
```

The `x`

axis is hours and the `y`

axis is scaled so that 1 = a capsule of dried Pm mushroom powder.

```
plot(function(x) pm(x), 0,6, xlab="hours",ylab="capsules")
```

Using this function, you can estimate the effect of multiple capsules taken at different times. For example, suppose you take 2 at 10am, 1 at 10:45am, 1 at 11:15am, and 1 at 11:45am. If you space out 5 capsules like this then the peak is only about 3.5 capsules worth at 1pm:

```
plot(function(x) 2*pm(x) + pm(x-.75) + pm(x-1.25) + pm(x-1.75), 0,7, xlab="hours",ylab="capsules")
```

For background, see Psilomethoxin, My Case Report.