在Stata的相关矩阵中添加mean和sd列
，可以通过以下步骤实现：
- 首先，确保你已经加载了你要分析的数据集。
- 使用
correlate命令计算相关矩阵。例如，假设你要计算变量x和y的相关系数，可以使用以下命令：
correlate x y- 默认情况下，
correlate命令只会显示相关系数矩阵。要添加mean和sd列，可以使用pwcorr命令。该命令会在相关系数矩阵的右侧添加mean和sd列。例如，使用以下命令计算变量x和y的相关系数，并添加mean和sd列：
pwcorr x y, star(0.05) meansd其中，star(0.05)表示在相关系数小于0.05的情况下显示星号。
- 运行命令后，Stata会显示相关系数矩阵，其中包括mean和sd列。
在Stata中，相关矩阵的mean列显示每个变量的平均值，sd列显示每个变量的标准差。这些统计量可以帮助你了解变量的中心趋势和离散程度。
推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云云服务器（CVM）：https://cloud.tencent.com/product/cvm
- 腾讯云云数据库 MySQL 版：https://cloud.tencent.com/product/cdb_mysql
- 腾讯云人工智能平台（AI Lab）：https://cloud.tencent.com/product/ailab
- 腾讯云物联网平台（IoT Explorer）：https://cloud.tencent.com/product/iotexplorer
- 腾讯云移动开发平台（移动推送）：https://cloud.tencent.com/product/umeng_push
- 腾讯云对象存储（COS）：https://cloud.tencent.com/product/cos
- 腾讯云区块链服务（BCS）：https://cloud.tencent.com/product/bcs
- 腾讯云元宇宙平台（QCloud Metaverse）：https://cloud.tencent.com/product/qcloud_metaverse
请注意，以上链接仅供参考，具体产品选择应根据实际需求进行评估。
相关·内容
我正在尝试创建相关矩阵，它还包括每个变量的均值和标准差。 ** Set variables used in Summary and Correlation
local variables relationship commission anxiety enjoyment negotiation_efficacy similarity_values similarity_behaviors SPT_confidence own_SPT_effort
** Descriptive statistics
estpost summarize `variables'
matrix ta
浏览 80提问于2020-04-28得票数 0
1回答
我想重现一项研究的结果，该研究比较了Stata中多个组的平均值。我没有原始数据，但有关于均值、sd和n的特定于组的信息。
在Stata中，可以使用以下代码比较两组：
ttesti n1 mean1 sd1 n2 mean2 sd2
有相当多的组和相当多的变量来比较它们，这就是为什么我在寻找一个循环函数，它允许我安全地编写代码。
关键是一个循环函数，它使用列中相应的组特定信息(n，mean，sd)对所有行(比较组的变量)执行所描述的t测试。
浏览 2提问于2016-09-20得票数 0
1回答
我使用社区贡献的命令estout在Stata中创建一个汇总统计表。
我的代码在每个平均估计数下面的括号中显示标准偏差：
sysuse auto, clear
eststo clear
estpost summarize price mpg weight
esttab, cells(mean sd(par))
-------------------------
(1)
mean/sd
-------------------------
price 6165.257
浏览 3提问于2018-04-07得票数 0
我试着使用qwraps2()包获得汇总表，但是我的代码并没有生成一个想要的结果，而是一个奇怪的表达式，比如
\$\textbackslash pm\$，
没有显示变量名和
大多数数字没有逗号。想知道我能不能解决这三个问题。
我的数据是公司财务报表，公司分为大公司和小公司。我还把时间划分为危机前和危机后两个时期.所以我的分类小组有两个:小企业和大公司，危机前和危机后。因此，我的表有4列作为四个合并类别。
structure(list(firmid = structure(c("000020", "000020", "000020",
"0
浏览 8提问于2021-12-15得票数 2
2回答
我想知道是否有可能得到一个描述性的统计表(平均值，sd和n)，其中包括1个连续变量和2个分类变量。我在找我从Stata那里得到的东西： table with stata 我尝试使用dplyr，但没有成功： library(dplyr)
library(tidyverse)
data(mtcars)
mtcars%>%
group_by(carb, am)%>%
summarise_each(funs(mean, sd, n()), mpg) %>%
spread(key=am, value=mean)
# Groups: carb [6]
carb s
浏览 19提问于2020-05-02得票数 0
回答已采纳
我使用Stata 14中的tabstat命令和社区贡献的latabstat命令来生成相应的LaTeX代码：
tabstat ListOfVariables, stat(count mean sd min max)
当listOfVariables只包含一个变量时，变量将显示在行中，而stats则显示在列中。但是，当listOfVariables中有多个变量时，这些变量出现在列中，而stats则显示在行中。当列中列出许多变量时，这特别奇怪。
如何将变量从列切换到行？
浏览 7提问于2016-12-24得票数 0
回答已采纳
1回答
我试图在论文中报告描述性的统计数据。
因为我正在处理面板数据，所以我需要根据年份报告平均值和SD。
所以我希望年份放在第一行，变量放在第一列。
bys wave: asdoc variable1 variable2, replace stat(mean sd)
然而，这段代码产生了相反的结果。我该如何处理这个问题？(我使用的是Stata 14)
浏览 1提问于2021-06-12得票数 0
1回答
我在Stata中有下面的代码，我正在尝试将这些代码转换为python
foreach var in income_pc_r income_village_pc_r income_migration exp_total_pc_r sav_cash12_any_r sav_cash12_amtcon_pc_r anygorl_r_ffv gorl_r_ffv_condamt anygorl_g_ffv gorl_g_ffv_condamt hhsize_r migrant_r { if inlist("`var'","income_pc_r",&
浏览 6提问于2022-03-25得票数 0
回答已采纳
我有一些关于Stata的数据，其中有一些变量，如logTA和class。我有一千多个观察结果，logTA没有任何缺失的值。数据如下：
logTA class
-------- --------
.
.
21.26871
.
现在，我要做的是根据以下简单规则为class变量赋值：
replace class = 1 if logTA < mean(logTA) - sd(logTA)
replace class = 2 if logTA >= mean(logTA) - sd(logTA) &
浏览 6提问于2015-05-07得票数 0
回答已采纳
我想在Stata中生成一个表，其中包含4个不同组的平均值、差值和t值。 特别是，假设您有2x2研究设计，并希望显示结果变量的平均值和标准差，并添加另一列来测试处理一的差异，以及另一行包含处理二的差异和t值。 我有以下代码： clear
sysuse auto
gen large_trunk = (trunk > 14)
gen price_large_trunk = price if large_trunk == 1
gen price_small_trunk = price if large_trunk == 0
eststo price_domestic: qui estpost
浏览 14提问于2019-01-12得票数 0
回答已采纳
我希望在使用esttab时能够将统计数据集中在多个列上。
在我的玩具示例中，我希望N跨越两列：
sysuse auto, clear
est clear
qui estpost sum if foreign == 1
qui est store sum_foreign
qui estpost sum
qui est store sum_all
esttab sum_foreign sum_all, ///
replace ///
cells("mean(fmt(3)) sd(fmt(3))") ///
nonum ///
collabels(
浏览 0提问于2018-11-07得票数 1
回答已采纳
2回答
我试图在R中模仿table Stata命令，它执行摘要统计表。该命令允许您在生成的单元格中创建具有不同统计信息的交叉表。例如，在下面的示例中，我跨越了三个变量(category1、category2和category3)，并将metric1的均值和标准差以及均值和标准差metric2作为列向量。
在Stata上用下列单行获得所述行为。
table category1 category2 category3 ,c(mean metric1 sd metric1 mean metric2 sd metric2)
期望输出:对表的解释。
这里，产生的交叉表的每个列向量，假设交叉表的X包含X = [
浏览 3提问于2020-12-18得票数 0
回答已采纳
最近我一直在研究统计模拟，在生成随机数据以匹配相关矩阵的规范后，我希望将每一列转换为具有特定平均值和标准差的列。我在下面的代码中成功地做到了这一点，但它非常混乱，我想知道是否有更有效的方法来做到这一点。
#Input Correlation Matrix
sigma <- matrix(c(1.00, -0.03, 0.39, -0.05, -0.08,
-0.03, 1.00, 0.07, -0.23, -0.16,
0.39, 0.07, 1.00, -0.13, -0.29,
浏览 0提问于2013-06-16得票数 2
回答已采纳
3回答
我刚刚遇到了一个.do文件，我需要转换成R，因为我没有Stata许可证；我的Stata是生锈的，所以有人能确认代码正在执行我认为的那样吗？
为了重现性，我将把它转化为我在网上找到的数据集，特别是，它是Chatterjee，Hadi和Price教科书的一部分。
以下是Stata代码：
collapse (min) min_protein = protein ///
(mean) avg_protein = protein ///
(median) median_protein = protein ///
(sd) sd_protein
浏览 5提问于2015-06-22得票数 6
回答已采纳
我希望使用社区贡献的命令LaTeX将汇总统计表(均值、中位数、最小值、最大值、sd)导出到esttab，但将报告的统计数据舍入为两个小数。
下面是我的Stata代码：
eststo sumstats1: quietly estpost sum var1 var2 var3, detail
esttab sumstats1, cells("mean p50 min max sd") nonumbers label
esttab sumstats1 using "$repodir/output/tables/sumstats.tex", booktabs ///
浏览 0提问于2019-07-19得票数 0
回答已采纳
2回答
我在Stata中使用tabstat，并使用estpost和esttab将其输出到LaTeX。我有过
tabstat
若要按组显示统计信息，请执行以下操作。例如,
tabstat assets, by(industry) missing statistics(count mean sd p25 p50 p75)
我的问题是，是否有一种方法可以让tabstat (或其他Stata命令)显示按平均值排序的输出，以便那些具有较高均值的类别位于首位。默认情况下，当我使用industry时，Stata按tabstat的字母顺序显示。
浏览 4提问于2014-09-23得票数 4
回答已采纳
我想将以下Stata代码转换为R：
forvalues j=1/17 {
forvalues i=1870(1)2015 {
qui sum 'var' if country==`j' & year<=`i'
scalar mean_x = r(mean)
scalar std_x = r(sd)
replace stdup_`var'=(`var'-mean_xt)/std_xt if country==`j' & year==`i'
浏览 10提问于2022-08-08得票数 0
我想在我的ggplot中使用theme_stata。使用以下代码，除了主题之外，所有内容都可以正常工作：
df3 <- data_summary(panel_data,
varname = "invest_amnt",
groupnames = c("trend_id", "path_id")
)
df3$path_id <- as.factor(df3$path_id)
df3$trend_id <- as.factor(df3$trend_
浏览 0提问于2019-07-22得票数 0
1回答
我希望根据相关矩阵而不是原始数据进行回归。我已经看过了，但搞不懂。我在R里怎么做？
以下是一些代码：
#Correlation matrix.
MyMatrix <- matrix(
c(1.0, 0.1, 0.5, 0.4,
0.1, 1.0, 0.9, 0.3,
0.5, 0.9, 1.0, 0.3,
0.4, 0.3, 0.3, 1.0),
nrow=4,
ncol=4)
df <- as.data.frame(My
浏览 4提问于2020-12-30得票数 0
回答已采纳
在R中添加差分变量的正确语法是什么？类似于Stata中的'D.variable‘命令。我尝试在R中生成一个差分变量。我尝试了几种方法，但无法获得与Stata相同的输出。R代码： DF <- DF%>%
group_by(id) %>%
mutate(D_MDR = F_MDR - dplyr::lag(F_MDR))
summary(DF$D_MDR) 输出： Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
-0.7724 -0.1456 -0.0098 -0.0113 0.12
浏览 35提问于2021-07-17得票数 0
我正在尝试学习如何使用R中的MetRology包。我一直在使用GUM手册附录H.2中的数据作为一个简单的例子，尝试蒙特卡洛不确定性传播，包括相关性。我有两个变量V和I，因此我有一个2x2的相关矩阵，I，V在非对角线元素中的相关性。如果我把这个矩阵提供给uncertMC函数，我会得到这样的错误："Error in eigen(Sigma，symmetric = TRUE)：0x0矩阵“。如果我不包括相关矩阵，我就得不到错误。为什么在包含相关性的情况下会出现错误？metRology手册中有一个包含四个变量和一个4x4相关矩阵的示例，适用于我。有没有什么明显的地方我做错了？ 提前谢谢你！
浏览 60提问于2020-04-01得票数 0
回答已采纳
我已经使用psych::cortest.mat生成了相关矩阵的比较。我想把输出放在一个Sweave文件中，以便用knitr进行生产。当我使用Hmisc::latex()函数时，它可以工作，但它也为每个结果生成大约7位数，这使得它非常不吸引人。我可以在knitr中使用markup参数生成输出，但是我的文档中的所有其他表都使用latex输出更高效地生成(results='asis')。
有什么想法？
#Sample data
variable1<-rnorm(100, mean=7, sd=1)
variable2<-rnorm(100, mean=4, sd=1)
v
浏览 1提问于2015-07-30得票数 1
我希望在Stata的相关矩阵中，p<.05有一颗星，p<.001有两颗星。这是我目前使用的代码。代码仍然会生成一个相关矩阵，但是在它们应该出现的地方没有出现任何星星。感谢您的帮助！ asdoc corr RELATIONSHIP anxiety BEH_SIM SIM_VALUES sptconf NEG_EFFICACY spteffort SPTEFFORT_OTHER COOP_MOTIV COMP_MOTIV, star(0.5), replace
浏览 100提问于2020-05-03得票数 0
回答已采纳
1回答
我正在使用esttab + tabstat生成一个.tex文件，以便在LaTeX中打开。我快要得到我想要的了，但有一个问题：
在平均值之后，我如何在同一条线上得到我的标准差？它当前显示在均数后面的行上。
一辆MWE跟着。请注意，我实际上是在创建两个表并将它们附加到另一个表上。这只是在Stata中显示为两个单独的表，但是在我稍微修改代码以便保存到一个文件而不是输出到屏幕上之后，它在LaTeX中工作。如果有一种方法不像我所做的那样追加，并且一次做每件事，那将是超级的，但我没有意识到一个。还请注意，我正在遵循的代码，以便在这两个程序之间切换。
sysuse auto, replace
*create
浏览 3提问于2014-07-23得票数 1
回答已采纳
1回答
我试图使用Stata中的estout按子组创建一个汇总统计表。我的努力已经停止了。
我希望每个子组的计数出现在最后一个表中。我知道最后一行中的noobs删除了观察计数。但是，如果我只包括noobs，则只包括总观察计数，而不包括单个子组的计数。
例如：
. sysuse auto
(1978 Automobile Data)
. estpost tabstat price mpg rep78, by(foreign) statistics(mean sd) ///
columns(statistics) listwise
Summary statistics: mean sd
fo
浏览 0提问于2018-05-09得票数 0
回答已采纳
我想用平均值来绘制每一个数字列，而标准差是一条横列。如何对iris数据集执行此操作？
我正在尝试转换我的数据集，使在ggplot2中绘图变得容易。
我试过的
iris %>%
dplyr::select_if(is.numeric) %>%
dplyr::summarise(avg_sepal_length = mean(Sepal.Length),
avg_sepal_width = mean(Sepal.Width),
avg_petal_length = mean(Petal.Length),
浏览 7提问于2022-07-05得票数 0
回答已采纳
3回答
我在这里有点困惑，找不到一个好的答案。
我有一个数据收集，我正在试图聚合：
dt <- data.frame(age=rchisq(20,10),group=sample(1:2,20,rep=T))
当我聚合这个dataframe并将它保存到一个新的dataframe时，它只保存了2个观测数据和2个变量到全局环境中：
ag<-aggregate(age ~ group, dt, function(x) c(mean = mean(x), sd = sd(x)))
group age
1 1 9.119008
2 2 9.740361
即列组和
浏览 2提问于2021-01-21得票数 1
回答已采纳
1回答
我试图模拟变量，了解它们的边际分布及其相关矩阵。我知道我们可以使用像copula这样的软件包，但我不熟悉如何实现它。有人能帮上忙吗
#mean(w1)=0.6, sd(w1)=0.38; w1 is normally distributed
#mean(w2)=0.31; w2 is binary
#mean(w3)=0.226; w3 is binary
cor
w1 w2 w3
w1 1.0000000 -0.3555066 -0.1986376
w2 -0.3555066 1.0000000 0.1030849
w3 -0.1
浏览 1提问于2017-12-06得票数 2
回答已采纳
我试图使用ttest命令存储的结果在Stata中进行功率分析。
我的具体做法如下：
. ttest r, by(predictor)
Two-sample t test with equal variances
------------------------------------------------------------------------------
Group | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval]
---------+--------------------------
浏览 0提问于2018-09-08得票数 0
回答已采纳
我有以下data.frame：
dat <- data.frame(bifpar = c(0.5, 0.75, 1.0), sd_distfeed = rep(2,3), sd_nodist = rep(3,3), mean_distfeed = rep(10,3), mean_nodist = rep(20,3))
dat
bifpar sd_distfeed sd_nodist mean_distfeed mean_nodist
1 0.50 2 3 10 20
2 0.75 2
浏览 1提问于2020-12-03得票数 0
2回答
我试图使用列中静态变量和值的组合来运行模拟，对输出进行求和，并将各个输出存储在向量或数据帧中。
mean1 <- 2.4
sd1 <- 0.5
df <- data.frame(x = c(2, 3, 4), y = c(5, 6, 7))
我想做的是：
将x列中的每一行除以y列中的每一行，除以正态分布，使用mean1和sd1求和得到的行值，因此我每次模拟时只有一个值。
我想我理解如果不是逐行的话，我将如何得到值，所以对于第1行，应该是：
v1 <- replicate(n = 1, expr = rnorm(n = 100, mean = mean1, sd = sd
浏览 1提问于2021-03-04得票数 0
回答已采纳
我正在处理一个数据集，其中隐藏了我感兴趣的变量的值。我有这个变量的范围(min，max)，平均值和sd，并且对于每个观察值，我都有关于观察值所在的十进制数的信息。有没有什么方法可以使用Stata中的随机数生成器或rnormal()命令套件来计算这个变量的一些值？大致是这样的：
set seed 1
gen imputed_var=rnormal(mean,sd,decile) if decile==1
感谢大家在这方面的帮助，谢谢！
浏览 0提问于2020-06-15得票数 0
我目前正在写一篇论文的复制品。作者使用Stata，我对它知之甚少，所以我必须将它翻译成R。我有一个关于以下代码的问题： use "${directory_data}/income_dataset.dta", clear
reghdfe log_income post61_sc post65_sc male std_score mean_score_class privatista non_bocciato if tipo_scuola2==1 & laureato==1, vce(cluster liceo_anno) absorb(liceo anno_matur
浏览 27提问于2020-11-23得票数 1
回答已采纳
1回答
我正在学习R。我想使用简单、可读的R代码为出版物生成一个摘要统计表。该表应该将变量作为行，交替的平均值和SD作为列，两个分组变量也作为列。所有值应四舍五入为两位数，包括零(必要时添加零)。
以mtcars数据集为例，我希望表看起来如下所示，比较4、6和8辆汽缸车(自动或手动)：
| |4 0 | |4 1 | |6 0 | |6 1 | |8 0 | |8 1 | |
|:----|:---------|:-------|:----
浏览 2提问于2017-02-02得票数 0
回答已采纳
我在data.frame中有一个6的列表
它有3列：
id、T_C、Sales
T_C是测试还是控制
有人在这里帮助了我，我学会了如何通过循环找到mean()和sd()，而不是做单独的语句。
现在我的目标是从6个列表中删除异常值，并生成一个包含6个的列表(在删除异常值之后)。
data.frames (Str)#这是str中的6的列表
我能够获得每个列表的mean()和sd()，如下所示：
list_mean_sd <- lapply(dfList,
function(df)
{
浏览 2提问于2016-12-23得票数 0
回答已采纳
1回答
表中的子组和
、、、、
下面是我尝试过的Stata代码：
eststo clear
sysuse auto, clear
eststo Dom: estpost sum rep78 mpg turn trunk weight length if foreign==0
eststo For: estpost sum rep78 mpg turn trunk weight length if foreign==1
esttab Dom For, cells("mean(fmt(2))" "sd") ///
nonumber nodepvars noobs se collabels(
浏览 0提问于2018-09-13得票数 0
回答已采纳
1回答
假设我有一个非常大的这种形式的相关矩阵：
t1.rep1 = rnorm(n=100,mean=10,sd=)
t2.rep1 = t1.rep1 + rnorm(n=100,mean=3,sd=2)
t3.rep1 = t1.rep1 + rnorm(n=100,mean=2,sd=2)
t1.rep2 = rnorm(n=100,mean=2,sd=1)
t2.rep2 = t1.rep2 + rnorm(n=100,mean=0.5,sd=0.5)
t3.rep2 = t1.rep2 + rnorm(n=100,mean=0.7,sd=0.9)
t1.rep3 = rnorm(n=100
浏览 2提问于2014-02-25得票数 0
回答已采纳
我有一个表示组成员资格的Group列，以及许多其他包含数值的列。对于包含数值的每一列，我希望获得每个子组的平均值、标准差和样本大小。 以内置虹膜数据集为例，我编写了以下代码： lapply(names(iris)[1:4], function(x) {
iris %>%
group_by(Species) %>%
summarise(mean = mean(noquote(x), na.rm = T),
sd = sd(noquote(x), na.rm = T),
n = n())
}) 但是，每组数字
浏览 21提问于2021-07-06得票数 1
回答已采纳
我有一个数据表，列有三栏:准确性、特异性和敏感性。对于每一列，我想要计算每一列的平均值和sd。我希望结果是一个有六列(Accuracy.mean、Accuracy.sd、Specificity.mean、Specificity.sd、.)的数据表。到目前为止我这么做了：
test<-function(x){
c(mean = mean(x),sd = sd(x))
}
b<-a[ ,lapply(.SD, test)]
有了这个，我得到了每列两行。我该如何解决这个问题？
浏览 0提问于2018-03-18得票数 3
回答已采纳
我有以下Python代码和输出：
>>> import numpy as np
>>> s = [12.40265325, -1.3362417499999921, 6.8768662500000062, 25.673127166666703, 19.733372250000002, 21.649556250000003, 7.1676752500000021, -0.85349583333329804, 23.130314250000012, 20.074925250000007, -0.29701574999999281, 17.078694250000
浏览 3提问于2013-10-04得票数 9
回答已采纳
1回答
在完成以下代码后，我想在list表中的每个数据框中添加一列：
#list of my dataframes
df <- list(df1,df2,df3,df4)
#compute stats
stats <- function(d) do.call(rbind, lapply(split(d, d[,2]), function(x) data.frame(Nb= length(x$Year), Mean=mean(x$A), SD=sd(x$A) )))
#Apply to list of dataframes
table <- lapply(df, stats)
例
浏览 3提问于2014-06-23得票数 0
3回答
R循环问题-遍历数据
、、、
我是新的R和需要一些帮助解决一个问题。总之，我有一个具有不同值的dataframe : 10行6列。每一列代表一个变量:列1- n1，列2- mean1，列3- varaince1，列4- n2，column 5- mean2，列6- variance2。每一行都是这些变量的不同组合。我希望遍历每一行并生成两个样本--样本1- -带有n1、mean1和sd1的随机正态变量(variance1 sqrt)和样本2--带有n1、mean1和sd1的随机正态变量(variance1 sqrt)。有人能告诉我什么是最好的方法吗？谢谢你的帮助。
下面是使用dput()函数的示例数据：
structure
浏览 6提问于2022-04-24得票数 1
1回答
创建函数外的矩阵
、、、
到目前为止我还找不到解决方案。我有一个7列的数据帧X。我希望从矩阵中收集的第6列和第7列获得统计数据。我想用一行代码来解决这个问题。我需要一个函数，它取数据帧X的第6列和第7列，并计算平均值，变量，sd，... 我使用的代码是这样的： matrix(c(mean(X$6), mean(X$6, trim = 0.05), var(X$6), sd(X$6),
mean(X$7), mean(X$7, trim = 0.05), var(X$7), sd(X$7)),
nr = 4, dimnames = list(c("mean",
浏览 17提问于2019-10-17得票数 1
回答已采纳
我有一个数据集df，其中我希望显示组间的均值和st偏差比较，同时在输出中保留id、start、end和value列。
id start end value
a 5/1/2020 6/1/2020 2
a 6/1/2020 7/1/2020 3
a 7/1/2020 8/1/2020 4
a 8/1/2020 9/1/2020 20
b 5/1/2020 6/1/2020 15
b 6/1/2020 7/1/2020 2
b 7/1/2020 8/1/2020
浏览 1提问于2020-11-17得票数 1
回答已采纳
3回答
我在R中有一个数据表，如下所示
DT = data.table(a = c(1,2,3,4,5), a_mean = c(1,1,2,2,2), b = c(6,7,8,9,10), b_mean = c(3,2,1,1,2))
我想再创建两个列a_final和b_final，分别定义为a_final =(a-a_mean)和b_final =(b-b_mean)。在我的实际用例中，可以有大量这样的列对，我想要一个符合R的数据表精神的可伸缩的解决方案。
我尝试了一些类似的东西
DT[,paste0(c('a','b'),'_final') :=
浏览 2提问于2019-12-30得票数 3
回答已采纳
我如何按列分组，然后计算R中其他列的平均偏差和标准差？
例如，考虑著名的Iris数据集。我想做一些类似于按物种分组的事情，然后计算花瓣/萼片长度/宽度的平均值和sd值。我知道这与拆分-申请-合并有关，但我不知道如何从那里着手。
我能想到的是：
require(plyr)
x <- ddply(iris, .(Species), summarise,
Sepal.Length.Mean = mean(Sepal.Length),
Sepal.Length.Sd = sd(Sepal.Length),
Sepal.Width.Mean = mean(Sepal.Wi
浏览 1提问于2016-05-26得票数 3
回答已采纳
1回答
R:多组的独立t检验
、、、
我有一个名为'dat‘的数据集，列有5列: month；mean0；sd0；mean1；sd1。它看起来如下(但有数字)：
月mean0 sd0 mean1 sd1
1
2
3.
。。
48
我想用一个独立的(而不是配对的)t检验来比较mean0和mean1在1到48之间的每个月。理想情况下，输出将放在另一个名为‘dat1 1’的数据中，列用于: t-statisitc、自由度(DF)和p值。就像这样：
月t统计DF值
1
2
3.
。。
48
我试过使用dplyr和broom包，但似乎找不出答案。任何帮助都将不胜感激。
浏览 0提问于2018-09-20得票数 0
回答已采纳
请考虑以下几点：
library(tidyverse)
df <- tibble(x = rnorm(100), y = rnorm(100, 10, 2), z = x * y)
df %>%
mutate_all(funs(avg = mean(.), dev = sd(.), scaled = (. - mean(.)) / sd(.)))
是否有方法通过引用mean和dev列来避免两次调用avg和dev。我想的是
df %>%
mutate_all(funs(avg = mean(.), dev = sd(.), scaled = (. - avg) / de
浏览 0提问于2018-11-04得票数 8
回答已采纳
2回答
summarize命令在Stata中创建各种标量。例如，您可以随后通过gen mean=r(mean)存储平均值或最小/最大值。还可以通过summarize varname, detail选项获得更复杂的度量。通过此，还可以得到50%百分位数的中位数。
我的目标是存储中位数。有相应的标量吗？
在像summarize这样的标准操作之后，我在哪里可以获得关于存储标量的信息？据我所见，它们没有被列在Stata手册中。
浏览 6提问于2017-06-22得票数 0
回答已采纳
我对新的dplyr::summarise()函数有一些问题。
这是数据
df <- data.frame(id = factor(1:10),
group = factor(rep(letters[1:2],each = 5)),
w1 = rnorm(10),
w2 = rnorm(10),
w3 = rnorm(10),
dummy = as.character(LETTERS[1:10]),
浏览 1提问于2020-09-02得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
