一亩三分地

 找回密码 注册账号

扫描二维码登录本站


北美版丁香园
美国和加拿大
疫情地图实时动态追踪

热门职场讲座
Career in Tech
职场晋升之路

Salarytics=Salary Analytics
查询工资数据
系统自动计算每年收入

科技公司如何
用数据分析驱动产品开发
coupon code 250off 立减$250
游戏初创公司招聘工程师、UIUX Designer和游戏策划
坐标湾区
DreamCraft创始团队
招聘游戏开发工程师
查看: 701|回复: 3
收起左侧

[Leetcode] R请教,R大牛们请进

[复制链接] |试试Instant~ |刷题, leetcode
论坛匿名账号
论坛匿名账号  发表于 2020-2-14 09:05:52 |阅读模式
本楼: 👍   0% (0)
 
 
0% (0)   👎

注册一亩三分地论坛,查看更多干货!

您需要 登录 才可以下载或查看,没有帐号?注册账号

x
哪位大神能帮我把这个run一下,用R

install.packages("DelayedEffect.Design")
library("DelayedEffect.Design")
pow.SEPPLE

这个是package 里面还有package,我看不到里面的了。帮我看看,谢谢大牛们!

上一篇:到哪里可以看看别人写的优秀的代码
下一篇:请教Restore IP Addresses题
我的人缘0
vcyy 2020-2-14 09:20:33 | 显示全部楼层
本楼: 👍   0% (0)
 
 
0% (0)   👎
全局: 👍   0% (0)
 
 
0% (0)    👎
Run完后里面显示是个function?
function (lambda1, t1, p, N, HR, tao, A, ap = 0.5, alpha = 0.05,
    nsim = 10000)  ...

你的问题是什么?
回复

使用道具 举报

我的人缘0
wx2123 2020-2-14 12:18:19 | 显示全部楼层
本楼: 👍   0% (0)
 
 
0% (0)   👎
全局: 👍   81% (18)
 
 
18% (4)    👎
function (lambda1, t1, p, N, HR, tao, A, ap = 0.5, alpha = 0.05,
    nsim = 10000)
{
    tmp <- get_lambda_t1_p(lambda1, t1, p)
    lambda1 <- tmp$lambda1
    t1 <- tmp$t1
    p <- tmp$p
    checkParms(1, t1, N, HR, tao, A, ap, alpha, nsim, NULL)
    Nupper <- round(N + qnorm(0.99) * sqrt(N))
    a <- N/A
    tvec <- c(0, t1)
    ratevec1 <- c(lambda1, HR * lambda1)
    ratevec2 <- c(lambda1, lambda1)
    p.val <- rep(NA, nsim)
    which <- 0
    for (i in 1:nsim) {
        t <- 0
        e <- rep(0, Nupper)
        tmp <- log(runif(Nupper))/a
        for (k in 1:Nupper) {
            t <- t - tmp[k]
            if (t > A)
                break
            e[k] <- t
        }
        ef <- e[e != 0]
        ns <- length(ef)
        if (!ns)
            next
        Z <- rbinom(ns, 1, ap)
        tmp <- Z == 1
        ef.trt <- ef[tmp]
        Z.trt <- Z[tmp]
        tmp <- !tmp
        ef.ctr <- ef[tmp]
        Z.ctr <- Z[tmp]
        n1 <- length(Z.trt)
        n0 <- length(Z.ctr)
        if (!n1 || !n0)
            next
        Tt <- rpexp(n1, rate = ratevec1, t = tvec)
        Tc <- rpexp(n0, rate = ratevec2, t = tvec)
        T <- c(Tt, Tc)
        eff <- c(ef.trt, ef.ctr)
        Z <- c(Z.trt, Z.ctr)
        tmp <- T > t1
        data2.1 <- T[tmp] - t1
        data2.2 <- tao - t1 - eff[tmp]
        trt <- Z[tmp]
        X <- pmin(data2.1, data2.2)
        evt <- as.numeric(data2.1 <= data2.2)
        p.val[i] <- getPval(which, X, evt, trt)
    }
    p.val <- p.val[is.finite(p.val)]
    m <- length(p.val)
    if (!m) {
        stop("ERROR: p-value could not be computed")
    }
    else if (m < nsim) {
        warning(paste("WARNING: only ", m, " simulations used in computing the p-value",
            sep = ""))
    }
    X <- mean(p.val <= alpha)
    X
}
<bytecode: 0x000001b645812748>
<environment: namespace:DelayedEffect.Design>
回复

使用道具 举报

我的人缘0
solarboy18 2020-2-15 04:17:08 | 显示全部楼层
本楼: 👍   0% (0)
 
 
0% (0)   👎
全局: 👍   50% (1)
 
 
50% (1)    👎
Thanks! I want to see the function of "getPval",can you help me?
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册账号

本版积分规则

隐私提醒:
■为防止被骚扰甚至人肉,不要公开留微信等联系方式,请以论坛私信方式发送。
■特定版块可以超级匿名:https://pay.1point3acres.com/tools/thread
■其他版块匿名方法:http://www.1point3acres.com/bbs/thread-405991-1-1.html

手机版|||一亩三分地

GMT+8, 2020-4-2 01:30

Powered by Discuz! X3

© 2001-2013 Comsenz Inc. Design By HUXTeam

快速回复 返回顶部 返回列表