Displaying 3 results from an estimated 3 matches for "fctk1".
Did you mean:
fct1
2006 Jun 02
1
doubt with integrate ()
...ink this may have something to do with the function (it is not
returning a vector but a number), but I don't see exactly what.
The function to integrate was defined with a for loop first:
for (i in 1:(nbcol-1))
{tdon1<-tdon[,i]
tdon11<-subset(tdon1,tdon1!="NA")
fctk1<-function(x)
{density(tdon11, bandwidth=h[i], kernel="gaussian")$y}
for (j in (i+1):nbcol)
{tdon2<-tdon[,j]
tdon21<-subset(tdon2,tdon2!="NA")
fctk2<-function(x)
{density(tdon21,bandwidth=h[j],kernel="gaussian")$y}
diffctk<-f...
2006 May 17
0
script S-plus -> R
...tdon,na.rm=T,unbiased=T)
h<-1.06*ect*(nlc^-0.2)
cpt<-0
nlgn<-choose(nbcol,2)
Results<-matrix("*",nrow=nlgn,ncol=3,dimnames=list(NULL,c("Ech.1","Ech.2","Overlap")))
for (i in 1:(nbcol-1))
{tdon1<-tdon[,i]
tdon11<-tdon1[tdon1!="NA"]
fctk1<-function(x)
{ksmooth(tdon11,kernel="normal",bandwidth=h[i],x.points=x)$y}
for (j in (i+1):nbcol)
{tdon2<-tdon[,j]
tdon21<-tdon2[tdon2!="NA"]
fctk2<-function(x)
{ksmooth(tdon21,kernel="normal",bandwidth=h[j],x.points=x)$y}
diffctk<-function(x)
{abs(fctk1...
2006 Jul 11
0
script problem to obtain pairs of overlap values
...5, 11) , sp.2=c (4, 2, 4, NA, 11,
),sp.3=c(5, 4, 2, 6, 13), sp.4=c(3 , 11, NA, 5, 3), sp.5=c(2 ,5 ,2, 9, 9))
> h
[1] 1.047 2.973 0.887 1.520 2.955
Here is the code:
for (i in 1:(nbcol-1)) # nbcol<-ncol(tdon)
{tdon1<-tdon[,i]
tdon11<- subset(tdon1,tdon1!="NA")
fctk1<-function(x)
{density (tdon11, bw=h[i], kernel="gaussian")$y}
for (j in (i+1):nbcol)
{tdon2<-tdon[,j]
tdon21<- subset(tdon2,tdon2!="NA")
fctk2<-function(x)
{density (tdon21, bw=h[j], kernel="gaussian")$y}
diffctk<-function(x...