如何修复套索中的“storage.mode(y) <-“double”错误:更改因子的存储模式无效”?

2024-04-05

使用岭回归和套索回归时出现以下错误。

storage.mode(y)

我的数据集如下所示:

> as_data_frame(blca.only)
# A tibble: 315 x 435
   Tags  TCGA.BL.A5ZZ.01 TCGA.XF.AAN1.01 TCGA.DK.A2I4.01 TCGA.XF.A9T4.01
   <fct>           <dbl>           <dbl>           <dbl>           <dbl>
 1 ENSG~           2001.           2170.             0                0 
 2 ENSG~              0               0             64.4              0 
 3 ENSG~          68078.          93147.        108379.           78136.
 4 ENSG~              0               0              0                0 
 5 ENSG~         160562.         176357.        153682.          204135.
 6 ENSG~         217109.         399909.        215158.          253031.
 7 ENSG~              0               0              0                0 
 8 ENSG~           1583.           1236.          1982.            4618.
 9 ENSG~        3535627.        7037139.       2725637.         3036243.
10 ENSG~           1843.           3997.          1282.            1607.
# ... with 305 more rows, and 430 more variables: TCGA.DK.AA6Q.01 <dbl>,
#   TCGA.XF.A8HB.01 <dbl>, TCGA.2F.A9KR.01 <dbl>, TCGA.FD.A3N6.01 <dbl>,
#   TCGA.GV.A3QK.01 <dbl>, TCGA.ZF.A9RE.01 <dbl>, TCGA.FD.A43S.01 <dbl>,
#   TCGA.K4.A5RI.11 <dbl>, TCGA.ZF.AA4U.01 <dbl>, TCGA.GV.A3QF.01 <dbl>,
#   TCGA.UY.A9PE.01 <dbl>, TCGA.FD.A3B4.01 <dbl>, TCGA.GC.A3YS.01 <dbl>,
#   TCGA.G2.A2EF.01 <dbl>, TCGA.BT.A42C.01 <dbl>, TCGA.FD.A62O.01 <dbl>,
#   TCGA.BL.A13J.01 <dbl>, TCGA.E5.A4U1.01 <dbl>, TCGA.E7.A6MD.01 <dbl>,
#   TCGA.CF.A47Y.01 <dbl>, TCGA.5N.A9KI.01 <dbl>, TCGA.DK.A3IT.01 <dbl>,
#   TCGA.DK.AA77.01 <dbl>, TCGA.ZF.AA4R.01 <dbl>, TCGA.GU.A763.01 <dbl>,
#   TCGA.DK.A3IL.01 <dbl>, TCGA.FD.A6TC.01 <dbl>, TCGA.FD.A3SP.01 <dbl>,
#   TCGA.GC.A6I1.01 <dbl>, TCGA.FD.A6TK.01 <dbl>, TCGA.BL.A0C8.01 <dbl>,
#   TCGA.ZF.AA51.01 <dbl>, TCGA.ZF.AA5H.01 <dbl>, TCGA.CF.A9FM.01 <dbl>,
#   TCGA.H4.A2HQ.01 <dbl>, TCGA.DK.AA74.01 <dbl>, TCGA.XF.AAMR.01 <dbl>,
#   TCGA.BT.A20N.01 <dbl>, TCGA.DK.A2HX.01 <dbl>, TCGA.CF.A27C.01 <dbl>,
#   TCGA.DK.AA6M.01 <dbl>, TCGA.XF.A9SH.01 <dbl>, TCGA.CF.A7I0.01 <dbl>,
#   TCGA.KQ.A41P.01 <dbl>, TCGA.CF.A3MG.01 <dbl>, TCGA.CF.A9FF.01 <dbl>,
#   TCGA.FJ.A3ZF.01 <dbl>, TCGA.K4.A6MB.01 <dbl>, TCGA.HQ.A5NE.01 <dbl>,
#   TCGA.XF.A9SV.01 <dbl>, TCGA.FD.A6TB.01 <dbl>, TCGA.H4.A2HO.01 <dbl>,
#   TCGA.ZF.AA4N.01 <dbl>, TCGA.E7.A7XN.01 <dbl>, TCGA.G2.AA3B.01 <dbl>,
#   TCGA.CF.A47T.01 <dbl>, TCGA.K4.A3WU.01 <dbl>, TCGA.GU.A766.01 <dbl>,
#   TCGA.GV.A3JV.01 <dbl>, TCGA.GC.A3WC.11 <dbl>, TCGA.BT.A20R.01 <dbl>,
#   TCGA.DK.A2I6.01 <dbl>, TCGA.ZF.AA53.01 <dbl>, TCGA.ZF.A9R1.01 <dbl>,
#   TCGA.E7.A97P.01 <dbl>, TCGA.C4.A0EZ.01 <dbl>, TCGA.DK.A3X2.01 <dbl>,
#   TCGA.FD.A5BU.01 <dbl>, TCGA.K4.A6FZ.01 <dbl>, TCGA.KQ.A41Q.01 <dbl>,
#   TCGA.CF.A9FH.01 <dbl>, TCGA.BL.A13I.01 <dbl>, TCGA.ZF.A9R4.01 <dbl>,
#   TCGA.DK.A3IU.01 <dbl>, TCGA.CU.A0YN.01 <dbl>, TCGA.2F.A9KT.01 <dbl>,
#   TCGA.XF.A9SL.01 <dbl>, TCGA.DK.AA75.01 <dbl>, TCGA.CF.A9FL.01 <dbl>,
#   TCGA.DK.A6B1.01 <dbl>, TCGA.K4.A3WV.11 <dbl>,
#   TCGA.BL.A13I.01_Rep85 <dbl>, TCGA.GV.A3QH.01 <dbl>,
#   TCGA.XF.A9SM.01 <dbl>, TCGA.ZF.AA5P.01 <dbl>, TCGA.DK.A6B0.01 <dbl>,
#   TCGA.G2.A2ES.01 <dbl>, TCGA.CF.A1HS.01 <dbl>, TCGA.FJ.A3Z7.01 <dbl>,
#   TCGA.HQ.A2OF.01 <dbl>, TCGA.C4.A0F6.01 <dbl>, TCGA.FD.A5C1.01 <dbl>,
#   TCGA.ZF.A9R9.01 <dbl>, TCGA.GU.A42P.01 <dbl>, TCGA.FD.A43U.01 <dbl>,
#   TCGA.S5.A6DX.01 <dbl>, TCGA.GD.A3OP.11 <dbl>, TCGA.K4.A4AC.01 <dbl>,
#   TCGA.GD.A2C5.01 <dbl>, TCGA.4Z.AA84.01 <dbl>, ..

I tried:

> lasso.1 <- glmnet(x = as.matrix(blca.only[,-1]), y = factor(blca.only$Tags), alpha = 1)

Error in storage.mode(y) <- "double" : 
  invalid to change the storage mode of a factor
> Y <- factor(blca.only$Tags)
> lasso.1 <- glmnet(x = as.matrix(blca.only[,-1]), y = blca.only$Tags, alpha = 1)
Error in storage.mode(y) <- "double" : 
  invalid to change the storage mode of a factor

> Y <- factor(blca.only$Tags)
> lasso.1 <- glmnet(x = as.matrix(blca.only[,-1]), y = Y, alpha = 1)

Error in storage.mode(y) <- "double" : 
  invalid to change the storage mode of a factor

我的数据类型如下所示:

    > typeof(blca.only)

[1] "list"

我不知道如何解决这个问题。先感谢您!希望你能帮我解决这个问题。


在你的情况下,我假设你想要一个逻辑回归。因此您必须指定family争论。尝试这个:

lasso.1 <- glmnet(x = as.matrix(blca.only[,-1]), y = factor(blca.only$Tags), alpha = 1,
family = "binomial")

也可以看看这个问题 https://stats.stackexchange.com/questions/72251/an-example-lasso-regression-using-glmnet-for-binary-outcome更多细节。

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

如何修复套索中的“storage.mode(y) <-“double”错误:更改因子的存储模式无效”? 的相关文章

  • R:使用管道将单个参数传递到函数中的多个位置

    我试图专门使用管道来重写以下代码 使用babynames包中的babynames数据 library babynames library dplyr myDF lt babynames gt group by year gt summari
  • 如何在R中逐行写入文件

    我正在尝试逐行读取 csv 文件 并且仅选择左侧的第二个和第三个单元格以及右侧的第三个单元格 例如 如果这一行中有 17 个单元格 我将选取第 15 个单元格 然后我想合并这 3 个单元格 用逗号分隔 然后将此行写入一个新的 csv 文件
  • R:在 Shiny 中,如何修复应用于“反应性”类对象的“xtable”没有适用的方法

    我收到此错误 Error in UseMethod xtable no applicable method for xtable applied to an object of class reactive UI R library shi
  • 当子集长度为零时,如何简洁地处理子集?

    从向量中排除元素x x lt c 1 4 3 2 我们可以减去位置向量 excl lt c 2 3 x excl 1 1 2 这也是动态工作的 excl lt which x which max x gt quantile x 25 1 2
  • 使用 R 中的 tidyverse 重新调整因子和重新排序因子

    我想使用这些功能重新调平 and 重新排序 在我的数据框中 我了解重新调整级别的工作原理 但我不明白为什么我在 data frame 中看不到级别的变化 例如 假设我有鸢尾花数据集 library tidyverse head iris g
  • 如何在 R 中查找平衡面板数据(又名,如何查找面板中的哪些条目在给定窗口内完整)

    我有来自 Compustat 的大量数据 我向其中添加了一些手工收集的数据 认真地从一堆旧书中手工收集 但我不想手工收集整个面板 只想随机选择一个子集 为了找到更大的集合 我从中随机选择 我想从 Compustat 的平衡面板开始 我看到p
  • for 循环与 cor.test 在许多类别上

    我正在尝试在 R 中编写一个循环 它将循环遍历 3 个不同的物种 以计算两个连续变量 Redness 和 VarNormAbund 之间的相关性 我的循环正在运行 但 3 个物种中每一个的输出都是相同的 这让我认为循环卡在第一个物种上 co
  • 如何将 mcmc.list 转换为 bugs 对象?

    我正在使用rjagsR 库 功能coda samples产生一个mcmc list 例如 来自example coda samples library rjags data LINE LINE recompile LINE out lt c
  • 将 csv 文件上传到shinyApps.io

    我的应用程序在本地运行良好 并且我能够成功地将应用程序部署到shinyapps io 服务器 但是当我尝试使用shinyapps URL 在浏览器中加载应用程序时 收到以下错误消息 错误对象 数据 不是成立 我认为这是因为 data 变量从
  • 在 R 中收集多组列[重复]

    这个问题在这里已经有答案了 我有一个宽数据框 需要将其收集或融化成一个高数据框 我遇到的问题是我有几组列需要保持关联 分组 我每个表单提交有 2 个用户 每个用户有 3 列数据 我想将这 6 列基本上以 3 组的形式堆叠起来 以便每个用户都
  • 一次评论多个对象

    假设您有一个变量 a 到 j 的列表 for x in 1 10 assign letters x x 您将如何评论那些最近创建的对象 我尝试过类似的事情 for x in 1 10 comment get letters x lt pas
  • 单击 R Shiny 中的按钮后将输入字段重置为 null

    我正在构建一个应用程序 用户可以在其中按列输入表的数据值 单击 添加 按钮后 输入的值将按列附加到现有值 例如 如果输入 col1 2 3 并单击 ADD 我们将在显示屏中看到 col1 2 3 如果输入 col2 4 7 并单击 ADD
  • R 中的线性模型 - 乘法表达式

    我有 3 个数值变量A B and C 我正在尝试创建一个能够预测的线性模型A 我使用的表达式是B C为了预测A 然而 当查看输出时 我无法得到我的方程 因为我得到了额外的变量 但我不知道它是什么 这是我的代码 MyData lt read
  • 如何更新 R 中的包?

    我想将一个 R 包升级到已经可用的新版本 我试过 update packages c R2jags 但它什么也没做 控制台上没有输出 没有错误 什么也没有 我使用了相同的语法install packages但也许我做错了什么 我一直在看 u
  • R 的 read.csv 在第一列名称前面加上垃圾文本[重复]

    这个问题在这里已经有答案了 我已将数据从 SQL Server Management Studio 中的结果网格导出到 csv 文件 csv 文件看起来正确 但是当我使用 read csv 将数据读入 R 数据帧时 第一个列名称前面带有 如
  • ggplot2 - 在绘图顶部添加辅助 y 轴

    对于出版物 我需要向现有绘图添加第二个 y 轴 我遇到了一种方法来做到这一点 https rpubs com kohske dual axis in ggplot2 https rpubs com kohske dual axis in g
  • 将其他数据集的点添加到ggplot2

    关于这个主题已经有很多问题 但我找不到能回答我的具体问题的问题 我有一个barplot see testplot1 and testplot3如下 绘制数据集 bardata如下 并希望从另一个数据集向其添加点 pointdata 请参阅简
  • geom_polygon 的渐变填充

    此代码生成一个包含 3 个多边形的图表 我正在创建一个显示 3 个多边形的图表 如果有更好的方法来绘制多边形 我不太感兴趣 实际上这些多边形代表事件 并且这些事件有一个持续时间 首先 我感兴趣的是使用渐变填充每个多边形的可能性 librar
  • 安装 gplots 时出错

    我正在 OSX v 10 9 2 上运行 R v 3 0 3 当尝试使用以下命令在 R studio 中安装 gplots 包时 出现错误 gt library gplots Error in library gplots there is
  • 为什么在 data.frame 中预先指定类型会比较慢?

    我预先分配了一个大 data frame 以便稍后填写 我通常这样做NA是这样的 n lt 1e6 a lt data frame c1 1 n c2 NA c3 NA 我想知道如果我预先指定数据类型是否会让事情变得更快 所以我测试了 f1

随机推荐

  • 使用硬件加速内容截取 WKWebview 的屏幕截图

    我在截屏时遇到严重问题WKWebview内容当有硬件加速内容 一些在 iframe 内运行的特定赌场游戏 到目前为止 我使用了像大家建议的标准截图方式 UIGraphicsBeginImageContextWithOptions conta
  • PHP 中可选包含

    我有一个包含常规配置的配置文件 在 git 存储库中 以及一个覆盖配置属性的本地配置文件 在存储库中被忽略 现在本地配置文件包含在配置文件的开头 include once local config php 但我希望包含是有条件的 仅当文件
  • 防止 PR 完成后删除分支

    在 Azure Devops 中 我有一个带有开发分支的 git 存储库 我们从此分支创建多个功能分支 并通过拉取请求将代码合并到开发中 一旦我们完成拉取请求 功能分支就会被删除 我想阻止这种情况发生 我想保留这些功能分支 我怎样才能做到这
  • 使用 Celery 创建动态队列

    这是我的场景 当用户登录我的网站时 我会为给定用户排队一堆任务 通常每个任务需要 100 毫秒 每个用户有 100 多个任务 这些任务排队到默认的 Celery 队列中 并且我有数百个工作线程正在运行 当任务在后端完成时 我使用 webso
  • 将新的 Date() 格式设置为 EEE MMM dd HH:mm:ss zzz yyyy

    我使用 new date 显示日期 时间 目前显示的是 Thu May 31 2012 13 04 29 GMT 0500 CDT 我需要这个 Thu May 31 13 04 29 CDT 2012 我该如何格式化它 您可以使用正则表达式
  • 两个 .Net 应用程序之间的高效通信

    我目前正在用 c 编写一个 Net 应用程序 它有两个主要组件 数据生成器 生成大量数据的组件 Viewer 能够可视化生成器创建的数据的 WPF 应用程序 这两个组件目前是我的解决方案中的两个单独的项目 此外 我正在使用棱镜4 0框架 以
  • 安排 Web Api 方法按设定的时间间隔运行

    在我当前的项目中 需要安排一个方法以设定的时间间隔运行 例如每周一次 目前这是通过 Windows 服务创建 HttpClient 并点击所需的控制器方法来完成的 我想知道这是否可以在 Web Api 项目本身中实现自动化 而不是使用外部服
  • JsonIgnore 在 System.Web.Mvc.Controller 中不起作用

    我有一个 Web API 项目和一个带有一些属性的简单类 其中一些已标记
  • Android SQLite 数据库,为什么删除表并在升级时重新创建

    在我正在关注的教程以及更多地方我看到了这一点 onUpgrade gt 删除表 如果存在 然后重新创建表 这样做的目的是什么 private static class DbHelper extends SQLiteOpenHelper pu
  • 使用 ggplot R 处理多图

    我有一个大数据框 我正在使用 ggplot ggplot geom line data DATA aes logl PercPos group name col blue geom line data DATA aes logl PercN
  • 法新社的“find_theorems”

    我怎样才能使用find theorems搜索整个正式证据档案馆 AFP 的机制 我已将存档下载到我的计算机上 并且可以从中导入理论 例如 如果我写imports Kleene Algebra Kleene Algebra Models那么该
  • 在 SQL 中组合连续行中特定列的值

    我正在将日志导出到 SQL 问题在于 如果消息列的长度超过 800 个字符 则日志记录 API 会对其进行修剪 并使用消息的其余部分创建新条目 我使用 logparser 将这些数据直接导出到 SQL 中 因此 对于其余的讨论 我们假设我必
  • MySQL 使用 Grandtotal 进行行小计

    我需要创建一个 MariaDB SQL 它允许我对两列 借方和贷方列 求和以获得差额 并返回每种不同帐户类型的小计 表格如下 Account Debit Credit acc1 1 2 acc1 1 4 acc2 3 2 acc2 2 1
  • 组合列并删除 NaNs Pandas

    如果我有一个像这样的 pandas 数据框 a1 0 116667 NaN NaN a2 NaN 0 516667 NaN a3 NaN 0 006667 NaN a4 NaN NaN 0 426667 a5 NaN 0 506667 Na
  • 导航到页面会增加内存使用量 Windows Universal 8.1

    我正在创建一个 Windows Universal 8 1 应用程序 每次我导航到一个页面 然后导航回来 然后再次导航到该页面 该页面的新实例就会保存在内存中 显然 垃圾收集器会在一段时间后释放内存 但是如果不需要 我宁愿不使用内存 有没有
  • Shutil make_archive 生成​​嵌套的 .zip 文件 [重复]

    这个问题在这里已经有答案了 我一直在尝试使用make archive https docs python org 3 library shutil html shutil make archive from shutil图书馆 这是代码 带
  • 使用相同的值更新 mysql 表并仍然获得时间戳更新

    所以我有这个stamp timestamp DEFAULT NOW ON UPDATE NOW 行在我的表上 即使我正在执行的更新基本上是所有字段上的相同数据 我也需要它进行更新 在表的声明中是否有任何方法可以做到这一点 就像除此之外的其他
  • 如何向 woocommerce 编辑订单页面添加选项?

    如何在此处添加另一个选项 例如订单操作 订单总计 我知道如何创建选项页面 但不知道如何编辑 woocommerce 选项 有什么办法吗 在 woocommerce 中自定义订单和优惠券非常简单 因为它们只是另一种自定义帖子类型 您可以在其中
  • 是否可以使用 JMS 进行 IBM MQ 消息分段?

    是否可以使用 JMS 实现消息分段 就像使用本机 IBM API 一样此处显示 http www 01 ibm com support docview wss uid swg21405730 我读过的一种可能的解决方案是 JMS 的消息分组
  • 如何修复套索中的“storage.mode(y) <-“double”错误:更改因子的存储模式无效”?

    使用岭回归和套索回归时出现以下错误 storage mode y 我的数据集如下所示 gt as data frame blca only A tibble 315 x 435 Tags TCGA BL A5ZZ 01 TCGA XF AA