- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在尝试在泰勒图中绘制多个模型,并且在代码方面遇到了一些困难。我已经设法生成图表(参见图片),但无法弄清楚如何减少轴,因为它们太大,标准化轴以标记为 1,2,3,4 并在相关性上添加刻度线 - 用刻度线标记 我希望每 0.1 有一次主要刻度,每 0.05 到 0.9 有一次次要刻度,之后我尝试在此时将主要刻度设置为 0.95,将次要刻度设置为每 0.01(如果这有意义的话)。任何有关上述内容的帮助/建议都会有帮助。我已经在“plotrix”包中使用了“taylor.diagram”(并阅读了它的指南 - 但我对 R 相对缺乏经验)并附加了我的(有些基本的)代码到目前为止,但我的情节看起来相当困惑。谢谢
all.models <- as.data.frame(cbind(Sy.One, Sy.Two, Sy.Three, Sy.Four, Sy.Five, Sy.Six, Sy.Seven, Sy.Eight, Sy.Nine, Sy.Ten))
taylor.diagram(CSR, Sy.One, sd.arcs=T, ref.sd=T, pcex=1.5, main=NULL, pos.cor=F,
xlab="Standard Deviation (cm)", ylab="Standard Deviation (cm)")
for (i in 1:dim(all.models)[2]) {
model.wanted <- all.models[,i]
taylor.diagram(CSR, model.wanted, sd.arcs=T, ref.sd=T, pcex=1.5, col=i, add=T, pos.cor=F)}
# Add legend
model.names <- c("Sy=1%","Sy=2%","Sy=3%","Sy=4%","Sy=5%","Sy=6%","Sy=7%","Sy=8%","Sy=9%","Sy=10%")
legend("top", model.names, pch=19, col=i, cex=1.0, bty="n", ncol=5)
最佳答案
一种选择是使用不同的软件包,例如 openair
,它可能更灵活。由于您的特定要求,使用为您的要求设计的代码可能会更容易。我编写了一些代码来生成以下图,该图接近您想要的图。您可以破解代码以将绘图调整为您想要的格式。
# code to make a Taylor diagram
# formulas found in http://www-pcmdi.llnl.gov/about/staff/Taylor/CV/Taylor_diagram_primer.pdf
# and http://rainbow.llnl.gov/publications/pdf/55.pdf
# correlations and tick marks (only major will have a line to the center)
# minor will have a tick mark
correlation_major <- c(seq(-1,1,0.1),-0.95,0.95)
correlation_minor <- c(seq(-1,-0.95,0.01),seq(-0.9,9,0.05),seq(0.95,1,0.01))
# test standard deviation tick marks (only major will have a line)
sigma_test_major <- seq(1,4,1)
sigma_test_minor <- seq(0.5,4,0.5)
# rms lines locations
rms_major <- seq(1,6,1)
# reference standard deviation (observed)
sigma_reference <- 2.9
# color schemes for the liens
correlation_color <- 'black'
sigma_test_color <- 'blue'
rms_color <- 'green'
# line types
correlation_type <- 1
sigma_test_type <- 1
rms_type <- 1
# plot parameters
par(pty='s')
par(mar=c(3,3,3,3)+0.1)
# creating plot with correct space based on the sigma_test limits
plot(NA
,NA
,xlim=c(-1*max(sigma_test_major),max(sigma_test_major))
,ylim=c(-1*max(sigma_test_major),max(sigma_test_major))
,xaxt='n'
,yaxt='n'
,xlab=''
,ylab=''
,bty='n')
#### adding sigma_test (standard deviation)
# adding semicircles
for(i in 1:length(sigma_test_major)){
lines(sigma_test_major[i]*cos(seq(0,pi,pi/1000))
,sigma_test_major[i]*sin(seq(0,pi,pi/1000))
,col=sigma_test_color
,lty=sigma_test_type
,lwd=1
)
}
# adding horizontal axis
lines(c(-1*max(sigma_test_major),max(sigma_test_major))
,c(0,0)
,col=sigma_test_color
,lty=sigma_test_type
,lwd=1)
# adding labels
text(c(-1*sigma_test_major,0,sigma_test_major)
,-0.2
,as.character(c(-1*sigma_test_major,0,sigma_test_major))
,col=sigma_test_color
,cex=0.7)
# adding title
text(0
,-0.6
,"Standard Deviation"
,col=sigma_test_color
,cex=1)
#### adding correlation lines, tick marks, and lables
# adding lines
for(i in 1:length(correlation_major)){
lines(c(0,1.02*max(sigma_test_major)*cos(acos(correlation_major[i])))
,c(0,1.02*max(sigma_test_major)*sin(acos(correlation_major[i])))
,lwd=2
,lty=correlation_type
,col=correlation_color
)
}
# adding minor tick marks for correlation
for(i in 1:length(correlation_minor)){
lines(max(sigma_test_major)*cos(acos(correlation_minor[i]))*c(1,1.01)
,max(sigma_test_major)*sin(acos(correlation_minor[i]))*c(1,1.01)
,lwd=2
,lty=correlation_type
,col=correlation_color
)
}
# adding labels for correlation
text(1.05*max(sigma_test_major)*cos(acos(correlation_major))
,1.05*max(sigma_test_major)*sin(acos(correlation_major))
,as.character(correlation_major)
,col=correlation_color
,cex=0.5)
# adding correlation title
text(0
,max(sigma_test_major)+0.5
,"Correlation"
,col=correlation_color
,cex=1)
#### adding rms difference lines
# adding rms semicircles
for(i in 1:length(rms_major)){
inds <- which((rms_major[i]*cos(seq(0,pi,pi/1000))+sigma_reference)^2 + (rms_major[i]*sin(seq(0,pi,pi/1000)))^2 < max(sigma_test_major)^2)
lines(rms_major[i]*cos(seq(0,pi,pi/1000))[inds]+sigma_reference
,rms_major[i]*sin(seq(0,pi,pi/1000))[inds]
,col=rms_color
,lty=rms_type
,lwd=1
)
}
# adding observed point
points(sigma_reference
,0
,pch=19
,col=rms_color
,cex=1)
# adding labels for the rms lines
text(-1*rms_major*cos(pi*rms_major/40)+sigma_reference
, rms_major*sin(pi*rms_major/40)
,as.character(rms_major)
,col=rms_color
,cex=0.7
,adj=1)
# adding title
text(0
,-1.5
,'Centered RMS Difference'
,col=rms_color
,cex=1
,adj=0.5)
###################### adding points #####################
names <- paste("model",seq(1,8),sep='')
correl_names <- seq(-0.6,0.8,by=0.2)
std_names <- seq(2,4,by=0.26)
color_names <- topo.colors(length(names))
points(std_names*cos(acos(correl_names))
,std_names*sin(acos(correl_names))
,col=color_names
,pch=19
,cex=1.5)
# making legend
par(xpd=TRUE)
legend(0,-2
,names
,pc=19
,col=color_names
,ncol=3
,bty='n'
,xjust=0.5)
关于R - 泰勒图绘制,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/24999338/
我学习 SDL 二维编程已有一段时间了,现在我想创建一个结合使用 SDL 和 OpenGL 的程序。我是这样设置的: SDL_Init(SDL_INIT_VIDEO); window = SDL_Cr
尝试查找可在地块中使用的不同类型项目的列表 来自不同样本的投影类型: projection = list(type = "equirectangular") projection = list(typ
我正在尝试使用 Java Graphics API 绘制 GIF,但无法使用下面的代码成功绘制 GIF。仅绘制 GIF 的第一张图像或缩略图,但不播放。 public void paintCompon
我目前正在使用 JFrame 并尝试绘制一个矩形,但我不知道如何执行代码 paint(Graphics g),如何获取 Graphics 对象? package com.raggaer.frame;
这个领域的新手,希望得到一些帮助。 我有一个"Missile.java" 类,我在那里画东西。我想绘制一个 ImageView,我正在使用以下代码: ImageView v = (ImageView)
下面列出了圆形的例子 这是我的 JavaScript 代码。 最佳答案 假设您的 randomColor 是正确的,您只需要: 从 canvas.onclick 中移除 context.clearR
我在绘制和缩放 ImageView 时遇到问题。请帮帮我.. 当我画一些东西然后拖动或缩放图像时 - 绘图保留在原处,如您在屏幕截图中所见。而且我只需要简单地在图片上绘图,并且可以缩放和拖动这张图片。
我们可以在形式之外绘制图像和文本...我的意思是在字面上... 我知道问这个问题很愚蠢但是我们能不能... 最佳答案 您可以通过创建表单并将其 TransparentColor 属性设置为背景色来“作
我在绘制/布局期间收到 3 个对象分配警告 super.onDraw(canvas); canvas.drawColor(Color.WHITE); Paint textPaint = new Pai
我有一个示例时间序列数据框: df = pd.DataFrame({'year':'1990','1991','1992','1993','1994','1995','1996',
我试图想出一种简洁的方法来绘制 R 数据框中所有列的 GridView 。问题是我的数据框中既有离散值又有数值。为简单起见,我们可以使用 R 提供的名为 iris 的示例数据集。我会使用 par(mf
我有一个由 10 列和 50 行组成的 data.frame。我使用 apply 函数逐列计算密度函数。现在我想绘制我一次计算的密度。 换句话说,而不是绘图... plot(den[[1]]) plo
我想知道我们如何才能在第一个和第二个组件之外绘制个人,如下所示: 最佳答案 这可能有效: pc.cr <- princomp(USArrests, cor = TRUE) pairs(pc.cr$lo
我是Pandas和matplotlib的新手,想绘制此DataFrame season won team matches pct_won 0 20
我正在尝试为 distplot 子图编写一个 for 循环。 我有一个包含许多不同长度列的数据框。 (不包括 NaN 值) fig = make_subplots( rows=len(asse
我想创建一个具有密度的 3d 图。 我使用函数 density 首先为特定的 x 值创建一个二维图,然后该函数创建密度并将它们放入 y 变量中。现在我有第二组 x 值并将其再次放入密度函数中,然后我得
全部, 我一直在研究全局所有 MTB 步道的索引。我是 Python 人,所以对于所有涉及的步骤,我都尝试使用 Python 模块。 我能够像这样从 OSM 立交桥 API 中获取关系: from O
我正在使用 e1071 包中的支持向量机对我的数据进行分类,并希望可视化机器实际如何进行分类。但是,在使用 plot.svm 函数时,出现无法解决的错误。 脚本: library("e1071") d
我制作了以下图表,它是使用 xts 对象创建的。 我使用的代码很简单 plot(graphTS1$CCLL, type = "l", las = 2, ylab = "(c)\nCC for I
在绘制状态图时,您如何知道哪些状态放在框中,哪些状态用于转换箭头?我注意到转换也是状态。 我正在查看 this page 上的图 1 : 最佳答案 转换不是状态。转换是将对象从一种状态移动到下一种状态
我是一名优秀的程序员,十分优秀!