如何使用Metafor/R去除Forest Plot中的预打印文本？
Metafor/R是一个用于进行元分析的R语言包。Forest Plot是元分析中常用的一种图形展示方式，用于展示各个研究的效应估计值及其置信区间。预打印文本是指在Forest Plot中显示的每个研究的标签或其他额外信息。
要去除Forest Plot中的预打印文本，可以使用Metafor/R提供的函数和参数进行操作。具体步骤如下：
- 导入Metafor/R包：首先需要在R环境中安装并加载Metafor包，可以使用以下命令完成：install.packages("metafor") library(metafor)
- 创建Forest Plot：使用Metafor/R提供的函数创建Forest Plot，例如使用
forest()函数。该函数可以接受多个参数，用于设置Forest Plot的各种属性，例如效应估计值、置信区间、研究标签等。 - 去除预打印文本：在创建Forest Plot时，可以通过设置相应的参数来去除预打印文本。具体取决于Forest Plot的具体实现方式和使用的参数，可以参考Metafor/R的官方文档或函数说明来了解可用的参数和选项。
- 绘制Forest Plot：完成上述设置后，使用Metafor/R提供的函数绘制Forest Plot，例如使用
plot()函数。该函数会将Forest Plot绘制在R环境的图形设备中，可以进一步调整和保存。
需要注意的是，Metafor/R是一个专门用于元分析的R语言包，其主要功能是进行元分析相关的统计计算和图形展示。对于Forest Plot中的预打印文本的具体操作和去除方式，可能需要根据具体的需求和Forest Plot的实现方式进行调整和修改。
推荐的腾讯云相关产品：腾讯云服务器（https://cloud.tencent.com/product/cvm）可以提供云计算资源，满足元分析等计算需求。
相关·内容
1回答
我对R很陌生，我正在使用metafor包的forest.default实现创建一个带有数据集的森林地块(我只有森林地块的数据中的x、ci.ub和ci.ub )
我几乎完成了，我的绘图看起来类似于，只是我的数据在图形顶部的黑色水平边框上。我做过这方面的工作，但我有一个森林地块总是在同一行的水平边框，并与白线叠加也没有帮助，正如提到的。
我的问题是:如何在数据集中添加空白，以便森林地块本身将黑色边框向上移动？
谢谢!
编辑:下面是代码片段
require(metafor)
metafor::forest(
x = as.vector(t(data_set[, forest_plot$
浏览 2提问于2017-01-04得票数 3
回答已采纳
我正在进行一个元分析在R与元the包。我遵循了关于的教程
如何在“所有研究的RE模型”中添加τ^2？
library(metafor)
### load BCG vaccine dataset
data(dat.bcg)
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, measure="OR",
slab=paste(author, year, sep=", "), method="REML")
###Plot Forest Plot
forest
浏览 2提问于2018-01-11得票数 1
回答已采纳
我正在使用metafor包的函数forest()为我的元分析创建一个森林图。除了在底部预先打印的文本“对数均值比”和“对数光盘”，我想添加文本“对数响应率”和"LnR“。我知道这两个名称意味着相同的效果大小，但为了连续性，我想在整篇论文(和绘图)中使用相同的名称。对于摘要菱形，可以通过mlab="Summary"定义文本。对于我所说的两个位置，有没有类似的代码？或者，有没有办法摆脱预先打印好的文本，手动添加我的文本？ 这是我在我的数据子集上使用的代码： Modsub <- rma.mv (yi = yi, V = vi,random= list(~ 1|ID_E
浏览 42提问于2021-11-12得票数 0
1回答
R区并排森林地块
、、
这是一些通用数据。目前，我使用元计算包和par()函数来并排绘制图。但是，我想在理想的情况下去掉第二张图中的学习栏和这两幅图中的日志RR估计。我希望这两片森林也能更紧密地联系在一起。抱歉之前没有提供数据-我是新来的。
library(metafor)
par(mfrow=c(1,2))
### load BCG vaccine data
data(dat.bcg)
### calculate log relative risks and corresponding sampling variances
dat <- escalc(measure="RR", ai=tpo
浏览 3提问于2014-06-15得票数 4
回答已采纳
2回答
目前，我正在做一项流行率的元分析。在R中有很多用于meta分析的包(即meta，metafor，MAd等)。然而，我找不到一个合适的流行率命令。你有任何关于meta分析和R的流行率的经验吗？您会使用哪个命令？
library(metafor)
res <- rma(p(logit), SE(logit), data=data, method="DL")
forest(res)
有什么想法吗？
向托蒂致以最好的问候
浏览 0提问于2017-05-05得票数 2
1回答
我正在使用metafor包在R中创建Forest图。我想将文本添加到我的图中，以使用text()函数创建标签。我想知道最简单的方法是确定我想要的文本的x，y坐标。目前，我只是猜测并查看它的外观，然后根据需要进行编辑。有没有一种方法可以在我的图上覆盖一个网格，或者一些引导我的东西(然后删除它)？谢谢!
浏览 2提问于2018-07-26得票数 2
我试图缩小R中由xlab包生成的森林地块中metafor (标题为："Log Risk Ratio")和x轴(见下图)之间的空间，但没有成功。
有什么办法可以做到这一点吗？
library(metafor)
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
forest(dat$yi, dat$vi, mgp = c(2, .3, 0)) # Tried setting `mgp` without success!
浏览 3提问于2019-11-20得票数 1
回答已采纳
1回答
我目前正在使用"metafor“包在R中进行一次分解。在做我的研究时，我发现了R中的一个不同的元分析包，即"meta“。我更喜欢后一个包创建的森林图(从设计上讲)，但不幸的是，有些数据与我用metafor创建的图不同。 具体地说，只有I^2和合并估计的数据是不同的。 meta_1 <- rma(yi=yi, vi=vi, measure="SMD", method="ML", slab=Citation, data=dat)
forest(meta_1)
meta_2 <- metagen(yi,vi^.5,data = dat,
浏览 110提问于2020-08-25得票数 1
回答已采纳
1回答
我会在一片森林中加两支箭，作为一份去日记的论文。下面是metafor::forest的演示图
require(metafor)
data(dat.bcg)
res <- rma(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg,
slab=paste(author, year, sep=", "))
forest(res, cex=.8, order=order(dat.bcg$ablat), addfit=F)
我想要的是沿着x轴添加两个箭头，比如(在红色方
浏览 2提问于2014-11-03得票数 2
回答已采纳
我在R包metafor中使用了metafor。我希望使用xlim参数来指定绘图的宽度。但是，当x值之一超出该范围时，此值将被重写。以下是两个例子：
library(metafor)
par(mfrow=c(2,1))
rmaObject <- rma.uni(yi=c(1.3, -16), sei=c(0.8, 0.7))
forest(rmaObject, atransf=exp, xlim=c(-5,12), alim=c(-5,12))$xlim
forest(rmaObject, atransf=exp, xlim=c(-5,12))$xlim
在这两种情况下，我认为函数应该按
浏览 3提问于2016-02-26得票数 0
回答已采纳
1回答
如何保存R中的非图形图？
、、、
我有一个来自metafor包的森林地块。因为这不是一个图形，所以我不能使用我最喜欢的ggsave()。我尝试过像png(filename="forest_plot_bmj.png", res=315, width=3312, height=1228)这样的函数，但是我得到的结果与我在控制台输出中看到的不同。你的建议是什么？
浏览 8提问于2022-05-09得票数 0
回答已采纳
我使用metafor软件包在R中拟合了一个稳健的随机效应元回归模型.
我的完整以及可复制的R代码如下所示。
问题：
(1)灰色钻石的含义和解释是什么？
(2)当我有主持人的时候，我不会得到一个整体的平均效果，对吗？
library(metafor)
d <- read.csv("https://raw.githubusercontent.com/izeh/m/master/d.csv", h = T) ## DATA
res <- robust(rma.uni(yi = dint, sei = SD, mods = ~es.type, data = d, slab
浏览 1提问于2019-08-15得票数 1
2回答
我用下面的代码在r中使用元‘ve包做了一个预告：
res <- metafor::rma(cohens_d, variance, data = my_data)
par(mar=c(3.4,0,0,0))
par(cex=2.5, font=4)
metafor::forest.rma(res, alim=c(-3.75, 3.75), xlab = "Cohen's D with 95% CI",
slab = my_data$Paper)
这给了我一个形象：
我想翻转x轴，意思是右边有负号。对怎么做有什么想法吗？
谢
浏览 4提问于2017-04-05得票数 1
回答已采纳
1回答
我正在使用Metafor包来绘制森林图，我想知道是否有人知道从文本注释中省略上下配置项的方法，例如，其中一个组是引用，结果固定为1。
下面是一个示例：
library(metafor)
par(mar=c(5,4,1,2))
forest(x = c(1, 0.9, 1.1),
ci.ub = c(1, 0.98, 1.18),
ci.lb = c(1, 0.82, 1.02),
refline = 1)
我只想从研究1行中删除1.00，1.00，只留下1.00。
我能想到的唯一可能是在上面画一些白色的东西，但这会很麻烦，而且我有
浏览 1提问于2018-09-23得票数 1
实际上，我希望使用元计算包来生成我的森林地块。它目前的工作非常好，除了它产生日志赔率比，而我想要纯粹的优势比本身。在元元code代码中是否有这样做的简单方法？
#Metafor library
library(metafor)
#ReadXL library to import excel sheet
library(readxl)
#Name the data sheet from the excel file
abcd<- read_excel("analysis.xlsx")
#View the data sheet with view(abcd)
dpi=
浏览 0提问于2018-03-10得票数 1
回答已采纳
我有一个关于在森林地块中显示稳健估计的快速问题。 我使用metafor包中的rma.mv()来拟合多级随机效果模型。为了获得可靠的估计，我使用clubSandwich包中的conf_int()和coef_test()函数。但是，来自metafor的forest()函数只接受rma类的对象，这些对象不是由用于健壮估计的函数生成的。 有没有办法使用forest()函数创建漂亮的森林图，同时以某种方式集成来自clubSandwich函数的健壮的总体估计？ 感谢任何建议或替代解决方案！
浏览 31提问于2021-08-13得票数 3
回答已采纳
我正在用metafor包做一个森林地块。在(http://www.metafor-project.org/doku.php/plots:forest_plot_with_subgroups)和几个论坛文章的帮助下，我写了一个代码，但不幸的是，字母是模糊的。我试着用windows ()把图做大一点，这确实有点帮助，但还没准备好把它放到纸上。它看起来是这样的：also blurry with zooming in，当我只使用效应大小的子集进行子组分析时，森林图的可读性更好，但直到大约15个效应大小为止)。由于行间距太窄，我甚至不能考虑在其间放置子组的汇总多边形。 这是我用于我的数据的代码(多级模
浏览 77提问于2021-09-30得票数 0
2回答
在metafor库的文档中，它说，如果未指定可选参数psize，则点大小与估计的精度成正比。据我所知，精度是方差的倒数。
元库如何从置信区间计算点大小？
如何只放大点大小？cex也会影响文本大小。
示例
library(metafor)
model_data <- read.table(dec = ",", text="
OR lower upper
age 0,9
浏览 5提问于2016-09-05得票数 1
回答已采纳
1回答
我正在尝试一种方法来隐藏我们在森林地块右侧看到的信息(ES(95% CI)。我正在使用R中的metafor包。
代码如下
library(metafor)
par(mar=c(4,4,1,2))
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, measure="RR",
slab=paste(author, year, sep=", "), method="REML")
forest(res, xlim=c(-16, 6), at=log(
浏览 0提问于2017-11-26得票数 2
我的数据：
dat <- structure(list(study = structure(c(2L, 6L, 7L, 4L, 1L, 3L, 5L), .Label = c("Ansari et al.", "De Santis et al.", "Dewan et al.", "Garbossa et al.", "lockney et al", "Sughrue et al. ", "Wu et al."), class = "factor"), year
浏览 1提问于2018-01-10得票数 1
回答已采纳
1回答
我正在创建一个包含3个部分的图形，这3个部分是用metafor包创建的森林地块。如何减少子图之间的空白，使它们更接近？
以下是示例数据和代码
library(metafor)
dat <- read.csv(("https://raw.githubusercontent.com/aelhak/data/main/dat2.csv"))
graphics.off()
grDevices::cairo_pdf("plot.pdf", width = 4.5, height = 6)
layout(mat = matrix(c(1, 2, 3), nrow =
浏览 31提问于2021-03-12得票数 0
回答已采纳
1回答
我需要获取R字符向量中的小于或等于符号(不是<=，即≤)。字符向量将为metafor提供"slab“标签。
所以我想
txt<-c("<=2", "3-5", "6-7",">=8")
但使用适当的小于/大于或等于符号
以及以后的
forest( ..... slab=txt ...)
所以
txt<-c(expression(1 <= 2), "3-5", "6-7",">=8")
运行正常，因为"1 <= 2
浏览 0提问于2014-01-03得票数 1
1回答
上下文:我遵循一个现有的R脚本，使用元from包中的林函数在R中创建一个林地块。
问题:情节显示，然而，研究的名字并没有显示在上面，我不知道为什么。
请看这里的数据框架。我希望author和year显示在下面的森林地块中，而不是“学习X"：
📷
见由此产生的森林地块：
📷
代码如下：
library("here") # detect root directory
library("tidyverse") # data science
library("robumeta")
library("metafor") # for
浏览 0提问于2022-02-01得票数 0
回答已采纳
1回答
我在r中使用metafor包做了一个前图：
#forestplot
par(mar=c(2.75,0,0,0))
par(font=2,mgp=c(2.3, 0.8, 0), xaxs = "i")
metafor::forest.rma(res, alim=c(-1, 3.75), xlim = c(-5, 3.5), at =
c(-1, -.5, 0 , .5, 1, 1.5, 2, 2.5, 3), xlab = "Cohen's D with 95% CI",
slab = data$Author, annotate = FALSE,
浏览 17提问于2018-03-03得票数 1
回答已采纳
当使用军官和rvg图形时，我希望我的图像没有背景白方。这样做的原因是有些函数在顶部留下了大量的空空间，这通常会与标题重叠。我可以进去删除每一个，但当制作20个或更多的情节，这变得乏味，快速。
具体来说，我想将它与元数据中的林函数一起使用，因为这里在示例中使用了这个函数。
删除背景前的.
**
在删除背景之后。
**
library(officer)
doc <- read_pptx()
doc <- add_slide(doc, "Title and Content", "Office Theme")
doc <- ph_wit
浏览 0提问于2018-08-08得票数 1
1回答
我正在研究R中的土地利用/覆被变化模拟，我有两个栅格图，一个是“土地使用/覆盖类”，另一个是“毁林风险信息”。我使用风险栅格识别森林像素(土地利用/覆盖类别之一)更有可能被砍伐森林。到目前为止，我有一个可以工作的R代码，下面是一个可以复制的例子：
#create land-use/cover raster
real <- raster(nrows=25, ncols=25, vals=round(rnorm(625, 3), 0))
real[ real > 3 ] <- 3 #just so we end with 3 land-use/cover classes
real
浏览 0提问于2018-04-27得票数 0
回答已采纳
问题陈述
我正在编写一个latex文档，并使用R的pdf()函数创建我的绘图。为了匹配我的文档的字体，我在pdf函数中指定了字体族：pdf(family = "LM Roman 12")。但是，在图(林图)中，我使用par(font = ..)更改为粗体字体，以突出显示标题。不幸的是，这已经不起作用了，因为我已经在pdf()中指定了字体系列。
示例代码
### Required packages
library(metafor)
library(extrafont)
### My dataset reproduced with: dput()
data.rct.forest
浏览 0提问于2019-06-18得票数 10
回答已采纳
1回答
我正在尝试在一个网格中绘制几条Kaplan Meijer曲线，在R中有一个共享的图例。这里(common legend in arrange_ggsurvplots) ggsurvplot_facet是建议的。然而，在我的例子中，每个图都是基于使用不同变量的不同拟合，所以ggsurvplot_facet (至少在我的理解中)在这里是不适用的。请看我下面的数据和脚本。有没有一种方法可以在分享一个图例的同时安排这些情节？我尝试过使用ggarrange和grid_arrange_shared_legend，但是我不能让它与ggsurvplot对象一起工作。 library(survival)
li
浏览 18提问于2020-01-28得票数 3
回答已采纳
1回答
我试图使用meta:: forest ()函数在R中创建一个森林图，并且在试图更改左列时存在格式化问题。总体效果统计数据消失，“作者”栏成为中心。
下面是一个例子，说明当我不更改左边列名时，森林地块是什么样子的：
您可以看到每个“子组”的最底部，为了整体效果，统计数据是shown.However，当我试图将左边的列更改为'Author‘和’度量值‘时，会将a行添加到我的代码'leftcols( )’中，这些统计数据将被移除，并且列变得中心对齐：
具体来说，我想要帮助的是：
将左列名更改为“作者”和“度量值”，而不删除整个统计信息
保持列左对齐
增加每
浏览 10提问于2020-03-29得票数 5
如何在置信区间列(3d列)上方添加标题"OR CI 95%“，并在第一列上方添加标题"Cohort”？
这是我的代码：
library(metafor)
y <- scan(what=numeric())
0.9158
0.6304
se <- scan(what=numeric())
-0.08461987
-0.165519
meta_analysis <- rma(y,se^2,method="REML")
forest.rma(meta_analysis,xlab="Odds Ratio", slab=c(&
浏览 2提问于2018-12-18得票数 0
我正在为多因素meta分析绘制一个总结森林图，并希望强调治疗方法(网络腿)之间的分类差异。虽然我喜欢metafor生成的干净的数字，但我希望能够增加用于绘制汇总效果和置信区间的空间，在不损失权重的情况下增加点的大小。 ? 如何增加绘图大小(即按比例“放大”汇总效应点和置信度区间)，使其在图形中心占据更多空间并更易于阅读？(这样做不会影响相邻的文本(注释和标签)，也不会消除权重。) 我使用了cex选项，它确实增加了绘图和汇总点的大小，但也增加了因子描述和数字注释的文本(Forest 1 cex = 1.55)。 ? 我还尝试通过psize (Forest 1 psize = 2)增加点
浏览 16提问于2019-04-18得票数 0
我有一个带有一些重要指标的RandomForestRegressor()图(RMSE、MAE和MAPE)。我已经有了图，但我需要将这些指标放在图外的文本框中。下面是它看起来的样子： sn.set_style('dark')
# plot the model
plot = sn.scatterplot(plant_target, predicted_values_plant_1)
plot.set(xlabel='Given', ylabel='Prediction',title='Usine Random Forest Cross Va
浏览 49提问于2021-06-17得票数 0
回答已采纳
假设下面的代码(如中给出的)
library(metafor)
### decrease margins so the full space is used
par(mar=c(4,4,1,2))
### fit random-effects model (use slab argument to define study labels)
res <- rma(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg, measure="RR",
slab=paste(author, year, sep=&
浏览 19提问于2016-08-25得票数 0
回答已采纳
我想用metafor来绘制ggplot (或类似的包)的元回归图。
我发现了这个网站()，它解释了如何在ggplot2和metafor的散点图上绘制报告的值和拟合模型的回归线。
然而，作者在网站上设置的代码本身并没有正常运行。似乎我必须以某种方式计算预测，但由于我是一个新手R用户，我不知道我如何能做到这一点。我不能像对待rma()的结果一样对待lm()的结果。
Model2<-rma.uni(yi,vi,mods=~Intensity+Method,method="ML",data=ROM)
summary(Model2)
in.plot<-ggplot(data
浏览 2提问于2018-01-19得票数 0
回答已采纳
我正在使用metafor包处理R中的森林图，并尝试将整个x轴(alim)向右移动以容纳ilab列。
我仍然不允许发布图像，所以我的当前绘图现在显示为文本和x轴重叠的地方：
|ilab text| |mean [ci.lb, ci.ub]|
|---measure values + ci---|
我想要这样的东西
|ilab text| |mean [ci.lb, ci.ub]|
|---measure values + ci---|
尽管
浏览 3提问于2016-04-28得票数 1
我想使用R中的forestmodel包将多个逻辑回归模型(OR，具有95%的特定预测因子CI)绘制成一个图。 使用这个包，我可以生成单独的森林地块，但我不知道如何合并它们。 # Load packages and data
library(tidyverse)
library(forestmodel)
mydata <- read.csv("https://stats.idre.ucla.edu/stat/data/binary.csv")
str(mydata)
> 'data.frame': 400 obs. of 4 variables
浏览 93提问于2021-04-19得票数 0
我想知道是否可以(以某种形式)访问元数据包中的-林命令中显示的信息。
我正在检查/验证结果，我希望生成值的输出。到目前为止，计算结果都已确认，但我希望它们可以用于打印、保存等，而不是手工打印出来。
示例代码如下：
es <- read.table(header=TRUE, text = "
b se_b
0.083 0.011
0.114 0.011
0.081 0.013
0.527 0.017
" )
library(metafor)
es.est <- rma(yi=b, sei=se_b, dat=es, method="DL")
浏览 2提问于2017-01-25得票数 1
回答已采纳
? from sklearn.metrics import confusion_matrix
from sklearn.metrics import plot_confusion_matrix
print('*** Test Accuracy ***',forest.score(X_test,y_test))
disp = plot_confusion_matrix(forest, X_test, y_test,
display_labels=[0,1],
浏览 42提问于2020-01-31得票数 4
回答已采纳
在使用R中的metafor创建森林地块时，由于对这个包的了解有限，我很难对这个图进行一些美学上的更改。谢谢!
，我想删除图底部的总体估计，用红色框表示。这是我正在使用的代码，res是一个rma.uni对象。
forest(res, at=c(-0.5, -1, 0, 0.5, 1), xlim=c(-16,6), atransf = mytransf,
ilab=cbind(Z0, Z1, ZD0, ZD1), ilab.xpos=c(-9.5,-8,-6,-4.5),
cex=.75,
header="Subgroups",
浏览 6提问于2022-10-03得票数 0
1回答
我正在尝试在R中绘制时间序列图，一切似乎都正常(我所有的标签都正确显示，轴按我想要的方式进行索引，我就是不能得到任何点来显示。我的图表看起来完全是空白的。
为什么会发生这种情况？
我的代码如下：
forest <- c(0.096999565, 0.099238207, 0.101582401, 0.103624695, 0.105185874, 0.106321613, 0.107100639, 0.107698571, 0.108253954)
armenia.forest <- data.frame(forest)
plot.ts(armenia.forest, xlim
浏览 14提问于2018-05-31得票数 0
回答已采纳
1回答
? 我正在学习AI/ML，并试图从这个示例表单中获取文本。 import cv2
import pytesseract
pytesseract.pytesseract.tesseract_cmd = r'C:\Users\Pranav\AppData\Local\Programs\Tesseract-OCR\tesseract.exe'
image = cv2.imread('image2.png')
gray = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY)
blur = cv2.GaussianBlur(gra
浏览 34提问于2021-04-15得票数 0
回答已采纳
1回答
我已经编写了一个SQL查询，尝试按森林类型获取一组树的一些信息。但是，当我运行此查询时，Count(T_Plots.ID_plot) AS [Plots/Forest Type]返回与Count(T_Trees.ID_Tree) AS [Trees/Forest Type]相同的数字。这些数字太高了，因为整个T_Plots表只有175条记录，而在结果表中，它返回的结果高达4290条。 我忽略了什么导致这些错误的数字，以及如何获得每种森林类型的正确地块数量？ SELECT
T_Plots.Forest_type,
Count(T_Plots.ID_plot) AS [Plots/Forest
浏览 11提问于2021-09-01得票数 0
我们在R中使用包"meta“和"metafor”执行元分析。这工作得很好，但是当我们制作一个图形时，它不适合屏幕。文章的名称正在下降，上面和最下面的行是看不见的。我们已经尝试过减小字体大小等，但这并不能解决问题。
代码：
subganal.study.qual=rma(yi, vi, data=ies.logit.qual.nomis, mods=~study.qual, method="DL")
pes.summary.study.qual.nomis=metaprop(Positive, Total, study, data=db.nomis.qual, s
浏览 5提问于2020-08-17得票数 2
我想知道如何定义，如果多边形有一个开边或没有。我认为，如果多边形完全被邻居包围，它就没有开边。
使用奇妙的poly2nb(fc)，我可以得到邻居的列表:但是，从这个列表中，我不知道有多少邻居必须被邻居完全包围？情况如下：
在这两种情况下，我的中心red多边形都有3个邻居，但是有开边(左)，或者被邻居完全包围(右)。如果使用raster格式和皇后case，则完全包围单元需要8个邻居。如果它少了，它就是开放的细胞。但是，我可以从poly2nb(fc)和nb对象获得类似的信息吗？当然，我的数据可以包含单个多边形之间的条条和空白，所以我不想完全依赖于重叠的边缘或其他东西。
我的真实数据可以在
浏览 0提问于2020-01-27得票数 2
回答已采纳
2回答
假设以下代码(如Viechtbauer，2010年所示)：
library(metafor)
data("dat.bcg", package = "metafor")
dat <- escalc(measure = "RR", ai = tpos, bi = tneg, ci = cpos, di = cneg, data = dat.bcg, append = TRUE)
res <- rma(ai = tpos, bi = tneg, ci = cpos, di = cneg, data = dat, measure = "
浏览 5提问于2016-08-10得票数 5
回答已采纳
1回答
我想使用ggplot2包绘制一个森林图，并且我对我的输出很满意(请参阅下面的森林图)。
这幅图的特征是回归模型中给定变量的水平(优势比和置信度内部)，以及参考水平。
问题是需要大量的体力劳动才能生成这个地块。
First problem，我希望引用级别跟随图中给定变量的其他级别，所以我手动输入了每个这样的引用级别(参见下面的foresttable )。为了使ggplot2工作，我输入了参考水平的任意负比特率和置信区间值，然后将地块限制设置为从零到大的正数。
第二个问题，因为我的原始变量在一个列中，所以我手动输入了颜色，这很费时。
还有更简单的方法来生成这样的情节吗？任何帮助都将不胜感激。
#
浏览 1提问于2018-04-06得票数 2
回答已采纳
我正在努力处理一个查询，它试图在“select”语句中从列中的值的一个column....based创建3个单独的列。
我的例子：
我有两张桌子
forest.plots
columns = id(PK), a, b
forest.plot_measurements
columns = id(PK), plot_id (FK), measure_type (int), value_real
当plot_measurement具有measure_type =1时，这是斜率测量；如果measure_type =2，则它是wet_weight测量
期望的结果是有一个具有标题的表：
plot_id,
浏览 1提问于2014-12-04得票数 0
2回答
使用R制作两个组的自定义森林地块
、、、、
我真的很想用R制作一个类似这样的情节，我已经尝试过来自metafor包的metafor命令，但似乎无法创建任何类似的东西。我也有两个组，多个分类变量和百分比-基本上，完全相同的数据显示在这幅图。有人能给我建议如何用R来解决这个问题吗？
浏览 4提问于2015-02-10得票数 1
回答已采纳
我是元分析流行的药物使用rma.glmm命令在metafor包。由于它们遵循二项分布，所以使用log-变换。
所以，我的模型是：
med_A<-rma.glmm(xi=A, ni=Sample_size, measure="PLO")
森林地块由以下代码组成，其结果附在下面。
forest(med_A, atransf = transf.ilogit, ilab = cbind(A, Sample_size),
ilab.xpos = c(-7,-5), slab = Study, refine=1)
比例尺看上去很奇怪，因为1.00出现了两次，而且分布也
浏览 2提问于2016-01-27得票数 0
回答已采纳
我想用RandomForest模型的Dalex创建解释器，但是当我尝试的时候，我得到了这个错误： Error: $ operator is invalid for atomic vectors 有人知道如何解决这个问题吗？ 下面是我的代码： library(DALEX)
library(ranger)
model <- ranger(Species ~ .,data = iris,probability =T,classification = T)
explanation <- explain(model,
data=iris,
浏览 19提问于2021-01-15得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
