我目前正在对几项临床试验中的生存数据进行荟萃分析.
为此,我使用相同的方法从已发布的分析中获取代码.但是,使用已发布分析中的数据运行此代码时,我无法复制其结果.事实上,结果未能收敛到任何合理的估计.
代码本身(不包括数据)应该是正确的,因为它直接来自作者.我假设问题必须是w /初始值或如何运行采样的参数,但是在玩了很多初始值之后,燃烧,稀释等的长度......我没有得到有意义的结果.
我将不胜感激任何关于如何运行它(初始值等等)的建议,以使其正常运行.或者,如果代码中存在问题或者数据的设置方式与代码不匹配,那么知道这将是有用的.
作为旁注,我正在使用R2WinBUG进行分析,尽管我单独使用WinBUG遇到了同样的问题.
关于该方法的一些额外背景:
这种方法的工作方式是通过使用随机效应估计多个研究中处理之间重新参数化的Weibull分布的形状和尺度参数的差异.
威布尔分布被重新参数化,使得危险率的对数是a + b*log(t),其中a是比例参数,b是形状参数.由此,您可以计算一定间隔内给定数量的患者中给定数量的失败的似然函数.
不幸的是,这篇文章是公开的,但如果你可以访问这里是链接:http: //onlinelibrary.wiley.com/doi/10.1002/jrsm.25/abstract;jsessionid=2BA8F0D9BEF9A33F84975618D33F8DD9.f03t03?userIsAuthenticated=false&deniedAccessCustomisedMessage=
输入模型的变量的快速摘要:
NT:包括的单独治疗数量.
N:主数据集中的行数.NS:研究数量
s:研究数据线对应(编号为1:6)
r:该治疗/研究的间隔失败的患者数量
n:该治疗/研究的间隔开始时有风险的患者数
t:此数据对应的处理(编号1:3)
b:指示哪个治疗是与其他治疗比较的基线治疗(每条线设置为1).
bs:治疗是本研究的对照组
bt:治疗是本研究的研究部门
WinBUGS代码(包括数据):
#Winbugs code for random effects networks meta-analysis model
Model
{
for (i in 1:N)
{ # N=number of data points in dataset
#likelihood
r[i]~ dbin(p[i],n[i])
p[i]<-1-exp(-h[i]*dt[i]) # hazard h over interval [t,t+dt] # expressed as deaths per unit person-time (e.g. months)
#random effects model
log(h[i])<-nu[i]+log(time[i])*theta[i]
nu[i]<-mu[s[i],1]+delta[s[i],1]*(1-equals(t[i],b[i]))
theta[i]<-mu[s[i],2]+ delta[s[i],2]*(1-equals(t[i],b[i]))
}
for(k in 1 :NS)
{ …Run Code Online (Sandbox Code Playgroud) 我正在开发一个Angular2最终应用程序,它目前有两个模块:
的AppModule:
/**
* Created by jamdahl on 9/21/16.
*/
// Angular Imports
import {NgModule} from '@angular/core';
import {BrowserModule} from '@angular/platform-browser';
import {HttpModule} from '@angular/http';
import {FormsModule, ReactiveFormsModule} from '@angular/forms';
import {CoreModule} from '../core-module/core.module';
import {UserService, AuthService, AuthComponent} from '../core-module/core.module';
// Components
import {HomePageComponent} from './components/home-page.component';
//import {enableProdMode} from '@angular/core';
//enableProdMode();
@NgModule({
imports: [
BrowserModule,
HttpModule,
FormsModule,
ReactiveFormsModule,
CoreModule
],
declarations: [
AuthComponent,
HomePageComponent
],
providers: [
AuthService,
UserService
],
bootstrap: [
HomePageComponent
]
})
export class AppModule …Run Code Online (Sandbox Code Playgroud) 我想在使用d3和SVG制作的绘图中有阴影,但是阴影与相邻元素重叠会给我带来麻烦。有关当前外观,请参见下图。请注意,中间的六边形似乎具有变化的标高,因为在其中一些阴影上渲染了阴影。我想做的是设置阴影,使其仅在背景上渲染,而不在其他相邻的十六进制上渲染。
这是当前如何定义阴影的代码:
var filter = defs.append("filter")
.attr("id", "drop-shadow")
.attr("height", "130%");
// SourceAlpha refers to opacity of graphic that this filter will be applied to
// convolve that with a Gaussian with standard deviation 3 and store result
// in blur
filter.append("feGaussianBlur")
.attr("in", "SourceAlpha")
.attr("stdDeviation", 1)
.attr("result", "blur");
// translate output of Gaussian blur to the right and downwards with 2px
// store result in offsetBlur
filter.append("feOffset")
.attr("in", "blur")
.attr("dx", 1)
.attr("dy", 1)
.attr("result", "offsetBlur");
// overlay original SourceGraphic over translated …Run Code Online (Sandbox Code Playgroud) 我用下面的虚拟表制作了一个简单的闪亮仪表板。当我运行它并且窗口太薄时,表的列越过框的边缘。如果我把窗户弄得足够宽,那么桌子就可以装在盒子里了。
谁能提供更一致的格式帮助?解决方案可能是使框和表具有固定的宽度,使表列的宽度随框的宽度而变化,或者其他原因。
问题的底部是每个状态的屏幕快照,用于说明。
用户界面
library(shinydashboard)
header <- dashboardHeader(title="Dashboard")
sidebar <- dashboardSidebar(sidebarMenuOutput("menu"))
body <- dashboardBody(
tabItems(
tabItem(tabName = "Home",
fluidRow(
column(width = 4,
box(title = "Column 1", width = NULL, solidHeader = TRUE,
tableOutput("Table1")),
box(
title = "Title 1", width = NULL, solidHeader = TRUE,
"Box content")
),
column(width = 4,
box(
title = "Column 2", width = NULL, solidHeader = TRUE,
tableOutput("Table2"))
),
column(width = 4,
box(title = "Column 2", width = NULL, solidHeader = TRUE),
box(title = "Title 3", …Run Code Online (Sandbox Code Playgroud) css ×2
angular ×1
bayesian ×1
d3.js ×1
javascript ×1
r ×1
r2winbugs ×1
rstudio ×1
shiny ×1
svg ×1
typescript ×1
web-frontend ×1
winbugs ×1
winbugs14 ×1
xtable ×1