如何解决函数中的嵌套向量应用于R中的数据帧
我有以下基本数据(问题末尾提供了dput):
P <- 72
E_c <- 80000
head (deflections)
L I_comp
1 60 17299.21
2 70 25760.94
3 80 32734.69
4 90 51343.59
5 100 60167.30
6 110 64887.87
其中b
被定义为b = L - a
,而a
是可以定义为向量的矢量
a <- seq(10,L,10)
和x
是a
这样的嵌套序列,
x <- seq(1,a,1)
我已经编写了下面的函数,以计算所有x
值的所有a
值的最大偏转量。
delta_less <- function(x){
a_seq <- seq(10,10)
L <- x[1]
I_comp <- x[2]
X_seq <- seq(a_seq)
delta <- ((P*(L-a_seq)*X_seq)/(6*E_c*I_comp*L))*(L^2-(L-a_seq)^2-X_seq^2)
c(val = max(delta),pos = which.max(delta))
}
test <- cbind(deflections,vars = t(apply(deflections,1,delta_less)))
问题是我的代码无法识别嵌套序列。仅当我给a_seq
一个奇异的数值时,它才起作用。
提供的数据:
dput(deflections)
structure(list(L = c(60,70,80,90,100,110,120,130,140,60,140),I_comp = c(17299.2063000114,25760.9420686703,32734.6949444858,51343.5889683982,60167.3001695375,64887.8729936444,83451.9473750139,103000.852590734,112631.744898846,23539.3443398425,23544.8558575609,32808.739059987,40880.4964809276,60171.1172692296,64894.6680546358,83469.4777437875,103114.132264558,120223.534252571,25772.7904955165,32810.8231970971,51345.5936366056,74673.7079705815,80752.3694583712,147386.853127916,25767.8092758621,40881.3376639665,55608.9342154021,80694.6568665257,119068.355471205,119402.817753462,147386.853127916)),class = "data.frame",row.names = c(NA,-36L))
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。