[R] Help with implementing Whittaker-Henderson graduation for raw-data
percival_bueser at cocolife.com
Fri Feb 3 01:42:19 CET 2017
I eventually was able to make the code anyway. Here it is:
datafile <- read.table("test_file.txt", header = TRUE)
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)
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)
which will reveal the results using the whittaker function itself.
However, when I use this command
I get the following error:
Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck =
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
before library(pracma), I get the following error:
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.
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
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.
> On 2 Feb 2017, at 03:44, Percival Bueser <percival_bueser at cocolife.com>
> 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:
> On the attached .txt file, The x's are the independent variables and the
> 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
> 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.
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> PLEASE do read the posting guide
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help