[R] KhmaladzeTest
raul sanchez
r_o_sanchez2003 at yahoo.com.ar
Sat Jul 8 02:39:13 CEST 2006
Hello. I am a beginer in R and I can not implement the KhmaladzeTest in the following command. Please help me!!!!!!!!!!!
PD: I attach thw results and the messages of the R program
R : Copyright 2006, The R Foundation for Statistical Computing
Version 2.3.1 (2006-06-01)
ISBN 3-900051-07-0
R es un software libre y viene sin GARANTIA ALGUNA.
Usted puede redistribuirlo bajo ciertas circunstancias.
Escriba 'license()' o 'licence()' para detalles de distribucion.
R es un proyecto colaborativo con muchos contribuyentes.
Escriba 'contributors()' para obtener mas informacion y
'citation()' para saber como citar R o paquetes de R en publicaciones.
Escriba 'demo()' para demostraciones, 'help()' para el sistema on-line
de ayuda,
o 'help.start()' para abrir el sistema de ayuda HTML con su navegador.
Escriba 'q()' para salir de R.
> utils:::menuInstallLocal()
package 'quantreg' successfully unpacked and MD5 sums checked
updating HTML package descriptions
> utils:::menuInstallLocal()
package 'foreign' successfully unpacked and MD5 sums checked
updating HTML package descriptions
> utils:::menuInstallLocal()
package 'Rcmdr' successfully unpacked and MD5 sums checked
updating HTML package descriptions
> local({pkg <- select.list(sort(.packages(all.available = TRUE)))
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
> local({pkg <- select.list(sort(.packages(all.available = TRUE)))
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
quantreg package loaded: To cite see citation("quantreg")
> local({pkg <- select.list(sort(.packages(all.available = TRUE)))
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
> local({pkg <- select.list(sort(.packages(all.available = TRUE)))
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
Loading required package: tcltk
Loading Tcl/Tk interface ... done
--- Please select a CRAN mirror for use in this session ---
also installing the dependencies 'acepack', 'scatterplot3d', 'fBasics',
'Hmisc', 'quadprog', 'oz', 'mlbench', 'randomForest', 'SparseM',
'xtable', 'chron', 'fCalendar', 'its', 'tseries', 'DAAG', 'e1071', 'mvtnorm',
'zoo', 'strucchange', 'sandwich', 'dynlm', 'leaps'
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/acepack_1.3-2.2.zip'
Content
type 'application/zip' length 55667 bytes
URL abierta
downloaded 54Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/scatterplot3d_0.3-24.zip'
Content
type 'application/zip' length 540318 bytes
URL abierta
downloaded 527Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/fBasics_221.10065.zip'
Content
type 'application/zip' length 3327499 bytes
URL abierta
downloaded 3249Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/Hmisc_3.0-12.zip'
Content
type 'application/zip' length 1993038 bytes
URL abierta
downloaded 1946Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/quadprog_1.4-8.zip'
Content
type 'application/zip' length 38626 bytes
URL abierta
downloaded 37Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/oz_1.0-13.zip'
Content
type 'application/zip' length 39640 bytes
URL abierta
downloaded 38Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/mlbench_1.1-1.zip'
Content
type 'application/zip' length 1324913 bytes
URL abierta
downloaded 1293Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/randomForest_4.5-16.zip'
Content
type 'application/zip' length 209710 bytes
URL abierta
downloaded 204Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/SparseM_0.68.zip'
Content
type 'application/zip' length 728594 bytes
URL abierta
downloaded 711Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/xtable_1.3-2.zip'
Content
type 'application/zip' length 56703 bytes
URL abierta
downloaded 55Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/chron_2.3-4.zip'
Content
type 'application/zip' length 101287 bytes
URL abierta
downloaded 98Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/fCalendar_221.10065.zip'
Content
type 'application/zip' length 754551 bytes
URL abierta
downloaded 736Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/its_1.1.1.zip'
Content
type 'application/zip' length 194287 bytes
URL abierta
downloaded 189Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/tseries_0.10-3.zip'
Content
type 'application/zip' length 392748 bytes
URL abierta
downloaded 383Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/DAAG_0.79.zip'
Content
type 'application/zip' length 796328 bytes
URL abierta
downloaded 777Kb
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/e1071_1.5-13.zip'
Content
type 'application/zip' length 627972 bytes
URL abierta
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/mvtnorm_0.7-2.zip'
Content
type 'application/zip' length 214609 bytes
URL abierta
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/zoo_1.1-0.zip'
Content
type 'application/zip' length 695387 bytes
URL abierta
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/strucchange_1.2-13.zip'
Content
type 'application/zip' length 866668 bytes
URL abierta
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/sandwich_1.9-0.zip'
Content
type 'application/zip' length 675607 bytes
URL abierta
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/dynlm_0.1-2.zip'
Content
type 'application/zip' length 43966 bytes
URL abierta
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/leaps_2.7.zip'
Content
type 'application/zip' length 73301 bytes
URL abierta
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/relimp_0.9-6.zip'
Content
type 'application/zip' length 70677 bytes
URL abierta
probando la URL
'http://cran.au.r-project.org/bin/windows/contrib/2.3/multcomp_0.4-8.zip'
Content
type 'application/zip' length 310081 bytes
URL abierta
Versión del Rcmdr 1.1-7
> Datos <- read.dta("C:/Mis documentos/r1-r1.dta", convert.dates=TRUE,
convert.factors=TRUE, missing.type=TRUE, convert.underscore=TRUE,
warn.missing.labels=TRUE)
Aviso en read.dta("C:/Mis documentos/r1-r1.dta", convert.dates = TRUE,
:
'missing.type' only applicable to version 8 files
> fit1 <- rq(yyy ~ lll, tau = 0.5, data = Datos)
> fit1
Call:
rq(formula = yyy ~ lll, tau = 0.5, data = Datos)
Coefficients:
(Intercept) lll
4.249606 -0.929411
Degrees of freedom: 36 total; 34 residual
> summary(fit1)
Call: rq(formula = yyy ~ lll, tau = 0.5, data = Datos)
tau: [1] 0.5
Coefficients:
coefficients lower bd upper bd
(Intercept) 4.24961 2.36774 4.95617
lll -0.92941 -1.08243 -0.52124
> r1 <- resid(fit1)
> c1 <- coef(fit1)
> summary(fit1, se = "nid")
Call: rq(formula = yyy ~ lll, tau = 0.5, data = Datos)
tau: [1] 0.5
Coefficients:
Value Std. Error t value Pr(>|t|)
(Intercept) 4.24961 1.05525 4.02710 0.00046
lll -0.92941 0.22934 -4.05246 0.00043
> fit <- rqProcess( y.net ~ lgdp2 + fse2 + gedy2 + Iy2 + gcony2,
+ data = barro, taus = seq(.1,.9,by = .05))
Erro en inherits(x, "data.frame") : objeto "barro" no encontrado
> fit <- rqProcess(yyy ~ lll, data = Datos, taus = seq(.1,.9, by =
.05))
taus: 0.1 Erro en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
tau - h < 0: error in summary.rq
> f <- rqProcess(yyy ~ lll, data = Datos, taus = seq(.1,.9, by = .05))
taus: 0.1 Erro en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
tau - h < 0: error in summary.rq
> fit1 <- rqProcess(yyy ~ lll, data = Datos, taus = seq(.1,.9, by =
.05))
taus: 0.1 Erro en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
tau - h < 0: error in summary.rq
> T1 <- KhmaladzeTest(yyy ~ lll, taus = -1, nullH = "location")
Erro en eval(expr, envir, enclos) : objeto "yyy" no encontrado
> Datos <- read.dta("C:/Mis documentos/r1-r1.dta", convert.dates=TRUE,
convert.factors=TRUE, missing.type=TRUE, convert.underscore=TRUE,
warn.missing.labels=TRUE)
Aviso en read.dta("C:/Mis documentos/r1-r1.dta", convert.dates = TRUE,
:
'missing.type' only applicable to version 8 files
> T1 <- KhmaladzeTest(yyy ~ lll, taus = -1, nullH = "location")
Erro en eval(expr, envir, enclos) : objeto "yyy" no encontrado
> T1 <- KhmaladzeTest( yyy ~ lll, taus = -1, nullH = "location")
Erro en eval(expr, envir, enclos) : objeto "yyy" no encontrado
> data(barro)
> invisible(edit(barro))
> rqProcess(yyy ~ lll, data = Datos, taus = -1, nullH = "location")
$taus
[1] 0.3161186 0.3543310 0.4007676 0.4150787 0.4268406 0.4542215
0.5008570 0.5304535 0.5690997
[10] 0.6206522 0.6522563
$qtaus
[1] -0.037757189 -0.033227447 -0.032656681 -0.016152701 -0.015090004
-0.014413058
[7] -0.013738239 -0.013471008 -0.012253429 -0.010675005 -0.009866834
$Vhat
[,1] [,2] [,3] [,4] [,5] [,6]
[,7] [,8]
[1,] 0.3064789 0.5104205 0.6209815 -0.0850221 -0.2053440 -0.1631580
-0.1766553 -0.1768334
[,9] [,10] [,11]
[1,] -0.1483747 -0.1191922 -0.1026903
$vhat
[,1] [,2] [,3] [,4] [,5] [,6]
[,7] [,8]
[1,] 0.3064789 0.5104205 0.6209815 -0.0850221 -0.2053440 -0.1631580
-0.1766553 -0.1768334
[,9] [,10] [,11]
[1,] -0.1483747 -0.1191922 -0.1026903
attr(,"class")
[1] "rqProcess"
> rqProcess(yyy ~ lll, data = Datos,taus = seq(.1,.9, by = .05), nullH
= "location")
taus: 0.1 Erro en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
tau - h < 0: error in summary.rq
> KhmaladzeTest( yyy ~ lll, data = Datos, taus = seq(.1,.9, by = .05))
taus: 0.1 Erro en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
tau - h < 0: error in summary.rq
> T1 <- KhmaladzeTest(yyy ~ lll, taus = -1, nullH = "location")
Erro en eval(expr, envir, enclos) : objeto "yyy" no encontrado
> fit1 <- rqProcess(yyy ~ lll, data = Datos, taus = seq(.1,.9, by =
.05))
taus: 0.1 Erro en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
tau - h < 0: error in summary.rq
> KhmaladzeTest( yyy ~ lll, data = Datos, taus = -1)
Erro en qr(X[1:p, ]) : NA/NaN/Inf en llamada a una función externa (arg
1)
> KhmaladzeTest( yyy ~ lll, data = Datos, taus = -1, nullH =
"location")
Erro en qr(X[1:p, ]) : NA/NaN/Inf en llamada a una función externa (arg
1)
> data(barro)
> KhmaladzeTest( y.net ~ lgdp2 + fse2 + gedy2 + Iy2 + gcony2,
+ data = barro, taus = seq(.05,.95,by = .01))
taus: 0.05 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.06 0.07 0.08 0.09 Aviso en summary.rq(rq(formula, data = data,
tau = taus[i], method = "fn"), :
1 non-positive fis
0.1 Aviso en summary.rq(rq(formula, data = data, tau = taus[i], method
= "fn"), :
1 non-positive fis
0.11 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
2 non-positive fis
0.12 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
3 non-positive fis
0.13 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
3 non-positive fis
0.14 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
4 non-positive fis
0.15 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
3 non-positive fis
0.16 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
5 non-positive fis
0.17 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
4 non-positive fis
0.18 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
5 non-positive fis
0.19 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
4 non-positive fis
0.2 Aviso en summary.rq(rq(formula, data = data, tau = taus[i], method
= "fn"), :
4 non-positive fis
0.21 0.22 0.23 0.24 0.25 Aviso en summary.rq(rq(formula, data =
data, tau = taus[i], method = "fn"), :
1 non-positive fis
0.26 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.27 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.28 0.29 0.3 0.31 0.32 0.33 Aviso en summary.rq(rq(formula, data
= data, tau = taus[i], method = "fn"), :
1 non-positive fis
0.34 0.35 0.36 0.37 0.38 0.39 0.4 0.41 0.42 0.43 0.44 0.45
0.46 0.47 0.48 0.49 0.5 0.51 0.52 0.53 0.54 0.55 0.56 0.57
Aviso en summary.rq(rq(formula, data = data, tau = taus[i], method =
"fn"), :
1 non-positive fis
0.58 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
2 non-positive fis
0.59 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.6 Aviso en summary.rq(rq(formula, data = data, tau = taus[i], method
= "fn"), :
1 non-positive fis
0.61 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.62 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.63 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.64 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.65 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.66 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.67 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.68 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.69 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
2 non-positive fis
0.7 Aviso en summary.rq(rq(formula, data = data, tau = taus[i], method
= "fn"), :
2 non-positive fis
0.71 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
2 non-positive fis
0.72 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
1 non-positive fis
0.73 0.74 0.75 0.76 0.77 0.78 0.79 Aviso en
summary.rq(rq(formula, data = data, tau = taus[i], method = "fn"), :
2 non-positive fis
0.8 0.81 0.82 0.83 0.84 Aviso en summary.rq(rq(formula, data =
data, tau = taus[i], method = "fn"), :
2 non-positive fis
0.85 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
3 non-positive fis
0.86 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
3 non-positive fis
0.87 0.88 0.89 Aviso en summary.rq(rq(formula, data = data, tau =
taus[i], method = "fn"), :
1 non-positive fis
0.9 0.91 0.92 Aviso en summary.rq(rq(formula, data = data, tau =
taus[i], method = "fn"), :
2 non-positive fis
0.93 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
9 non-positive fis
0.94 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
13 non-positive fis
0.95 Aviso en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
14 non-positive fis
$nullH
[1] "location"
$Tn
[1] 6.029988
$THn
lgdp2 fse2 gedy2 Iy2 gcony2
1.3198411 1.1505847 1.2480321 0.8466474 1.0036528
attr(,"class")
[1] "KhmaladzeTest"
> data(Datos)
Aviso en data(Datos) : data set 'Datos' not found
> Datos <- read.dta("C:/Mis documentos/r1-r1.dta", convert.dates=TRUE,
convert.factors=TRUE, missing.type=TRUE, convert.underscore=TRUE,
warn.missing.labels=TRUE)
Aviso en read.dta("C:/Mis documentos/r1-r1.dta", convert.dates = TRUE,
:
'missing.type' only applicable to version 8 files
> data(Datos)
Aviso en data(Datos) : data set 'Datos' not found
> Datos <- read.dta("C:/Mis documentos/r1-r1.dta", convert.dates=TRUE,
convert.factors=TRUE, missing.type=TRUE, convert.underscore=TRUE,
warn.missing.labels=TRUE)
Aviso en read.dta("C:/Mis documentos/r1-r1.dta", convert.dates = TRUE,
:
'missing.type' only applicable to version 8 files
> KhmaladzeTest( yyy ~ lll, data = Datos, taus = -1, nullH =
"location")
Erro en qr(X[1:p, ]) : NA/NaN/Inf en llamada a una función externa (arg
1)
> invisible(edit(Datos))
> local({pkg <- select.list(sort(.packages(all.available = TRUE)))
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
> fit <- rqProcess( yyy ~ lll, data = Datos, taus = seq(.1,.9,by =
.05))
> fff <- rqProcess( yyy ~ lll, data = Datos, taus = seq(.1,.9,by =
.05))
taus: 0.1 Erro en summary.rq(rq(formula, data = data, tau = taus[i],
method = "fn"), :
tau - h < 0: error in summary.rq
> fff <- rqProcess(yyy ~ lll, data = Datos, taus = -1, nullH =
"location")
> khmaladze.test(fit, nullH = "location")
Error: no se pudo encontrar la función "khmaladze.test"
> khmaladze.test(fff, nullH = "location")
Error: no se pudo encontrar la función "khmaladze.test"
> Khmaladze.test(fff, nullH = "location")
Error: no se pudo encontrar la función "Khmaladze.test"
> KhmaladzeTest(fff, nullH = "location")
Erro en terms.default(formula, data = data) :
no terms component
> KhmaladzeTest(fff,data = Datos)
Erro en terms.default(formula, data = data) :
no terms component
> KhmaladzeTest(fff, data = Datos)
Erro en terms.default(formula, data = data) :
no terms component
> KhmaladzeTest(fff, data = Datos, nullH = "location")
Erro en terms.default(formula, data = data) :
no terms component
> khmaladze.test(fit, nullH = "location")
Error: no se pudo encontrar la función "khmaladze.test"
> KhmaladzeTest(fff, data = Datos, taus = -1, nullH = "location")
Erro en terms.default(formula, data = data) :
no terms component
> help.search("KhmaladzeTest")
> KhmaladzeTest( yyy ~ lll, data = Datos, taus = -1, nullH =
"location" , trim = c(0.05, 0.95))
Erro en qr(X[1:p, ]) : NA/NaN/Inf en llamada a una función externa (arg
1)
> KhmaladzeTest( fff, data = Datos, taus = -1, nullH = "location" ,
trim = c(0.05, 0.95))
Erro en terms.default(formula, data = data) :
no terms component
> ls()
[1] "barro" "c1" "Datos" "fff" "fit1" "GLM.1" "r1"
> summary(fff)
Length Class Mode
taus 11 -none- numeric
qtaus 11 -none- numeric
Vhat 11 -none- numeric
vhat 11 -none- numeric
> KhmaladzeTest(formula, data = NUL, taus = -1, nullH = "location"
+ >
> KhmaladzeTest(fff, data = Datos, taus = -1, nullH = "location")
Erro en terms.default(formula, data = data) :
no terms component
> KhmaladzeTest(fff, data = Datos)
Erro en terms.default(formula, data = data) :
no terms component
> KhmaladzeTest( fff, data = Datos, taus = -1, nullH = "location" ,
trim = c(0.05, 0.95))
Erro en terms.default(formula, data = data) :
no terms component
> KhmaladzeTest( fit1, data = Datos, taus = -1, nullH = "location" ,
trim = c(0.05, 0.95))
Erro en qr(X[1:p, ]) : NA/NaN/Inf en llamada a una función externa (arg
1)
> fix(Datos)
---------------------------------
1GB gratis, Antivirus y Antispam
Correo Yahoo!, el mejor correo web del mundo
Abrí tu cuenta aquí
More information about the R-help
mailing list