使用函数创建列

时间:2015-03-02 11:44:59

标签: r apply

我在尝试使用基于函数的条件计算创建新列时遇到问题。

我有一些小数据集用于根据海拔高度(CalcAlt)插入参考温度(Tref)。

当我尝试进行单个计算时,该功能有效,但当我尝试将该函数应用于数据集时,我会遇到问题,目的是创建一个新的Tref列。

代码摘录如下。

欢迎任何建议!

我收到错误"申请出错(FUN = FUN_Tref,Airlines $ AC_MODEL,Airlines $ CalcAlt):   昏暗(X)必须有正长度"。

这是我第一次使用apply函数。

我哪里错了?我是R的亲戚,所以可以承认迷路了!

史蒂夫

    CalcAlt <- c(200,200,400,400,600,600,800,800,1000,1000)
    AC_MODEL <-c("320-232","321-231","320-232","321-231","320-232","321-231","320-232","321-231","320-232","321-231" )
    Airlines <- data.frame(AC_MODEL , CalcAlt) 

    #Create dataframe showing Tref 
    V2533_alt <- c(-2000,-1000,0,1000,2000,3000,4000,5000,6000,7000,8000,9000,10000,11000,12000,13000,14000,14500)
    V2533_Tref <- c(19.2, 17.1,15,13.8,11.4,11.9,12,11.9,13.9,15.9,17.8,19.5,21.2,21.8,22.8,25.1,28.5,30.7)

    V2533 <- data.frame(V2533_alt, V2533_Tref) 
    V2533

    V2527_alt <- c(0,5200,14500)
    V2527_Tref <- c(31, 25,25)

    V2527 <- data.frame(V2527_alt, V2527_Tref) 
    V2527

    #Create function to calculate Tref based on aircraft type and calculated altitude
    FUN_Tref <- function(AC_type, CalcAlt){
    if(AC_type =="320-232"){
            #systematic calculation
            Tref <- approx(V2527_alt, V2527_Tref, xout = CalcAlt)

            }
    if(AC_type =="321-231"){
            #systematic calculation
            Tref <- approx(V2533_alt, V2533_Tref, xout = CalcAlt)
            }

    Tref <- as.numeric(Tref[2])
    return(Tref)
    }
    #END-OF_FUNCTION
    ############################
    #Apply function to create new column Tref
    Airlines$Tref <- apply( FUN = FUN_Tref, Airlines$AC_MODEL, Airlines$CalcAlt)

1 个答案:

答案 0 :(得分:1)

使用mapply

Airlines$Tref <- mapply( FUN = FUN_Tref, Airlines$AC_MODEL, Airlines$CalcAlt)
#   AC_MODEL CalcAlt     Tref
#1   320-232     200 30.76923
#2   321-231     200 14.76000
#3   320-232     400 30.53846
#4   321-231     400 14.52000
#5   320-232     600 30.30769
#6   321-231     600 14.28000
#7   320-232     800 30.07692
#8   321-231     800 14.04000
#9   320-232    1000 29.84615
#10  321-231    1000 13.80000

来自帮助:

  

mapply将FUN应用于每个...参数的第一个元素,第二个元素,第三个元素,等等。

PS:我还没有检查过这是否可以进行矢量化。