Skip to content

Commit

Permalink
Merge pull request #106 from dominikl/add_examples_2
Browse files Browse the repository at this point in the history
Add formular for the model (yeast notebook)
  • Loading branch information
dominikl authored Nov 3, 2023
2 parents 478529c + d67b267 commit 7b5d64d
Showing 1 changed file with 4 additions and 2 deletions.
6 changes: 4 additions & 2 deletions examples/Yeast_Replication.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -79,8 +79,8 @@ countCells <- function (t) {
### Analyse all time points

```{r}
timepoints <- (1:tmax) # analyze every time point
#timepoints <- seq(1,tmax,7) # only use every seventh time point (faster)
#timepoints <- (1:tmax) # analyze every time point
timepoints <- seq(1,tmax,4) # only use every fourth time point (faster)
df <- data.frame(TimePoint=timepoints)
df$CellCount <- sapply(df$TimePoint, countCells)
df$TimePoint <- sapply(df$TimePoint, function(x) x * tInMinutes)
Expand All @@ -100,6 +100,8 @@ model <- lm(log(df$CellCount)~ df$TimePoint)
print(summary(model))
```

This gives us a model for the number of yeast cells N at time point t (in min): $N = e^{4.72710 + 0.00374 * t}$ .

### Lets see how well it fits visually

```{r}
Expand Down

0 comments on commit 7b5d64d

Please sign in to comment.