r - 将列表中的字符串值应用于默认的 drc 函数
问题描述
我正在运行几个模型和模型的功能。我需要将列表中的值应用到默认的 R 函数,例如 modelFit() 我的想法是对多个模型运行相同的函数,然后编译结果。
我正在尝试在 R 中使用循环和“应用”函数系列,但没有成功。
#package drc is necessary
library(drc)
#my data
rates <- c(.1,.1,.1,1,1,1,10,10,10,100,100,100,1000,1000,1000,.1,.1,.1,1,1,1,10,10,10,100,100,100,1000,1000,1000)
prod <- c("A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","B","B","B","B","B","B","B","B","B","B","B","B","B","B","B")
resp <- c(.295,.32,.301,.155,.1501,.148,.05,.03,.044,.002,.001,.0015,.001,.0005,.0003,.312,.337,.299,.265,.2501,.248,.102,.103,.114,.02,.01,.015,.003,.0002,.0007)
data.test <- data.frame(rates,prod,resp) #my data frame
#my models
m1 <- drm(resp~rates, fct=LL.4(), data=data.test[data.test$prod=="A",])
m2 <- drm(resp~rates, fct=LL.4(), data=data.test[data.test$prod=="B",])
#lack of fit test
modelFit(m1)
modelFit(m2)
#I can get the modelFit p-values this way:
modelFit(m1)$"p value"[2]
modelFit(m2)$"p value"[2]
#I have several models. I want to create a loop to give me only the p value for each model fit. I want to use that for other information given by the summaries and function
#list of models
modelsList <- c("m1","m2")
#I can print the strings with the loop
for(i in modelsList){
print(i)
}
我的想法是使用字符串为默认的“drc”函数 modelFit() 添加信息。结果将是一个列表,其中包含循环中所有结果中的所有 p 值,但我收到错误消息:
错误:$ 运算符对原子向量无效
#Not working
for(i in modelsList){
modelFit(i)$"p value"[2]
}
#Error
# Error: $ operator is invalid for atomic vectors
#Trying to use lapply
#Following this logic
lapply(modelsList, function(x) print(x))
#I could not get the results
lapply(modelsList, function(x) modelFit(x)$"p value"[2])
有了结果,我将继续创建一个包含所有模型和相应 p 值的 data.frame。
解决方案
很高兴见到你,达尼洛
我认为你错过了drm
功能。
modelsList <- c("m1","m2")
prod = c("A","B")
type = data.frame(modelsList,
prod)
#I can print the strings with the loop
for(i in modelsList){
print(i)
}
for(i in modelsList){
model_Prod = type[modelsList == i,]$prod
drm_result<-drm(resp~rates, fct=LL.4(), data=data.test[data.test$prod==model_Prod,])
print(modelFit(drm_result)$"p value"[2])
}
# lapply
lapply(modelsList, function(x) {
model_Prod = type[modelsList == x,]$prod
x<-drm(resp~rates, fct=LL.4(), data=data.test[data.test$prod==model_Prod,])
modelFit(x)$"p value"[2]
}
)
推荐阅读
- oop - 如何在代码块 IDE 中组织低级设计的类
- regex - Laravel 从正则表达式输入产生错误
- google-cloud-platform - 如何在 Google Cloud CDN 上设置默认首页 index.html
- node.js - Puppeteer 连续获取单元格
- javascript - 我试图在我的 Wix 网站上的几个页面上隐藏标题。似乎没有办法做到这一点。想法?
- r - 具有多个范围的 if 语句 (R)
- python - 如何在这段代码中使用 pylint(not-an-iterable)?
- regex - 正则表达式:二进制字符串至少包含某个整数的 3 个
- ruby - 如何在 jekyll 网站中添加谷歌翻译插件?
- django - 为什么 django 迁移会暂时占用磁盘空间?