问题标签 [r-grid]

我有英格兰肯特郡 17,306 个池塘的北向和东向坐标，以及几乎所有池塘的面积（以平方米为单位）。我正在尝试创建一个 1 公里的网格，它给出每个网格正方形的平均池塘面积，同时为没有池塘的网格正方形提供 0 值。我一直在寻找与我类似的问题，并找到了一个使用薄板样条算法生成英国上空的网格化降雨数据，然后将表面绘制到该网格上并将数据写入表格的问题（如何生成网格化输出在R中并消除不在陆地上的网格方块？）。

我已经能够在少量数据上使用此代码来产生类似的结果。以下是我的少量数据的示例。

使用所有 17,306 个池塘它会使我的计算机崩溃，但更重要的是我希望能够以某种方式对其进行调整，而不是通过 tps 产生的预测值，我会得到我想要的，即每个网格正方形的平均池塘面积。我一直觉得这非常困难，所以如果有人能提供解决方案或为我指明正确的方向，我将不胜感激。

亲切的问候，

艾丹

我想最大化以下相关性：

其中 qhat 是长度为 90 的时间序列；mu_1,mu_2 ,mu_3 , andmu_4 are essentially the parameters I would like to run the grid search over;I is the 4 X 4 identity matrix;A is a 4 X 4 matrix of coefficients; hereinv represents the matrix inverse operator;D` 是一个维度为 4 X 90 的矩阵（4 个变量的时间序列）。我搜索并发现了一个类似的问题，并试图破解解决方案以满足我的需求。我在 R 中尝试了以下内容（为冗余代码道歉，我想分解所有内容，所以我尝试的内容很清楚）：

该对象的param_space最终尺寸为 14641 X 4。我使用as.vector上面的函数来避免 R 中出现的关于使用非数字对象执行矩阵运算的错误。t(varset)上面广义等式中的 D 在实际代码中已替换为此处。物体的尺寸t(varset)是 4 X 90。

这是一条可能/可能没有用的信息：当我简单地复制和粘贴时

在我的控制台中，我得到与运行上面编写的代码时相同的错误。即，

另外：警告信息：

所以我想我不恰当地设置了网格搜索。任何帮助，将不胜感激。这实际上是我的第一篇文章，所以我也很欣赏这个资源。不要挑剔，但我试图避免在 R 中循环。

-克里斯

我想为美国和加拿大的网格单元着色。我的目标与这个问题非常相似：R Plot Filled Longitude-Latitude Grid Cells on Map 然而，这个问题只涉及美国，我不知道如何添加加拿大。

通过修改此处的代码，我能够绘制美国和加拿大的地图：https ://groups.google.com/forum/#!topic/ggplot2/KAKhoE0GO4U

这是为网格单元创建假属性数据的代码，可在此处找到：http ://www.numbertheory.nl/2011/11/08/drawing-polar-centered-spatial-maps-using-ggplot2/

这是ggplot2早期 Stack Overflow 帖子中使用的代码，用于在 Lower 48 的地图上覆盖属性网格：

如何结合这两个ggplot语句将虚假属性数据网格覆盖到美国和加拿大的地图上？谢谢你的任何建议。

我正在尝试将 ASCII TOMS 网格格式的文件读入 R。我已经能够以在 R 中打开的方式读取它。但是，我以线性矩阵的形式打开。该文件包含的内容摘要可在此处获得：

可以在此处下载数据集的示例：

该数据集是 2006 年 1 月的数据集，我只是将其重命名以方便访问，因为我需要使用很多数据集。我在使用中阅读它：

如何将其读入 R，以便每个纬度的数据都在一行上？我觉得这将有助于建立一个适当的数据结构。注意：让我尝试一下，简单地说，我理解：
这意味着结构是这样的，一行表示标题，例如 lat=-89.9，接下来的 144 行有 20 个元素，每个元素属于行 lat=-89.9；所以我现在的问题是在下一个“lat = ...”之前将所有这些元素读入一行。

此外，我只是尝试使用以下方法遍历一组文件：

我想boxplot在直方图下方放置一个。我已经想出了如何做到这一点，但是boxplot和直方图大小相同，我想缩小boxplot. 当我减小宽度时，边缘的空间保持不变。但是，我想减小整个事物的宽度。

这是我到目前为止所拥有的：

编辑：
如果我像这样使用 grid.arrange() ：

比例与我想要的完全一样，但我无法弄清楚如何调整上面的第一个示例，以使轴再次对齐。

任何人？

如何更改使用 grid.curve 函数绘制的曲线的颜色和粗细？我尝试向函数参数添加参数，但不支持它们。我可以以某种方式更改默认的 par() 吗？

我想为我的数据绘制一个回归平面：

结构（列表（L = C（96.4155，76.803，71.5615，68.193，65.6975，74.627，67.82，67.82，64.26，62.06，62.06，60.35，68.35，68.284 58.27, 56.42, 54.91, 53.65), a = c(-0.8375, -20.47, -22.1875, -22.6125, -22.5845, -2.2415, -12.91, -16.16, -17.9, -18.77, -2.4, -10.28, - 13.64, -15.33, -16.37, -2.4495, -8.66, -11.74, -13.46, -14.55, -2.4645, -7.58, -10.38, -12.06, -13.09), b = c(1.437, -24.5915, -29.3275) , -31.892, -33.524, -5.151, -15.22, -19.61, -22.65, -24.8, -6.239, -13.24, -17.11, -19.64, -21.55, -6.7845, -12.18, -15.45, -17.68, - 19.47, -7.124, -11.53, -14.38, -16.41, -17.89)), .Names = c("L", "a", "b"), row.names = c(NA, 25L), class = “数据帧”）

我从这里尝试了一个例子。这是我的尝试：

3d-Plot 看起来不错，但对于网格，我收到此错误：

我的代码中的问题在哪里，是否有另一种方法可以解决我的问题？谢谢你的帮助！

更新：

我已经用 plane3d() 试过了，但得到了这个错误：

但是，在寻找其他解决方案后，我发现我的意图是错误的。我正在寻找这样的东西......

数据
意图

...但是不知道如何实现它，或者根据我对 R 的实际了解是否有可能。有人知道吗？

更新2：

我已经使用以下代码尝试了 scatter3d() 函数：

但收到此错误：

我已经发现，我必须使用 choose.k 来调整自由度，但它总是针对 gam 函数进行解释，我不知道如何在 scatter3d 函数中使用它。任何想法？

我写了一个保存ggplots的函数：

我曾经grid.arrange从另外两个创建一个情节：

我可以这样保存它：

但我想用我的功能来做到这一点。我调用我的函数：

不幸的是，它不起作用。知道如何重新排列函数以使其也适用于grid绘图吗？

使用我称之为的低级graphics封装工具绘制一个高度比宽度大 2 倍的“交叉”矩形：

该图有一个很好的功能：保留了纵横比，因此如果我更改绘图窗口的大小，我会得到相同的高宽比例。

如何获得与grid图形等效的结果？

我是 R 新手，但我需要制作一个带有网格的图，网格必须在 1000 步中具有 0.004 到 0.08 的间隔。我已经有了间隔的序列，但我不知道如何让网格和轴正确？随着plot(data)我得到很好的数据作为点和线我得到一个很好的线。我试过grid(NA,1000,lwd=2)了，但这只会让我得到一个情节中的网格线。axis(1, xaxp = seq(0.004, 0.08,by = (0.08-0.004)/1000), las=2)当我尝试这个时，我得到错误图形参数xaxp的长度错误，所以我很茫然。