`[<-`(`*tmp*` 中出错,下标越界 下标越界

2023-12-31

在下面的代码中,我尝试创建一个矩阵,列出每个城市的 opt.lam。运行循环后,前两个城市始终有效,然后对于此后的任何城市都会出现错误。

这是我得到的错误。 (coefmatrix 工作正常,只是 lambdamatrix 产生此错误)。

错误于[<-(*tmp*, , i, 值 = c(0.577199381062121, 0.577199381062121, : 下标越界

这是我的代码:

lambdamatrix <- matrix(nrow=n,ncol=2)
rownames(lambdamatrix) <- cityIDs
colnames(lambdamatrix) <- c("lambda.min","lambda.1se")
for (i in 1:n) {
  data <- subset(simdata, city==cityIDs[i])
  x <- as.matrix(data.frame(data[,3:24]))
cvfit <- cv.glmnet(x, data$Y, family="poisson", offset=log(data$population))
opt.lam <- c(cvfit$lambda.min, cvfit$lambda.1se)
fit <- glmnet(x, data$Y, family= "poisson", offset=log(data$population))
abline(plot(fit, "lambda", label= TRUE, 
            main = cityIDs[i]), v=log(opt.lam), lty=2, lwd=3, 
                                    col=c("red","dark green"))
coefmatrix[,i] <- coef(fit, s=opt.lam[1])[1:23]
lambdamatrix[,i] <- c(cvfit$lambda.min, cvfit$lambda.1se)[1:n]
}`

In [,i], i是列索引器(而[i,]将是一个行索引器)。

既然你定义了lambdamatrix as matrix(nrow = n, ncol = 2),一旦你过去了i=2您要求的列不存在。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

`[<-`(`*tmp*` 中出错,下标越界 下标越界 的相关文章

  • R: pi[[j]] 中的错误:下标越界——数据帧列表上的 rbind

    我正在尝试重新绑定一个大的数据帧列表 outputDfList 它是通过将一个复杂的函数应用于一个大表而生成的 您可以通过以下方式重新创建outputDfList df1 data frame randomseq chr15q22 1 tr
  • ggplot2错误:美学必须是长度一,或者与数据长度相同问题:颜色、字母

    我收到此错误 错误 美学必须是长度一 或者与数据长度相同问题 颜色 字母 当我将 ggplot 与数据框一起使用时Z如图所示 Z lt data frame Name c A G C T T T AG AG GC GC CT CT AT A
  • 如何与 R 包 sf 进行“完整”联合

    我尝试使用三个多边形之间的并集sf st union 下图中显示了 ArcGIS Overlay Union All 的结果 我希望通过使用 R 中的 sf 包获得与 OUTPUT 中五个不同多边形类似的结果 library sf a1 l
  • for 循环和 into_iter() [重复]

    这个问题在这里已经有答案了 据我了解 for循环隐式调用into iter 在向量上循环遍历它 但是 在下面的代码中 为什么vec在 Question1 1 的第一个 for 循环中没有消耗 但在第二个 for 循环中 vec 消耗了 fn
  • 根据R中的行差异对行进行分组[重复]

    这个问题在这里已经有答案了 我有一组具有不同采样间隔的动物位置 我想要做的是采样间隔符合特定标准 例如低于特定值 的组和序列 让我用一些虚拟数据来说明 start lt Sys time timediff lt c rep 5 3 20 r
  • tm 包本身是否提供了组合文档术语矩阵的内置方法?

    tm 包本身是否提供了组合文档术语矩阵的内置方法 我在同一语料库上生成了 4 个文档术语矩阵 每个矩阵为 1 2 3 4 克 它们都非常大 200k 10k 因此将它们转换为数据帧然后绑定它们是毫无疑问的 我知道我可以编写一个程序来记录每个
  • gridExtra 2.0.0 更改标题大小

    我知道 gridExtra 已更新 因此 我想知道如何更改标题大小 这不再有效 grid arrange a b c d ncol 2 nrow 2 main textGrob Title gp gpar fontsize 15 font
  • 当我创建新变量时出了什么问题?

    我想根据原始变量施加的条件创建一个新变量 比方说 原始变量 var 是由 1 20 中的随机样本组成的向量 并且 当原来的 var 大于10时 新变量 newvar 被设置为缺失 当 var 小于10时 新变量 newvar 被设置为等于
  • 如何更改 ESS 中的智能分配键(“_”到“<-”)绑定

    在 emacs ESS 中 如何正确更改 ess smart S assign 的键绑定 我尝试的是添加 custom set variables ess smart S assign key to my emacs 但这让奇怪的事情发生了
  • for 循环范围内的初始值设定项列表

    我有从单个超类型派生的不同类型的对象 我想知道使用中有没有什么缺点std initializer在 for 循环范围内列出 如下所示 for auto object std initializer list object1 object2
  • ggplot 中的分层轴?

    我想知道是否可以在 GGLPOT2 或其他图形包 我只是更喜欢 ggplot 中制作分层 分段轴 我想要做的是获取下面的数据 制作一个堆积条形图 其中 x 轴上有周期 但在每个周期内 还有每种动物 那么每只动物内的条形颜色将是 颜色 变量
  • 使用 data.table 而不是 data.frame 进行子集化

    我正在处理一个包含 300 万行和 10 列的数据框 并且正在对其进行一些子集化 我下面有一些玩具代码 当我子集化时 需要很长时间 如果我使用 data table 和 data table 上的子集会更快吗 这是一些玩具代码 s lt c
  • 不同元素的ggplot字体大小

    我知道在创建 ggplot 图后我可以使用theme get 返回所有主题元素的详细信息 这对于弄清楚诸如此类的事情非常有帮助strip text x等等 但我有两件事无法弄清楚 1 在下面的ggplot图形中 代表短语 被土拨鼠卡住的木头
  • 在 RGUI 中自动保存 R 源

    是否可以告诉 RGui 自动保存 R 源代码 以便它们在系统崩溃时不会丢失 我非常怀念这个功能 因为我的系统时不时地崩溃 而且我不想仅仅因为需要这个功能而使用 R Studio 或其他一些 IDE 我喜欢 RGui 提供的简单性 None
  • 如何为每个条形图制作具有定义水平边框的堆叠条形图

    我有一些数据想以一种我不知道如何在条形图中显示的方式 希望你能帮我解决这个问题 我的表由 4 列组成 簇 0 6 IgG Status mild high mild low Severe High 患者 1 16 和值 每个簇的标准化值 这
  • 通过另一个函数将数据和列名称传递给 ggplot

    我将直接跳到一个示例并在后进行评论 cont lt data frame value c 1 20 variable c 1 20 1 20 1 5 1 20 2 group rep c 1 2 3 each 20 value variab
  • 使用 kableExtra 增加行/行间距

    有没有办法在 r markdown 或 bookdown 中使用 kableExtra 增加 pdf 输出的行间距 library knitr library kableExtra kable head iris 5 caption Iri
  • 如何在 R 中使用别名运行系统可执行文件?

    假设我正在 R 中运行系统命令来运行executable inputfile lt path myfile txt 我该如何更换 path myfile txt在下面的命令中inputfile如下面命令所示 system executabl
  • 实例着色器矩阵的设置

    我想绘制实例立方体 我可以打电话GL DrawArraysInstanced PrimitiveType Triangles 0 36 2 成功地 我的问题是所有立方体都绘制在相同的位置和相同的旋转 我如何为每个立方体单独更改它 要创建不同
  • R Markdown 文档标题中的希腊字母

    R markdown 文档的标题中是否可以包含希腊字母 我试过这个 title Amylase author author date 8 March 2017 output pdf document keep tex true toc ye

随机推荐

  • 为什么用户可以设置新密码,而无需输入与我的批处理代码正确的旧密码?

    我最近写了一个文件夹锁定程序 我唯一的问题是当我需要更改密码时 例如我去部分newpass 该程序似乎不会将旧密码与新密码进行比较 而只是使用给定的任何输入 这违背了拥有密码的目的 因为每个人都可以更改密码 而不是用它解锁文件夹 有人可以仔
  • 在 Chrome 中重新加载动画 GIF 时出现问题

    我有一个可以重新加载 gif 动画的应用程序 它在 Safari 中始终有效 但在 Chrome 中时断时续 我相信这个问题与提到的问题类似here https groups google com forum topic shiny dis
  • Excel - 基于ID的工作表值查找[关闭]

    很难说出这里问的是什么 这个问题是含糊的 模糊的 不完整的 过于宽泛的或修辞性的 无法以目前的形式得到合理的回答 如需帮助澄清此问题以便重新打开 访问帮助中心 help reopen questions 我有两张工作表 第一个列出了我的客户
  • 如何获取处理器和硬盘的制造序列号和 ID?

    如何使用 Matlab 获取以下硬件属性 主板制造序列号 处理器 ID 处理器制造序列号 硬盘ID 硬盘制造序列号 是否有任何函数或类负责检测其他机器硬件组件属性的属性 我知道可以使用系统或控制台命令来完成 但我不知道如何完成 不过 我更喜
  • 将 1 亿个文件写入 s3

    我的主要目标是根据每条记录的 id 将记录拆分为文件 目前有超过 150 亿条记录 而且肯定还会增加 我需要一个使用 Amazon EMR 的可扩展解决方案 我已经为一个包含大约 9 亿条记录的较小数据集完成了这项工作 输入文件采用 csv
  • Gensim LDA 中的主题明智文档分布

    python 有没有办法映射属于某个主题的文档 例如 主要是 主题 0 的文档列表 我知道有多种方法可以列出每个文档的主题 但我该如何反过来呢 Edit 我正在使用以下 LDA 脚本 doc set for file in files ne
  • 如何在 Visual Studio 2010 中构建这个简单的 C++/SWIG/C# 项目?

    我需要帮助设置一个简单的 C C SWIG 项目 我很难将使用 SWIG 绑定的 C 项目组合在一起 我正在使用 Visual Studio 2010 和最新版本的 SWIG 我的代码很简单 cpp 文件 h pragma once cla
  • 是否有任何方便的代码覆盖率工具可与 JUnit 一起使用?

    是否有方便的代码覆盖率工具可与 JUnit 一起使用 我想大多数代码覆盖率工具都可以很容易地与 JUnit 一起使用 以前我用过EMMA http emma sourceforge net 这很好 如果您是 Eclipse 用户 有一个 E
  • 如何以编程方式清空浏览器缓存?

    我正在寻找一种以编程方式清空浏览器缓存的方法 我这样做是因为应用程序会缓存机密数据 并且我想在您按 注销 时删除这些数据 这可以通过服务器或 JavaScript 发生 当然 仍然不鼓励在外国 公共计算机上使用该软件 因为存在更多危险 例如
  • 如何在 Mongoose 中验证数组及其元素

    我有这个模式 我在其中验证了数组的元素book 但我不知道如何验证数组本身 var DictionarySchema new Schema book 1 type String required true 2 String 3 String
  • 我可以从运行时加载的共享对象访问主机进程的符号吗?还有其他选择吗?

    在我的场景中 我想要一个插件 它是在运行时加载的共享对象 用于访问 主机应用程序 中的符号 以便我可以向我的应用程序添加任何功能 我已经尝试过 但没有找到任何方法来做到这一点 而且我不知道这是否可能 那么 我可以以某种方式做到这一点 或者使
  • 为什么 `add cl, 2` 以 x86 十六进制打印 `80 c1 02`?

    试图缩小我的问题范围 以便我能够克服这个障碍 This https stackoverflow com questions 28664856 how to interpret x86 opcode map noredirect 1 comm
  • 在 C# 中通过 HTTP 侦听器获取发布的文件

    我使用 C 制作了一个简单的 http 服务器 我知道如何获取发布的数据并输出它们 这是我的 C 代码 public static void start HttpListener listener new HttpListener list
  • 用于图像创建的 Scala 库 [关闭]

    Closed 这个问题正在寻求书籍 工具 软件库等的推荐 不满足堆栈溢出指南 help closed questions 目前不接受答案 有没有好的 Scala 或 Java 图像处理库 对于简单的任务 例如与其他一些图像合成图像然后生成拇
  • Spark 2.0中DataFrame的GroupBy操作需要花费大量时间

    在我的一项 Spark 作业 EMR 5 0 0 上的 2 0 中 我有大约 5GB 的数据与 30 行交叉连接 数据大小几 MB 我还需要按它分组 我注意到我花费了大量时间 使用 1 个 m3 xlarge 主节点和 6 个 m3 2xl
  • 是否可以将 requirejs 模块转换为 commonjs ?

    It s 已经可以将 commonjs 模块转换为 requirejs http requirejs org docs commonjs html autoconversion 但我仍然想知道是否可以做相反的事情 有没有办法将 requir
  • 依赖注入和工作单元模式

    我有一个两难的处境 我使用 DI 即 工厂 为自制 ORM 提供核心组件 容器根据请求提供数据库连接 DAO 映射器及其结果域对象 这是映射器和域对象类的基本概述 class Mapper public function construct
  • 我如何在 xcode iphone 中创建一个滑动菜单,就像主 android 菜单滑动菜单一样?

    我一直在网上寻找标题中描述的滑动菜单的示例 我所需要的只是知道我应该从 iphone 库中查找哪些项目来制作这个 我不想占用别人的时间让他们写出代码 但如果有一点指导将不胜感激 我们在 iPhone 应用程序中创建了一个滑动抽屉 我们使用以
  • PHP 中仅用一个空格替换多个空格和换行符

    我有一个带有多个换行符的字符串 字符串 This is a dummy text I need to format this 期望的输出 This is a dummy text I need to format this 我正在使用这个
  • `[<-`(`*tmp*` 中出错,下标越界 下标越界

    在下面的代码中 我尝试创建一个矩阵 列出每个城市的 opt lam 运行循环后 前两个城市始终有效 然后对于此后的任何城市都会出现错误 这是我得到的错误 coefmatrix 工作正常 只是 lambdamatrix 产生此错误 错误于 l