如何将刻度长度设置为1英寸？

Question

 x<-seq(-3,3,by=0.1)
 y=0.5*x^3 
 ## Set png size to be slightly bigger than plot to take account of one line margin
 png( "~/myplot.png" , width = 600 , height = 2600 ,units="px",res=100)
 ## Set graphical parameters to control the size of the plot area and the margins
 ## Set size of plot area to 6 inches wide by 15 inches tall
 par( pin = c(6,26) )
 ## Remove margins around plot area
 par( mai = c(0,0,0,0) )
 ## Remove outer margins around plot area
 par( omi = c(0,0,0,0) )
 ## Set y-axis take use entire width of plot area (6 inches) and to have 7 tick marks (-3,-2,-1,0,1,2,3)
 par( xaxp = c(-3,3,7) )
 ## Set x-axis take use entire height of plot area (26 inches) and to have 27 tick marks (-13,-12,...,0,...11,12,13)
 par( yaxp = c(-13,13,27) )
 ## Create the plot
 plot( x , y , type = "l" , frame.plot = FALSE , axes = FALSE )
 axis( 1 , pos = 0 , at = seq( -3 , 3 , by = 1 ) , labels = seq( -3 , 3 , by = 1 ) )
 axis( 2 , pos = 0 , at = seq( -13 , 13 , by = 1 ) , labels = seq( -13 , 13 , by = 1 ) )
 text(0.5,5,expression(f(x)==frac(1,2)*x^3) )
 ## Turn PNG device off
 dev.off()

我们可以得到6英寸宽，26英寸高的图形，但x轴或y轴上的刻度单位不是1英寸;请参阅附件g1.png：

原因是R不会使用整个空间来设置x轴和y轴，R保持很小的空间;请参阅附件g2.png：

我可以png( "~/myplot.png" , width = 610 , height = 2700 ,units="px",res=100)，但我如何才能将刻度单位的物理长度设为1英寸？

Answer 1

要指定所需的确切设置，最好输出为PDF。屏幕图可以在重新调整大小等时更改。您需要设置各种图形参数来处理边距。 base图形中的绘图区域和适当的par名称可以很好地勾勒出here。此代码应该准确地提供您所要求的内容：

修改

看起来PDF设备在外面保持一条线的边界，无论如何。默认线高约0.2英寸。因此，我们将宽度和高度增加0.5英寸，我们得到的间距正好一英寸。

这是我电脑上的测量结果！

## Set pdf size to be slightly bigger than plot to take account of one line margin
pdf( "~/myplot.pdf" , width = 6.5 , height = 15.5 )

## Set graphical parameters to control the size of the plot area and the margins
## Set size of plot area to 6 inches wide by 15 inches tall
par( pin = c(6,15) )

## Remove margins around plot area
par( mai = c(0,0,0,0) )

## Remove outer margins around plot area
par( omi = c(0,0,0,0) )

## Set y-axis take use entire width of plot area (6 inches) and to have 7 tick marks (-3,-2,-1,0,1,2,3)
par( yaxp = c(0,1,7) )

## Set x-axis take use entire height of plot area (15 inches) and to have 27 tick marks (-13,-12,...,0,...11,12,13)
par( xaxp = c(0,1,27) )

## Create the plot
plot( x , y , type = "l" , frame.plot = FALSE , axes = FALSE )
axis( 1 , pos = 0 , at = seq( -3 , 3 , by = 1 ) , labels = seq( -3 , 3 , by = 1 ) )
axis( 2 , pos = 0 , at = seq( -13 , 13 , by = 1 ) , labels = seq( -13 , 13 , by = 1 ) )
text(0.5,5,expression(f(x)==frac(1,2)*x^3) )

## Turn PDF device off
dev.off()

结果：