HI peter
there is the function that i want to programm (joint in pdf folder).
my data is > dataexp v Ti1 26 5.792 26 1579.523 26 2323.704 28
68.855 28 426.076 28 110.297 28 108.298 28 1067.609 30 17.0510 30
22.6611 30 21.0212 30 175.8813 30 139.0714 30 144.1215 30 20.4616 30
43.4017 30 194.9018 30 47.3019 30 7.7420 32 0.4021 32 82.8522 32
9.8823 32 89.2924 32 215.1025 32 1.7526 32 0.7927 32 15.9328 32
3.9129 32 0.2730 32 0.6931 32 100.5832 32 27.8033 32 13.9534 32
53.2435 34 0.9636 34 4.1537 34 0.1938 34 0.7839 34 8.0140 34
31.7541 34 7.3542 34 6.5043 34 8.2744 34 33.9145 34 32.5246 34
3.1647 34 4.8548 34 2.7849 34 4.6750 34 1.3151 34 12.0652 34
36.7153 34 72.8954 36 1.9755 36 0.5956 36 2.5857 36 1.6958 36
2.7159 36 25.5060 36 0.3561 36 0.9962 36 3.9963 36 3.6764 36
2.0765 36 0.9666 36 5.3567 36 2.9068 36 13.7769 38 0.4770 38
0.7371 38 1.4072 38 0.7473 38 0.3974 38 !
1.1375 38 0.0976 38 2.38
NB> X<-Ti
thanks for helping mehafida
Date: Wed, 8 Aug 2012 03:13:28 -0700
From: ml-node+s789695n463956...@n4.nabble.com
To: hafida...@hotmail.fr
Subject: Re: help to program my function
Hi
Maybe it is time for you to read some basic stuff like R intro. It seems
to me that you expect R to behave like some other language you know but
probably your expectation is wrong.
See inline
>
> HI
>
> >i have a problem please help me to solve it:
> http://r.789695.n4.nabble.com/file/n4639434/aj.pdf aj.pdf
>
> >i want to calculate the vecteur a[j] where j: 1...8
>
> >this is the code in R:
>
> >aj.fun <- function(j, i, X, z, E, beta0, beta1){
> + n <- length(X)
> + iX <- order(X)
> + iz <- order(z)
> + e1 <- -(beta)*z[ iz[1:(i - 1)] ]
where do you get beta
> + numer <- E[j] - sum( X[ iX[1:(i - 1)] ] * exp(e1) )
> + e2 <- -(beta)*z[ iz[i:n] ]
> + denom <- sum( exp(e2) )
> + numer/denom
> + }
>
> > iX<-order(X)
> > iX
> [1] 75 37 29 60 73 20 69 55 30 70 72 38 26 35 65 61 74 50 71 57 25 54
64 76
> 56
> [26] 58 48 67 46 63 28 62 36 49 47 66 1 42 41 19 39 43 22 51 68 33 27 9
15
> 11
> [51] 10 59 32 40 45 44 52 16 18 34 4 53 21 23 31 7 6 13 14 12 17 24 5
8
> 2
> [76] 3
>
> > iZ<-order(Z)
> > iZ
> [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
23 24
> 25
> [26] 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47
48 49
> 50
> [51] 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72
73 74
> 75
> [76] 76
>
> >e1 <- -(beta)*Z[ iZ[1:(i - 1)] ]
> Warning message:
> In 1:(i - 1) : numerical expression has 76 elements: only the first used
As somebody already mentioned i is probably vector and in this case only
first value is taken. i seems to have the firs value 3.
> > e1
> [1] -442 -1664
>
> >> numer <- E[j] - sum( X[ iX[1:(i - 1)] ] * exp(e1))
> Warning message:
> In 1:(i - 1) : numerical expression has 76 elements: only the first used
> > numer
> [1] 9.5 9.5 9.5 9.5 9.5 9.5 9.5 9.5
Here j is vector of 8 values therefore 8 values
>
> >> e2 <- -(beta)*Z[ iZ[i:n] ]
> Warning message:
> In i:n : numerical expression has 76 elements: only the first used
> > e2
> [1] -442 -1664 -442 -1792 -476 -1792 -476 -1792 -510 -1920 -510
-1920
> [13] -510 -1920 -510 -1920 -510 -1920 -510 -2048 -544 -2048 -544
-2048
> [25] -544 -2048 -544 -2048 -544 -2048 -544 -2048 -544 -2048 -578
-2176
> [37] -578 -2176 -578 -2176 -578 -2176 -578 -2176 -578 -2176 -578
-2176
> [49] -578 -2176 -578 -2176 -578 -2304 -612 -2304 -612 -2304 -612
-2304
> [61] -612 -2304 -612 -2304 -612 -2304 -612 -2304 -646 -2432 -646
-2432
> [73] -646 -2432 -646 -2432
Strange, here first value of i seems to be 1 as n shall be 76 and final e2
length is 76.
> > denom <- sum( exp(e2) )
> > numer/denom
> [1] 4.313746e+192 4.313746e+192 4.313746e+192 4.313746e+192
4.313746e+192
> [6] 4.313746e+192 4.313746e+192 4.313746e+192
>
> >my problem that the vecteur a[j] could not have the same number!!!
I do not understand. Your numer is 9.5 repeted 8 times. If you divide it
by one number you will get nine times the same number.
You send us a code but no data so it is difficult to understand what is
your goal. It would be better to send input data
j, i, X, z, E, beta0, beta1
and assumed result in whole not in chunks scattered in several mails.
Regards
Petr
>
>
> thank you in advance
> hafida
>
>
>
> --
> View this message in context: http://r.789695.n4.nabble.com/help-to-
> program-my-function-tp4639434.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> [hidden email] mailing list
> 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.
______________________________________________
[hidden email] mailing list
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.
If you reply to this email, your message will be added to the
discussion below:
http://r.789695.n4.nabble.com/help-to-program-my-function-tp4639434p4639565.html
To unsubscribe from help to program my function, click here.
NAML
--
View this message in context:
http://r.789695.n4.nabble.com/help-to-program-my-function-tp4639434p4639605.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
