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: jlu...@ria.buffalo.edu [mailto:jlu...@ria.buffalo.edu] Sent: 30 March 2015 16:01 To: T.Riedle Cc: r-help@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" <tr...@kent.ac.uk<mailto:tr...@kent.ac.uk>> Sent by: "R-help" <r-help-boun...@r-project.org<mailto:r-help-boun...@r-project.org>> 03/29/2015 08:59 AM To "r-help@r-project.org<mailto:r-help@r-project.org>" <r-help@r-project.org<mailto:r-help@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@r-project.org<mailto:R-help@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.
______________________________________________ R-help@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 and provide commented, minimal, self-contained, reproducible code.