[R-br] RES: smooth.bibasis

David Matta davidhmatta em gmail.com
Terça Julho 21 11:32:49 BRT 2015


Utilizo a fda....

Em 20 de julho de 2015 23:40, Mauro Sznelwar <sznelwar em uol.com.br> escreveu:

>  Que biblioteca precisa para isto?
>
>
>
>
>
> Caros,
>
>
>
> Estou utilizando da função smooth.bibasis para suavizar uma superficie a
> partir do produto tensorial de bases b-splines (codigo abaixo).
>
>
>
> Ocorre que quando altero o número de pontos observados (ns e nt) e
> mantenho o número de bases o algoritmo trava.  Alguém já fez algo parecido?
>
>
>
> Grato
>
>
>
> David
>
>
>
>
>
>
>
> snbasis = 15
>
> sbasis  = create.bspline.basis(c(0,1),snbasis)
>
> tnbasis = 15
>
> tbasis  = create.bspline.basis(c(0,1),tnbasis)
>
> lambdas = 1e-3
>
> fdPars  = fdPar(sbasis, 2, lambdas)
>
> lambdat = 1e-3
>
> fdPart  = fdPar(tbasis, 2, lambdat)
>
>
>
> ns = 15
>
> nt = 15
>
>
>
> sarg = seq(0,1,len=ns)
>
> targ = seq(0,1,len=nt)
>
>
>
> sigma = 0.2
>
> y0 = outer(sin(2*pi*sarg),cos(2*pi*targ))
>
> y  = y0 + rnorm(matrix(0, ns,nt))*sigma
>
>
>
> result = smooth.bibasis(sarg, targ, y, fdPars, fdPart)
>
>
>
>
>
>
> ------------------------------
>   [image: Avast logo] <https://www.avast.com/antivirus>
>
> Este email foi escaneado pelo Avast antivírus.
> www.avast.com <https://www.avast.com/antivirus>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150721/a7a7e7d6/attachment.html>


Mais detalhes sobre a lista de discussão R-br