Sorry this is the code
library(mvtnorm) #Load mvtnorm package
N<-200;P<-9
W=matrix(NA, nrow=N, ncol=2)
xi=matrix(NA, nrow=N, ncol=2)
Eta=numeric(N)
phi<-matrix(data=c(1.0,0.3,0.3,1.0),ncol=2)
yo<-matrix(data=NA,nrow=N1,ncol=P); p<-numeric(P); v<- numeric(P)
for (t in 1:25) {
for (i in 1:N) {
xi1[i,]=rmvnorm(1, c(0,0), phi)
delta<-rnorm(1,0,sqrt(0.3))
eps<-rnorm(3,0,1)
Eta[i]=0.6*W[i,1]+0.6*W[i,1]*xi[i,1]+0.6*W[i,2]*xi[i,2]+0.6*W[i,2]*xi[i,2]*xi[i,2]+0.6*xi[i,1]+0.6*xi[i,2]+0.6*xi[i,1]*xi[i,1]+delta
v[1]<-1.0+xi[i,1]; v1[2]<-1.0+0.8*xi1[i,1]; v1[3]<- 1.0+0.8*xi1[i,1]
v[4]<-1.0+xi[i,2]; v1[5]<-1.0+0.8*xi1[i,2]; v1[6]<- 1.0+0.8*xi1[i,2]
v[7]<-1.0+Eta[i]+eps[1]; v[8]<-1.0+0.8*Eta[i]+eps[2];
v[9]<-1.0+0.8*Eta[i]+eps[3]
for (j in 1:9) { if (v[j]>=0) yo[i,j]<-1 else yo[i,j]<-0 } }
}#end
On 18 September 2015 at 13:54, PIKAL Petr <petr.pi...@precheza.cz> wrote:
> Hi
>
> You do not say anything what the code shall do and after getting third
> error of similar type, I gave up further investigation of your code.
>
> Error: object 'delta' not found
>
> Your error can have something to do with v[j] beeing NA
>
> Cheers
> Petr
>
> > -----Original Message-----
> > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of thanoon
> > younis
> > Sent: Friday, September 18, 2015 12:04 PM
> > To: r-help@r-project.org
> > Subject: [R] help
> >
> > Hi,
> >
> >
> > I need you help to correct the code below please
> > I have this error
> > Error in if (v[j] >= 0) yo[i, j] <- 1 else yo[i, j] <- 0 :
> > missing value where TRUE/FALSE needed
> >
> >
> >
> >
> >
> > library(mvtnorm) #Load mvtnorm package
> > N<-200;P<-9
> > W=matrix(NA, nrow=N, ncol=2)
> > xi=matrix(NA, nrow=N, ncol=2)
> > Eta=numeric(N)
> > phi<-matrix(data=c(1.0,0.3,0.3,1.0),ncol=2)
> > yo<-matrix(data=NA,nrow=N1,ncol=P); p<-numeric(P); v<- numeric(P)
> > for (t in 1:25) {
> > for (i in 1:N) {
> > xi1[i,]=rmvnorm(1, c(0,0), phi)
> > delta1<-rnorm(1,0,sqrt(0.3))
> > eps<-rnorm(3,0,1)
> >
> > Eta[i]=0.6*W[i,1]+0.6*W[i,1]*xi[i,1]+0.6*W[i,2]*xi[i,2]+0.6*W[i,2]*xi[i
> > ,2]*xi[i,2]+0.6*xi[i,1]+0.6*xi[i,2]+0.6*xi[i,1]*xi[i,1]+delta
> > v[1]<-1.0+xi[i,1]; v1[2]<-1.0+0.8*xi1[i,1]; v1[3]<-
> > 1.0+0.8*xi1[i,1]
> > v[4]<-1.0+xi[i,2]; v1[5]<-1.0+0.8*xi1[i,2]; v1[6]<-
> > 1.0+0.8*xi1[i,2]
> > v[7]<-1.0+Eta[i]+eps[1]; v[8]<-1.0+0.8*Eta[i]+eps[2];
> > v[9]<-1.0+0.8*Eta[i]+eps[3]
> >
> > for (j in 1:9) { if (v[j]>=0) yo[i,j]<-1 else yo[i,j]<-0 } }
> >
> > }#end
> >
> >
> >
> >
> > Any help would be appreciated.
> >
> >
> > Regards
> >
> >
> > --
> > Thanoon Y. Thanoon
> > PhD Candidate
> > Department of Mathematical Sciences
> > Faculty of Science
> > University Technology Malaysia, UTM
> > E.Mail: thanoon.youni...@gmail.com
> > E.Mail: dawn_praye...@yahoo.com
> > Facebook:Thanoon Younis AL-Shakerchy
> > Twitter: Thanoon Alshakerchy
> > H.P:00601127550205
> >
> > [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide http://www.R-project.org/posting-
> > guide.html
> > and provide commented, minimal, self-contained, reproducible code.
>
> ________________________________
> Tento e-mail a jakékoliv k němu připojené dokumenty jsou důvěrné a jsou
> určeny pouze jeho adresátům.
> Jestliže jste obdržel(a) tento e-mail omylem, informujte laskavě
> neprodleně jeho odesílatele. Obsah tohoto emailu i s přílohami a jeho kopie
> vymažte ze svého systému.
> Nejste-li zamýšleným adresátem tohoto emailu, nejste oprávněni tento email
> jakkoliv užívat, rozšiřovat, kopírovat či zveřejňovat.
> Odesílatel e-mailu neodpovídá za eventuální škodu způsobenou modifikacemi
> či zpožděním přenosu e-mailu.
>
> V případě, že je tento e-mail součástí obchodního jednání:
> - vyhrazuje si odesílatel právo ukončit kdykoliv jednání o uzavření
> smlouvy, a to z jakéhokoliv důvodu i bez uvedení důvodu.
> - a obsahuje-li nabídku, je adresát oprávněn nabídku bezodkladně přijmout;
> Odesílatel tohoto e-mailu (nabídky) vylučuje přijetí nabídky ze strany
> příjemce s dodatkem či odchylkou.
> - trvá odesílatel na tom, že příslušná smlouva je uzavřena teprve
> výslovným dosažením shody na všech jejích náležitostech.
> - odesílatel tohoto emailu informuje, že není oprávněn uzavírat za
> společnost žádné smlouvy s výjimkou případů, kdy k tomu byl písemně zmocněn
> nebo písemně pověřen a takové pověření nebo plná moc byly adresátovi tohoto
> emailu případně osobě, kterou adresát zastupuje, předloženy nebo jejich
> existence je adresátovi či osobě jím zastoupené známá.
>
> This e-mail and any documents attached to it may be confidential and are
> intended only for its intended recipients.
> If you received this e-mail by mistake, please immediately inform its
> sender. Delete the contents of this e-mail with all attachments and its
> copies from your system.
> If you are not the intended recipient of this e-mail, you are not
> authorized to use, disseminate, copy or disclose this e-mail in any manner.
> The sender of this e-mail shall not be liable for any possible damage
> caused by modifications of the e-mail or by delay with transfer of the
> email.
>
> In case that this e-mail forms part of business dealings:
> - the sender reserves the right to end negotiations about entering into a
> contract in any time, for any reason, and without stating any reasoning.
> - if the e-mail contains an offer, the recipient is entitled to
> immediately accept such offer; The sender of this e-mail (offer) excludes
> any acceptance of the offer on the part of the recipient containing any
> amendment or variation.
> - the sender insists on that the respective contract is concluded only
> upon an express mutual agreement on all its aspects.
> - the sender of this e-mail informs that he/she is not authorized to enter
> into any contracts on behalf of the company except for cases in which
> he/she is expressly authorized to do so in writing, and such authorization
> or power of attorney is submitted to the recipient or the person
> represented by the recipient, or the existence of such authorization is
> known to the recipient of the person represented by the recipient.
>
--
Thanoon Y. Thanoon
PhD Candidate
Department of Mathematical Sciences
Faculty of Science
University Technology Malaysia, UTM
E.Mail: thanoon.youni...@gmail.com
E.Mail: dawn_praye...@yahoo.com
Facebook:Thanoon Younis AL-Shakerchy
Twitter: Thanoon Alshakerchy
H.P:00601127550205
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
