小编sco*_*324的帖子

禁止在Rmarkdown HTML输出中从RStan运行消息

以下是我的代码.

```{r message=FALSE, warning=FALSE, cache=0,eval=TRUE, error=FALSE}
stan_m1 <- rethinking::map2stan(
  alist(
    y ~ dbinom(n, p),
    logit(p) <- alpha + bP*P + bA*A + bV*V,
    alpha ~ dnorm(0, 10),
    bP ~ dnorm(0, 5),
    bA ~ dnorm(0, 5),
    bV ~ dnorm(0, 5)
  ),
  data = d1, chains=2 , iter=2500 , warmup=500, debug=FALSE , 
  verbose=FALSE,refresh=-1
)
Run Code Online (Sandbox Code Playgroud)

knit我的rmarkdown文件作为html输出时,我得到以下输出.有没有办法关掉这些消息?

R-输出

r knitr r-markdown stan

8
推荐指数
1
解决办法
871
查看次数

威布尔分布的更新函数

Weibull 分布的更新函数m(t)如下t = 10所示。在此输入图像描述

我想找到 的值m(t)。我写了下面的r代码来计算m(t)

last_term = NULL
gamma_k = NULL
n = 50
for(k in 1:n){
  gamma_k[k] = gamma(2*k + 1)/factorial(k)
}

for(j in 1: (n-1)){
  prev = gamma_k[n-j]
  last_term[j] = gamma(2*j + 1)/factorial(j)*prev
}

final_term = NULL
find_value = function(n){
  for(i in 2:n){
  final_term[i] = gamma_k[i] - sum(last_term[1:(i-1)])
  }
  return(final_term)
}
all_k = find_value(n)

af_sum = NULL
m_t = function(t){
for(k in 1:n){
af_sum[k] = (-1)^(k-1) * all_k[k] …
Run Code Online (Sandbox Code Playgroud)

statistics r distribution gamma-function weibull

5
推荐指数
1
解决办法
343
查看次数

R中如何检查列表是否包含某个元素

我有以下清单

A = list(c(1,2,3,4), c(5,6,7,8), c(4,6,2,3,1), c(6,2,1,7,12, 15, 16, 10))
A
[[1]]
[1] 1 2 3 4

[[2]]
[1] 5 6 7 8

[[3]]
[1] 4 6 2 3 1

[[4]]
[1]  6  2  1  7 12 15 16 10
Run Code Online (Sandbox Code Playgroud)

我想检查该元素2是否存在于每个列表中。如果存在,那么我需要分配1给相应的列表。

先感谢您。

r list

4
推荐指数
1
解决办法
6952
查看次数

标签 统计

r ×3

distribution ×1

gamma-function ×1

knitr ×1

list ×1

r-markdown ×1

stan ×1

statistics ×1

weibull ×1