[R] Help with implementing Whittaker-Henderson graduation for raw-data

[Percival Bueser]

Good day!

I eventually was able to make the code anyway. Here it is:

getwd()
setwd("C:/R")
datafile <- read.table("test_file.txt", header = TRUE)
library(SparseM)
var2 = datafile$y
d = 2
lambda = 1600
    require(SparseM, warn.conflicts = FALSE)
    m <- length(var2)
    E <- as(m, "matrix.diag.csr")
    class(E) <- "matrix.csr"
    Dmat <- diff(E, differences = d)
    B <- E + (lambda * t(Dmat) %*% Dmat)
    z <- solve(B, var2)
    z
write.table(z, "C:/R/test.txt", sep="\t")  

(then transfer the results to Excel and that's it)

Now, my problem is: 

I'd like to use the following command:
    z2 <- whittaker(var2, lambda, d)
    z2
which will reveal the results using the whittaker function itself. 

However, when I use this command
    library(pracma)
I get the following error:

> library(pracma)
Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck =
vI[[j]]) : 
  there is no package called 'quadprog'
Error: package or namespace load failed for 'pracma'

So this must mean that I should have quadprog, but when I tried
    library(quadprog)
before library(pracma), I get the following error:

> library(quadprog)
Error in library(quadprog) : there is no package called 'quadprog'

I'd like to get around these so I can use the whittaker function itself. How
can I do this? Thank you very much.

Regards,

Percy

-----Original Message-----
From: Berend Hasselman [mailto:bhh at xs4all.nl] 
Sent: Friday, February 3, 2017 2:08 AM
To: Percival Bueser
Cc: r-help at r-project.org
Subject: Re: [R] Help with implementing Whittaker-Henderson graduation for
raw-data


Why don't you just try the function whittaker from the pracma R package?
There is an example which should be adequate for finding out what to do.

Berend Hasselman

> On 2 Feb 2017, at 03:44, Percival Bueser <percival_bueser at cocolife.com>
wrote:
> 
> Good day everyone!
> 
> I would appreciate if anyone can help me regarding the following: I would
> like to implement the Whittaker-Henderson smoothing to the raw data on the
> attached .txt file, based on the description on this link:
> 
> https://artax.karlin.mff.cuni.cz/r-help/library/pracma/html/whittaker.html
> 
> On the attached .txt file, The x's are the independent variables and the
y's
> are the dependent variables. The signal to be smoothed is y, lambda = 1600
> and d = 2.
> 
> Can anyone please send me a sample R script, or a link to an R script
which
> I can adapt, where I can get both the (1) smoothed graph and the (2)
> smoothed values of y for each x given, and then import both the smoothed
> graph and the smoothed values of y to Microsoft Excel?
> 
> Thank you very much.
> 
> Regards,
> 
> Percy
> 
> 
> <Sample_data.txt>______________________________________________
> 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
> and provide commented, minimal, self-contained, reproducible code.