ddply的加权平均值错误(R，ddply)
ddply是R语言中的一个函数，用于对数据进行分组操作。它可以根据指定的变量对数据进行分组，并对每个分组应用指定的函数进行计算。
加权平均值错误是指在使用ddply函数计算加权平均值时出现的错误。加权平均值是一种统计方法，它考虑了每个数据点的权重，以更准确地计算平均值。然而，在使用ddply函数计算加权平均值时，可能会出现错误，这可能是由于数据中存在缺失值、权重值不合理或其他数据处理问题导致的。
为了解决加权平均值错误，可以采取以下几个步骤：
- 检查数据：首先，需要检查数据是否存在缺失值或异常值。可以使用R语言中的函数如is.na()、complete.cases()等来检查数据的完整性。
- 调整权重：如果发现权重值不合理或存在问题，可以对权重进行调整或重新计算。可以根据具体情况，使用合适的方法来计算权重，例如根据数据的重要性、样本的大小等进行调整。
- 数据处理：在使用ddply函数之前，可以对数据进行预处理，例如删除缺失值、处理异常值等。可以使用R语言中的函数如na.omit()、na.exclude()等来处理缺失值。
- 使用其他函数：如果ddply函数在计算加权平均值时仍然存在问题，可以尝试使用其他函数来进行计算。R语言中还有其他函数如aggregate()、tapply()等可以用于数据分组和计算。
在腾讯云的产品中，与数据处理和分析相关的产品有腾讯云数据仓库（Tencent Cloud Data Warehouse，CDW）、腾讯云数据湖（Tencent Cloud Data Lake，CDL）等。这些产品提供了强大的数据处理和分析能力，可以帮助用户高效地进行数据处理和计算。
更多关于腾讯云数据仓库的信息，请访问：腾讯云数据仓库产品介绍
更多关于腾讯云数据湖的信息，请访问：腾讯云数据湖产品介绍
相关·内容
我有一个数据框，看起来像这样：
site date var dil
1 A 7.4 2
2 A 6.5 2
1 A 7.3 3
2 A 7.3 3
1 B 7.1 1
2 B 7.7 2
1 B 7.7 3
2 B 7.4 3
我需要向此数据框添加一个名为wt的列，该列包含计算加权平均值所需的加权因子。对于site和date的每个组合，都必须推导出该加权因子。
我使用的方法是首先构建一个计算权重因子的函数：
浏览 1提问于2011-09-28得票数 11
回答已采纳
1回答
我想对一些时间序列数据应用分析权重，但不知道如何在R中实现这一点。我正在转录一些Stata代码，代码使用collapse和[aweight='weightVar']。
Stata码
collapse temp [aweight='weightVar], by(year);
如何将分析权重应用于数据，使用下面的croparea作为对每年每个id的temp的加权变量？
样本数据
df <- structure(list(id = c(1, 1, 1, 1, 2, 2, 2, 2), year = c(1900,
1900, 1900, 1900, 1901, 190
浏览 4提问于2016-01-03得票数 2
回答已采纳
我正尝试在一个名为DFA_CMO的表上计算加权平均值(每视图数的AverageViewingTime)的聚合。DFA_CMO有5个维度，Site就是其中之一。
代码如下：
ddply(DFA_CMO,.(Site),summarize, wm = weighted.mean(DFA_CMO$AverageViewingTime, DFA_CMO$Views, ,na.rm=TRUE))
但结果令人失望，因为它为每个字段显示了相同的值：
Site wm
1 Advideum 21.17633
2 bbc.com 21.17
浏览 2提问于2013-06-21得票数 1
1回答
在寻找解决方案和摆弄之后，我正在寻求帮助，试图在箱线图上显示加权平均值(我也曾尝试将其交叉发布到ggplot2邮件列表中)。
下面我提供一个玩具示例。
#data
value <- c(5, 7, 8, 6, 7, 9, 10, 6, 7, 10)
category <- c("one", "one", "one", "two", "two", "two",
"three", "three", "three",
浏览 0提问于2013-04-25得票数 3
回答已采纳
1回答
带加权中间值的枢轴表
、、、
这个任务看起来很简单，但我仍然找不到答案。我在枢轴表中看到了加权平均值(即均值)的解决方案，但没有看到加权中间值。
我需要创建一个带有加权中间值的枢轴表(使用两个列，一个列具有值，另一个列具有权重--每个值的数量要按照中值计算)。换句话说，…
city value weight
Phoenix 100 40
Phoenix 140 12
Phoenix 500 15
San Francisco 700 5
San Francisco 300 5
San Francisco
浏览 0提问于2016-02-29得票数 2
回答已采纳
2回答
我想知道是否有任何方法可以跨矩阵向量化加权算术。例如，用于计算给定weights以下的values的加权平均值 set.seed(12321)
values = matrix(sample(10:40, 6, replace = TRUE), ncol = 2)
weights = matrix(sample(1:5, 6, replace = TRUE), ncol = 2) 我们可以这样做： w.mean = numeric()
for (i in 1:dim(values)[1]){
w.mean[i] = stats::weighted.mean(x=values[i,], w
浏览 41提问于2019-02-04得票数 2
我想在回归表的底部添加一行，用于列出因变量的加权平均值。通常，我会跑
reg y x1 x2 x3
estadd ysumm, mean
eststo r1
esttab r1 using results.tex, replace label title("Title") long nomtitles cells("b(fmt(a3) star)" t(par fmt(2))) stats(r2 N ymean, labels("R-squared" "Observations" "Mean of Y"))
然而，
浏览 10提问于2013-09-17得票数 1
我能够成功地在因子水平内获得给定向量的简单平均值，但在尝试将其带到下一步加权观察值时，我无法使其工作。这是可行的：
> tapply(exp.f,part.f.p.d,mean)
1 2 3 4 5 6 7 8 9 10
0.8535996 1.1256058 0.6968142 1.4346451 0.8136110 1.2006801 1.6112160 1.9168835 1.5135006 3.0312460
但这
浏览 2提问于2011-02-02得票数 4
回答已采纳
1回答
我有一个数据框，如下所示(8列- myPOSIX列采用'y-m-d h:s‘格式)
head(new)
Date.and.Time..UTC. Receiver Transmitter Latitude Longitude ndiffs29912 flag
1 07/10/2010 15:53 VR2W-107619 A69-1303-29912 48.56225 -53.89144 NA FALSE 2010-10-07
2 07/10/2010 15:56 VR2W-107619 A69-1303-29912
浏览 1提问于2012-02-20得票数 2
回答已采纳
我正在用数据帧InsectSprays (在R库中)练习R，使用plyr包：
ddply(InsectSprays,.(spray), summarize, sum = sum(count))
Error in .fun(piece, ...) : argument "by" is missing, with no default
我不明白这个错误是什么意思，教程执行了完全相同的命令，输出如下：
spray sum
1 A 174
2 B 184
3 C 25
我还想创建一个新的变量(仍然是错误的)：
sp
浏览 0提问于2018-11-27得票数 0
2回答
我对plyr还不熟悉，我想取类内值的加权平均值来重新构造多变量的数据。使用以下代码，我知道如何对一个变量(如x2)执行此操作：
set.seed(123)
frame <- data.frame(class=sample(LETTERS[1:5], replace = TRUE),
x=rnorm(20), x2 = rnorm(20), weights=rnorm(20))
ddply(frame, .(class),function(x) data.frame(weighted.mean(x$x2, x$weights)))
但是
浏览 3提问于2013-08-22得票数 2
回答已采纳
我试图显示不同因素水平(样本)和加权平均值(weight=coverage)的变量(等位基因特定表达式)。
我做了一些样本数据：
set.seed(2)
x <- sample(c("A","B","C"), 100, replace=T)
y <- rnorm(100)
w <- ceiling(rnorm(100,200,200))
df <- data.frame(x, y, w)
library(ggplot2)
ggplot(df, aes(x=factor(x), y=y, weight=w)) +
ge
浏览 0提问于2013-04-05得票数 1
回答已采纳
2回答
我在使用ddply生成表时遇到了困难：
ddply(Data, .(Data$Condition), summarise, n = length(Data$Condition),
mean = mean(Data$Composite), std = sd(Data$Composite), var = var(Data$Composite))
出现以下错误：
Error in `[.data.frame`(col, i) : undefined columns selected
但是，我试图在ddply中使用的df中的列是明确定义的。我也在R里测试过这个，很明显是R：
可以检索它们
浏览 1提问于2021-04-13得票数 0
2回答
总结和lm模型误差
、、、
在尝试创建统计表时，我遇到了summary.lm r.squared值的问题。
首先，我从csv中读取数据。
df <- as.data.frame(read.csv("BCO.csv", header = TRUE, stringsAsFactors = FALSE))
df <- df[,2:4]
然后我开始观察趋势..。
CLDD_trend <- ddply(df, .(STATION_NAME), function(z)coef(lm(CLDD_yr ~ year, data = z)))
这是我遇到问题的地方..。
CLDD_rsq <- d
浏览 3提问于2016-04-18得票数 1
回答已采纳
1回答
我试图使用approx()和dplyr在现有数组中插值值。我的初始代码是这样的..。
p = c(1,1,1,2,2,2)
q = c(1,2,3,1,2,3)
r = c(1,2,3,4,5,6)
Inputs<- data.frame(p,q,r)
new.inputs= as.numeric(c(1.5,2.5))
library(dplyr)
Interpolated <- Inputs %>%
group_by(p) %>%
arrange(p, q) %>%
mutate(new.output=a
浏览 0提问于2016-09-27得票数 1
回答已采纳
4回答
我是一个非常新的R，并正在与数据集的支付细节从我运行的网站。我有一个数据框架，其中包括以下列：
¦ Date ¦ Amount ¦ Type
每一行都是金额不同的单个付款。我想找一种方法，不循环通过25,000左右的付款，在数据集中的每一天的总支付金额。
我尝试过sum(sub$Amount[sub$Date == unique(sub$Date)])，但这会引发以下错误：
Warning message:
In sub$Date == unique(sub$Date) :
longer object length is not a multiple of shorter object
浏览 0提问于2014-09-04得票数 1
回答已采纳
3回答
我正在尝试弄清楚为什么我在使用ddply时收到一条错误消息。
示例数据：
data<-data.frame(area=rep(c("VA","OC","ES"),each=4),
sex=rep(c("Male","Female"),each=2,times=3),
year=rep(c(2009,2010),times=6),
bin=c(110,120,125,125,110,130,125,80,90,90,80,140),
shell_length=c(.4,4,1,2
浏览 1提问于2015-11-19得票数 14
下面的几乎最小代码的要点是在函数f中应用ddply，其中ddply的聚合函数(助手)是在f的定义中自定义的。
不幸的是，我不明白为什么整个代码片段的获取会在eval(expr，envir，expr)中产生一个错误:无法找到函数“助手”。当助手函数独立于函数f运行时，代码工作。当我用的未注释调用替换为时，代码将无错误地运行。
您能解释错误并提供解决方案或解决方法吗？用plyr 1.8.1和R 3.0.3测试
rm (list = ls())
library(plyr)
f <- function() {
dfx <- data.frame(
group = c(rep
浏览 2提问于2014-04-22得票数 0
回答已采纳
1回答
我有一个data.frame，它看起来像这样(实际上是100万行)：
`> df
R.DMA.NAMES quarter daypart allpersons.imp rate station spot.id
1 Wilkes.Barre.Scranton.Hztn Q22014 afternoon 0.0 30 WSWB 13048713
2 Nashville Q12014 primetime
浏览 4提问于2014-11-08得票数 0
回答已采纳
我试图从一个基于同一DateTime组的贸易数据中计算数量加权平均价格(VWAP)。来自一个小数据框架(20,000项)的样本数据如下：
testdata[,c(5,8,10,11)]
transactiontime price volume totalEquity
334 2014-12-01 01:30:00.110000 19330 500000 966500
335 2014-12-01 01:30:00.830000 19340 8000 15472
336 2014-12-01 01:30:00.830000 19340
浏览 1提问于2015-08-17得票数 0
回答已采纳
1回答
我想在用ggplot2创建的堆叠面积图上添加错误条。
我的csv文件如下所示：
Day Cat Val Error
0 A 0 0.00
0 B 44.77 1.16
0 C 54.64 0.88
13 A 1.34 0.32
13 B 22.78 0.45
13 C 38.33 2.12
19 A 1.95 0.35
19 B 24.00 2.25
19 C 40.30 3.86
我试过这个：
ggplot(data=mydata, aes(
浏览 3提问于2013-04-17得票数 2
回答已采纳
1回答
我想取一个组的加权平均值，它使用r中的聚合函数。
以下是我的数据：
set.seed(1980)
Group_1 <- sample(letters[1:4], 50, TRUE)
Group_2 <- sample(letters[8:13], 50, TRUE)
Weight <- sample(seq(1,50), 50, TRUE)
Value <- sample(seq(1,50), 50, TRUE)
DF <- data.frame(Group_1, Group_2, Weight, Value)
head(DF)
我想取Value列的加权平均值
浏览 4提问于2020-10-18得票数 3
回答已采纳
我有一个类似的问题，这里已经讨论过了：
在我的例子中，我在自己的包中使用了ddply函数，当我检查R-package时，我得到了一个提示，“全局变量'VARX‘没有可见的绑定。”
地址代码行是：
subsample <- ddply(my_data, .(VARX), function(x){x[sample(nrow(x), 1), ]})
变量VARX是包(sysdata)提供的数据帧中的一列。什么是改进代码以获得合适R-package的好实践？
浏览 1提问于2013-03-05得票数 0
回答已采纳
3回答
在函数中循环以填充数据帧
、、、
我有三个气象站，为每个已注册的变量创建一个气象数据寄存器和数据帧。例如：
HR是包含三个站点的已登记相对湿度数据的数据帧，而feed是每个变量应该保存我加权平均值的数据帧。在这种情况下，我只对相对湿度应用这个函数来填充feed数据帧中的第二列。
HR<-data.frame(M5125=c(70,75,NA,NA),
M5126=c(73,NA,NA,71.5),
PRAA=c(NA,NA,NA,NA))
precip<-data.frame(M5125=c(16,NA,11,NA),
M5
浏览 3提问于2020-05-19得票数 0
回答已采纳
1回答
:)我想问你一些关于R的plot的事情，如果有人能帮我，我会很高兴的！
我写了一段Heston Model in R代码。它产生了一个期权价格的向量，让我们称之为H1。
每个期权价格(向量H1的每个元素)对应于一天。向量非常长( 3556个元素)。我想把它画出来，并对图表进行分析。我使用了函数plot(.....)。然后我希望在axis x上有dates，在axis y上有我可选的prices。因此，我使用函数axis(1，z) (其中z是包含所有3556日期的向量)和axis(2，H1) (其中H1包含所有3556期权价格)。
关键是所有的日期和期权价格都包含在我的图中:/它看起来非常糟糕，
浏览 0提问于2014-01-16得票数 0
1回答
应用函数而不是使用循环
、、、、
我已经在文档和论坛中搜索了很长一段时间，但是对于更复杂的函数，我仍然很难理解如何使用apply函数而不是R中的循环。(对于像apply(data，1，sum)这样的函数，它还好)
例如，我有以下函数
AOV_GxT=function(Trial_group,trait,df){
sub_table=df[which(df$Trial.group == Trial_group),]
aov_GxT = anova(aov(sub_table[,trait] ~ Genotype + Treatment + Treatment/Rep.number + Ge
浏览 3提问于2013-12-10得票数 1
回答已采纳
当比较两个整数和用户定义函数中的数据子集时，我在R中得到了一个“更长的对象长度而不是更短的对象长度的倍数”警告。
用户定义的函数只返回从dataframe获取的整数子集的中位数：
function(s){
return(median((subset(EDB,as.integer(validSession) == as.integer(s)))$absStudentDeviation))
}
(我最初没有as.integer矫顽器在那里。我把它们放在那里是为了调试，文本，而且我仍然收到一个错误。)
我遇到的具体错误是：
In as.integer(validSession) == a
浏览 1提问于2011-12-18得票数 1
1回答
分别处理R中变量的不同值
、、、
我在R上相当业余，所以我希望它听起来不会很傻。这是:我有一个我正在处理的数据集。其中一个变量(x1)有三个分类值(countrya，countryb，countryc)。数据集有许多变量和观察结果，但我想为每个国家单独分析它。我应该准备一个数据帧吗?我如何用给定的信息来做这件事？假设数据集称为data；变量称为x1；此变量采用的值是countrya、countryb和countryc，我想分别分析它们。我希望这对代码有帮助。谢谢..。
浏览 1提问于2012-04-10得票数 0
4回答
我有一段聚合代码，它可以很好地工作，但是对于一个有10e6行的数据帧，运行速度有点慢。我在R方面并不是很有经验，所以请为我那令人畏缩的代码道歉！
我只想对一个公共键的值做一个基本的汇总和求和...
例如，从……
key val
1 a 5
2 b 7
3 a 6
为了..。
key val
1 a 11
2 b 7
我能做的最多就是..。
keys = unique(inp$key)
vals = sapply(keys, function(x) { sum(inp[inp$key==x,]$val) })
out = data.frame(ke
浏览 0提问于2011-07-25得票数 9
回答已采纳
1回答
我想我需要你的帮助。
我有一张这样的桌子：
Product Price Quantity Quality
1 10 100 100
2 15 1200 88
1 12 88 95
现在，我想要计算依赖于数量的每个产品的平均值。我在Excel中可以做到这一点，但是我不知道我应该如何在R中解决这个问题。
如果我使用这样的东西：
weighted.mean(Test_weighted_avg$Price, Test_weighted_avg$Qua
浏览 3提问于2018-02-18得票数 1
回答已采纳
1回答
我试图在R中按组计算加权平均值，但它只是返回整个数据集的加权平均值，并且无法确定我的问题发生在哪里。下面是我的密码。注意，在weighted.mean函数中，如果不为列类型指定数据帧名称，则不会返回任何内容，因此不确定引用数据的方式是否导致了问题。
unit_averages = selected_units %>%
group_by(`Length x Width`,Date) %>%
summarise(index_mean = weighted.mean(selected_units$"Wtd Avg Price",w=selected_units
浏览 3提问于2022-06-10得票数 0
回答已采纳
1回答
在每个非NA之前分组拆分数据
、、、、
我希望根据列“高度”将我的数据分割成子集，每个子集都有一个带有值的行和带有NAs的0-Inf行。
这是为了以后能够将函数应用到子集，特别是根据它们的“直径”值、每个子集中的行数、sqrt(和(直径^2))等对行进行排序。
1.
我的数据看起来是这样的：
> dput(df[1:300, c("IDD", "Height", "Diameter")])
structure(list(IDD = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,
13, 14, 15, 16, 17, 18, 1
浏览 1提问于2015-03-13得票数 0
回答已采纳
我在这个问题上挣扎了好几个小时，这似乎适合于plyr或*申请。有人能给我指出一个比我下面列出的方案更不笨重的R方案吗？
Bkdg:我看了很多与R金融相关的套餐，我找不到一个流行的软件包，它能很好地处理在时间序列中间消失的证券，以及按计划变化的权重。我正在为这个特定的问题构建自己的解决方案，但我宁愿使用现有的解决方案。
问题的症结在于，我想使用plyr来“循环”一个日期子集内的证券列表。一些证券在日期范围内消失。(我使用的是来自数据的前向rtns，tha没有生存偏差。)我希望每个日期范围的输出是所选证券累计回报的数据框架。我可以使用这个(连同初始权重)与其他日期范围相结合来计算各种组合度量。
我
浏览 3提问于2013-10-01得票数 1
回答已采纳
我在R中使用下面的函数将受试者/样本分割成训练和测试集，它工作得很好。然而，在我的数据集中，受试者被分为两组(病人和对照组)，因此，我希望将数据分成两组，同时保持每个培训和测试集中患者和对照组所占比例与完整数据集的比例相同。我怎样才能在R中做到这一点？如何修改以下功能，使其在将数据拆分为培训和测试集时考虑到组间的隶属关系？
# splitdf function will return a list of training and testing sets#
splitdf <- function(dataframe, seed=NULL) {
if (!is.null(seed))
浏览 3提问于2013-09-22得票数 1
回答已采纳
2回答
R的ave()函数比它的名字要有用得多-它基本上是tapply()的一个版本，它允许您返回与输入长度相同的向量，并将这些值重新排列为与输入相同的顺序。
> x <- 1:10
> ave(x, x %% 2, FUN=function(d) d-mean(d))
[1] -4 -4 -2 -2 0 0 2 2 4 4
您可以使用ddply()实现类似的效果，但它需要两份额外的数据副本，以及几个辅助变量：
> x <- 1:10
> val <- ddply(data.frame(x=x, id=1:10), .(x %% 2),
浏览 0提问于2011-12-28得票数 5
回答已采纳
我已经使用lm()对R中的数据集执行了多元回归分析，并且我能够使用下面的函数提取一年中每一天的系数。我还想提取一年中每一天的R2，但这似乎不是以相同的方式工作的。 这个问题和：Print R-squared for all of the models fit with lmList几乎是一样的，但是当我尝试这个的时候，我得到了“错误：$ operator对于原子向量是无效的”。如果可能的话，我也想把它包含在同一个函数中。如何以这种方式提取每个doy的R2？ #Create MR function for extracting coefficients
getCoef <- functi
浏览 47提问于2019-07-01得票数 1
回答已采纳
1回答
在R中处理长格式数据
、、、
我创建了一个长格式的数据框，它有5个变量(公司名称，关键项目，原产国，年份和价值)我想做的是创建一个表，其中包含每年一个关键项目的值的平均值。到目前为止我做到了以下几点：
for(i in dat$Year) {
dat[dat$KeyItem == "Total Assets" & dat$Year == i, ]
FUN1 <- mean(dat$Value, na.rm = TRUE)
}
我在R方面几乎没有任何经验，如果这个问题太简单了，我很抱歉。提前感谢你们的帮助。
浏览 4提问于2012-10-25得票数 0
回答已采纳
1回答
我是R的新手，我正在学习如何编写循环。我想知道我尝试用这个循环做的事情是否有意义：
require(plyr)
year = c("1993", "2009")
for (i in 1:2) {
temp <- subset(census.long,trt=="RIL" & dbh > 50 & date==year[i])
tab.res[year[i]] <- ddply(.data=temp, .(spcode),
.fun=summarise,
浏览 1提问于2013-06-12得票数 0
1回答
我是R的新手，所以就当这篇文章很有价值吧。
我已经写了一个函数，它接受4个参数并返回一个数据帧。压缩版本如下所示。
Advantage <- function(tRos, tTat, cRos, cTat){
#case 1.1 tRos is lower
if((tRos > 0 | cRos > 0) & cRos > tRos & cTat < tTat){
tRosAd <- (cRos - tRos) * cTat * -1
tTatAd <- (tTat - cTat) * tRos
r &l
浏览 1提问于2013-07-15得票数 0
回答已采纳
2回答
动物园对象时间序列聚合
、、、
我有一个R动物园的对象。动物园对象(z)按日期编制索引，并有多个列：
V1 (聚合值是“所选”行中所有值的总和)V2(聚合值是“所选”行中所有值的q1第一四分位数)V3(聚合值是“所选”行中所有值的最小值)V4(聚合值是“所选”行中所有值的第一个值)v5(聚合值是所有值中所有值的最后一个值)<code>H 210</code></code>F 211</code>
我希望以不同的方式(即使用不同的函数)聚合每个“列”中的数据，但是在相同的行数上进行聚合。
我希望使用一个函数进行聚合，该函数允许我指定要聚合的行数。例如：
my_aggregate
浏览 1提问于2012-02-05得票数 1
一些sql数据库将with cube修饰符到group by操作。我的没有这个功能。
基本上，如果我有一个数据集，比如：
+------+-----------+---------+---------+
| sum | source_id | type_id | variety |
+------+-----------+---------+---------+
| 491 | 1 | 1 | 1 |
| 2008 | 1 | 2 | 1 |
| 33 | 1 | 3 |
浏览 0提问于2011-04-20得票数 2
回答已采纳
我正在尝试计算一些简单的比率，并使用R的括号表示比率的基线。
现在，我正在努力定义一个允许我参数化基线的函数。我不想硬编码它，因为我有几个这样的代码。而且，我真的不明白R到底在做什么，我非常好奇如何实现期望的行为。
下面是一些基于示例数据的代码：
data("singer", package = "lattice")
# this is what I want, but what currently doesn't work
my_ratio <- function(voice) {
ddply(singer, ~ voice.part,
浏览 0提问于2012-09-01得票数 2
回答已采纳
1回答
请解释在培训Spacy NER模型时列的含义：
E # LOSS TOK2VEC LOSS NER ENTS_F ENTS_P ENTS_R SCORE
--- ------ ------------ -------- ------ ------ ------ ------
0 0 0.00 78.11 26.82 22.88 32.41 0.27
26 200 82.40 3935.97 94.44 94.44 94.44 0.94
59
浏览 26提问于2022-09-14得票数 0
回答已采纳
我正在尝试在R中做一个"group by“风格的加权均值。使用一些基本的均值，下面的代码(使用Hadley的plyr包)工作得很好。
ddply(mydf,.(period),mean)
如果我对weighted.mean使用相同的方法，我会得到以下错误"'x‘和'w’必须具有相同的长度“，这是我不能理解的，因为weighted.mean部分在ddply之外工作。
weighted.mean(mydf$mycol,mydf$myweight) # works just fine
ddply(mydf,.(period),weighted.mean,mydf$my
浏览 0提问于2010-07-19得票数 9
回答已采纳
1回答
ID<-c("R1","R2","R2","R3","R3","R4","R4","R4","R4","R3","R3","R3","R3","R2","R2","R2","R5","R6")
event<-c("a","b","b","M&#
浏览 1提问于2016-01-08得票数 4
回答已采纳
我得到的是ROC曲线下具有不同区域的向量(不同的研究)，例如，
a <- c(.91, .85, .76, .89)
我也有每项研究的绝对参与人数，例如，
n <- c(50, 34, 26, 47)
我计算了以下区域的加权平均值
weighted.mean(a, n)
在R中，有没有一种方法可以根据我掌握的信息来计算加权平均值的95%置信区间？我研究了pROC，但据我所知，你需要每条ROC曲线的原始数据(我没有)。我将非常感谢您的任何建议！
浏览 0提问于2016-06-23得票数 4
2回答
我有以下数据
df <- data.frame(
Type=rep(LETTERS[1:6],3),
Level=rep(1:3,each=6),
Value=1:18)
我想增加两栏，
添加“R1”，即同一级别的值与类型A值的比率
添加'r2'，值与同一类型的1级值的比率
到最后
> df
Type Level Value r1 r2
1 A 1 1 1.000000 1.0
2 B 1 2 2.000000 1.0
3 C 1 3 3.0
浏览 4提问于2012-11-18得票数 3
回答已采纳
我是R的新手，有一个足够简单的问题可以回答；任何帮助都将不胜感激。
情况是，我有一个制表符分隔的数据矩阵(data matrix.txt)，如下图所示，最后一列包含组信息。
sampleA sampleB sampleC Group
obs11 23.2 52.5 -86.3 1
obs12 -86.3 32.5 -84.7 1
obs41 -76.2 35.8 -16.3 2
o
浏览 0提问于2012-11-22得票数 1
2回答
由于某种原因，自从升级到R-2.13.0和升级到plyr_1.5.1.tar.gz后，我得到了比我预期的更多的结果……我在一个旧版本的plyr上尝试过(不幸的是，版本不确定，因为我刚刚覆盖了它……)
library(plyr)
dd <-data.frame(matrix(rnorm(216),72,3),c(rep("A",24),rep("B",24),
rep("C",24)),c(rep("J",36),rep("K",36)))
colnames(dd) <- c("v1"
浏览 0提问于2011-04-15得票数 2
回答已采纳
1回答
我是R的新手，所以我的头衔可能没有任何意义。基本上，我有这个数据集，这是在两个不同地点收集的昆虫的长度。当我使用'summary‘函数时，我得到了两个站点的值，我想计算每个站点的平均值：
> ls()
[1] "simulies"
> ls.str()
simulies : 'data.frame': 615 obs. of 2 variables:
$ Site : Factor w/ 2 levels "Lovering","Orford": 2 2 2 2 2 2 2 2 2 2 ..
浏览 0提问于2017-10-05得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
