小编Dom*_*nik的帖子

我有一堆4个栅格.我想要一个像素和它的8个邻居之间的平均时间相关性.

一些数据:

library(raster)  

r1=raster(matrix(runif(25),nrow=5))
r2=raster(matrix(runif(25),nrow=5))
r3=raster(matrix(runif(25),nrow=5))
r4=raster(matrix(runif(25),nrow=5))
s=stack(r1,r2,r3,r4)

所以对于位置x的像素,在NE,E,SE,S等位置有8个邻居,我想要平均值

cor(x,NE)
cor(x,E)
cor(x,SE)
cor(x,S)
cor(x,SW)
cor(x,W)
cor(x,NW)
cor(x,N)

以及在结果栅格中保存在位置x的平均值.边缘单元将是NA,或者如果可能的话,标记用于计算与其接触的单元(3或5个单元)的平均相关性.谢谢!

如何区分4种不同的因素(不使用大小)？是否可以使用空心点和实心点来区分ggplot2中的变量？

test=data.frame(x=runif(12,0,1),
     y=runif(12,0,1),
     siteloc=as.factor(c('a','b','a','b','a','b','a','b','a','b','a','b')),
     modeltype=as.factor(c('q','r','s','q','r','s','q','r','s','q','r','s')),
     mth=c('Mar','Apr','May','Mar','Apr','May','Mar','Apr','May','Mar','Apr','May'),
     yr=c(2010,2011,2010,2011,2010,2011,2010,2011,2010,2011,2010,2011))

其中x是观测值,y是建模结果,我想比较几个因素的不同模型版本.谢谢!

如何使用相对路径包含外部图形？

我开始我的 rmd 文件：

```{r setup,include=F}
knitr::opts_chunk$set(cache=F,echo=F)
knitr::opts_knit$set(root.dir = '/Users/dosc3612/odrive/CUDrive/real-time_interpolation',
child.path='/Users/dosc3612/odrive/CUDrive/real-time_interpolation')
```

然后包括：

# Heading 1
some text referring to figure 1. ![Figure 1. a map from qgis](data/gis/fig1alt-location_map2.png)

我的 rmd 文件在 /Users/dosc3612/odrive/CUDrive/real-time_interpolation/reports/real-time 中，所以我的理解是它将从该目录中查找数字。但基于块选项帮助我认为我可以使用 child.path 在文档中设置相对路径，但事实并非如此。

“data/gis”是“/Users/dosc3612/odrive/CUDrive/real-time_interpolation”中的子文件夹

谢谢

我有一个数组列表,其中每个数组只有2个带数字的元素,否则为NA.所有阵列中的任何给定元素只有1个数字.

这是一些示例数据.

ar1=array(NA,dim=c(2,3,4))
ar1[1,1,1]=100
ar2=array(NA,dim=c(2,3,4))
ar2[2,3,4]=200
ar3=array(NA,dim=c(2,3,4))
ar3[2,1,4]=300
ar=list(ar1,ar2,ar3)

谢谢!

我有一个数据集,我已经提取了事件发生的日期.日期的格式为MMDDYY,虽然MatLab没有显示前导零,所以它经常是MDDYY.

有没有找到平均值或中位数(我可以使用)日期的方法？median当有奇数天时工作正常但是对于偶数,我相信它是平均两个不产生合理值的中间数.我一直在尝试将日期转换为MatLab格式regexp并将其重新组合在一起,但我还没有让它工作.谢谢

dates=[32381 41081  40581  32381  32981 41081   40981  40581];

我正在尝试找到一个适合我的数据的方程式.我认识到形状,y=-exp(x)但nls(y~-a*exp(x*b))各种参数起始值失败.y是负的,因此log(y)~log(a)+ b x的"简单"拟合效果不佳.我尝试了log(y + 2)~a + b x来使所有东西都正面但是没有产生正确的形状.

有人能为这些数据提供一些帮助吗？谢谢!

Y = swediff

X = avgdate

dat2=structure(list(swediff = c(0.0379635202678687, 0.0845477936160927, 
0.146010217481196, 0.0416237104326292, 0.0659140490644253, 0.134535534695029, 
0.0095147654468483, 0.238456044233877, 0.276025694437364, 0.29435448415394, 
0.00301157777812485, 0.19171002685605, 0.277759059448242, 0.00400780564144798, 
0.342605838471721, 0.236804884903432, 0.151048712082562, 0.188620966368049, 
-0.0615972418208484, -0.00184933102124457, -0.0163171325413688, 
0.00370250929658511, 0.30014673206306, 0.135354035472228, 0.00699671782210069, 
0.0174510674253347, -0.0145499677497698, 0.0113155610814752, 
-0.0683884523999768, 0.20157093417998, 0.186320361855075, -0.115609443650563, 
0.069177592825418, -0.0161221161393796, 0.150181081582068, 0.0632121126749741, 
0.0769960292118834, 0.061783685314432, 0.0442014176783082, -0.00990798027657931, 
-0.00186219548019918, 0.0274216740478325, 0.118878480695049, 
0.0592089915185285, -0.00823096478874009, 0.120750948230554, 
0.278594307094423, -0.0111994006625954, 0.0379360193757585, 6.29460162030332e-05, 
0.0602068958909111, 0.173755367986025, 0.135902420389977, 0.124863098282806, …