> require(MASS)
> y=c(49,-67,8,16,6,23,28,41,14,29,56,24,75,60,-48)
> a <- fitdistr(y,"t")
> a
m s df
26.847104 22.705444 2.502808
( 7.253774) ( 9.024268) ( 2.189144)
> x <- read.table("/home/matheus/Downloads/abbey-shares.txt", header=TRUE,
+ sep="", na.strings="NA", dec=".", strip.white=TRUE)
> x
X0.0000
1 0.0135
2 0.0067
3 -0.0066
4 0.0133
5 -0.0033
6 -0.0132
7 -0.0201
8 -0.0034
9 0.0000
10 -0.0034
11 0.0068
12 -0.0271
13 0.0035
14 0.0312
15 0.0269
16 0.0066
17 0.0000
18 -0.0098
19 -0.0033
20 0.0033
21 0.0000
22 0.0164
23 0.0000
24 0.0000
25 -0.0065
26 -0.0033
27 -0.0065
28 -0.0132
29 -0.0133
30 0.0169
32 -0.0101
33 0.0000
34 -0.0068
35 -0.0034
36 0.0171
37 -0.0101
38 -0.0034
39 0.0424
40 -0.0098
41 -0.0066
42 0.0066
43 0.0165
44 -0.0097
45 -0.0098
46 -0.0033
47 -0.0133
48 0.0067
> b <- fitdistr(x,"t")
> b pressiono o enter e nada acontece, não me dá a estimação
Em 8 de novembro de 2012 09:48, Leonard Mendonça de Assis
<assis.leonard@gmail.com > escreveu:
Matheus,
sem um dump deste seu suposto X, nem um CMR, fica difícil saber o
que está acontecendo
[]s
Leonard de Assis
http://about.me/ldeassis
Em 08/11/2012 03:54, Matheus Monteiro escreveu:
é o seguinte:y=c(49,-67,8,16,6,23,28,41,14, 29,56,24,75,60,-48) fitdistr
(y, "t")
_______________________________________________
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.
R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br http://www.leg.ufpr.br/r-br-guia ) e forneça código mínimo reproduzível.