我正在使用R包“ pt”来计算累积前景理论值。
第一个输入如下:
choice_ids <- c(1, 1, 1, 1, 2, 2, 2, 2)
gamble_ids <- c(1, 1, 1, 2, 1, 1, 2, 2)
outcome_ids <- c(1, 2, 3, 1, 1, 2, 1, 2)
objective_consequences <- c(2500, 2400, 0, 2400,2500, 0, 2400, 0)
probability_strings <- c("0.33", "0.66", "0.01", "1.0","0.33", "0.67", "0.34", "0.66")
my_choices <- Choices(choice_ids=choice_ids,gamble_ids=gamble_ids,outcome_ids=outcome_ids,objective_consequences=objective_consequences,probability_strings=probability_strings)
之后
tk_1992_utility <- Utility(fun="power", par=c(alpha=0.88, beta=0.88, lambda=2.25))
linear_in_log_odds_prob_weight <- ProbWeight(fun="linear_in_log_odds", par=c(alpha=0.61, beta=0.724))
comparePT(my_choices,prob_weight_for_positive_outcomes=linear_in_log_odds_prob_weight,prob_weight_for_negative_outcomes=linear_in_log_odds_prob_weight,utility=tk_1992_utility, digits=4)
## cid gid ev pt ce rp
## 1 1 1 2409 881.3 2222 187
## 2 1 2 2400 943.2 2400 -0.000000000001819
## 3 2 1 825 312.6 684.2 140.8
## 4 2 2 816 307.2 670.9 145.1
comparePT命令将pt值作为输出,但还有很多其他值。但是,我只想将pt值作为输出,这有可能吗?不幸的是,我看了看包裹,但找不到里面的公式。
答案 0 :(得分:1)
似乎ProjectManagers
已从CRAN上移除,但可以从github档案库中安装:
pt
library(devtools)
install_github("cran/pt")
是S4函数。检查这些与常规S3类型有些不同。首先,使用comparePT()
查看可用的方法,然后再将showMethods()
用于感兴趣的方法。
getMethod()
但是,showMethods("comparePT")
# Function: comparePT (package pt)
# object="Choices"
getMethod("comparePT", "Choices")
# Method Definition:
#
# function ...
的输出只是一个常规的data.frame,因此您可以照常使用comparePT()
对其进行子集化。并将$
包裹起来,以其编码为字符。
as.numeric()