小编Cir*_*rus的帖子

我正在尝试使用ggplot2制作温度趋势图.我不能用给定的斜率和截距绘制abline(它们是Sen Estimate用于趋势计算的结果).

df是我数据的子集.

df<- structure(list(Mean_Tmin = c(-6.85888797419416, -6.61640313272608, 
-4.10521557474695, -3.33048575952967, -4.85158854541956, -6.90958783702363, 
-5.4738307882025, -5.50757613208475, -3.34384859553843, -4.79079906097484, 
-4.64407818570133, -4.17247075814977, -4.00324802670715, -2.92271555321199, 
-4.26302362769172, -4.05299360203852, -3.65186880762282, -2.6784292470159, 
-3.55928476504212, -3.26383609036539, -3.75455004498325, -4.01641178860892, 
-4.95576810887431, -3.85043378628423, -3.39114296653151, -4.62589589895125, 
-2.62611368991667, -3.60484580346817, -3.55386283157491, -3.76902382151618, 
-3.05841501472943, -1.77644594829943, -4.01934355211525, -2.35357265559614, 
-3.63012525456073, -4.61818077798782, -2.78669513101182, -1.87421430381448, 
-2.99847245406934, -3.57236610167573, -3.56963065732644, -2.8243000106372, 
-1.99565301030409, -3.3461888162997, -2.58940007000144, -3.19945356820737, 
-5.37830757225912, -2.73135885451205, -1.88970245530541, -2.54034752481066, 
-3.38038340627931, -3.13416288370415, -2.8610910675591, -2.89723228973215, 
-2.3992604730445, -2.68507391337318, -2.92441273949878, -3.33097198455173, 
-4.43334081889723, -4.53411741435393, -2.96651491236555), ID = c(1, 
 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, …

我是相对较新的R.我正在尝试从netCDF文件中获取温度数据的不同点(lat,lon)的时间序列.我的示例数据文件在这里,这里是小文件.我已经尝试过netCDF包和我到目前为止使用的代码

library(ncdf)
obsdata = open.ncdf("obs.nc")

print.ncdf(obsdata) 

obsdatadates = obsdata$dim$time$vals
obsdatadates = as.Date(obsdatadates,origin = '1950-01-01') 
obsdatadates
obsoutput = get.var.ncdf(obsdata, varid = 'tasmin', start = c(1,1,1),
                         count = c(1,1,22280))
dim(obsoutput)
datafinal=merge(obsdatadates,obsoutput)

任何人都可以帮我获取该数据的特定点(lat,lon)的时间序列数据帧(第一列)和另一个数据值.在这种情况下,我正在寻找特定纬度点(并重复多个兴趣点)和给定变量(在这种情况下)的时间序列(1950-01-01至2010-12-31,其中数据为) tasmin).非常感谢您的帮助.谢谢,拜访

如何在奇数编号的刻面ggplot中的空白区域中注释一些文本.让我们有一个刻面的ggplot,数据如下,有2行和2列.所以有2行,第2列的空白区域.

  df<- data.frame(Ara = rep("XTX", each = 3),
              Len = c(744, 750, 755),
              Mon = c("Sep", "Oct","Nov"),
              Value=c(11.224,10.15,4.23))
  df
  facetplot<-ggplot(df, aes(x=Value, y=Len, shape=Ara))+
      geom_point(size=5.0)+
      theme(legend.position = c(.7, .4), legend.direction="vertical")+
      facet_wrap(~Mon,scales="free_x", nrow=2)
  facetplot

现在我试图在空间中注释一些文本但不能(在图像中用红色书写).我正在寻找类似于legend.position的注释文本.有没有人对此有任何想法.或者可能的工作是什么.谢谢.

嗨,我有从2010年到2050年的年度数据.我正在尝试制作线图,以便我的x轴从2010年开始到2050年结束,显示5年间隔iei希望x轴断裂为(2010年,2015年,2020年......) .,2045,2050).不幸的是,它从2013年开始到2048年结束.为什么？怎么解决这个问题？我的数据,代码和结果图片如下.

plot1<- ggplot(test, aes(yr, y=value)) +
geom_line()+
scale_x_date(breaks = date_breaks("5 year"),labels=date_format("%Y")) +
scale_y_continuous(breaks=seq(-4,6, by=1))
plot1

数据

 dput(test)
structure(list(value = c(2.47099989218436, 3.09640452148661, 
1.32121989082519, 0.742309399667898, 0.180070229458727, 2.2475619117108, 
0.606470664265897, 2.12742331755353, 2.73894680327422, 2.22368873465667, 
1.58381022102847, 2.10136510397371, 1.74582199030396, 2.21689521610787, 
2.51618709804907, 1.87243814589322, 1.92029935267449, 1.79383249007138, 
0.537680017904451, 1.2415782984683, 3.62075008273724, 4.50975793125965, 
3.70660640492563, 4.16317150909305, 2.24008439109747, 2.24587596633027, 
3.63019754286973, 4.28513572439197, 3.61013179034863, 4.20010027834161, 
2.06766292535187, 4.34833637648799, 5.71460894423653, 4.12185659615561, 
3.93305702163007, 3.29384139246081, 3.2915580598453, 4.21009646693621, 
4.32889796119913, 4.99213117815761), yr = structure(c(14610, 
14975, 15340, 15706, 16071, 16436, 16801, 17167, 17532, 17897, 
18262, 18628, 18993, 19358, 19723, 20089, 20454, 20819, …

我正在尝试合并两个数据帧,如下所示.但是这给了我合并错误,如下所示.我探索并无法弄清楚我的代码中出了什么问题.为什么会出现这个错误？示例数据,代码和错误如下

dput(head(dataframe, 10))
structure(list(Year = c(1979L, 1979L, 1979L, 1979L, 1979L, 1979L, 
1979L, 1979L, 1979L, 1979L), Days = 1:10, Months = structure(c(5L, 
5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L), .Label = c("April", "Aug", 
"Dec", "Feb", "Jan", "July", "June", "March", "May", "Nov", "Oct", 
"Sep"), class = "factor"), SWE = c(201, 200, 199, 198, 197, 196, 
194, 192, 191, 190), Date = structure(c(3287, 3288, 3289, 3290, 
3291, 3292, 3293, 3294, 3295, 3296), class = "Date")), .Names = c("Year", 
"Days", "Months", "SWE", …

我有一个数据框df,想制作一个点图.我无法将所有值绘制(或显示)为点,; 20分中只显示了8分.我认为,这是因为Lat和Lon的重叠？有没有办法可以显示所有绘制的值,即20点而不仅仅是8？谢谢.

dput(df)
structure(list(Lat = c(55.5, 49.5, 49.5, 60, 58.5, 55.5, 54, 
49.5, 55.5, 52.5, 58.5, 55.5, 55.5, 55.5, 55.5, 52.5, 55.5, 49.5, 
49.5, 48, 58.5), Lon = c(-132, -126, -124.5, -139.5, -136.5, 
-132, -129, -126, -132, -127.5, -136.5, -132, -132, -132, -132, 
-127.5, -132, -126, -126, -123, -136.5), Value = c(169, 308, 
178, 110, 224, 280, 212, 365, 224, 276, 121, 159, 166, 342, 96, 
283, 51, 375, 198, 284, 201)), .Names = c("Lat", "Lon", "Value"
), …