[R] Sem error - subscript out of bounds
What´s happening with this following code: require(sem) Celpe.Mod.RAM - matrix(c( # path parametroInicio Produção- T1, gamma.11, NA, Produção- T2, gamma.12, NA, Produção- T3, gamma.13, NA, Produção- T4, gamma.14, NA, Compreensão - T1, gamma.21, NA, Compreensão - T2, gamma.22, NA, Compreensão - PI, gamma.23, NA, Certificado - Produção,alpha.11, NA, Certificado - Compreensão, alpha.12, NA, Compreensão - Produção,alpha.21, NA ), ncol=3, byrow=TRUE) Celpe.cov - cov(Celpe.Dados.Fatorial[3:7]) Celpe.cfa.sem - sem(ram=Celpe.Mod.RAM, S=Celpe.cov, N=1443,refit=TRUE,debug=TRUE) Here comes the Results: Celpe.cfa.sem - sem(ram=Celpe.Mod.RAM, S=Celpe.cov, N=1443,refit=TRUE,debug=TRUE) observed variables: [1] 1:T1 2:T2 3:T3 4:T4 5:PI latent variables: [1] 6:Produção7:Compreensão 8:Certificado parameters: [1] 1:gamma.11 2:gamma.12 3:gamma.13 4:gamma.14 5:gamma.21 6:gamma.22 7:gamma.23 8:alpha.11 9:alpha.12 [10] 10:alpha.21 RAM: heads to from parameter start [1,] 1 16 1NA [2,] 1 26 2NA [3,] 1 36 3NA [4,] 1 46 4NA [5,] 1 17 5NA [6,] 1 27 6NA [7,] 1 57 7NA [8,] 1 68 8NA [9,] 1 78 9NA [10,] 2 6710NA Error in startvalues(S, ram, debug = debug, tol = start.tol) : subscript out of bounds Celpe.cov T1T2T3T4PI T1 1.2459212 0.6072653 0.6410553 0.4030267 0.7001731 T2 0.6072653 0.9369251 0.5762526 0.3293337 0.6282211 T3 0.6410553 0.5762526 1.5346904 0.3888291 0.6064140 T4 0.4030267 0.3293337 0.3888291 0.5377340 0.3840199 PI 0.7001731 0.6282211 0.6064140 0.3840199 1.9750506 version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status Patched major1 minor9.0 year 2004 month04 day 28 language R []s Leonard Assis Estatístico / CONFE 7439 UIN : 41-764-523 / Skype : lmassis / msn : [EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] sporadic errors with nlrq() / optim()
In my computer, there's no error
WIN XP
PENTIUM4 1.67GHz
1GB RAM
R 1.9.0 Patched
[]s
Leonard Assis
Estatístico / CONFE 7439
UIN : 41-764-523 / Skype : lmassis / msn : [EMAIL PROTECTED]
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Peter Dalgaard
Sent: Thursday, May 06, 2004 6:21 PM
To: Rogers, James A [PGRD Groton]
Cc: '[EMAIL PROTECTED]'
Subject: Re: [R] sporadic errors with nlrq() / optim()
Rogers, James A [PGRD Groton] [EMAIL PROTECTED] writes:
Dear List,
Apologies if this is a known problem ... I wasn't able to find it on
the bug list, but it is a problem that does not seem to occur with a
MAC build of R 2.0, so perhaps this problem has already been addressed for
the future.
I am getting *sporadic* errors when refitting the same model to the
same data set, using nlrq() in the nlrq package. The algorithm is not
stochastic, so I would expect to get errors either every time, or never.
###
library(stats) # or library(nls) if using R 1.9.0
library(nlrq)
test - data.frame(x = c(7.60090245954208, 6.90775527898214,
6.21460809842219, 5.52146091786225, 4.60517018598809,
3.91202300542815,
3.2188758248682 , 2.52572864430826, 1.83258146374831,
7.60090245954208, 6.90775527898214, 6.21460809842219,
5.52146091786225, 4.60517018598809, 3.91202300542815, 3.2188758248682
, 2.52572864430826, 1.83258146374831, 6.21460809842219, 6.21460809842219),
y = c( 11.0161506644269, 9.84267541313937,
8.66146668057266, 7.48099216286952, 6.50578406012823,
6.24027584517077, 5.63121178182137, 5.71702770140622,
5.64190707093811, 10.8983676287705, 9.91857318995417,
8.74608021735751, 7.58120982619635, 6.361302477573 , 5.91889385427315,
5.63835466933375 , 5.80211837537706, 5.64897423816121,
8.6195692580331 , 8.70367275835886)
)
i - 1
while(i 500) {# I usually hit an error within 50 iterations
cat(i, \n)
nlrq(y ~ SSfpl(x, A, B, xmid, scal), data = test)
i - i + 1
}
###
Errors occur with version 1.8.1 and 1.9.0 on both Windows (two
different
machines) and UNIX, but not on version 2.0 on a MAC (these are the
only R version - OS permutations I was able to get reports on easily).
Anyone understand what is happening here?
There's no version 2.0, only 2.0.0 Under development (unstable).
If you look at the NEWS file for that version (or for 1.9.0-patched) under
1.9.1 changes, you'll find
o The L-BFGS-B option of optim() apparently needs part of its
workspace zeroed. (PR#6720)
which gave problems of the sort you experiences.
--
O__ Peter Dalgaard Blegdamsvej 3
c/ /'_ --- Dept. of Biostatistics 2200 Cph. N
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] sporadic errors with nlrq() / optim()
Only theese: :(
495
496
497
498
499
There were 50 or more warnings (use warnings() to see the first 50)
warnings()
Warning messages:
1: multi-argument returns are deprecated in: return(coef, w)
2: multi-argument returns are deprecated in: return(coef, w)
3: multi-argument returns are deprecated in: return(coef, w)
4: multi-argument returns are deprecated in: return(coef, w)
5: multi-argument returns are deprecated in: return(coef, w)
6: multi-argument returns are deprecated in: return(coef, w)
7: multi-argument returns are deprecated in: return(coef, w)
8: multi-argument returns are deprecated in: return(coef, w)
9: multi-argument returns are deprecated in: return(coef, w)
10: multi-argument returns are deprecated in: return(coef, w)
...
50: multi-argument returns are deprecated in: return(coef, w)
[]s
Leonard Assis
Estatístico / CONFE 7439
UIN : 41-764-523 / Skype : lmassis / msn : [EMAIL PROTECTED]
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Leonard assis
Sent: Thursday, May 06, 2004 6:42 PM
To: 'Peter Dalgaard'; 'Rogers, James A [PGRD Groton]'
Cc: [EMAIL PROTECTED]
Subject: RE: [R] sporadic errors with nlrq() / optim()
In my computer, there's no error
WIN XP
PENTIUM4 1.67GHz
1GB RAM
R 1.9.0 Patched
[]s
Leonard Assis
Estatístico / CONFE 7439
UIN : 41-764-523 / Skype : lmassis / msn : [EMAIL PROTECTED]
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Peter Dalgaard
Sent: Thursday, May 06, 2004 6:21 PM
To: Rogers, James A [PGRD Groton]
Cc: '[EMAIL PROTECTED]'
Subject: Re: [R] sporadic errors with nlrq() / optim()
Rogers, James A [PGRD Groton] [EMAIL PROTECTED] writes:
Dear List,
Apologies if this is a known problem ... I wasn't able to find it on
the bug list, but it is a problem that does not seem to occur with a
MAC build of R 2.0, so perhaps this problem has already been addressed
for
the future.
I am getting *sporadic* errors when refitting the same model to the
same data set, using nlrq() in the nlrq package. The algorithm is not
stochastic, so I would expect to get errors either every time, or never.
###
library(stats) # or library(nls) if using R 1.9.0
library(nlrq)
test - data.frame(x = c(7.60090245954208, 6.90775527898214,
6.21460809842219, 5.52146091786225, 4.60517018598809,
3.91202300542815,
3.2188758248682 , 2.52572864430826, 1.83258146374831,
7.60090245954208, 6.90775527898214, 6.21460809842219,
5.52146091786225, 4.60517018598809, 3.91202300542815, 3.2188758248682
, 2.52572864430826, 1.83258146374831, 6.21460809842219, 6.21460809842219),
y = c( 11.0161506644269, 9.84267541313937,
8.66146668057266, 7.48099216286952, 6.50578406012823,
6.24027584517077, 5.63121178182137, 5.71702770140622,
5.64190707093811, 10.8983676287705, 9.91857318995417,
8.74608021735751, 7.58120982619635, 6.361302477573 , 5.91889385427315,
5.63835466933375 , 5.80211837537706, 5.64897423816121,
8.6195692580331 , 8.70367275835886)
)
i - 1
while(i 500) {# I usually hit an error within 50 iterations
cat(i, \n)
nlrq(y ~ SSfpl(x, A, B, xmid, scal), data = test)
i - i + 1
}
###
Errors occur with version 1.8.1 and 1.9.0 on both Windows (two
different
machines) and UNIX, but not on version 2.0 on a MAC (these are the
only R version - OS permutations I was able to get reports on easily).
Anyone understand what is happening here?
There's no version 2.0, only 2.0.0 Under development (unstable).
If you look at the NEWS file for that version (or for 1.9.0-patched) under
1.9.1 changes, you'll find
o The L-BFGS-B option of optim() apparently needs part of its
workspace zeroed. (PR#6720)
which gave problems of the sort you experiences.
--
O__ Peter Dalgaard Blegdamsvej 3
c/ /'_ --- Dept. of Biostatistics 2200 Cph. N
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Lme4 error
What's happening? Celpe.Mixed - lme(Nota~Tarefa + Corretor %in% Tarefa,random=~1|ID,data=Celpe.Dados11) Error in EMsteps-(`*tmp*`, value = controlvals) : error code 6 from Lapack routine dtrtri I've trying to fit a Mixed model using lme from package lme4 and the error above has occourred. What Should I Do to solve this? Best regards []s Leonard Assis Estatístico / CONFE 7439 UIN : 41-764-523 / Skype : lmassis / msn : [EMAIL PROTECTED] [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] dlls
Is there any way to compile R 1.9 Sources in VC++ ? []s Leonard Assis Estatístico / CONFE 7439 UIN : 41-764-523 / Skype : lmassis / msn : [EMAIL PROTECTED] -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley Sent: Friday, April 30, 2004 6:33 AM To: Ben.Flood Cc: '[EMAIL PROTECTED]' Subject: Re: [R] dlls On Fri, 30 Apr 2004, Ben.Flood wrote: I am trying to run an algorithm in a large dataset and want to speed this process up by using C/C++ functions. I have been attempting to create a dll library using Microsoft Visual C++ in Windows and call that using dyn.load(file.dll) but when I use is.loaded(function) I invariably get FALSE. If anybody has done this and could give me more detailed instructions it would be very much appreciated, It's a `Microsoft Visual C++ in Windows' issue, not an R issue. I expect you are not exporting any entry points (the default under that `system', amazingly): readme.packages tells you how to use VC++ in detail, and there are worked examples in the on-line complements to `S Programming'. Do you have a very good reason to make things difficult for yourself by not using the recommended compilers? Just a thought -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Formulae Construction
Where Can I Have More information about Formulae in R (Nesti8ng, Crossing Factors, Etc). I´ve tried to obtain this information in R Help and that source were a kind of inconclusive for my Doubts. []s Leonard Assis Estatístico CONFE 7439 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] How to purchase R on CD
Why Dont U Download and Burn into a CD at your Home? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Barry Rowlingson Sent: Monday, March 29, 2004 3:08 PM To: [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Subject: Re: [R] How to purchase R on CD [EMAIL PROTECTED] wrote: My organization has a policy against installing software downloaded from the Internet. There is a waiver procedure, but it is difficult and lengthy. It is much easier to get approval to purchase software, even if it is a CD copy of open-source. However, a search on r-project.org didn't seem to show any way of obtaining the distribution in this way. If you don't want a Linux version, and you dont want a Mac version, but you do want a Windows version, GNU-Win sell CDs with masses of goodies on, including R (version 1.8.0 at the moment though) http://gnuwin.epfl.ch/apps/en/bestlist.html Just dont waste all your time playing Tux Racer. A whole $3.95 from osdepot.com. Baz __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Ordered Logistic Mixed Model
What is the best solution to fit an Ordered Logistic (Or Probit) Mixed Model in R? []s Leonard Assis Estatístico - CONFE 7439 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] Ordered logit/probit
Anyone Knows how Shoud I Fit a Mixed Effects Ordered Logistic Model in R? (lme?) []s Leonard Assis Estatístico - CONFE 7439 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Valentin Stanescu Sent: quarta-feira, 24 de março de 2004 15:15 To: [EMAIL PROTECTED] Subject: [R] Ordered logit/probit Hello everyone I am trying to fit an ordered probit/logit model for bank rating prediction. Besides polr() in MASS package which is not written especially for this as far as I know, do you know how else I can do this? I already found the modified polr () version on the Valentin STANESCU Enrst and Young Tel. 402 4000 -- The information contained in this communication is intended solely for the use of the individual or entity to whom it is addressed and others authorized to receive it. It may contain confidential or legally privileged information. If you are not the intended recipient you are hereby notified that any disclosure, copying, distribution or taking any action in reliance on the contents of this information is strictly prohibited and may be unlawful. If you have received this communication in error, please notify us immediately by responding to this email and then delete it from your system. Ernst Young is neither liable for the proper and complete transmission of the information contained in this communication nor for any delay in its receipt. Note: If you have received a delivery failure report, it may be due to the change in the Ernst Young e-mail domain from eyi.com to ey.com. Could you please make the necessary amendment, if required, and resend the message. [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] substitute question
Anyone Knows How Shoud I Build Rasch Models (Partial Credit Models) In R? []s Leonard Assis Estatístico - CONFE 7439 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
