2014-12-24 32 views

cevap

13

Ayrıca plotrix paketinde arctext deneyebilirsiniz: daha yüksek özelleştirme fırsatları için

library(plotrix) 

# set up a plot with a circle 
plot(x = 0, y = 0, xlim = c(-2, 2), ylim = c(-2, 2)) 
draw.circle(x = 0, y = 0, radius = 1) 

# add text 
arctext(x = "wrap some text", center = c(0, 0), radius = 1.1, middle = pi/2) 
arctext(x = "counterclockwise", center = c(0, 0), radius = 1.1, middle = 5*pi/4, 
     clockwise = FALSE, cex = 1.5) 
arctext(x = "smaller & stretched", center = c(0, 0), radius = 1.1, middle = 2*pi , 
     cex = 0.8, stretch = 1.2) 

enter image description here

(bir ifade; güzel vinyet bakınız), bir görünüme sahip olabilir circlize paketinde. facing = "bending"'u circos.text'da ayarlayarak, metin bir daire etrafında sarar.

library(circlize) 

# create some angles, labels and their corresponding factors 
# which determine the sectors 
deg <- seq(from = 0, to = 300, by = 60) 
lab <- paste("some text", deg, "-", deg + 60) 
factors <- factor(lab, levels = lab) 

# initialize plot 
circos.par(gap.degree = 10) 
circos.initialize(factors = factors, xlim = c(0, 1)) 
circos.trackPlotRegion(ylim = c(0, 1)) 

# add text to each sector 
lapply(factors, function(deg){ 
    circos.updatePlotRegion(sector.index = deg, bg.col = "red") 
circos.text(x = 0.5, y = 0.5, labels = as.character(deg), facing = "bending") 
}) 
circos.clear() 

enter image description here

Güncelleme

:
circlize version 0.2.1 yılında circos.text iki yeni seçenek vardır: bending.inside orijinal bending ve bending.outside ile aynıdır ( the vignette incir 11 bakınız). Böylece, bending.outside kullanılarak arsa alt yarısında metni çevirmek kolaydır:

circos.par(gap.degree = 10) 
circos.initialize(factors = factors, xlim = c(0, 1)) 
circos.trackPlotRegion(ylim = c(0, 1)) 

lapply(factors[1:3], function(deg){ 
    circos.updatePlotRegion(sector.index = deg, bg.col = "red") 
    circos.text(x = 0.5, y = 0.5, labels = as.character(deg), facing = "bending.outside") 
}) 

lapply(factors[4:6], function(deg){ 
    circos.updatePlotRegion(sector.index = deg, bg.col = "red") 
    circos.text(x = 0.5, y = 0.5, labels = as.character(deg), facing = "bending.inside") 
}) 
circos.clear() 

enter image description here

+2

Keşke bunu iki kez yapabilirim çünkü 'circlize' düzenlemesi kafasındaki çiviye çarptı. –

+0

Teşekkürler @RichardScriven. 'circlize' gerçekten güzel bir paket. Yazar vinyetlerle harika bir iş çıkardı. Şerefe. – Henrik

+0

Bunu bir süredir bakmadım, fakat bunun genel bir eğriyle (cevabımdaki “günah” örneğim gibi) çalışmayacağını düzelttim mi? –

25

Evet, ve burada kod ücretsizdir, bedavadır :-). Bunu bir süre önce yazdım ama hiç bir CRAN paketinde yayınlamadığını düşünmüyorum.

# Plot symbols oriented to local slope. 
# Interesting problem: if underlying plot has some arbitrary aspect ratio, 
# retrieve by doing: Josh O'B via SO: 
# myasp <- with(par(),(pin[2]/pin[1])/(diff(usr[3:4])/diff(usr[1:2]))) 
# so make that the default value of argument 'asp' 
# Default is 'plotx' is vector of indices at which to 
# plot symbols. If is_indices=FALSE, only then turn to splinefun to 
# calculate y-values and slopes; and user beware. 
# 
# 6 Feb 2014: added default col arg so can stick in a color vector if desired 
# TODO 
# 
slopetext<-function(x,y,plotx, mytext, is_indices=TRUE, asp=with(par(), (pin[1]/pin[2])*(diff(usr[3:4])/diff(usr[1:2]))),offsetit= 0, col='black', ...) { 
if (length(x) != length(y)) stop('data length mismatch') 
if (!is.numeric(c(x,y,plotx))) stop('data not numeric') 
if(is_indices) { 
    # plotting at existing points. 
    if(any(plotx<=1) | any(plotx>= length(x))) { 
     warning("can't plot endpoint; will remove") 
     plotx<-plotx[(plotx>1 & plotx<length(x))] 
    } 
    lows<-plotx-1 
    highs<-plotx+1 
    # then interpolate low[j],high[j] to get slope at x2[j] 
    slopes <- (y[highs]-y[lows])/(x[highs]-x[lows]) #local slopes 
    # sign(highlow) fix the rotation problem 
    angles <- 180/pi*atan(slopes/asp) + 180*(x[lows] > x[highs]) 
    intcpts <- y[highs]-slopes*x[highs] 
    ploty <- intcpts + x[plotx]*slopes 
    # change name, so to speak, to simplify common plotting code 
    plotx<-x[plotx] 
    }else{ 
    #interpolating at plotx values 
     if (any(plotx<min(x)) | any(plotx>max(x))) { 
      warning("can't plot extrapolated point; will remove") 
      plotx<-plotx[(plotx>min(x) & plotx<max(x))] 
     } 
     spf<-splinefun(x,y) 
     ploty<-spf(plotx) 
     angles <- 180/pi * atan(spf(plotx,1)/asp) #getting first deriv, i.e. slope 
    } #end of else 
xlen<-length(plotx) # better match y and mytext 
# The trouble is: srt rotates about some non-centered value in the text cell 
# Dunno what to do about that. 
dely <- offsetit*cos(angles) 
delx <- offsetit*sin(angles) 
# srt must be scalar 
mytext<-rep(mytext,length=xlen) 
col <- rep(col,length=xlen) 
for (j in 1:xlen) text(plotx[j], ploty[j], labels=mytext[j], srt= angles[j], adj=c(delx,dely),col=col[j], ...) 
} 

Düzenleme:

x <- 1:100 
y <- x/20 + sin(x/10) 
plot(x,y,t='l') 
slopetext(x=x,y=y,plotx=seq(10,70,by=10),mytext=letters[1:8]) 

Bu örnekte, üçüncü bağımsız değişken metin yerleştirilmesi için (x, y) her onuncu değeri seçilir: David mükemmel öneri, bir örnek olay başına. is_indices=FALSE vakasına sadık kalmadığımı ve spline fit'in aşırı durumlarda metinlerinizi komik yollara yerleştirebileceğimi uyarmalıyım.

enter image description here

+2

Belki kullanımı ve üretiminin bazı örnek eklemek? Yakında daha kolay olurdu :) –

+1

@DavidArenburg Sonuçların ilanı için teşekkürler! –

İlgili konular