Displaying 2 results from an estimated 2 matches for "spval".
Did you mean:
pval
2006 Oct 16
2
help with plot()
Hi,
I am new to R and I have been trying it out. I ran into a problem with the
plot() function. Below is my code:
> d <- read.table("c:/test/sp.txt",header=0)
> spval <- d[,2]
> plot(spval,type="l")
Warning messages:
1: graphical parameter "type" is obsolete in: plot.window(xlim, ylim, log,
asp, ...)
2: graphical parameter "type" is obsolete in: axis(side, at, labels, tick,
line, pos, outer, font, lty, lwd,
3: graphical parame...
2009 Feb 01
0
possible memory leak involving looping, optimization, and gam
...r on a Windows XP machine with R 2.6.1 (32-bit), nor on Mac OSX (Leopard) with R2.6.2.
mgcv versions are 1.4-1 for Linux and 1.3-29 for Windows and Mac.
n=700
x=runif(n)
dist=abs(matrix(x,n,n,byrow=TRUE)-matrix(x,n,n,byrow=FALSE))
Sigma=exp(-dist/0.3)
y=t(chol(Sigma))%*%rnorm(n)
psFun=function(spVal,dfWanted){
return((summary(gam(y~s(x,k=400),sp=spVal))$edf-dfWanted)^2)
}
for(i in 1:10000){
spVal=optimize(psFun,c(.00001,500),50)$minimum
print(i)
}
Note that the same issue seems to arise regardless of whether I use uniroot, optimize, nlminb, or nlm to do the optimization.
-chris
-----...