Функция вызова, которая суммирует фрейм данных на основе параметра и применяется к фрейму данных и nls

Я работаю с данными о погоде и счетами за коммунальные услуги и пытаюсь оценить модель нелинейной регрессии.

Я прихожу с проблемой. Функция, которую я вызываю для расчета статистики погоды, градусо-дней отопления и охлаждения (HDD и CDD), не может быть применена к фрейму данных, и nls не может ее использовать. Ясно, что я упускаю что-то довольно очевидное в аргументах функций.

Может ли кто-нибудь указать, где я ошибаюсь с функциями HDD и CDD ниже?

Вот простой пример проблемы с некоторым кодом, который генерирует поддельные данные о погоде и счетах.

# Generate Fake Weather Data
CZ<-c(1,2)
Date<-c('2001-01-01','2001-01-02','2001-01-03','2001-01-04')
Weather<-expand.grid(CZ,Date)
names(Weather)<-c("CZ","Date")
Weather$AvgTemp<-rnorm(8,mean= 60,sd=20)

#Generate Fake Billing Data
ID<-as.numeric(1:10)
CZ<-c(1,2)
StartDate<-'2001-01-01'
EndDate<-'2001-02-01'
FakeBilling<-data.frame(cbind(ID,CZ,StartDate,EndDate))
FakeBilling$KWH<-rnorm(10,mean=1000, sd=200)

#Heating and cooling degree functions
HDD<- function(b,CZ,StartDate,EndDate) {
    Temps<-Weather$AvgTemp[Weather$CZ==CZ&as.Date(Weather$Date) >=as.Date(StartDate) &     as.Date(Weather$Date) < as.Date(EndDate)];

    sum((b-Temps)/(1+exp(-5*(b-Temps))))
}


CDD <- function(b,CZ,StartDate,EndDate) {
Temps<- Weather$AvgTemp[as.character(Weather$CZ)==as.character(CZ) &     as.Date(Weather$Date) >=as.Date(StartDate)& as.Date(Weather$Date) < as.Date(EndDate)]

    sum((Temps-b)/(1+exp(-5*(Temps-b))))
}

#these work
HDD(60,1,'2001-01-01','2001-02-01')
# [1] 29.34333
CDD(60,1,'2001-01-01','2001-02-01')
# [1] 53.49393

# This does not. Lots of warnings about length
HDD(60,FakeBilling$CZ,FakeBilling$StartDate,FakeBilling$EndDate)
# [1] NA
# Warning messages:
#   1: In is.na(e1) | is.na(e2) :
#   longer object length is not a multiple of shorter object length
# 2: In `==.default`(Weather$CZ, CZ) :
#   longer object length is not a multiple of shorter object length
# 3: In `>=.default`(as.Date(Weather$Date), as.Date(StartDate)) :
#   longer object length is not a multiple of shorter object length
# 4: In `<.default`(as.Date(Weather$Date), as.Date(EndDate)) :
#   longer object length is not a multiple of shorter object length

# Would like to run this but get similar error.
nls(KWH~load + heatload*(HDD(base,CZ,StartDate,EndDate)) ,start=c(load=200,     heatload=.1,base=65), data=FakeBilling, na.action=na.omit)
# Error in numericDeriv(form[[3L]], names(ind), env) : 
#   Missing value or an infinity produced when evaluating the model
# In addition: Warning messages:
#   1: In is.na(e1) | is.na(e2) :
#   longer object length is not a multiple of shorter object length
# 2: In `==.default`(Weather$CZ, CZ) :
#   longer object length is not a multiple of shorter object length
# 3: In `>=.default`(as.Date(Weather$Date), as.Date(StartDate)) :
#   longer object length is not a multiple of shorter object length
# 4: In `<.default`(as.Date(Weather$Date), as.Date(EndDate)) :
#   longer object length is not a multiple of shorter object length
# 5: In is.na(e1) | is.na(e2) :
#   longer object length is not a multiple of shorter object length
# 6: In `==.default`(Weather$CZ, CZ) :
#   longer object length is not a multiple of shorter object length
# 7: In `>=.default`(as.Date(Weather$Date), as.Date(StartDate)) :
#   longer object length is not a multiple of shorter object length
# 8: In `<.default`(as.Date(Weather$Date), as.Date(EndDate)) :
#   longer object length is not a multiple of shorter object length