Tag: 情节

从Python中的数据点中find移动平均数

我再次玩Python,并且find了一个带有例子的整洁的书。 其中一个例子是绘制一些数据。 我有一个两列的.txt文件,我有数据。 我绘制的数据很好,但在练习中说:进一步修改程序来计算和绘制数据的运行平均值,定义如下: $Y_k=\frac{1}{2r}\sum_{m=-r}^r y_{k+m}$ 在这种情况下r=5 (而y_k是数据文件中的第二列)。 让程序在同一个图上绘制原始数据和运行平均值。 到目前为止,我有这样的: from pylab import plot, ylim, xlim, show, xlabel, ylabel from numpy import linspace, loadtxt data = loadtxt("sunspots.txt", float) r=5.0 x = data[:,0] y = data[:,1] plot(x,y) xlim(0,1000) xlabel("Months since Jan 1749.") ylabel("No. of Sun spots") show() 那么如何计算总和? 在Mathematica中它很简单,因为它是符号操作(例如Sum [i,{i,0,10})),但是如何在python中计算每一个十点的数据并对其求平均,直到结束点数? 我看了一下这本书,却发现什么都不能解释这个: heltonbiker的代码做了诡计^^:D from __future__ import division […]

将多个alignment的图块放在一个页面上时,避免浪费空间

我想把四个地块放在一个页面上。 轴标签只能在边框上打印,也就是只有底部的X轴标签,左侧的Y轴标签。 这既适用于整个轴的名称,也适用于单独的刻度标记。 我可以使用下面的代码沿着这些行生成一些东西: pdf(file = "ExampleOutput.pdf", width = 6.61, height = 6.61, pointsize = 10 ) set.seed(42) catA <- factor(c("m100", "m500", "m1000", "m2000", "m3000", "m5000")) catB <- factor(20:28) samples <- 100 rsample <- function(v) v[ceiling(runif(samples, max=length(v)))] Tab <- data.frame(catA = rsample(catA), catB = rsample(catB), valA = rnorm(samples, 150, 8), valB = pmin(1,pmax(0,rnorm(samples, 0.5, 0.3)))) […]

情节传说没有边界和白色背景

我在一个情节中有一个传说(从一个基本语句)经过它。 我怎样才能实现在传说附近隐藏的abline? 这应该可以通过设置图例背景白色,没有边界,但我怎么能做到这一点? 假设图表应该是这样的: windows.options(width=30, height=12) plot(1:10) abline(v=seq(1,10,1), col='grey', lty='dotted') legend(4.8, 3, "This legend text should not be disturbed by the dotted grey lines") 如果图例干扰点图的点:我怎样才能实现在图例附近的图例变得不可见(如上),但点仍然可见? windows.options(width=30, height=12) plot(1:10) abline(v=seq(1,10,1), col='grey', lty='dotted') legend(1, 5, "This legend text should not be disturbed by the dotted grey lines, but the plotted dots should still be visible") 最后:有没有办法在图例中引入换行符?

在R中绘制3D数据

我有一个3D数据集: data = data.frame( x = rep( c(0.1, 0.2, 0.3, 0.4, 0.5), each=5), y = rep( c(1, 2, 3, 4, 5), 5) ) data$z = runif( 25, min = (data$x*data$y – 0.1 * (data$x*data$y)), max = (data$x*data$y + 0.1 * (data$x*data$y)) ) data str(data) 我想绘制它,但R alwyas的内置函数给出了错误 预期增加“x”和“y”值 # ### 3D Plots ###################################################### # built-in function […]

网格在R图中

有一个命令,可以很容易地添加一个网格到R的情节?

八度不会阴谋

当我尝试在GNU Octave上绘制一个图表,并尝试使用plot时,它给了我下面的输出 set terminal aqua enhanced title "Figure 1" size 560 420 font "*,6" dashlength 1 ^ line 0: unknown or ambiguous terminal type; type just 'set terminal' for a list 我正在使用Mac OS X 10.9.2。 我曾尝试使用 octave:79> setenv("GNUTERM","X11") 但我仍然得到相同的错误。

如何在matplotlib中closuresimshow()的模糊效果?

我想做一个概率的颜色图,然而imshow会产生有零概率的点的模糊值。 我怎样才能摆脱围绕真实网格点的模糊周边? 例: import numpy as np import matplotlib.pyplot as plt a=np.asarray([[ 0.00000000e+00 , 1.05824446e-01 , 2.05086136e-04, 0.00000000e+00], [ 1.05824446e-01 , 3.15012305e-01 , 1.31255127e-01 , 1.05209188e-01], [ 2.05086136e-04 , 1.31255127e-01 , 0.00000000e+00 , 0.00000000e+00], [ 0.00000000e+00 ,1.05209188e-01 , 0.00000000e+00 , 0.00000000e+00]]) im=plt.imshow(a,extent=[0,4,0,4],origin='lower',alpha=1,aspect='auto') plt.show()

R图:大小和分辨率

我已经堆积到这个问题:我需要用DPI = 1200和具体的打印尺寸来绘制图像。 默认情况下,PNG看起来不错… png("test.png",width=3.25,height=3.25,units="in",res=1200) par(mar=c(5,5,2,2),xaxs = "i",yaxs = "i",cex.axis=1.3,cex.lab=1.4) plot(perf,avg="vertical",spread.estimate="stddev",col="black",lty=3, lwd=3) dev.off() 但是当我应用这个代码时,图像变得非常糟糕,它不能缩放(适合)所需的大小。 我错过了什么? 如何“适合”形象的情节? ,

如何指定实际的x轴值作为R中的x轴刻度

我正在创build一个R的情节,我不喜欢由R绘制的X轴值。 例如: x <- seq(10,200,10) y <- runif(x) plot(x,y) 这在X轴上绘制了具有以下值的图表: 50, 100, 150, 200 但是,我想绘制variablesx中存储的10,20, 30 … 200个值作为X轴值。 我已经搜遍了无数的博客和简短的手册 – 经过几个小时的search,最接近我find有用的是以下(总结)的指示: 调用plot()或par() ,指定参数xaxt='n' call axis()例如axis(side = 1, at = seq(0, 10, by = 0.1), labels = FALSE, tcl = -0.2) 我试了一下,结果绘图根本没有X轴值。 有人可能知道如何做到这一点? 我不敢相信以前从来没有人尝试过这样做。

创build特定大小的绘图窗口

我怎样才能创build一个新的屏幕R绘图窗口与特定的宽度和高度(像素等)?