Re: [R] Need hep for converting date data in POSIXct

2013-07-11 Thread laila
d/%m/%Y %H:%M") location<- subset(Geo,select=c(lat.comp,long)) time1<- time[,1] lat<- location[,1] long<- location[,2] library(argosfilter) vmask(lat,long,time1,25) #[1] "end_location" "end_location" "not" "not"

Re: [R] Need hep for converting date data in POSIXct

2013-07-11 Thread laila
nvert 'time1' to class "POSIXct" I send you the example again: > library(argosfilter) > setwd("C:/Users/Usuario/Dropbox/Laila Aranda/PUFGRA") > Geo = read.table("2370001_PUFGRA_2009_Gough_000_retarded10_both.trj",header=FALSE,sep = "

Re: [R] Need hep for converting date data in POSIXct

2013-07-09 Thread laila
Muchas gracias Thank very much I have converted my data!😊😊 Laila Date: Tue, 9 Jul 2013 01:56:18 -0700 From: ml-node+s789695n4671137...@n4.nabble.com To: laila_...@hotmail.com Subject: Re: Need hep for converting date data in POSIXct Hi Laila, There is only one

Re: [R] Need hep for converting date data in POSIXct

2013-07-08 Thread laila
_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real,NA_real_), class = c("POSIXct", "POSIXt"), tzone = "") Script I have done is: library(argosfilter) setwd("C:/Users/Usuario/Dropbox/Laila Aranda/PUFGRA") Geo = read.table("2370001_PUFGRA_2009_Gough_

[R] need help in using Hessian matrix

2009-11-02 Thread Laila Alkhalfan
Hi I need to find the Hessian matrix for a complicated function from a certain kind of data but i keep getting this error Error in f1 - f2 : non-numeric argument to binary operator the data is given by U<-runif(n) Us<-sort(U) tau1<- 2 F1tau<- pgamma((tau1/t

[R] intigrate function and absolute error

2009-11-01 Thread Laila Alkhalfan
14 withou the other part(with absolute error <7.2e-05) how can we do that? thank you and take care Laila [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the p

[R] Bootstrapping

2008-07-28 Thread laila khalfan
use the boot function?what should be the ran.gen and the mle?Thank you in advance Laila [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

[R] parametric bootstrap

2008-07-28 Thread laila khalfan
;- pgamma((tau/theta1),shape,1) n1<-sum(Us<- boot(data=c(Ti,Tj) , statistic= func1 , R=1000,sim='parametric', ran.gen=func2, mle= m0(Ti,Tj) ) boot.ci(bootobject, conf=.90, type= 'bca') Thank you and take care. Laila __ R-help@r

[R] parametric bootstrap

2008-07-26 Thread laila khalfan
HiI am using the boot function in boot package. I am facing a problem in getting the values of bias and st.error in the output. r<-36n<-40shape<-2theta11<-exp(1) # (=2.718282)theta21<-exp(.5) #( =1.648721)m0<- function(Ti,Tj) #a function that generates the MLestimates{ loglik<-function(ti,t

[R] parametric bootstrap

2008-07-25 Thread laila khalfan
Hi I am using the boot function in boot package. I am facing a problem in getting the values of bias and st.error in the output. r<-36n<-40shape<-2theta11<-exp(1) # (=2.718282)theta21<-exp(.5) #( =1.648721) m0<- function(Ti,Tj) #a function that generates the estimates{ loglik<-function(ti,tj,

[R] parametric bootstrap

2008-07-25 Thread laila khalfan
Hi I am trying to find a parametric bootstrap confidence interval and when I used the boot function I get zero bias and zero st.error? What could be my mistake? Thank you and take care. Laila [[alternative HTML version deleted]] __ R-help@r

[R] need help in parametric bootstrap

2008-07-25 Thread laila khalfan
bootstrap Confidence intervals. I know that I should use first the boot function with sim='parametric' and ran.gen=?? but I dont understand what is the difference between the ran.gen and the statistics. Thank you and take care. Laila [[alternative HTML versi

[R] For loop

2008-07-03 Thread laila khalfan
HiI have a specific sample coming from a gamma(alpha,theta1) distribution and then divided into two parts first part follows a gamma(alpha,theta1) the second is gamma(alpha,theta2) then I would like to find the mle`s for theta1 and theta2 which I found. Now I would like to simulate those estimat

[R] FW: For loop

2008-07-03 Thread laila khalfan
HiI have a specific sample coming from a gamma(alpha,theta1) distribution and then divided into two parts first part follows a gamma(alpha,theta1) the second is gamma(alpha,theta2) then I would like to find the mle`s for theta1 and theta2 which I found. Now I would like to simulate those es

[R] For loop

2008-07-03 Thread laila khalfan
matrix(Tj,ncol=1) cc<-data.frame(rbind(c1,c2))[,1] cc # the special sample that I need to find the! mle`s for estimates<- as.data.frame(t(m0(cc))) }estimates Thanks in advance Laila _