[R] generating phi using function()

[T.Riedle]

Hi,
I am struggling with following function
> phi <- function(w1, w2, j, k, K){
+   zaehler <- (k/K)^(w1-1)*(1-k/K)^(w2-1)
+   nenner <- sum( ((1:K)/K)^(w1-1)*(1-(1:K)/K)^(w2-1))
+   return( zaehler/nenner )
+ }
> phi(c(1, 1), 44L, 1)
Error in phi(c(1, 1), 44L, 1) : argument "k" is missing, with no default


Hence, I have changed the function to

phi <- function(w, k, K){
+ w1 <- w[1]
+ w2 <- w[2]
+ zaehler <- (k/K)^(w1-1)*(1-k/K)^(w2-1)
+ nenner <- sum( ((1:K)/K)^(w1-1)*(1-(1:K)/K)^(w2-1))
+ return( zaehler/nenner )
+ }

 Unfortunately, when running the midas regression I get the following error message

 m22.phi<- midas_r(rv~mls(rvh,1:max.lag+h1,1,phi), start = list(rvh=c(1,1)))

Error in X[, inds] %*% fun(st) : non-conformable arguments

I guess the problem is w but I do not find a solution how to produce the formula shown in the attached file where the exponents are w1 and w2, respectively.

Thanks for your help


From: JLucke at ria.buffalo.edu [mailto:JLucke at ria.buffalo.edu]
Sent: 30 March 2015 16:01
To: T.Riedle
Cc: r-help at r-project.org; R-help
Subject: Re: [R] generating phi using function()

Your function phi has 5 arguments with no defaults.  Your call only has 3 arguments.  Hence the error message.
> phi <- function(w1, w2, j, k, K){
+   zaehler <- (k/K)^(w1-1)*(1-k/K)^(w2-1)
+   nenner <- sum( ((1:K)/K)^(w1-1)*(1-(1:K)/K)^(w2-1))
+   return( zaehler/nenner )
+ }
> phi(c(1, 1), 44L, 1)
Error in phi(c(1, 1), 44L, 1) : argument "k" is missing, with no default



>





"T.Riedle" <tr206 at kent.ac.uk<mailto:tr206 at kent.ac.uk>>
Sent by: "R-help" <r-help-bounces at r-project.org<mailto:r-help-bounces at r-project.org>>

03/29/2015 08:59 AM

To

"r-help at r-project.org<mailto:r-help at r-project.org>" <r-help at r-project.org<mailto:r-help at r-project.org>>,

cc

Subject

[R] generating phi using function()







Hi everybody,
I am trying to generate the formula shown in the attachment. My formula so far looks as follows:

phi <- function(w1, w2, j, k, K){
zaehler <- (k/K)^(w1-1)*(1-k/K)^(w2-1)
nenner <- sum( ((1:K)/K)^(w1-1)*(1-(1:K)/K)^(w2-1))
return( zaehler/nenner )
}

Unfortunately something must be wrong here as I get the following message when running a midas regression

m22.phi<- midas_r(rv~mls(rvh,1:max.lag+h1,1,phi), start = list(rvh=c(1,1)))
Error in phi(c(1, 1), 44L, 1) : argument "K" is missing, with no default
Called from: .rs.breakOnError(TRUE)
Browse[1]> K<-125
Browse[1]> 125

Could anybody look into my phi formula and tell me what is wrong with it?

Thanks in advance.

______________________________________________
R-help at r-project.org<mailto:R-help at r-project.org> mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html>
and provide commented, minimal, self-contained, reproducible code.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: psi.png
Type: image/png
Size: 2199 bytes
Desc: psi.png
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20150330/69a1a9c5/attachment.png>