Atualização do R no Linux

Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave? André Oliveira Souza

O que eu uso e', antes de instalar a nova versao: writeLines(installed.packages()[,1], con="packages.dump") e depois de instalar (ou recompilar o R), assumindo que vc nao apagou o "packages.dump". source("http://www.bioconductor.org/biocLite.R") packs <- readLines("packages.dump") biocLite(packs) b ps: eu dependo de muitas coisas do BioC, se vc nao usa, basta ignorar a linha source() e substituir biocLite() por install.packages(). 2011/4/19 Andre Oliveira <andreolsouza@yahoo.com.br>:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave?
André Oliveira Souza
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)

Andre nao sei como0 atualizou, o sistema que usa etc etc Uma atuualizacao nao deveria remover os pacotres.. eles forma removidos? se nao removeu: update.packages(ask=F) se removeu: quer istalar o que? todos? se for...poderia usar options(repos="http://cran.r-project.org/") update.packages(ask=F) x <- packageStatus() st <- x$avai["Status"] install.packages(rownames(st)[which(st$Status=="not installed")]) Por gentileza procure ser mais claro nas mensagens On Tue, 19 Apr 2011, Andre Oliveira wrote:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave? André Oliveira Souza

Eu atualizei por sudo apt-get update && sudo apt-get upgrade no terminal. Os pacotes da versão 2.12 ainda estão na pasta /home/andre/R/i486-pc-linux-gnu-library/2.12 só que depois da atualização foi criada uma pasta nova /home/andre/R/i486-pc-linux-gnu-library/2.13. Seria só copiar os arquivos da pasta 2.12 para a pasta 2.13? abs obrigado André Oliveira Souza ________________________________ De: Paulo Justiniano <paulojus@leg.ufpr.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Terça-feira, 19 de Abril de 2011 9:15:00 Assunto: Re: [R-br] Atualização do R no Linux Andre nao sei como0 atualizou, o sistema que usa etc etc Uma atuualizacao nao deveria remover os pacotres.. eles forma removidos? se nao removeu: update.packages(ask=F) se removeu: quer istalar o que? todos? se for...poderia usar options(repos="http://cran.r-project.org/") update.packages(ask=F) x <- packageStatus() st <- x$avai["Status"] install.packages(rownames(st)[which(st$Status=="not installed")]) Por gentileza procure ser mais claro nas mensagens On Tue, 19 Apr 2011, Andre Oliveira wrote:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave?
André Oliveira Souza

vc pode ate' copiar, mas vai precisar definitivamente executar update.packages(ask=FALSE) o que vai baixar todos os pacotes do mesmo jeito... b 2011/4/19 Andre Oliveira <andreolsouza@yahoo.com.br>:
Eu atualizei por sudo apt-get update && sudo apt-get upgrade no terminal.
Os pacotes da versão 2.12 ainda estão na pasta /home/andre/R/i486-pc-linux-gnu-library/2.12 só que depois da atualização foi criada uma pasta nova /home/andre/R/i486-pc-linux-gnu-library/2.13.
Seria só copiar os arquivos da pasta 2.12 para a pasta 2.13?
abs obrigado
André Oliveira Souza
________________________________ De: Paulo Justiniano <paulojus@leg.ufpr.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Terça-feira, 19 de Abril de 2011 9:15:00 Assunto: Re: [R-br] Atualização do R no Linux
Andre
nao sei como0 atualizou, o sistema que usa etc etc
Uma atuualizacao nao deveria remover os pacotres.. eles forma removidos?
se nao removeu: update.packages(ask=F)
se removeu: quer istalar o que? todos?
se for...poderia usar
options(repos="http://cran.r-project.org/") update.packages(ask=F) x <- packageStatus() st <- x$avai["Status"] install.packages(rownames(st)[which(st$Status=="not installed")])
Por gentileza procure ser mais claro nas mensagens
On Tue, 19 Apr 2011, Andre Oliveira wrote:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave?
André Oliveira Souza
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)

Fiz o teste aqui somente fiz uma cópia dos arquivos da pasta 2.12 para 2.13 e esta tudo funcionando. obrigado a todos André Oliveira Souza ________________________________ De: Benilton Carvalho <beniltoncarvalho@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Terça-feira, 19 de Abril de 2011 12:02:02 Assunto: Re: [R-br] Res: Atualização do R no Linux vc pode ate' copiar, mas vai precisar definitivamente executar update.packages(ask=FALSE) o que vai baixar todos os pacotes do mesmo jeito... b 2011/4/19 Andre Oliveira <andreolsouza@yahoo.com.br>:
Eu atualizei por sudo apt-get update && sudo apt-get upgrade no terminal.
Os pacotes da versão 2.12 ainda estão na pasta /home/andre/R/i486-pc-linux-gnu-library/2.12 só que depois da atualização foi criada uma pasta nova /home/andre/R/i486-pc-linux-gnu-library/2.13.
Seria só copiar os arquivos da pasta 2.12 para a pasta 2.13?
abs obrigado
André Oliveira Souza
________________________________ De: Paulo Justiniano <paulojus@leg.ufpr.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Terça-feira, 19 de Abril de 2011 9:15:00 Assunto: Re: [R-br] Atualização do R no Linux
Andre
nao sei como0 atualizou, o sistema que usa etc etc
Uma atuualizacao nao deveria remover os pacotres.. eles forma removidos?
se nao removeu: update.packages(ask=F)
se removeu: quer istalar o que? todos?
se for...poderia usar
options(repos="http://cran.r-project.org/") update.packages(ask=F) x <- packageStatus() st <- x$avai["Status"] install.packages(rownames(st)[which(st$Status=="not installed")])
Por gentileza procure ser mais claro nas mensagens
On Tue, 19 Apr 2011, Andre Oliveira wrote:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave?
André Oliveira Souza
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Andre, muitos pacotes usam codigos compilados. Ao simplesmente copiar, voce esta' assumindo que toda a "estrutura interna" continua exatamente a mesma apos a mudanca.... e isso nao e' verdade (se fosse, vc nao precisaria atualizar o R). Os pacotes que nao usam codigo compilado tem grandes chances (nao 100%) de funcionar... mas os binarios compilados podem depender de porcoes que sofreram modificacoes, dai' tudo pode acontecer (inclusive o codigo rodar, dar resultados errados e vc nem ver erro). b 2011/4/19 Andre Oliveira <andreolsouza@yahoo.com.br>:
Fiz o teste aqui somente fiz uma cópia dos arquivos da pasta 2.12 para 2.13 e esta tudo funcionando.
obrigado a todos
André Oliveira Souza
________________________________ De: Benilton Carvalho <beniltoncarvalho@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Terça-feira, 19 de Abril de 2011 12:02:02 Assunto: Re: [R-br] Res: Atualização do R no Linux
vc pode ate' copiar, mas vai precisar definitivamente executar
update.packages(ask=FALSE)
o que vai baixar todos os pacotes do mesmo jeito...
b
2011/4/19 Andre Oliveira <andreolsouza@yahoo.com.br>:
Eu atualizei por sudo apt-get update && sudo apt-get upgrade no terminal.
Os pacotes da versão 2.12 ainda estão na pasta /home/andre/R/i486-pc-linux-gnu-library/2.12 só que depois da atualização foi criada uma pasta nova /home/andre/R/i486-pc-linux-gnu-library/2.13.
Seria só copiar os arquivos da pasta 2.12 para a pasta 2.13?
abs obrigado
André Oliveira Souza
________________________________ De: Paulo Justiniano <paulojus@leg.ufpr.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Terça-feira, 19 de Abril de 2011 9:15:00 Assunto: Re: [R-br] Atualização do R no Linux
Andre
nao sei como0 atualizou, o sistema que usa etc etc
Uma atuualizacao nao deveria remover os pacotres.. eles forma removidos?
se nao removeu: update.packages(ask=F)
se removeu: quer istalar o que? todos?
se for...poderia usar
options(repos="http://cran.r-project.org/") update.packages(ask=F) x <- packageStatus() st <- x$avai["Status"] install.packages(rownames(st)[which(st$Status=="not installed")])
Por gentileza procure ser mais claro nas mensagens
On Tue, 19 Apr 2011, Andre Oliveira wrote:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave?
André Oliveira Souza
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)

Pode ser sim e depois roda um update.packages() On Tue, 19 Apr 2011, Andre Oliveira wrote:
Eu atualizei por sudo apt-get update && sudo apt-get upgrade no terminal.
Os pacotes da versão 2.12 ainda estão na pasta /home/andre/R/i486-pc-linux-gnu-library/2.12 só que depois da atualização foi criada uma pasta nova /home/andre/R/i486-pc-linux-gnu-library/2.13.
Seria só copiar os arquivos da pasta 2.12 para a pasta 2.13?
abs obrigado
André Oliveira Souza
_______________________________________________________________________________________________________________________________________ De: Paulo Justiniano <paulojus@leg.ufpr.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Terça-feira, 19 de Abril de 2011 9:15:00 Assunto: Re: [R-br] Atualização do R no Linux
Andre
nao sei como0 atualizou, o sistema que usa etc etc
Uma atuualizacao nao deveria remover os pacotres.. eles forma removidos?
se nao removeu: update.packages(ask=F)
se removeu: quer istalar o que? todos?
se for...poderia usar
options(repos="http://cran.r-project.org/") update.packages(ask=F) x <- packageStatus() st <- x$avai["Status"] install.packages(rownames(st)[which(st$Status=="not installed")])
Por gentileza procure ser mais claro nas mensagens
On Tue, 19 Apr 2011, Andre Oliveira wrote:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave? André Oliveira Souza

On Tue, 2011-04-19 at 04:28 -0700, Andre Oliveira wrote:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave?
André Oliveira Souza
Andre Qual o sistema operacional que você usa?
-- []s Tura

Ubuntu LTS 10.04 obrigado André Oliveira Souza ________________________________ De: Bernardo Rangel Tura <tura@centroin.com.br> Para: r-br@listas.c3sl.ufpr.br Enviadas: Quarta-feira, 20 de Abril de 2011 6:36:55 Assunto: Re: [R-br] Atualização do R no Linux On Tue, 2011-04-19 at 04:28 -0700, Andre Oliveira wrote:
Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave?
André Oliveira Souza
Andre Qual o sistema operacional que você usa?
-- []s Tura _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

On Wed, 2011-04-20 at 11:11 -0700, Andre Oliveira wrote:
Ubuntu LTS 10.04
obrigado
André Oliveira Souza
Andre, Eu estou usando o a versão 10.10 mas acho que não há grande diferença. O detalhe é que quando: instalo o faço como root (sudo make install) deste forma os pacotes ficam sempre no mesmo diretório. De todo o jeito se você copiar o diretório depois basta update.packages(ask=FLASE,checkBuld=TRUE) -- []s Tura

Olá pessoal, Estou precisando de ajuda para entender e consertar o erro que surge na minha programação. Eu tenho um procedimento iterativo que em determinado ponto a função de otimização dá o seguinte erro: Error in optim(par = theta.old, fn = fun, gr = gradient, control = control, : initial value in 'vmmin' is not finite Quando eu testo a função no ponto inicial ela produz valor igual a -1035.843. Não sei bem o que está acontecendo e por enquanto estou sem idéias de como resolver. Alguém poderia me ajudar ? Segue parte da programação para quem quiser tentar reproduzir. #Inicio library(maxLik) library(fGarch) serie=c(-0.079365084,-0.069496156,0.751957634,0.589682298,0.137093637,-0.441 328623,-0.315022903,-0.712451064,-1.018074469,0.040120362,0.91835343,0.58453 6266,0.108604444,0.433242111,0.058933309,0.225589817,0.693192915,0.213820667 ,1.187100752,0.172546055,-0.258930857,0.115163161,0.229929206,0.601060818,0. 408183205,-0.056855872,0.359542444,-0.103945202,-0.378895973,-0.247054478,-0 .390830296,-0.095556625,0.448324307,0.12365056,0.227877045,0.501396443,0.348 547289,0.487805845,0.28967922,0.363247268,0.463779702,0.295694172,0.45109397 6,-0.018372221,0.21107703,0.557664745,0.245845787,0.181719112,-0.573537008,- 0.009130336,0.428227159,0.634404883,0.39671858,0.502648936,0.606982137,0.080 060494,-0.142374112,-0.849546804,0.376479469,-0.080554938,-0.55670432,-0.360 815834,-0.461936422,0.163280152,0.036248301,0.162954952,0.171720409,-1.44620 284,-0.689499399,-0.064596505,0.359364589,0.119501784,-0.395821385,-0.852172 626,0.027903083,0.065076933,0.602271915,0.571114301,0.63178345,0.600602406,0 .398406901,0.090326083,-0.670353086,-0.20927172,0.245622138,0.127122509,-0.5 7332823,-0.1918422,0.073126146,0,-1.176484158,-1.546662715,-0.952882346,0.47 2858835,0.357883222,0.178462453,-0.103281545,0.720740706,0.539486008,0.61959 7858,0.276192405,0.422018975,0.71155199,0.37200061,0.442778476,0.297096775,- 0.432433106,-0.379919046,0.343798403,-0.670596079,-0.556190173,0.38328209,-0 .045551862,-0.182415228,0.082128033,-0.255731256,-0.715862088,-1.036950316,- 1.538059017,-0.511655465,0.303519156,-1.123821352,-1.107537685,-0.145369992, 0.647689441,-0.579823812,-0.135777346,1.004257185,1.079851862,0.218371793,2. 269201915,0.987138419,-0.496964027,-0.35120184,-0.111162587,0.342355181,0.31 3566613,0.229938018,0.110182731,0.302378019,-0.541260897,-0.267121238,0.2119 13287,0.28491357,0.931768868,0.498844624,0.045226359,-0.027133361,-0.1176630 45,-0.18129084,-0.409110076,0.018218255,-0.009108712,0.291068061,-0.09995911 6,-0.154679072,-0.758652695,-0.174480048,0.073502392,-0.64498519,0.184706369 ,0.202783736,0.395166622,0.503180239,0.346178719,0.027278927,0.625370235,0.2 70684998,-0.21648935,-0.027094152,0.180489175,-0.198556022,-0.743565268,-0.2 4604742,-0.493918502,-0.985955925,-1.314275338,0.028143909,0.178061056,-0.47 8671898,-0.566092652,-1.343069196,0.210748234,0.572520648,-1.350659931,-0.14 476671,1.409757224,0.522392792,-0.094777753,-0.028450851,0.274686415,-0.8167 96259,-1.073835822,0.163737094,-0.45334054,0.260680814,-0.067518691,-0.19316 2123,-0.474830099,-1.133156601,-1.464647657,0.089681636,1.227495119,1.406822 733,0.859701021,0.985795417,1.089283799,1.040260541,-0.317490228,0.317490228 ,0.418663681,0.24109804,0.77498312,0.321145388,-0.08248179,0.347794603,0.510 344678,0.46252293,0.882809084,0.206065566,0.214573171,-0.026796481,-0.214630 739,0.205696983,0.570157447,0.354704637,1.012822961,-0.412697605,0.246067441 ,0.786854908,0.772070632,0.284765261,0.043075598,0.464037956,0.222660022,-0. 411452663,0.274489793,0.162622538,0.622364428,-0.178639833,-0.1107467,-0.093 804637,-0.616176532, -0.163237289,0.077356143,0.838471716,0.907516512,0.362395576,0.595514446,0.7 41515583,0.157591373,0.479537918,-0.082515064,0.173203062,0.665272874,0.4492 18708,0.122164775,-0.489557114,0.285866164,0.601725673,0.072936508,-0.032409 658,-0.08107013,0.251103773,0.315011771,0.667124606,0.08808104,0.351550378,0 .723172381,0.229366969,0.480751212,0.541029273,-0.039106801,0.444913508,0.34 9854584,-0.233100339,-0.945725912,0.305798647,0.320475516,-0.289164361,0.187 661326,-0.375675483,-1.421488043,0.499822533,0.803976901,0.133369964,0.00783 9755,0.12535257,-0.148873681,0.477178517,-0.062451212,0.560661802,0.94289283 5,0.705199918,0.669917244,-0.03035362,0.537411883,0.527030282,-0.08263532,0. 614372203,-0.104625972,0.268817366,-1.064557621,0.135572816,0.112837113,-0.8 38025364,-0.806885705,0.061124696,0.068720653,0.243958343,-0.894877669,-0.63 1986689,-0.441672873,0.116418966,0.394813751,-0.627787375,-2.232797159,0.784 007184,0.926514166,1.15778907,0.692896759,0.015343307,-0.46132636,0.81357421 7,1.019097724,0.166326492,0.256507119,0.56350876,0.246941421,0.11204483,-0.3 73972014,-0.691784096,0.451705955,0.292518491,-0.585895164,-0.559040513,-0.7 60459939,0.053420843,-0.382205021,-0.145621792,0.306325864,0.511002904,-0.12 1793423,-0.007617016,-0.396886493,0.236794978,0.236235583,-0.450090396,-0.43 6765569,-0.70127055,0.370513739,-0.161930868,-0.177654239,0.562773717,0.3912 99843,0.526980255,0.44083069,-0.273390207,-0.495522264,-0.29083137,-0.276306 876,-0.833531105,0.270887515,-0.760341936,-0.672035037,0.180173172,-0.533627 571,0.42399561,0.889241429,-0.233245324,0.140012468,-0.225672249,-0.16373634 3,-0.312622373,0.546449447,0.729137577,0.708512783,0.482297586,0.022908633,- 0.275250575,0.107131936,0.068809973,0.639953426,0.015187182,-0.167185995,-0. 015212596,-0.434534703,0.411711469,-0.106577355,0,-0.144822618,-0.535374128, -0.069042234,-0.415257441,-0.254619944,0.262325588,0.130905189,0.575485954,- 0.283492509,-0.538545693,-0.324525088,-0.061933888,-0.676019514,-0.845604788 ,0.259138706,0.054882592,0.070519102,-0.047007209,-0.149002106,0.070607621,- 0.164828733,0.125608432,-0.503383161,0.196935744,1.1579818,0.450171515,0.371 029258,-0.239465566,0.139114328,-0.146848579,-0.99495574,0.093698765,0.29613 488,0.396070879,-0.349392801,-0.545938027,-0.533209136,-1.043240316,0.443916 702,0.489191444,-0.110253594,0.85520901,0.514300206,0.194122047,-0.015515904 ,0.371718851,0.123599861,-0.054056142,0.600694914,1.031212591,0.401957236,0. 919165667,0.179829211,-0.390039498,-0.18053262,0.225614896,-0.263266917,0.09 7865781,-0.339175,0.128268022,-0.218909316,-0.158820226,0.052968107,0.226688 929,0.399232198,-0.338868505,0.105548864,0.015069319,0.173134147,0.734360737 ,0.149209219,0.104314144,0.297442218,0.281752985,0.066614858,-0.237054708,-0 .416202756,-0.4553436,-0.247200396,0.149891357,-0.744168504,-0.113237471,-0. 348063305,-0.646365507,0.152462295,-0.076202092,0.570105303,1.355136299,0.67 0743311,0.555536409,0.044309874,0.346416414,0.080903178,-0.022058013,0.05146 113,0.571639519,-0.608394766,-1.042099523,-0.424376172,0.029839612,0.2011848 23,-0.350472004,-0.539286754, -0.534659958,-1.48331172,-0.561000511,-1.006901854,-0.187003327,0.14028527,- 1.143585738,-0.363034098,-0.380228595,0.21405645,-1.018873735,-0.819413325,1 .513057725,0.807157979,0.244027358,-1.401392635,-0.792479633,0.096393292,0.3 68649239,-0.851547562,-0.218067364,-1.089353102,-0.878458915,-1.219577899,0. 507890022,-0.733520001,0.484060368,0.365661506,1.121846295,0.303042945,-0.13 9121918,-0.098320368,-0.106622934,0.188563286,0.008190344,-0.410341158,-0.31 2989301,-1.061545878,-0.568563405,-1.155399739,-1.883778626,-0.190344408,0.0 60603439,1.161653127,0.213647904,0.892369836,-0.619564893,-1.431964447,-1.14 6542569,-1.389888666,-0.773646649,1.109692243,0.651181463,0.699242426,0.4172 10514,0.259875406,-0.3206103,-0.58319354,-0.797585191,0.061579065,0.38620253 4,-0.994857877,-1.59646646,-1.686350026,0.319503672,-1.210357706,-0.50867161 7,0.101941531,1.00457045,0.694510696,0.978711936,-0.370420986,0.171829118,-0 .117535387,-0.199221292,-1.112034065,-2.449556509,-0.301035034,-0.037693178, -0.861380687,-2.210778236,-2.09190049,0.820365362,0.490967208,2.112891162,0. 859439973,-0.304733118,-1.595880822,-0.534422882,-0.331804735,-1.506898067,0 .928493435,-0.009832358,0.851818926,1.114946835,-0.657769831,1.684108513,-0. 439939501,-0.74077976,-1.263256523,-1.467789736,0.18837065,0.227553895,-0.78 3773048,0.585929478,-1.556762305,0.531622697,0.259766356,0.933565227,-0.1780 94437,-0.825341532,-0.701405681,-0.060350032,0.090511396,0.531356205,1.16309 2788,1.481363235,0.46633719,-0.135790516,0.908218803,0.642349172,1.074103763 ,0.706583561,1.083110495,0.064986309,0.499908466,0.827667039,-0.018318373,0. 173889212,-0.439883407,-0.220669453,-1.399661354,-0.140114917,-0.149672619,0 .364435296,0.083912177,-1.096389186,0.018843038,-0.784985352,-0.92527011,-0. 730983345,-0.231951394,-0.038710927,0.386436574,-1.007858129,-0.537136894,-0 .766739233,1.167423863,-0.126872608,0.496810724,0.077707629,0.773773084,0.19 2511369,0.842674717,0.598946507,0.501112002,0.310749344,0.618616651,-0.72209 2508,-0.452788246,-0.416864649,0.256009242,0.444088225,0.188377187,0.0282259 96,1.503226103,0.784607339,0.641556595,0.209901971,0.844269914,0.58587756,-0 .062929835,-0.20704873,0.709382044,0.330519752,-0.55446395,0.937213983,0.540 470017,-0.088401703,-0.799151829,1.231575395,0.474392526,-0.853765693,0.0088 38998,-0.425155645,0.141919485,0.380413605,-0.851295378,-0.428380843,0.28579 1031,0.16039925,0.426477768,-0.622556262,-1.311081999,-0.416667269,-0.500478 773,-0.631436552,-1.591594623,-0.861832032,0.403737462,0.551325977,-1.445491 66,-0.654496052,-0.324087597,-2.670437322,-0.343794848,-0.246293411,-2.37559 8291,-2.474816969,0.763283747,-0.020552872,-3.312218363,-1.132490736,1.98932 1493,-0.337624292,-1.180048776,-1.660594553,-0.446794183,-2.018933377,1.7016 98522,1.143865623,1.312834954,-1.442902174,-1.519071886,-1.878563552,0.57063 0919,-0.559410723,1.591663631,-0.353982671,0.033239156,-0.834311531,-0.95404 5061,1.699694279,1.168443628,0.47012605,0.554741184,0.368125145,1.767492855, -0.084979822,0.392427719,-0.477479819,-1.954739503, -0.849585766,0.882118434,0.18415214,0.722491955,1.51421106,1.533321366,-0.65 8757022,-0.885148034,0.853670147,0.492018789,-0.156764415,-0.219883864,-0.39 911827,-0.48528421,0.643021243,-0.252498819,-0.825052279,-1.34732577,-0.4639 87884,0.30237604,0.075451364,-0.659426204,-1.167307205,-0.69372078,-0.465168 631,0.144212128,0.849941647,-0.176017647,0.811496941,0.664020292,0.821182356 ,-0.226232254,0.419739087,1.015128662,0.30783953,-0.318471607,-0.266169973,0 .965473523,0.431966114,-0.622200795,0.674752165,0.7328341,-0.073042211,0.551 711179,0.186683315,0.020721094,-1.208977727,-0.652221963,0.21085933,0.430876 606,0.094334685,-0.104821813,-0.083936633,-1.118863266,-1.432566727,-0.64816 042,-0.401803223,0.141358149,0.108601227,-0.27172453,-1.29260207,-0.86369709 8,0.288728685,-0.344349132,-0.089057114,-0.412923973,1.012526036,0.430820882 ,0.19821612,-0.651180492,-0.044301695,-2.014166781,-1.73302707,-1.05209364,0 .232180276,0.254777208,0.461574787,0.882577796,0.38724422,-0.569996983,-0.18 3087361,-0.747258751,-1.067295514,0.081608865,-1.775368481,-1.337493069,-0.0 24047132,0.419992218,-0.011975331,-1.89814077,-0.612146867,-0.172074772,-0.7 03663171,-0.210826642,1.111122543,-0.479971385,-1.591677472,0.387815586,0.13 7251254,1.091280671,0.847410659,0.536652099,1.197039866,0.611109684,0.250551 942,-0.693533809,-0.288392413,0.168329966,-0.578314865,-1.485115549,-0.35628 7618,-0.728803176,-1.009291478,-0.931058563,-0.291157874,0.594073668,-1.3704 05916,-1.856919894,-1.138830294,-0.753224572,-2.118041644,-0.203431277,-0.12 2257707,-3.612099594,-0.522784659,-1.096798973,0.186019943,2.443106756,0.417 653398,-0.487432202,1.496494604,0.137438174,0.807062721,0.027244245,1.379705 371,0.254913937,-1.512111976,-0.929349316,-1.857037845,-0.266125236,-1.85439 5272,-0.515686575,-0.330864144,0.388294153,1.226834218,1.617361745,-1.433212 674,0.042447825,-0.467920031,-2.21305195,-2.307738189,-0.551704971,-1.536869 183,-1.114600439,0.901802061,-1.224755052,-2.069951737,-1.617029527,0.589596 163,-0.95786173,0.144265473,-1.712481267,1.616327414,-2.417967698,-3.2215086 57,2.131489607,0.793917501,-2.215464978,-0.676020827,-0.016959213,-2.1428779 16,-1.836197432,1.018446029,0.835077922,0.69061048,2.380255182,0.63641146,-0 .401472603,-1.383988335,-0.136054443,-1.526223572,-1.304479755,-2.571291584, -0.341880675,0.341880675,-1.447465112,0.527130185,1.743237603,2.533116973,2. 673681115,2.158176613,2.144986045,-0.536020125,-0.342605768,1.379556083,1.13 792397,1.690804284,0.500157341,-1.160433282,-1.269858334,0.446287001,-0.6061 59138,1.854693729,1.882591578,0.5380076,-0.260956476,-0.678596065,2.0979255, 0.181680595,0.798497651,-0.195268556,0.315244575,-2.05943075,-0.583092031,0. 0922935,-0.973356468,-3.668395518,-1.885928112,-0.295858204,1.210747774,0.93 8822046,-0.630305117,1.560395209,0.239139213,0.238568702,-2.070533803,-1.899 205742,-0.148871089,-2.497406256,-1.212343031,-0.516708005,1.457127007,0.864 192449,0.151732306,-0.236127619,-0.626748897,0.559181327,0.757515616,-0.2854 98557,-1.337529463,-1.390222339,0.774730569,0.972786919, -0.238054865,0.068073522,1.753839221,1.460364486,0.869356742,0.829475358,0.2 42659668,-0.632144094,1.404721573,2.094317385,-0.551052327,-0.363780551,1.19 7053198,1.445583578,-0.619197025,0.06209252,-0.404292263,-0.093530794,0.5753 84528,1.370622827,2.507305488,1.304108583) col1= c(0,0.212503855,0.250737898,0.256797484,0.250305492,0.160437456,0.173729721, 0.123175522,0.09614655,0.249782138,0.257425607,0.28244612,0.277128364,0.3014 64221,0.264538497,0.316169592,0.279915472,0.32318193,0.209468695,0.321998132 ,0.167048465,0.304321622,0.350053239,0.295919452,0.340083767,0.223468037,0.3 56917935,0.191030138,0.079131164,0.116083359,0.064752145,0.177265257,0.35035 1888,0.330194001,0.408532531,0.32879264,0.39448722,0.334371817,0.413588784,0 .393507989,0.346181261,0.417956695,0.3529334,0.212361621,0.424608731,0.28338 0135,0.436066263,0.403421217,0.006691429,0.216343473,0.367386109,0.231427536 ,0.38155202,0.322086381,0.261913289,0.309234648,0.130551397,0.001955164,0.38 3522771,0.186214213,0.017758575,0.054178017,0.030746966,0.380450106,0.267562 367,0.38690199,0.398667468, 1.19E-05,0.032934345,0.235022879,0.377782426,0.342504993,0.068042338,0.00900 6484,0.283150746,0.304963542,0.278664107,0.295053726,0.267917743,0.28397787, 0.367908899,0.322654732,0.018352448,0.135376843,0.412759123,0.349517953,0.01 8740501,0.122786618,0.303874769,0.23466185,0.000221982,0.000722467,0.0305249 38,0.33627405,0.364826229,0.366036427,0.227678547,0.248940571,0.315207689,0. 28610075,0.387866096,0.355389049,0.245795125,0.370571392,0.349415299,0.39263 5746,0.066619545,0.077093116,0.387988774,0.019406508,0.036910174,0.374063789 ,0.226357521,0.136000304,0.312894459,0.083211245,0.006196966,0.002058155,0.0 00977855,0.100005811,0.375955254,0.014629121,0.022460739,0.233374422,0.28586 3942,0.086964093,0.225399844,0.170467672,0.159886956,0.373260944,0.020882207 ,0.225964619,0.167337184,0.19768543,0.264477672,0.359078091,0.364938593,0.37 3032543,0.335797667,0.373983659,0.082876927,0.160284844,0.38138136,0.3874278 64,0.168765172,0.328808252,0.299489672,0.253213033,0.193895861,0.15112638,0. 058685371,0.265594597,0.238862685,0.41859832,0.158882011,0.114909687,0.00251 4528,0.120797495,0.301210314,0.006444846,0.399527031,0.417728711,0.385256591 ,0.319258312,0.407701129,0.254095421,0.234087177,0.432177226,0.082170421,0.1 9789821,0.412213161,0.065565793,0.00086183,0.069791562,0.013975821,0.0007435 12,0.00092757,0.292803616,0.373001563,0.064246412,0.044916602,0.002335589,0. 378036919,0.304722977,0.005828006,0.228434279,0.090712776,0.323383304,0.2544 70629,0.276711946,0.378937635,0.036908617,0.019266861,0.357550777,0.10962933 6,0.382523913,0.243843496,0.175231395,0.069982826,0.006154193,0.004595993,0. 328680229,0.130741199,0.111053995,0.241850153,0.211214738,0.189429621,0.2034 65983,0.201530584,0.361916577,0.346261124,0.372199841,0.251557246,0.36712878 3,0.251522219,0.367180722,0.325118366,0.339932285,0.197772389,0.376795561,0. 382884844,0.257241908,0.152425072,0.390691114,0.297166467,0.382871419,0.1197 69412,0.085571153,0.397689497,0.208753022,0.221582711,0.389667031,0.29470876 7,0.342660126,0.400612561,0.062037631,0.407567454,0.374676749,0.262174992,0. 135578275,0.165247199,0.166761835,0.009404793,0.125918593,0.304782589,0.1266 9013,0.132691842,0.381183344,0.284255859,0.220683861,0.365878666,0.336843576 ,0.20434712,0.381537898,0.241574299,0.350704931,0.345421674,0.031897959,0.41 3999156,0.270611642,0.306194852,0.21545625,0.1697397,0.439990338,0.426954068 ,0.206608901,0.315439825,0.407780794,0.180653167,0.426334825,0.335085039,0.3 00514497,0.209642482,0.355588185,0.403478562,0.081377202,0.000961091,0.40506 3963,0.404041454,0.078265039,0.399971004,0.040990753,8.97E-05,0.328911042,0. 214002044,0.349044316,0.275874971,0.345959198,0.170158158,0.337988479,0.2119 06008,0.296303683,0.129430508,0.242510403,0.260204279,0.251464192,0.31302003 6,0.317282057,0.212356715,0.276288088,0.194258057,0.409072725,0.002443105,0. 35159162,0.338110897,0.010499928,0.017768208,0.307464899,0.309936472,0.40564 1746,0.0074005,0.034981102,0.07242701,0.340255938,0.36997755,0.025907951, 4.51E-5,0.259803783,0.223232008,0.171254485,0.284585237,0.305278799,0.142701 062,0.245090694,0.192093594,0.354995261,0.373275299,0.310975064,0.377770184, 0.336597179,0.121705383,0.045237927,0.343600229,0.3842013,0.052436965,0.0568 01196,0.028160545,0.304190306,0.093653614,0.182981315,0.394534874,0.32310385 2,0.180290687,0.246652848,0.04988185,0.421081721,0.426135859,0.02747738,0.03 0338844,0.007084011,0.386970248,0.133917294,0.11606115,0.28513906,0.38242107 7,0.308808826,0.357002461,0.076743294,0.021810764,0.068005863,0.067751449,0. 00233889,0.415972954,0.007393737,0.017975005,0.382830248,0.032654726,0.36139 7361,0.150353268,0.135022892,0.356776268,0.119436242,0.141400081,0.068454245 ,0.299209865,0.204487315,0.226231484,0.335711445,0.273880589,0.094603375,0.3 33947805,0.300136274,0.240736886,0.255161477,0.11754837,0.216159618,0.019686 147,0.376487328,0.138449932,0.217810628,0.093382487,0.004547791,0.151210866, 0.013140017,0.045860054,0.461576202,0.361546706,0.238235634,0.039110948,0.00 5622823,0.035485223,0.160216312,0.001975524,0.001698168,0.419748637,0.289020 798,0.299763394,0.187894111,0.105957788,0.290585832,0.078214849,0.355135607, 0.003711854,0.437558122,0.013883285,0.350657536,0.384464422,0.10392071,0.360 856612,0.132295731,0.000923896,0.32492525,0.409097421,0.376596974,0.05913986 4,0.02062899,0.024557833,0.002085163,0.35022789,0.332812175,0.194783879,0.17 1294141,0.322798259,0.385587781,0.250107856,0.381628989,0.346529476,0.205527 078,0.268372527,0.091082806,0.365818364,0.155811419,0.37476512,0.08310035,0. 158401987,0.405780128,0.102929467,0.326518882,0.059997966,0.351870934,0.0895 78406,0.111701153,0.275289844,0.451241566,0.388016893,0.027404087,0.32998644 5,0.219542533,0.41723301,0.11404275,0.378110727,0.329000281,0.450274079,0.45 8453959,0.279030125,0.035494448,0.006879419,0.006356498,0.04193038,0.3865906 43,0.000424914,0.129862193,0.027738114,0.003308235,0.378072707,0.156549115,0 .264257194,0.014742536,0.2603758,0.307530893,0.297176279,0.383261095,0.31693 0991,0.24059952,0.290104384,0.287199074,0.018362436,0.00275424,0.077106992,0 .284815928,0.393231983,0.075075164,0.030717506,0.032766025,0.000468228,0.075 703414,0.01938443,0.201907823,0.350721292,0.010470583,0.137499019,0.12419527 6,0.381642951,0.014160406,0.03679393,0.071467162,0.247641155,0.372533756,0.0 13782399,0.078633017,0.332155138,0.357358472,0.052669347,0.206414656,0.02548 8742,0.051765675,0.02421981,0.32685458,0.079792326,0.332353547,0.358778115,0 .157347,0.368170357,0.232366916,0.24242739,0.231729988,0.371689468,0.2786002 21,0.08615924,0.111553781,0.006071406,0.060725837,0.008570821,0.002004432,0. 235399031,0.319405417,0.15805937,0.367086713,0.223608481,0.099321338,0.01453 9267,0.041004072,0.028283964,0.104583105,0.19417067,0.297305158,0.28504975,0 .344848499,0.367300818,0.183038557,0.103899917,0.061078213,0.316682454,0.356 360716,0.027596404,0.007049582,0.012865384,0.357848216,0.044385881,0.1589788 48,0.334212782,0.207634571,0.283509249,0.212770343,0.178713445,0.355489007,0 .247786783,0.209298611,0.021917889,0.00121846,0.228541972,0.296095842,0.0895 35675,0.007449034,0.020894582, 0.279450684,0.331525412,0.085777471,0.27432822,0.248817464,0.049885023,0.197 850729,0.235050703,0.047983182,0.253836305,0.309328632,0.262183328,0.2084373 45,0.145605676,0.112673679,0.200483997,0.129234051,0.054738346,0.042069215,0 .352752618,0.359898708,0.108703995,0.311148845,0.024008975,0.322819005,0.362 044367,0.230742687,0.244756661,0.079599398,0.099054221,0.275686297,0.3281734 07,0.319761996,0.149579856,0.101385137,0.335414357,0.249273705,0.22200805,0. 292614469,0.180149083,0.276761146,0.181022053,0.320654017,0.328255787,0.2388 48538,0.284848217,0.359781001,0.11583447,0.183237933,0.00568392,0.232671197, 0.221344424,0.364212862,0.322429308,0.01166156,0.292383105,0.036914798,0.026 365876,0.053277205,0.183151345,0.261630542,0.362917963,0.013772953,0.0791001 67,0.038575391,0.12690866,0.23257446,0.329473097,0.320387044,0.234194538,0.3 73067285,0.207676308,0.295550216,0.328183079,0.380838711,0.285122603,0.03454 8592,0.086243422,0.091116702,0.394113508,0.349375665,0.381893378,0.278934437 ,0.028955711,0.235911752,0.284031558,0.377366811,0.212122726,0.301228405,0.2 49158305,0.173562318,0.250038167,0.379196173,0.054596798,0.166683335,0.31483 0819,0.227998302,0.0235392,0.099157239,0.336046479,0.035054479,0.289276037,0 .107280282,0.351762174,0.366533421,0.019995716,0.093830129,0.388254327,0.364 287273,0.356289662,0.032271525,0.00233148,0.111194555,0.082873087,0.05323153 1,0.002501336,0.054263415,0.351091724,0.314902164,0.010257493,0.09981527,0.1 80047852,0.000860213,0.225962758,0.243834519,0.00804413,0.016235214,0.294460 034,0.314459327,0.006792528,0.148134138,0.141157106,0.268772816,0.13261779,0 .079513303,0.242585489,0.049167563,0.167589608,0.243718997,0.217058298,0.096 149251,0.088360242,0.06077714,0.32278981,0.216125048,0.168631913,0.250214899 ,0.323059417,0.144679296,0.120330266,0.13332072,0.218635488,0.333520227,0.32 043225,0.347200131,0.110135939,0.291608176,0.34466037,0.182030083,0.01590125 8,0.123859423,0.25814924,0.351821259,0.284982503,0.132290539,0.138405185,0.1 54803174,0.108213872,0.260049173,0.329974341,0.259126003,0.233896196,0.17356 8581,0.14260591,0.292536642,0.200479563,0.057014994,0.016835001,0.146269114, 0.365776917,0.32187074,0.075705802,0.02027075,0.080421432,0.128117968,0.3486 3564,0.222884486,0.211567576,0.230950937,0.278497474,0.229191022,0.188428273 ,0.35108201,0.16395546,0.373978568,0.14332266,0.154287328,0.167632066,0.3487 59225,0.056682156,0.268319532,0.248569858,0.241946834,0.310408663,0.37439026 8,0.281766014,0.003569226,0.053952744,0.380409899,0.351111595,0.326926387,0. 206785955,0.21032321,0.003132033,0.00302612,0.075200726,0.131592374,0.349961 501,0.335371463,0.147360462,0.004725828,0.038185421,0.375655011,0.139676943, 0.236353376,0.099755434,0.151339968,0.349303119,0.376626851,0.044720131,0.24 9999553,0.000315491,0.010234119,0.066062601,0.363857776,0.365428783,0.336898 591,0.233059599,0.351914908,0.112276905,0.222388932,0.064026306,0.029002901, 0.325303792,0.004767376,0.031476522,0.296773998,0.344544559,0.294536057,0.00 6290423,0.139952082,0.250173513,0.104745051,0.228193423,0.17773878,0.1465866 94,0.01355857,0.349251219, 0.344489828,0.184284223,0.245197361,0.320274283,0.160477386,0.303301102,0.36 8248783,0.088082028,0.188159913,0.35676455,0.089108836,0.007836597,0.1711410 7,0.07533458,0.038688934,0.051030199,0.190381973,0.303704517,0.014698449,0.0 07997487,0.060584107,0.117000421,0.010529261,0.262537842,0.278073634,0.00135 8825,0.229978021,0.132274342,0.349222243,0.077327909,0.340448581,0.226280223 ,0.177718376,0.342249045,0.282056218,0.320456693,0.176297034,0.355947419,0.0 48597565,0.117660935,0.030882227,0.255311102,0.03329361,0.204158056,0.237087 577,0.34499483,0.191458677,0.132100212,0.058373917,0.3227317,0.196344244,0.0 14010601,0.021576668,0.211984445,0.069851194,0.117148211,0.268329021,0.09506 4666,0.032443097,0.069610929,0.318982077,0.137715024,0.343153219,0.046704145 ,0.155636537,0.023912422,0.014422042,0.14324337,0.301188747,0.062398131,0.22 5838113,0.321523241,0.056391047,0.082629253,0.268926392,0.290788614,0.307903 252,0.097832898,0.315738249,0.251105114,0.099387329,0.294250525,0.075385972, 0.098939921,0.022527851,0.262535494,0.346558225,0.084930828,0.326607857,0.14 3581482,0.075536649,0.080944328,0.136003967,0.142069089,0.247814941,0.271958 48,0.21710296,0.247988629,0.172376891,0.331027264,0.128733151,0.111759878,0. 336415869,0.198011773,0.128065267,0.132293339,0.325619351,0.26532889,0.17599 6374,0.095933754,0.348065068,0.282867749,0.270193284,0.351937037,0.020248933 ,0.187736048,0.332831436,0.102708561,0.002101241,0.0711197,0.275292471,0.234 214614,0.271549145,0.203765108,0.17309274,0.352829304,0.353728266,0.03084317 9,0.04864611,0.291569504,0.023128358,0.12577941,0.226935452,0.18961115,0.278 591151,0.343880286,0.265331602,0.176376176,0.319428553,0.282604331,0.2343621 05,0.049078969,0.049333543,0.278421998,0.237272388,0.243629724,0.324483036,0 .093372359,0.148932826,0.260389466,0.266572572,0.359680764,0.139687899,0.144 559313,0.069254342,0.185492279,0.22037257,0.194515834,0.154170189,0.16351567 1,0.325083589,0.195484736,0.276115019,0.312188429,0.137690429,0.035475327,0. 191613658) col2=c(0,0.787496145,0.749262102,0.743202516,0.749694508,0.839562544,0.82627 0279,0.876824478,0.90385345,0.750217862,0.742574393,0.71755388,0.722871636,0 .698535779,0.735461503,0.683830408,0.720084528,0.67681807,0.790531305,0.6780 01868,0.832951535,0.695678378,0.649946761,0.704080548,0.659916233,0.77653196 3,0.643082065,0.808969862,0.920868836,0.883916641,0.935247855,0.822734743,0. 649648112,0.669805999,0.591467469,0.67120736,0.60551278,0.665628183,0.586411 216,0.606492011,0.653818739,0.582043305,0.6470666,0.787638379,0.575391269,0. 716619865,0.563933737,0.596578783,0.993308571,0.783656527,0.632613891,0.7685 72464,0.61844798,0.677913619,0.738086711,0.690765352,0.869448603,0.998044836 ,0.616477229,0.813785787,0.982241425,0.945821983,0.969253034,0.619549894,0.7 32437633,0.61309801,0.601332532,0.999988084,0.967065655,0.764977121,0.622217 574,0.657495007,0.931957662,0.990993516,0.716849254,0.695036458,0.721335893, 0.704946274,0.732082257,0.71602213,0.632091101,0.677345268,0.981647552,0.864 623157,0.587240877,0.650482047,0.981259499,0.877213382,0.696125231,0.7653381 5,0.999778018,0.999277533,0.969475062,0.66372595,0.635173771,0.633963573,0.7 72321453,0.751059429,0.684792311,0.71389925,0.612133904,0.644610951,0.754204 875,0.629428608,0.650584701,0.607364254,0.933380455,0.922906884,0.612011226, 0.980593492,0.963089826,0.625936211,0.773642479,0.863999696,0.687105541,0.91 6788755,0.993803034,0.997941845,0.999022145,0.899994189,0.624044746,0.985370 879,0.977539261,0.766625578,0.714136058,0.913035907,0.774600156,0.829532328, 0.840113044,0.626739056,0.979117793,0.774035381,0.832662816,0.80231457,0.735 522328,0.640921909,0.635061407,0.626967457,0.664202333,0.626016341,0.9171230 73,0.839715156,0.61861864,0.612572136,0.831234828,0.671191748,0.700510328,0. 746786967,0.806104139,0.84887362,0.941314629,0.734405403,0.761137315,0.58140 168,0.841117989,0.885090313,0.997485472,0.879202505,0.698789686,0.993555154, 0.600472969,0.582271289,0.614743409,0.680741688,0.592298871,0.745904579,0.76 5912823,0.567822774,0.917829579,0.80210179,0.587786839,0.934434207,0.9991381 7,0.930208438,0.986024179,0.999256488,0.99907243,0.707196384,0.626998437,0.9 35753588,0.955083398,0.997664411,0.621963081,0.695277023,0.994171994,0.77156 5721,0.909287224,0.676616696,0.745529371,0.723288054,0.621062365,0.963091383 ,0.980733139,0.642449223,0.890370664,0.617476087,0.756156504,0.824768605,0.9 30017174,0.993845807,0.995404007,0.671319771,0.869258801,0.888946005,0.75814 9847,0.788785262,0.810570379,0.796534017,0.798469416,0.638083423,0.653738876 ,0.627800159,0.748442754,0.632871217,0.748477781,0.632819278,0.674881634,0.6 60067715,0.802227611,0.623204439,0.617115156,0.742758092,0.847574928,0.60930 8886,0.702833533,0.617128581,0.880230588,0.914428847,0.602310503,0.791246978 ,0.778417289,0.610332969,0.705291233,0.657339874,0.599387439,0.937962369,0.5 92432546,0.625323251,0.737825008,0.864421725,0.834752801,0.833238165,0.99059 5207,0.874081407,0.695217411,0.87330987,0.867308158,0.618816656,0.715744141, 0.779316139,0.634121334,0.663156424,0.79565288,0.618462102,0.758425701,0.649 295069,0.654578326,0.968102041, 0.586000844,0.729388358,0.693805148,0.78454375,0.8302603,0.560009662,0.57304 5932,0.793391099,0.684560175,0.592219206,0.819346833,0.573665175,0.664914961 ,0.699485503,0.790357518,0.644411815,0.596521438,0.918622798,0.999038909,0.5 94936037,0.595958546,0.921734961,0.600028996,0.959009247,0.999910341,0.67108 8958,0.785997956,0.650955684,0.724125029,0.654040802,0.829841842,0.662011521 ,0.788093992,0.703696317,0.870569492,0.757489597,0.739795721,0.748535808,0.6 86979964,0.682717943,0.787643285,0.723711912,0.805741943,0.590927275,0.99755 6895,0.64840838,0.661889103,0.989500072,0.982231792,0.692535101,0.690063528, 0.594358254,0.9925995,0.965018898,0.92757299,0.659744062,0.63002245,0.974092 049,0.999954887,0.740196217,0.776767992,0.828745515,0.715414763,0.694721201, 0.857298938,0.754909306,0.807906406,0.645004739,0.626724701,0.689024936,0.62 2229816,0.663402821,0.878294617,0.954762073,0.656399771,0.6157987,0.94756303 5,0.943198804,0.971839455,0.695809694,0.906346386,0.817018685,0.605465126,0. 676896148,0.819709313,0.753347152,0.95011815,0.578918279,0.573864141,0.97252 262,0.969661156,0.992915989,0.613029752,0.866082706,0.88393885,0.71486094,0. 617578923,0.691191174,0.642997539,0.923256706,0.978189236,0.931994137,0.9322 48551,0.99766111,0.584027046,0.992606263,0.982024995,0.617169752,0.967345274 ,0.638602639,0.849646732,0.864977108,0.643223732,0.880563758,0.858599919,0.9 31545755,0.700790135,0.795512685,0.773768516,0.664288555,0.726119411,0.90539 6625,0.666052195,0.699863726,0.759263114,0.744838523,0.88245163,0.783840382, 0.980313853,0.623512672,0.861550068,0.782189372,0.906617513,0.995452209,0.84 8789134,0.986859983,0.954139946,0.538423798,0.638453294,0.761764366,0.960889 052,0.994377177,0.964514777,0.839783688,0.998024476,0.998301832,0.580251363, 0.710979202,0.700236606,0.812105889,0.894042212,0.709414168,0.921785151,0.64 4864393,0.996288146,0.562441878,0.986116715,0.649342464,0.615535578,0.896079 29,0.639143388,0.867704269,0.999076104,0.67507475,0.590902579,0.623403026,0. 940860136,0.97937101,0.975442167,0.997914837,0.64977211,0.667187825,0.805216 121,0.828705859,0.677201741,0.614412219,0.749892144,0.618371011,0.653470524, 0.794472922,0.731627473,0.908917194,0.634181636,0.844188581,0.62523488,0.916 89965,0.841598013,0.594219872,0.897070533,0.673481118,0.940002034,0.64812906 6,0.910421594,0.888298847,0.724710156,0.548758434,0.611983107,0.972595913,0. 670013555,0.780457467,0.58276699,0.88595725,0.621889273,0.670999719,0.549725 921,0.541546041,0.720969875,0.964505552,0.993120581,0.993643502,0.95806962,0 .613409357,0.999575086,0.870137807,0.972261886,0.996691765,0.621927293,0.843 450885,0.735742806,0.985257464,0.7396242,0.692469107,0.702823721,0.616738905 ,0.683069009,0.75940048,0.709895616,0.712800926,0.981637564,0.99724576,0.922 893008,0.715184072,0.606768017,0.924924836,0.969282494,0.967233975,0.9995317 72,0.924296586,0.98061557,0.798092177,0.649278708,0.989529417,0.862500981,0. 875804724,0.618357049,0.985839594,0.96320607,0.928532838,0.752358845,0.62746 6244,0.986217601,0.921366983,0.667844862,0.642641528,0.947330653,0.793585344 ,0.974511258,0.948234325, 0.97578019,0.67314542,0.920207674,0.667646453,0.641221885,0.842653,0.6318296 43,0.767633084,0.75757261,0.768270012,0.628310532,0.721399779,0.91384076,0.8 88446219,0.993928594,0.939274163,0.991429179,0.997995568,0.764600969,0.68059 4583,0.84194063,0.632913287,0.776391519,0.900678662,0.985460733,0.958995928, 0.971716036,0.895416895,0.80582933,0.702694842,0.71495025,0.655151501,0.6326 99182,0.816961443,0.896100083,0.938921787,0.683317546,0.643639284,0.97240359 6,0.992950418,0.987134616,0.642151784,0.955614119,0.841021152,0.665787218,0. 792365429,0.716490751,0.787229657,0.821286555,0.644510993,0.752213217,0.7907 01389,0.978082111,0.99878154,0.771458028,0.703904158,0.910464325,0.992550966 ,0.979105418,0.720549316,0.668474588,0.914222529,0.72567178,0.751182536,0.95 0114977,0.802149271,0.764949297,0.952016818,0.746163695,0.690671368,0.737816 672,0.791562655,0.854394324,0.887326321,0.799516003,0.870765949,0.945261654, 0.957930785,0.647247382,0.640101292,0.891296005,0.688851155,0.975991025,0.67 7180995,0.637955633,0.769257313,0.755243339,0.920400602,0.900945779,0.724313 703,0.671826593,0.680238004,0.850420144,0.898614863,0.664585643,0.750726295, 0.77799195,0.707385531,0.819850917,0.723238854,0.818977947,0.679345983,0.671 744213,0.761151462,0.715151783,0.640218999,0.88416553,0.816762067,0.99431608 ,0.767328803,0.778655576,0.635787138,0.677570692,0.98833844,0.707616895,0.96 3085202,0.973634124,0.946722795,0.816848655,0.738369458,0.637082037,0.986227 047,0.920899833,0.961424609,0.87309134,0.76742554,0.670526903,0.679612956,0. 765805462,0.626932715,0.792323692,0.704449784,0.671816921,0.619161289,0.7148 77397,0.965451408,0.913756578,0.908883298,0.605886492,0.650624335,0.61810662 2,0.721065563,0.971044289,0.764088248,0.715968442,0.622633189,0.787877274,0. 698771595,0.750841695,0.826437682,0.749961833,0.620803827,0.945403202,0.8333 16665,0.685169181,0.772001698,0.9764608,0.900842761,0.663953521,0.964945521, 0.710723963,0.892719718,0.648237826,0.633466579,0.980004284,0.906169871,0.61 1745673,0.635712727,0.643710338,0.967728475,0.99766852,0.888805445,0.9171269 13,0.946768469,0.997498664,0.945736585,0.648908276,0.685097836,0.989742507,0 .90018473,0.819952148,0.999139787,0.774037242,0.756165481,0.99195587,0.98376 4786,0.705539966,0.685540673,0.993207472,0.851865862,0.858842894,0.731227184 ,0.86738221,0.920486697,0.757414511,0.950832437,0.832410392,0.756281003,0.78 2941702,0.903850749,0.911639758,0.93922286,0.67721019,0.783874952,0.83136808 7,0.749785101,0.676940583,0.855320704,0.879669734,0.86667928,0.781364512,0.6 66479773,0.67956775,0.652799869,0.889864061,0.708391824,0.65533963,0.8179699 17,0.984098742,0.876140577,0.74185076,0.648178741,0.715017497,0.867709461,0. 861594815,0.845196826,0.891786128,0.739950827,0.670025659,0.740873997,0.7661 03804,0.826431419,0.85739409,0.707463358,0.799520437,0.942985006,0.983164999 ,0.853730886,0.634223083,0.67812926,0.924294198,0.97972925,0.919578568,0.871 882032,0.65136436,0.777115514,0.788432424,0.769049063,0.721502526,0.77080897 8,0.811571727,0.64891799,0.83604454,0.626021432,0.85667734,0.845712672,0.832 367934,0.651240775, 0.943317844,0.731680468,0.751430142,0.758053166,0.689591337,0.625609732,0.71 8233986,0.996430774,0.946047256,0.619590101,0.648888405,0.673073613,0.793214 045,0.78967679,0.996867967,0.99697388,0.924799274,0.868407626,0.650038499,0. 664628537,0.852639538,0.995274172,0.961814579,0.624344989,0.860323057,0.7636 46624,0.900244566,0.848660032,0.650696881,0.623373149,0.955279869,0.75000044 7,0.999684509,0.989765881,0.933937399,0.636142224,0.634571217,0.663101409,0. 766940401,0.648085092,0.887723095,0.777611068,0.935973694,0.970997099,0.6746 96208,0.995232624,0.968523478,0.703226002,0.655455441,0.705463943,0.99370957 7,0.860047918,0.749826487,0.895254949,0.771806577,0.82226122,0.853413306,0.9 8644143,0.650748781,0.655510172,0.815715777,0.754802639,0.679725717,0.839522 614,0.696698898,0.631751217,0.911917972,0.811840087,0.64323545,0.910891164,0 .992163403,0.82885893,0.92466542,0.961311066,0.948969801,0.809618027,0.69629 5483,0.985301551,0.992002513,0.939415893,0.882999579,0.989470739,0.737462158 ,0.721926366,0.998641175,0.770021979,0.867725658,0.650777757,0.922672091,0.6 59551419,0.773719777,0.822281624,0.657750955,0.717943782,0.679543307,0.82370 2966,0.644052581,0.951402435,0.882339065,0.969117773,0.744688898,0.96670639, 0.795841944,0.762912423,0.65500517,0.808541323,0.867899788,0.941626083,0.677 2683,0.803655756,0.985989399,0.978423332,0.788015555,0.930148806,0.882851789 ,0.731670979,0.904935334,0.967556903,0.930389071,0.681017923,0.862284976,0.6 56846781,0.953295855,0.844363463,0.976087578,0.985577958,0.85675663,0.698811 253,0.937601869,0.774161887,0.678476759,0.943608953,0.917370747,0.731073608, 0.709211386,0.692096748,0.902167102,0.684261751,0.748894886,0.900612671,0.70 5749475,0.924614028,0.901060079,0.977472149,0.737464506,0.653441775,0.915069 172,0.673392143,0.856418518,0.924463351,0.919055672,0.863996033,0.857930911, 0.752185059,0.72804152,0.78289704,0.752011371,0.827623109,0.668972736,0.8712 66849,0.888240122,0.663584131,0.801988227,0.871934733,0.867706661,0.67438064 9,0.73467111,0.824003626,0.904066246,0.651934932,0.717132251,0.729806716,0.6 48062963,0.979751067,0.812263952,0.667168564,0.897291439,0.997898759,0.92888 03,0.724707529,0.765785386,0.728450855,0.796234892,0.82690726,0.647170696,0. 646271734,0.969156821,0.95135389,0.708430496,0.976871642,0.87422059,0.773064 548,0.81038885,0.721408849,0.656119714,0.734668398,0.823623824,0.680571447,0 .717395669,0.765637895,0.950921031,0.950666457,0.721578002,0.762727612,0.756 370276,0.675516964,0.906627641,0.851067174,0.739610534,0.733427428,0.6403192 36,0.860312101,0.855440687,0.930745658,0.814507721,0.77962743,0.805484166,0. 845829811,0.836484329,0.674916411,0.804515264,0.723884981,0.687811571,0.8623 09571,0.964524673,0.808386342) matrizEsperanca=cbind(col1,col2) parametros=c(0.208464913,6.95E-20,0.038663299,0.895792587,0.009160388,0.1062 25595,0.893674651,1.429250593,1.910965368,0.230183544) parametros[1]=pi1=mean(matrizEsperancas[,1]) #Verossimilhança ponderada logLikelihood=function(parametrosLike) { matrizDensidade=matrix(0,length(vetor),2) sigma21=rep(1,length(serie)) sigma22=rep(1,length(serie)) nu1=parametrosLike[7] nu2=parametrosLike[8] mu1=parametrosLike[9] for(i in 2:length(vetor)) { erro2=serie[i-1]^2 sigma21[i]=parametrosLike[1]+parametrosLike[2]*erro2+parametrosLike[3]*sigma 21[i-1] sigma22[i]=parametrosLike[4]+parametrosLike[5]*erro2+parametrosLike[6]*sigma 22[i-1] matrizDensidade[i,1]=log(dged(vetor[i],mu1,sqrt(sigma21[i]),nu1)) matrizDensidade[i,2]=log(dged(vetor[i],-(pi1/(1-pi1))*mu1,sqrt(sigma22[i]),n u2)) } somaLogLik=0 for(i in 1:length(vetor)) { for(j in 1:ncol(matrizDensidade)) { somaLogLik=somaLogLik+matrizEsperancas[i,j]*matrizDensidade[i,j] } } return(somaLogLik) } A = matrix(0,10,9) A[1,1]=A[2,2]=A[3,3]=A[4,4]=A[5,5]=A[6,6]=A[7,7]=A[8,8]=1 A[9,2]=A[9,3]=-1 A[10,5]=A[10,6]=-1 B = matrix(0,10,1) B[9]=B[10]=1 parametrosLike=parametros[2:length(parametros)] logLikelihood(parametrosLike) res = maxBFGS(logLikelihood, start=parametrosLike, constraints=list(ineqA=A, ineqB=B)) #Muito Obrigado pela ajuda. #Pedro
participantes (5)
-
Andre Oliveira
-
Benilton Carvalho
-
Bernardo Rangel Tura
-
Paulo Justiniano
-
Pedro Henrique Melo Albuquerque