Utilizo a fda.... Em 20 de julho de 2015 23:40, Mauro Sznelwar <sznelwar@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@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.