私は生物学の大学院生で、特定の動作を R のモデルにコーディングしようとしていますが、「翻訳で失われた」問題がいくつかあります。私が持っているコードは投稿に従います。私はこのシステムをモデル化しようとしています:
大理石が振動するバスタブを想像してみてください。バスタブの表面は関数「bathtub」で与えられます。「温度」パラメータを指定して、浴槽のビー玉の振動をシミュレートします。低温/ゼロ温度では、底に置き、温度としてください。増加すると、側面をより高く探索する必要があります。これらをベクトルに保存し、この確率的プロセスの特定の実現がどのように見えるかを確認できるようにしたいと考えています。
私の問題は、大理石が動いている表面を記述する関数ではなく、開始値、時間、およびその他のパラメーター (曲率パラメーター a と b、温度パラメーター、およびおそらく他のパラメーター)。
私は基本的に正規分布のためにこのコードに類似したものを探しています:
bm <- function(x, x0, t, sigma) {
return(dnorm(x, mean=x0, sd=sigma*sqrt(t)))
}
アイデア、コード、または役立つリソースへのリンクをありがとうございます。
BATHTUB#
##Loads packages
library(ggplot2) #for graphing purposes
##Defines the composite distribution bathtub, which is composed of two Beta distributions
bathtub <- function(b, a){
dbeta(x, 1, b) + dbeta(x, a, 1)
}
##Initiates the parameters.
b = 5 #How sharp the higher bound is
a = 20 #How sharp the lower bound is
x <- seq(0, 1, length=101) #high density sample of (0,1)
##Plots
bathtub <- bathtub(a,b) #R likes certain kinds of variables.
qplot(x, bathtub, geom="line") #Plots smooth function pdf(bathtub)
qplot(x, 1-bathtub, geom="line")#this is the "likelihood" function