R和SPSS结果之间的不同置信区间
R和SPSS都是常用的统计分析软件，在计算置信区间方面有些许差异。
R是一种免费开源的统计计算和图形化软件，它提供了丰富的统计分析函数和包，可以广泛应用于数据分析、数据可视化、机器学习等领域。在R中，计算置信区间的方法有多种，比如使用t.test()函数可以计算均值的置信区间，使用confint()函数可以计算其他统计量的置信区间。R的灵活性和扩展性较高，可以根据具体需求灵活地进行统计分析和置信区间的计算。
SPSS是一种商业化的统计分析软件，提供了直观友好的界面，适合非专业人士使用。SPSS也支持计算置信区间，通过分析菜单中的描述统计功能可以计算均值、比例等统计量的置信区间。SPSS的优势在于其用户友好性和操作简单性。
在实际使用中，R更常用于统计学界和科研领域，由于其灵活性和可扩展性，研究人员可以根据具体需求编写代码进行分析。而SPSS更常用于商业和市场调研领域，由于其操作简单性，非专业人士可以轻松进行数据分析。
对于R和SPSS结果之间的不同置信区间，需要具体问题具体分析。在实际应用中，根据数据的特点和分析目的选择适合的统计软件和方法进行计算。对于置信区间的计算，可以参考以下腾讯云相关产品和文档：
- R语言相关产品介绍：https://cloud.tencent.com/product/crs
- SPSS相关产品介绍：https://cloud.tencent.com/product/sps
相关·内容
我正在使用SPSS来生成ROC曲线，但是ROC曲线不能给我灵敏度和特异性的置信区间。
因此，如果有人能帮助我在SPSS中产生敏感性和特异性的置信区间，那将是对我最大的帮助。
谢谢你
浏览 40提问于2016-06-05得票数 1
1回答
R中的倾向分数
、、、
我对R中的倾向性分数有一个一般性的问题，我过去使用过SPSS Modeler，我开始主要使用R来构建机器学习模型。在SPSS Modeler中，我们使用了一个名为“倾向性分数”的选项，这意味着这些分数显示了结果的可能性。
那么，是否可以在R中对每个客户的倾向得分进行编码？如果有的话，能告诉我怎么走吗？任何其他得分也可以给出结果的可能性，例如R中的置信度(而不是置信区间)？
例如:下周纽约、伊利诺伊州、佐治亚州、夏威夷和加利福尼亚州有雨、不下雨、多云和晴天的可能性有多大？结果:纽约可能会下雨的可能性是5%，IL可能是晴天的可能性是60%，CA....etc不会下雨的可能性是10%。
请从IBM
浏览 0提问于2015-01-31得票数 1
我有一个小的数据集，总病例数在20以下。当我用这个数据集在SPSS上使用Kaplan-Meier估计进行生存分析时，我能够得到特定的95%的上下限置信区间值，但在R(使用survfit，conf.type“对数”) 0.95UCL没有被估计。我有一种感觉，当我使用小数据集时，这种情况经常发生。我不确定它是否准确，但我读到过SPSS提供log-log作为默认设置。那么，谁能让我知道为什么结果是这样的？提前谢谢你。
浏览 71提问于2019-10-11得票数 0
我之前把这个作为一个问题贴到了我最初关于伸缩的问题上，但这被认为是不合适的，并被删除了，要求作为一个新的问题提交，所以我们来了！ 我和我的同事使用r (prcomp)、SAS和SPSS执行PCA。我们都使用了缩放和居中。所有三个都产生了相同的变量/因子加载(程序之间的术语不同)。这里没有问题。 然而，个体在R中产生的分数与SAS和SPSS不同。当我使用scale(to.my.pca)时，结果与SAS和SPSS相同。所以，看起来R中的scale.=TRUE并没有传递给个人分数。使用r和SAS或r和SPSS的人能证实这一点吗？ 非常感谢！
浏览 22提问于2018-12-27得票数 2
我正在尝试使用R运行重复测量方差分析，并将其与SPSS输出进行比较，结果差别很大！也许我在哪里犯了一个错误，但我不能弄清楚，所以一些样本数据: id是主题。每个受试者对三个项目(res_1，res_2和res_3)进行一个评级。我想比较一下物品的整体效果。
id<-c(1,2,3,4,5,6)
res_1<-c(1,1,1,2,2,1)
res_2<-c(4,5,2,4,4,3)
res_3<-c(4,5,6,3,6,6)
## wide format for spss
table<-as.data.frame(cbind(id, res_1, res_2, re
浏览 8提问于2015-03-31得票数 1
我得到了它用R和spss分析mancova。我想比较一下他们的输出。但在R中结果与.Which中的代码mancova有很大不同，结果是正确的R还是spss？代码R写得正确吗？请帮帮我
library(jmv)
Data=read.csv(file.choose())
Data
attach(Data)
names(Data)
> y2=mancova(data=Data,deps=vars(postGOcorrect,postNOgoerror,postRTgoerror
+ ,postRTNogoerror)
+ ,covs =vars(p
浏览 1提问于2020-02-24得票数 0
1回答
我正在使用Ward的方法对包含1000个观察值和37个变量(均为5点likert标度)的数据集进行层次聚类分析。
首先，我在SPSS中通过
CLUSTER Var01 to Var37
/METHOD WARD
/MEASURE=SEUCLID
/ID=ID
/PRINT CLUSTER(2,10) SCHEDULE
/PLOT DENDROGRAM
/SAVE CLUSTER(2,10).
FREQUENCIES CLU2_1.
我还在R中执行了分析：
datA <- subset(dat, select = Var01:Var37)
dist <- dist(dat
浏览 6提问于2020-11-27得票数 1
我有三个因素:单词、类型和语域。在SPSS中，很容易在SPSS中进行成对比较(或简单比较)，语法是：
/EMMEANS=TABLES(word*register*type) COMPARE(type) ADJ (BONFERRONI)
它会给我一个这样的结果：
但是我如何在R中使用Multcomp包实现这一点呢？假设我有一个这样的anova结果：
library (multcomp)
anova <- aov (data, min~word*register*type)
summary (anova)
如何进行类似于SPSS的成对比较？
浏览 3提问于2016-06-22得票数 0
我安装了SPSS modeler 18.1。我想在SPSS modeler中使用R。所以我安装了R版本3.4.2和R essentials for SPSS Modeller 18.1。但是当我打开SPSS建模器时，R节点没有出现在调色板中。可能的问题是什么？
谢谢!
浏览 3提问于2017-11-22得票数 1
我想知道SPSS在计算UICI和LICI (上、下个体置信区间)时所做的事情。我之所以问这个问题，是因为当我们用一个简单的线性回归的输出表计算一个给定个体的“手工”预测区间时，我们得到了一个稍微不同的区间(高达0,005个差)。
我无法在网上找到如何获得用于此命令的代码，以便更仔细地查看SPSS在“选中”复选框以确定平均值和个人预测间隔时所做的工作。
谢谢你的帮忙,
浏览 1提问于2016-05-07得票数 0
回答已采纳
我以前在spss.system.file包中使用过几次memisc函数，而且它总是没有错误地工作。但是，最近我尝试导入了几个包含字符串变量的SPSS文件。似乎每当我使用spss.system.file导入这些变量时，就会产生6列额外的空空间。例如，如果以下字符向量是SPSS中的字符串变量：
Age <- c("18 years old", "18 years old", "24", "19 years old", "18", "18")
使用spss.system.file将其读入R中将产生一
浏览 7提问于2014-11-04得票数 1
回答已采纳
1回答
我使用代码R中的表测试spss中的受试者间效应计算，但结果不同，真的不知道原因是什么？，请帮助我。感谢您 fit2222<-manova(cbind(postrunning,postbalance,postbilateral,poststrngth,postupperlimb,postresponse
,postvisumotor,postdexterity)
~ group+prerunning+ prebalance+prebilateralcoordination + prestrength + preupperlimb
浏览 37提问于2020-03-26得票数 0
作为一名R用户，我现在正在尝试学习SPSS语法。
我习惯于在R脚本的前面添加命令rm(list=ls())，以确保在我继续工作之前R是空的。
SPSS有没有类似的命令？谢谢。
浏览 8提问于2011-06-20得票数 3
回答已采纳
我对n=452参与者进行了重复测量。SPSS的弗里德曼检验给出了df=2，卡方为36.970，而下面是R和python的输出。我无法解释卡方差值(R和python的30和SPSS的36 )的差异，尽管它不是太多。有人能帮我理解这个吗？
#Friedman rank sum test in R
friedman.test(data.matrix(x))
data: data.matrix(x)
Friedman chi-squared = 30.389, df = 2, p-value = 2.518e-07
#--------------------------------
浏览 1提问于2020-05-22得票数 0
2回答
实际问题
鉴于以下数据：
library(dplyr)
df <- tibble(v1 = 1:6, cond = c(1, 0, 1, 1, 0, 1))
## # A tibble: 6 × 2
## v1 cond
## <int> <dbl>
## 1 1 1
## 2 2 0
## 3 3 1
## 4 4 1
## 5 5 0
## 6 6 1
我想要计算累积和累积积的混合物。如果cond =1，则计算当前v1和前面计算的结果。如果con
浏览 11提问于2022-07-25得票数 1
对于我的两个变量(ms和gar)，我都有一个包含11236个案例的非常大的样本。现在我想在SPSS中计算Spearman的rho相关性和bootstrapping。
我找出了SPSS中带偏差校正和加速置信区间的自举的标准语法：
DATASET ACTIVATE DataSet1.
BOOTSTRAP
/SAMPLING METHOD=SIMPLE
/VARIABLES INPUT=ms gar
/CRITERIA CILEVEL=95 CITYPE=BCA NSAMPLES=10000
/MISSING USERMISSING=EXCLUDE.
NONPAR CORR
浏览 5提问于2016-09-17得票数 0
回答已采纳
1回答
我试图从我的数据中计算贝叶斯因子，对于我的混合效果模型，我在R和SPSS中得到了非常不同的结果。这对于线性模型很好，但不适用于二项式。这是R代码：
``memory.model = glmer(correct ~ (1|ps) + (1|item), data=memorystudy, family=binomial, glmerControl(optimize = "bobyqa"))
memory.model2 = glmer (correct ~ encoding + retrieval + (1|ps) + (1|item), data=memorystudy, fam
浏览 3提问于2018-05-15得票数 1
我正在做一个需要我使用SPSS的项目。我对R比较满意，但我需要呆在SPSS环境中，我将无法使用R插件。
我的问题与财政年度有关，并根据它们引入的额外日期进行调整。例如：
* fake data to permit syntax to be run.
data list free/number.
begin data
1
12
123
1234
12345
end data.
* calculate the number of days between the start/end of a fiscal year.
* this period contains a leap-day that
浏览 3提问于2016-10-19得票数 2
回答已采纳
目的是在R中模拟SPSS重新编码过程。复制命令很难翻译。
在SPSS中，我的代码是
RECODE A (1,2 = 1) (3,4 = copy) (8 thru hi = 3) (else = 1) into B.
在A上应用，它看起来像
A <- c(1,2,3,4,5,NA,7,8,9)
我得到以下(SPSS)结果：
A = 1,1,3,4,1,1,1,3,3
在R中，类似的代码将如下所示：
B <- Recode(A, recodes = ("c(1,2) = 1; c(3,4) = c(3,4); c(8,9) = 3; else = 1"), as.nu
浏览 3提问于2013-02-28得票数 3
1回答
我不知道如何解释在使用CausalImpact函数时在CausalImpact R包中获得的置信区间。
我感到困惑，因为我认为有一个矛盾-模型返回一个非常低的p-值(0.009)，这表明有一个偶然的效果，但“实际”线(实线)似乎很好地在95%的置信范围内的反事实。如果有因果影响，你不认为这条线会超出蓝带吗？
以下是我的研究结果：
下面是模型总结结果(我对大文本表示歉意)。
这里发生了什么事？
浏览 1提问于2017-09-21得票数 0
回答已采纳
我正在研究某一行业的收入增长率在收入分配的不同部分是否存在显著差异，以确定收入差距是否正在显著拉大或缩小。
QUANTREG模型
我在SPSS中做了分位数回归(我刚开始编码，只有Python的基本知识，所以我需要你的帮助)。因变量是指数化收入，自变量是时间(这个数据集中的季度)，人口统计学组，职业部门。我还添加了每个虚拟人的交互条件和时间变量。
因此(至少在我看来)，这个模型允许在三个层次上比较收入的变化：
(例如，与数据输入工作相比:数据分析增加了100欧元，数据科学增加了200欧元)time*data_scientist系数在90%Q中比在10%Q中大得多，这表明高收入数据科学家的收入随着
浏览 1提问于2021-01-02得票数 0
回答已采纳
我在受试者之间进行2(性别)x2(建议)，R和SPSS都报告了相同的方差分析:建议:F= 372.012，效果df = 1，错误df = 661；性别x建议:F= 45.449，效果df = 1，df =2。
在计算偏eta平方时，R(跨多个R包: rstatix和DescTools)报告了0.221条建议和0.031条性别x建议。但是，无论是SPSS还是使用Lakens的影响大小电子表格()，都得到了0.360的建议和0.064的性别x建议。
R是否以与标准不同的方式计算部分eta平方值？
这里是一个样本数据集，它只包含测试问题所必需的变量：
这里是我在RStudio:中使用的代码
# l
浏览 4提问于2020-08-12得票数 1
1回答
我读过关于如何计算R中各种列的方法的各种文章，但我遇到的问题是，到目前为止，我在R中使用的函数给我的结果与其他程序不同(SAS和SPSS)。
到目前为止，我已经使用了这些功能：
newdata2$birdsMean <- rowMeans(newdata2[,7:15], na.rm=TRUE)
(The bird data appears in columns 7 to 15.)
以及：
newdata2$mean=rowMeans(newdata2[,c("bird1", "bird2", "bird3", "bird4&#
浏览 2提问于2017-10-31得票数 0
我试图在R.中获得Bonferroni同时的置信区间，我有以下数据集供我练习：
df2 <- read.table(textConnection(
'group value
1 25
2 36
3 42
4 50
1 27
2 35
3 49
4 57
1 22
2 37
3 45
4 51'), header = TRUE)
我试过了
aov(formula = value ~ group, data = df2)
然而，这并不能同时输出置信区间。使用SAS，计算结果应如下：
浏览 4提问于2018-02-01得票数 3
回答已采纳
1回答
我有22名学生的样本，我用两种不同的方法来识别他们的个性。我的问题是，我的统计能力不太好，我想用SPSS统计软件中的变量来计算这两种方法之间的一致性程度(即学习者1在方法1中有X作为他的个性，X在方法2中作为他的人格)。
此外，我还可以使用SPSS中的一些图表为每种方法生成图表。
这可以帮助读者轻松查看这两种方法中是否存在一致的结果。谢谢。
浏览 3提问于2015-04-21得票数 0
1回答
R中方差分析的重复测度
、、、
我进行了一项研究。这是数据的结构：重复测量(mood_1和mood_2)因变量，群体是主体因素之间的变量，反刍动物是协变量。
我在R：中准备了这样的数据，但是，我从R中得到的结果与从SPSS中得到的结果不一样。
在SPSS中有意义，在R中没有显着性，在第一张图片中进行同样的思考到底缺少什么？
model <- lm(depression ~ group * ruminati)
anova(model)
浏览 9提问于2021-12-21得票数 0
我有这样的代码：
#load library foreign and mass
library(foreign)
library(MASS)
#read spss file
read.spss("C:\\Users\\Morisato\\Desktop\\Modelação Estatistica e Análise de Dados\\1º Ano\\1º Semestre\\Análise de Dados com Software Estatistico\\Exercicios Regressão Linear\\imc_mest.sav")
#save data
浏览 2提问于2013-12-21得票数 0
回答已采纳
2回答
直到最近，我用SPSS进行统计，但由于我已经不在大学，我改到了R。事情进展顺利，但我似乎无法复制我在SPSS中重复使用LMM所获得的结果。我确实在这里找到了一些看似相关的脚印，但这些并没有解决我的问题。
这是我试图在R中复制的SPSS脚本
MIXED TriDen_L BY Campaign Watering Heating
/CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1)
SINGULAR(0.000000000001) HCONVERGE(0,
ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVER
浏览 1提问于2016-04-27得票数 0
1回答
我只是想弄清楚R -> Python的一个问题:为什么这两个项目不能给出相同的结果？ 计算均值= 65，s= 22，n= 121的样本数据的95%置信区间。 R：tsum.test(n.x=121, mean.x=65, s.x=22) 给出了61.04014 68.95986的95%置信区间 Python：stats.norm.interval(alpha=0.95, loc=65, scale=22/np.sqrt(121)) 给出了(61.080072030919894, 68.9199279690801)的95%置信区间 我认为这些结果应该是相同的，或者我没有为R的tsum.t
浏览 27提问于2019-03-11得票数 0
因此，我尝试了两种不同的方法来获得和注释我的相关散点图的直线方程--用R和我得到的结果，对我来说，这看起来很奇怪，并且与SPSS中看起来不那么奇怪的结果不同。
我显然没有足够的“声誉点”来插入图像，但这是我使用的代码.
model <- lm(hads ~ fatigue, mydata)
eqn <- as.character(as.expression(
substitute(italic(y) == a + b * italic(x) * "," ~~ italic(r)^2 ~ "=" ~ r2,
list(a = format(coef(m
浏览 3提问于2014-04-15得票数 0
回答已采纳
2回答
将R值标签导出到SPSS文件
、、、、
我正在和一些只使用SPSS的人合作，虽然我使用了R。我试图在R中创建一个带有变量标签和值标签的数据，以便将因素变量导出到SPSS。但是，一旦在SPSS中打开文件，就会得到变量名，而不是值名。
有没有人知道如何做到这一点？
到目前为止我已经尝试过了。mtcars例子。
mtcars = modify(mtcars,{
var_lab(mpg) = "Miles/(US) gallon"
var_lab(cyl) = "Number of cylinders"
var_lab(disp) = "Displacement (cu.in.
浏览 2提问于2019-08-22得票数 3
回答已采纳
我需要了解如何计算线性回归的置信区间。我自己得到的值与Matlab中的值不同。所以，我一直在努力理解如何在Matlab中完成这个任务。在polyval.m中，我不理解下面代码的部分。似乎E等价于E，R= qr(V,0)。但不确定到底是什么。我的问题是：
[E,R] = qr(V,0) different from [Q,R]=qr(V);怎么样？
如何使用E(=V/R)来计算置信区间(以下δ)？
%S是一个包含三个元素的结构:原始X的Vandermonde矩阵的三角因子、自由度和残差的范数。
E= V/R；
E= sqrt(1+sum(E.*E,2))；
..。
delta = nor
浏览 3提问于2012-11-13得票数 0
回答已采纳
我正在使用spss代码。我想把spss计算的RXX结果打印到控制台(输出文件)。在spss中有任何可能的方法来做到这一点吗？
代码如下所示：
COMPUTE R = CR. /* CR is default name of matrix */
COMPUTE N = NCOL(R). /* Number of variables */
COMPUTE RXX = R(2:N,2:N). /* Correlations between the X’s */
COMPUTE RXY = R(2:N,1). /* Correlations between Y and the X’s */
CALL
浏览 10提问于2019-10-25得票数 1
回答已采纳
2回答
在R中，predict.lm基于线性回归的结果计算预测，并提供计算这些预测的置信区间。根据手册，这些间隔是基于拟合的误差方差，而不是基于系数的误差间隔。
另一方面，基于逻辑回归和泊松回归(以及其他一些)计算预测的predict.glm没有置信区间的选项。我甚至很难想象如何计算这样的置信区间，以便为泊松和逻辑回归提供有意义的见解。
在某些情况下，为此类预测提供置信区间是有意义的吗？如何解释它们呢？这些情况下的假设是什么？
浏览 3提问于2013-01-20得票数 70
回答已采纳
我有一个SPSS文件，但没有SPSS。所以我想在R.
如果我打开它使用：
library(foreign)
dat <- read.spss("file.sav", to.data.frame=TRUE)
我收到警告了
re-encoding from CP1252
Warning message:
In `levels<-`(`*tmp*`, value = if (nl == nL) as.character(labels) else paste0(labels, :
duplicated levels in factors are deprecated
浏览 4提问于2014-01-20得票数 2
回答已采纳
2回答
我正在从SPSS跳到R，我想知道你是怎么处理系统丢失的.
例如，如果我想将以下SPSS代码重写为R：
RECODE income (1 THRU 6 = copy) (else = SYSMIS) INTO income2
我能够写以下内容：
income_2018$income2 <- dplyr::recode(income_2018$income, '1' = 1L, '2' = 2L, '3' = 3L, '4' = 4L, '5'
= 5L, '6' = 6L)
如何处理系统错误(
浏览 1提问于2020-03-24得票数 3
回答已采纳
我需要使用SPSS命令提示符将语法导出转换为sav (我需要从wcf服务以编程方式完成)。我尝试使用命令提示符为它提供sps、dat和spj文件，其中spj是一个带有到sps的语法路径的xml，它做了一些事情，但没有结果(没有文件)。
我在spss中打开了语法文件，即使在SPSS中，我也只能选择将其保存为语法，而不是保存为sav文件。
从这3个文件(我猜是dat或csv，sps和spj)创建一个sav文件，我缺少什么？
浏览 5提问于2013-07-23得票数 0
1回答
我最初在SPSS中运行我的数据，因为弄清楚lmer包对我来说需要一些时间来学习。我花了几个星期用R写了一个脚本，但我在R中的输出与我使用SPSS得到的结果不同。 我有3个固定的效果:组，会话和TrialType。 当我在SPSS中运行混合模型时，我得到了交互组*会话p=.08或p=.02，这取决于我使用的协方差结构。这是我想使用R的部分原因，因为我没有足够的信息来帮助我决定使用哪种结构。 这是我在R中的模型。我正在使用对数似然测试来获得这个Group*会话交互的p值。 Mod2 = lmer(accuracy ~ group*session*trialtype + (trialtype|su
浏览 65提问于2019-07-01得票数 1
5回答
在我的公司，我们正在考虑逐步淘汰SPSS，选择R。在过渡期间，我们仍然会使用SPSS数据文件格式(.sav)。
将此SPSS数据文件导入R时遇到问题。当我将SPSS文件导入R时，我希望同时保留变量的值和值标签。foreign包中的read.spss()函数让我可以选择保留变量的值或值标签，但不能同时保留两者。
AFAIK，R允许因子变量具有值(级别)和值标签(级别标签)。我只是想知道是否有可能以某种方式修改read.spss()函数来合并这一点。
或者，我遇到了来自memisc包的spss.system.file()函数，它理应允许这种情况发生，但它要求一个单独的语法文件(codes.file
浏览 10提问于2013-01-30得票数 5
回答已采纳
2回答
我为这个问题的广泛性表示歉意，但我想知道是否有人有任何建议/建议。
一位朋友的公司使用SPSS进行统计数据、调查等。他们希望能够通过web界面对信息进行切片和切分。我对SPSS一无所知，所以可能有一个简单的web前端，但作为一个PHP/MySQL的人，我想知道SPSS数据是否可以以一种MySQL可以导入或解析的方式导出。然后，我可以构建一个AJAX前端来实时返回查询。
这个想法是一个结果表，当你选择“性别:男性”，“年龄:20-25”时，它会发生变化……每次更改都会更新浏览器中的结果。
所以..。
1) SPSS可以通过某种前端来实现这一点吗？
2)数据是否可以方便地导出为通用格式？
在我和朋
浏览 2提问于2010-01-27得票数 8
回答已采纳
在这方面存在许多问题，但我不能用它们来解决我的问题。我有一个数据样本，我想为它的曲线创建置信区间。在这里，我提供一个简单的例子：
import numpy as np
%matplotlib inline
import matplotlib.pyplot as plt
mean, lower, upper = [],[],[]
ci = 0.2
for i in range (20):
a = np.random.rand(100) # this is the output
MEAN = np.mean(a)
mean.append(MEAN)
std =
浏览 14提问于2022-04-18得票数 1
回答已采纳
我想找出致命剂量(LD50)及其在R中的置信区间。其他软件线Minitab，SPSS，SAS提供了三个不同版本的这类置信区间。在R中的任何包中，我都找不到这样的间隔(我还从sos包中使用了findFn函数)。
我怎样才能找到这样的间隔？我基于Delta方法编写了一种类型的间隔(不确定其正确性)，但希望使用R包中的任何已建立的函数。谢谢
MWE:
dose <- c(10.2, 7.7, 5.1, 3.8, 2.6, 0)
total <- c(50, 49, 46, 48, 50, 49)
affected <- c(44, 42, 24, 16, 6, 0)
finney
浏览 2提问于2015-01-20得票数 4
回答已采纳
注意，标记为可能重复的前一个问题不是重复，因为前一个问题涉及向后消除，而这个问题涉及向前条目。
我目前正在进行一个模拟，我想展示逐步回归是一个有偏估计量。特别是，先前的研究人员似乎使用了SPSS (或与其相同的内容)中的逐步过程之一。这涉及到使用r-平方变化的F值的p值来确定是否应该向模型中添加额外的变量。因此，为了使我的模拟结果产生最大的影响，我需要在R中复制SPSS逐步回归过程，而R有许多逐步过程(例如基于AIC)，我发现的过程与SPSS不一样。
我找到了这个。它看起来很管用，但是函数的输入和输出有点奇怪。我已经开始调整它，使它(a)以一个公式作为输入，(b)返回最佳拟合模型。函数的逻辑就
浏览 1提问于2013-05-12得票数 2
1回答
我将根据R中包含的数据集"CES11“数据集确定该比例的置信区间。我用于检测的变量是堕胎，其值为”否“和”是“。我想找出“是”的比例的置信区间，但R做了一个假设检验，并给出了第一级的置信区间，即“否”。如何将其更改为结果对应于等于“是”的级别？
Frequency counts (test is for first level):
abortion
No Yes
1818 413
1-sample proportions test without continuity correction
data: rbind(.Table), null probability
浏览 7提问于2022-02-12得票数 1
回答已采纳
我正在从SPSS (sav文件)导入数据。在此数据集中，有几个日期列，其中几个列包含缺少值的单元格。所有日期列都作为日期变量导入到R中，到目前为止还不错。所有来自SPSS的数据条目都被正确地导入到R中，到目前为止仍然很好。然而，SPSS中的一些缺失值在输入R时被解释为1582-10-14，我知道这是公历的开始，所以我认为也许SPSS中的一些单元格是0而不是NA，但是在检查这些单元格时，它们确实是空的，并且看起来与其他在R中作为NA正确导入的空单元完全相同。
对不起，我想不出一个可重复的例子，我不能把我的SPSS文件放在这里，完整性的原因。
我尝试过使用里约(import_list)和港湾(re
浏览 13提问于2022-11-02得票数 1
我遇到了一个有趣的问题，试图使我的R结果与SPSS的一个同事的结果一致。似乎在SPSS中，在求和时，NA + NA = NA和0 + NA = 0，但在R中，我似乎无法复制这一点。在R中，要么和到NA，要么和到0。但是，在我的例子中，如果sum命令中的所有值都丢失了，那么它们应该是NA，但是至少有一个0，应该是0。
nana <- c(NA, NA)
na0 <- c(NA, 0)
sum(nana) # sums to NA
sum(na0) # sums to NA
sum(nana, na.rm=T) # sums to 0
sum(na0, na.rm=T) # sums
浏览 3提问于2021-04-30得票数 0
回答已采纳
我试图在R中重新创建Stata代码中的调查统计数据，但我无法得到同样的置信区间。我把数据按感兴趣的县减去，然后看看有多大比例的学生不戴自行车头盔，按等级划分，以及这些百分比的置信间隔是多少。在Stata代码中，我们设置了调查权重，并创建了一个带有置信区间的频率表：
svyset psu [pweight = fwt_cnty], strata(strata) singleunit(centered)
svy, subpop(if c_valencia==1): tab grade helmet, per ci format(%8.1f) nomarginals vertical
结果如下：
-
浏览 19提问于2022-10-10得票数 0
2回答
我对多重治疗数据进行了kruskal wallis测试，其中我比较了五种不同的方法。
一位朋友向我展示了spss中的计算，结果包括了每种方法的平均秩。
在R中，当将kruskal.test应用于我的数据集时，我只获得chi2、df value和p-value。这些值与spss中的值相同，但我没有得到任何排名。
我怎样才能打印出计算的等级？我的代码如下所示：
comparison <- kruskal.test(all,V3,p.adj="bon",group=FALSE, main="over")
如果我打印出比较结果，我会得到以下结果：
Kruskal
浏览 7提问于2015-03-02得票数 0
1回答
我希望运行一个混合效应模型在R的基础上，我如何运行统计在SPSS与重复的测量方差。是我在SPSS中建立重复测量方差的方法。我将如何在R中将其转换为lme4？
Key：EBT100.是任务的名称，基因型是我的IV，我的学科内因素是天(5级)和Cue (9级)。是我的DV。
在R中，这是我试图运行的代码:在R中，这是我的代码：
lmeModel <- lmer(Att ~ Genotype*Day*Cue + (1|Subject)
我的基因型效应在R和SPSS之间是相同的(p~0.12)，但我的所有交互作用都是不同的(基因型x日，基因型x Cue，基因型x型x日x Cue)。
R (lme
浏览 0提问于2018-04-15得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
