11:05 a.m.
Hi gretl users,
I am looking for a simple way to determine the column for which scalar B
is, let's say, larger or equal to the entry A[i] and store this
information in scalar C. If B is never greater than any entry in A[i], C
is zero.
Reading out each entry via a loop is possible, but I am wondering
whether there is a short-cut to this.
C = 0
A = seq(4,1)
B = 3
# check at which entry A[i] the value of B >= A[i]
Maybe somebody a nice suggestion :-)
Artur
11:45 a.m.
This is my experiment:
CC=(A .<= B)?1:0
CC
CC (1 x 4)
0 1 1 1
But that wont give you the indexes.
Good luck
On Thu, Nov 21, 2013 at 11:05 AM, Artur T. <artur.tarassow(a)googlemail.com>wrote:
Hi gretl users,
I am looking for a simple way to determine the column for which scalar B
is, let's say, larger or equal to the entry A[i] and store this
information in scalar C. If B is never greater than any entry in A[i], C
is zero.
Reading out each entry via a loop is possible, but I am wondering
whether there is a short-cut to this.
C = 0
A = seq(4,1)
B = 3
# check at which entry A[i] the value of B >= A[i]
Maybe somebody a nice suggestion :-)
Artur
_______________________________________________
Gretl-users mailing list
Gretl-users(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-users
11:58 a.m.
Hi there,
this is an ordinary if then or else loop in Fortran and Basic, but instead
of storing the found result in "C"
I'suggest creating a (4x1) matrix columm vector and store the "found values
in the column vector at each "i" location
and I suspect it is the same in Gretl but I confess to being a novice at
Gretl sorry to say right now
Best wishes to all,
Hope these few commnts are of small help
Good luck
Richard Hudson
Dr RJF Hudson
----- Original Message -----
From: "Artur T." <artur.tarassow(a)googlemail.com>
To: "Gretl list" <gretl-users(a)lists.wfu.edu>
Sent: Thursday, November 21, 2013 9:05 PM
Subject: [Gretl-users] small practical issue
Hi gretl users,
I am looking for a simple way to determine the column for which scalar B
is, let's say, larger or equal to the entry A[i] and store this
information in scalar C. If B is never greater than any entry in A[i], C
is zero.
Reading out each entry via a loop is possible, but I am wondering
whether there is a short-cut to this.
C = 0
A = seq(4,1)
B = 3
# check at which entry A[i] the value of B >= A[i]
Maybe somebody a nice suggestion :-)
Artur
_______________________________________________
Gretl-users mailing list
Gretl-users(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-users
12:20 p.m.
On Thu, 21 Nov 2013, Artur T. wrote:
I am looking for a simple way to determine the column for which
scalar B
is, let's say, larger or equal to the entry A[i] and store this
information in scalar C. If B is never greater than any entry in A[i], C
is zero.
You say "the" column: does that mean A is a row vector and at most one
element of A is <= B? If so:
scalar i = iminr(A)
scalar C = A[i] <= B ? i : 0
If more than one element of A can satisfy the condition and you want the
index of the first such element:
scalar i = imaxr(A .<= B)
scalar C = A[i] <= B ? i : 0
Allin Cottrell
12:20 p.m.
Thanks Helio,
I actually just forgot the "." for the element-wise iteration over A.
The solution is given by summing up the resulting C which gives the
position A[i] for which B>=A[i]. If the sum(C)+1 is greater cols(A), B
is never greater than any A[i]:
<hansl>
A = seq(4,1)
matrix B = 3
C = (B.>=A)?0:1
matrix pos = 0
pos = sum(C)+1
pos
pos<\hansl>
@ Richard: This is my first idea, but I was looking for a simpler way
avoiding any use of loops here.
Artur
Am 21.11.2013 12:58, schrieb Dr RJF Hudson:
Hi there,
this is an ordinary if then or else loop in Fortran and Basic, but instead
of storing the found result in "C"
I'suggest creating a (4x1) matrix columm vector and store the "found values
in the column vector at each "i" location
and I suspect it is the same in Gretl but I confess to being a novice at
Gretl sorry to say right now
Best wishes to all,
Hope these few commnts are of small help
Good luck
Richard Hudson
Dr RJF Hudson
----- Original Message -----
From: "Artur T." <artur.tarassow(a)googlemail.com>
To: "Gretl list" <gretl-users(a)lists.wfu.edu>
Sent: Thursday, November 21, 2013 9:05 PM
Subject: [Gretl-users] small practical issue
> Hi gretl users,
>
> I am looking for a simple way to determine the column for which scalar B
> is, let's say, larger or equal to the entry A[i] and store this
> information in scalar C. If B is never greater than any entry in A[i], C
> is zero.
> Reading out each entry via a loop is possible, but I am wondering
> whether there is a short-cut to this.
>
> C = 0
> A = seq(4,1)
> B = 3
> # check at which entry A[i] the value of B >= A[i]
>
> Maybe somebody a nice suggestion :-)
>
> Artur
>
> _______________________________________________
> Gretl-users mailing list
> Gretl-users(a)lists.wfu.edu
> http://lists.wfu.edu/mailman/listinfo/gretl-users
_______________________________________________
Gretl-users mailing list
Gretl-users(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-users
--
For global instant messaging use the free and transparent Jabber chat:
t.artur(a)cryptolab.net
Avoid Facebook and join me on Diaspora: https://despora.de/i/9fba3608d45a
If you want to use encrypted emailing, here is my public PGP key:
-----BEGIN PGP PUBLIC KEY BLOCK-----
Version: GnuPG v1.4.12 (GNU/Linux)
mQINBFIFGVsBEACc0QWXQ8fijGzIvQUWQTTYE/0O/qKQRRjNX7O3UwkASTOZM12H
aZ3YzpqZOvfFnK/BMUCtOYfkoZJJwoLwPDio+wlmQ1DCDjNVYJhkxbjThs+pBdk0
YdyZ3b+bwylb8StnexA0EOSpdV8tOpWo3JVw2WXAs9y9MxierKHIr3aAsHNyGI4V
j+46Cz9mmfyQBJAA7tlK7Th1pibb18ww0fxFnAOmBNBvexCTb+z353FtgDdIc/Wq
LFLmdtP73Ntui6DYviGT+YvV2FBYFxF6x8TGwnYnnvQ88VGl+aflyVjp+LRUnbo0
pOD5LTAA6f6b5ZSZq8GzufLDrZxMRL1qS/+8lp8AAI0/x7y7GtRdznB4KYOAKd1Z
a/GgcuMl9N6T7B+azwxQa1th+ucMdSZO4PMOHjhsp0xYJftgEOG44GHJPqchS5Zd
nKphboRVJjV9KA/uLUV/S6n8Qhh6tRBLiQ2RPMSGjKTdKObMxJ9YbDirHChSgi5x
77iGf+3F413s6XCbazLYv+03lylFCBT3p9mYDWIMqJS2jCzWVjXbea6iu5ZMN5Ji
ESNiTHjN1pW97c+HaLxzo9/PobUceHYVzwKvEfBH4wVIzCExMlKRI0xwIp6ACk3E
XJPKsp8+a9HMp1qdCvi0zye/Th7FhIOPhVYlRgCf157O18MxLj4dJu2XewARAQAB
tCBBcnR1ciBULiA8dC5hcnR1ckBjcnlwdG9sYWIubmV0PokCPgQTAQIAKAUCUgUZ
WwIbIwUJA8JnAAYLCQgHAwIGFQgCCQoLBBYCAwECHgECF4AACgkQlEChRi+nhULG
FQ/9Fcz/aI7fcKOd9NG1Wb0ziyKHihg5uoEMFsOTFGDboC+EddrQwkGdkOpWhQzN
vnPCji+WuWHliG6DiRYViG+HYN+yBhaFB6ASe+jJYCXaoxI1u7e3UOscwOmg6zb7
TdR27jI0mguu0ZT1bFLWh7hDzqAa33ANWwX++J8ncc08SMw3/cJoddzUrkc1PcUY
M+mB2AgURVCC+TCzpsEuZ7jqTJAH3+aLUnxGMWm3hTlt2eirYieon+a9PlRqlTZD
0wOyd6z5mU+uUQr4lqfHNl9XOH/vP9oqDNW6uCItdGTmKtOVioE0BA+7p06tnobv
iUxNCc8dz81x+ZLoxjS2H6zCTf00PVf8ZdBWtxM0S0/9sRkuyTE1exwZOpob552m
HbQIIRIgB3+PwyM2mXgRSWXciITyH6LOvJoPDTp+vrHqUtUCX3LoIRP6ySqo7rEi
/UnF2WAwjjEPIVByvfS9ZVXQaL5ASHUV06yMAmMp7h6GfLgo4JNOr6wxl+sYgNwy
UimKVGx2OdAQoUPGN2VQA9SPxFtbDPTsH472s4FoEVGhq2OBSMdAVeKJgDDkrsQm
u5lk0+Or2GewwjxoBu77kRgihkqBX2ND6Z2pHs95g1nngrqM6sg895fp5Lo+MHAJ
+8cSUxs0diDA7ybG1j8KTx5Vx6yPLcmgj9MlHVaNiMnAKuC5Ag0EUgUZWwEQANDG
ZbQoate2ul8zMDwlr+DmlzSJBUMue4ABFXW0dJSZ+E8awJkYRL3jcjJqVVOrk9cT
gqy0Rum16SGswYrFa4xEE/gb8h97A/ZuJ9wOSe/Os5hySI+62cZSJrd+2AYI8rY7
3iW/QlaO38mFcsulSMg7+SqZjHpFI3tfMj+5e3vx1FsYazAR/jI4U02gP6OdHj5h
6oAS7OqlCSvN0j0LatTuc7LvfjLkSWCWyZUfe2pujUCc1gZ95wINzSamYcxfd+o4
MjOizFULw5zmuioiHWlp35nGLu03uOaBJgR85dTvyQqa4aQXLCfwvSICBEl6aAZp
ODtfZUOpT/sQdqIpGrc4rytz763BpnuVl92RzNoQtjsLE8wnJRzElFjYeOI4JZLL
RhEa0iNVpbFvv4rmOqVafxL0dR8wkCjiPMgWii/qm6crpsOnhOJ3Wov8EGVL53hI
0OOGxkQ0YV9PF8L45vcvK8TMSm+V19l53yQYOvxPECcCksn9uFCdMNvPoef7FZLO
nXIORP4IBcMTKf2N07QFbfkZRVh5SZgUOFxMtLIqBmWQqpmEAdO7SCybJkAgwdpJ
VliEWQSxFeLwyV9kQUmN7oQsQ8jYxO1C0vTnstDjZYC4EBPf08kCPLQXI1gM3/v2
POKDZVk6EAFQelVia8OIFyHA9mLwc9HV7b7KYwGrABEBAAGJAiUEGAECAA8FAlIF
GVsCGwwFCQPCZwAACgkQlEChRi+nhUKh5A/7B75inrTA4PMWaY6fjA4UnhFi6BBy
bY/QIpkv7zle4H6y0tbDsAziB80i0GsMqDka+LEmCr9BBmZNyD6vUeFbE0uECiE8
XupOmw/B3+Czscuo29CqtMYpwlMX9oVXFiaWd/LHsIMOgI5BstDaE8pBga+IMyKz
HsPdSDz+MqdQwp/9pY5FfB1/h4VpFNXOIsEDjimALvwCHW9F34Nv/tFCd39Li2sd
jktLLgJKVurEHisEfws4vJRsmPtsFdgqGGHm0r/wxJYROY8ki+veZCNpYmwVWAcY
iZK7AFbqqsNMdMyKXOiMGawq9LZSnpBmgiBzKJmqAL6BNqEO4/M8E+zz8rlrjD3/
16iTljR+FPriXF/g0KWli1C1Fx6ntCYQUElO9Hu6X4CEA5zT7T6vXTapRkMA51PM
xepwgASu9YWRAig8qZI18X2uh0I3VdhERxgAMmDn8V1qjTn/gd0S7RPYtduUHgNW
/sCV5B6vZsmaO0fImKCJdJjVZTCViW5h9XHFSNH7hXjnr/+4SKrvfm4Du1DPiGKn
4oD7/QUHeSpdfDKJo/E9nM39aKxgN6+RqJUd2J2Db2aNBCQwa1o+SOtrriqyCMjp
K+4gDumAFG2lg/PD/ki/aNucMROO9t/uPFQFhheaUF9R1rXJh66KyxQFsxKC00uZ
kJ87csTGWSWmWns=
=tQPP
-----END PGP PUBLIC KEY BLOCK-----
12:25 p.m.
Thanks Allin,
this was exactly what I was looking for! Very nice and neat solution :-)
Artur
Am 21.11.2013 13:20, schrieb Allin Cottrell:
On Thu, 21 Nov 2013, Artur T. wrote:
> I am looking for a simple way to determine the column for which scalar B
> is, let's say, larger or equal to the entry A[i] and store this
> information in scalar C. If B is never greater than any entry in A[i], C
> is zero.
You say "the" column: does that mean A is a row vector and at most one
element of A is <= B? If so:
scalar i = iminr(A)
scalar C = A[i] <= B ? i : 0
If more than one element of A can satisfy the condition and you want the
index of the first such element:
scalar i = imaxr(A .<= B)
scalar C = A[i] <= B ? i : 0
Allin Cottrell
1 p.m.
You are welcome,
but I'd be pretty sure a true grtl wizard will pop by shortly o hlp you out
I atill think the be way o hnl our problem i by using a volumn vector ond
let the program do the "bubble Sort" looping work
fingers crossed for you
Seems to me rom ht ou write that you r mkin things oo hard for yourselfl
WHy Sum things ?
Set up 2 vectors of size to hold all the numbers say A[4] nd B[4]
iNitialize B[i]=0
THen loop for i =1 to 4
If a[i].> b[i] then b[i]= a[i]
b[i] will hold your result all sorted for you
if you get me ?
Richard Hudson
Kind Regards
Thanks Hello,
I actually just forgot the "." for the element-wise iteration over A.
The solution is given by summing up the resulting C which gives the
position A[i] for which B>=A[i]. If the sum(C)+1 is greater cols(A), B
is never greater than any A[i]:
<hansl>
A = seq(4,1)
matrix B = 3
C = (B.>=A)?0:1
matrix pos = 0
pos = sum(C)+1
pos
pos<\hansl>
@ Richard: This is my first idea, but I was looking for a simpler way
avoiding any use of loops here.
Artur
Am 21.11.2013 12:58, schrieb Dr RJF Hudson:
> Hi there,
> this is an ordinary if then or else loop in Fortran and Basic, but
> instead
> of storing the found result in "C"
> I'suggest creating a (4x1) matrix columm vector and store the "found
> values
> in the column vector at each "i" location
> and I suspect it is the same in Gretl but I confess to being a novice at
> Gretl sorry to say right now
> Best wishes to all,
> Hope these few commnts are of small help
> Good luck
> Richard Hudson
>
>
>
>
> Dr RJF Hudson
> ----- Original Message -----
> From: "Artur T." <artur.tarassow(a)googlemail.com>
> To: "Gretl list" <gretl-users(a)lists.wfu.edu>
> Sent: Thursday, November 21, 2013 9:05 PM
> Subject: [Gretl-users] small practical issue
>
>
>> Hi gretl users,
>>
>> I am looking for a simple way to determine the column for which scalar B
>> is, let's say, larger or equal to the entry A[i] and store this
>> information in scalar C. If B is never greater than any entry in A[i], C
>> is zero.
>> Reading out each entry via a loop is possible, but I am wondering
>> whether there is a short-cut to this.
>>
>> C = 0
>> A = seq(4,1)
>> B = 3
>> # check at which entry A[i] the value of B >= A[i]
>>
>> Maybe somebody a nice suggestion :-)
>>
>> Artur
>>
>> _______________________________________________
>> Gretl-users mailing list
>> Gretl-users(a)lists.wfu.edu
>> http://lists.wfu.edu/mailman/listinfo/gretl-users
>
> _______________________________________________
> Gretl-users mailing list
> Gretl-users(a)lists.wfu.edu
> http://lists.wfu.edu/mailman/listinfo/gretl-users
>
--
For global instant messaging use the free and transparent Jabber chat:
t.artur(a)cryptolab.net
Avoid Facebook and join me on Diaspora: https://despora.de/i/9fba3608d45a
If you want to use encrypted emailing, here is my public PGP key:
-----BEGIN PGP PUBLIC KEY BLOCK-----
Version: GnuPG v1.4.12 (GNU/Linux)
mQINBFIFGVsBEACc0QWXQ8fijGzIvQUWQTTYE/0O/qKQRRjNX7O3UwkASTOZM12H
aZ3YzpqZOvfFnK/BMUCtOYfkoZJJwoLwPDio+wlmQ1DCDjNVYJhkxbjThs+pBdk0
YdyZ3b+bwylb8StnexA0EOSpdV8tOpWo3JVw2WXAs9y9MxierKHIr3aAsHNyGI4V
j+46Cz9mmfyQBJAA7tlK7Th1pibb18ww0fxFnAOmBNBvexCTb+z353FtgDdIc/Wq
LFLmdtP73Ntui6DYviGT+YvV2FBYFxF6x8TGwnYnnvQ88VGl+aflyVjp+LRUnbo0
pOD5LTAA6f6b5ZSZq8GzufLDrZxMRL1qS/+8lp8AAI0/x7y7GtRdznB4KYOAKd1Z
a/GgcuMl9N6T7B+azwxQa1th+ucMdSZO4PMOHjhsp0xYJftgEOG44GHJPqchS5Zd
nKphboRVJjV9KA/uLUV/S6n8Qhh6tRBLiQ2RPMSGjKTdKObMxJ9YbDirHChSgi5x
77iGf+3F413s6XCbazLYv+03lylFCBT3p9mYDWIMqJS2jCzWVjXbea6iu5ZMN5Ji
ESNiTHjN1pW97c+HaLxzo9/PobUceHYVzwKvEfBH4wVIzCExMlKRI0xwIp6ACk3E
XJPKsp8+a9HMp1qdCvi0zye/Th7FhIOPhVYlRgCf157O18MxLj4dJu2XewARAQAB
tCBBcnR1ciBULiA8dC5hcnR1ckBjcnlwdG9sYWIubmV0PokCPgQTAQIAKAUCUgUZ
WwIbIwUJA8JnAAYLCQgHAwIGFQgCCQoLBBYCAwECHgECF4AACgkQlEChRi+nhULG
FQ/9Fcz/aI7fcKOd9NG1Wb0ziyKHihg5uoEMFsOTFGDboC+EddrQwkGdkOpWhQzN
vnPCji+WuWHliG6DiRYViG+HYN+yBhaFB6ASe+jJYCXaoxI1u7e3UOscwOmg6zb7
TdR27jI0mguu0ZT1bFLWh7hDzqAa33ANWwX++J8ncc08SMw3/cJoddzUrkc1PcUY
M+mB2AgURVCC+TCzpsEuZ7jqTJAH3+aLUnxGMWm3hTlt2eirYieon+a9PlRqlTZD
0wOyd6z5mU+uUQr4lqfHNl9XOH/vP9oqDNW6uCItdGTmKtOVioE0BA+7p06tnobv
iUxNCc8dz81x+ZLoxjS2H6zCTf00PVf8ZdBWtxM0S0/9sRkuyTE1exwZOpob552m
HbQIIRIgB3+PwyM2mXgRSWXciITyH6LOvJoPDTp+vrHqUtUCX3LoIRP6ySqo7rEi
/UnF2WAwjjEPIVByvfS9ZVXQaL5ASHUV06yMAmMp7h6GfLgo4JNOr6wxl+sYgNwy
UimKVGx2OdAQoUPGN2VQA9SPxFtbDPTsH472s4FoEVGhq2OBSMdAVeKJgDDkrsQm
u5lk0+Or2GewwjxoBu77kRgihkqBX2ND6Z2pHs95g1nngrqM6sg895fp5Lo+MHAJ
+8cSUxs0diDA7ybG1j8KTx5Vx6yPLcmgj9MlHVaNiMnAKuC5Ag0EUgUZWwEQANDG
ZbQoate2ul8zMDwlr+DmlzSJBUMue4ABFXW0dJSZ+E8awJkYRL3jcjJqVVOrk9cT
gqy0Rum16SGswYrFa4xEE/gb8h97A/ZuJ9wOSe/Os5hySI+62cZSJrd+2AYI8rY7
3iW/QlaO38mFcsulSMg7+SqZjHpFI3tfMj+5e3vx1FsYazAR/jI4U02gP6OdHj5h
6oAS7OqlCSvN0j0LatTuc7LvfjLkSWCWyZUfe2pujUCc1gZ95wINzSamYcxfd+o4
MjOizFULw5zmuioiHWlp35nGLu03uOaBJgR85dTvyQqa4aQXLCfwvSICBEl6aAZp
ODtfZUOpT/sQdqIpGrc4rytz763BpnuVl92RzNoQtjsLE8wnJRzElFjYeOI4JZLL
RhEa0iNVpbFvv4rmOqVafxL0dR8wkCjiPMgWii/qm6crpsOnhOJ3Wov8EGVL53hI
0OOGxkQ0YV9PF8L45vcvK8TMSm+V19l53yQYOvxPECcCksn9uFCdMNvPoef7FZLO
nXIORP4IBcMTKf2N07QFbfkZRVh5SZgUOFxMtLIqBmWQqpmEAdO7SCybJkAgwdpJ
VliEWQSxFeLwyV9kQUmN7oQsQ8jYxO1C0vTnstDjZYC4EBPf08kCPLQXI1gM3/v2
POKDZVk6EAFQelVia8OIFyHA9mLwc9HV7b7KYwGrABEBAAGJAiUEGAECAA8FAlIF
GVsCGwwFCQPCZwAACgkQlEChRi+nhUKh5A/7B75inrTA4PMWaY6fjA4UnhFi6BBy
bY/QIpkv7zle4H6y0tbDsAziB80i0GsMqDka+LEmCr9BBmZNyD6vUeFbE0uECiE8
XupOmw/B3+Czscuo29CqtMYpwlMX9oVXFiaWd/LHsIMOgI5BstDaE8pBga+IMyKz
HsPdSDz+MqdQwp/9pY5FfB1/h4VpFNXOIsEDjimALvwCHW9F34Nv/tFCd39Li2sd
jktLLgJKVurEHisEfws4vJRsmPtsFdgqGGHm0r/wxJYROY8ki+veZCNpYmwVWAcY
iZK7AFbqqsNMdMyKXOiMGawq9LZSnpBmgiBzKJmqAL6BNqEO4/M8E+zz8rlrjD3/
16iTljR+FPriXF/g0KWli1C1Fx6ntCYQUElO9Hu6X4CEA5zT7T6vXTapRkMA51PM
xepwgASu9YWRAig8qZI18X2uh0I3VdhERxgAMmDn8V1qjTn/gd0S7RPYtduUHgNW
/sCV5B6vZsmaO0fImKCJdJjVZTCViW5h9XHFSNH7hXjnr/+4SKrvfm4Du1DPiGKn
4oD7/QUHeSpdfDKJo/E9nM39aKxgN6+RqJUd2J2Db2aNBCQwa1o+SOtrriqyCMjp
K+4gDumAFG2lg/PD/ki/aNucMROO9t/uPFQFhheaUF9R1rXJh66KyxQFsxKC00uZ
kJ87csTGWSWmWns=
=tQPP
-----END PGP PUBLIC KEY BLOCK-----
_______________________________________________
Gretl-users mailing list
Gretl-users(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-users
4035
days inactive
4035
days old
6 comments
4 participants
participants (4)
-
Allin Cottrell -
Artur T. -
Dr RJF Hudson -
Hélio Guilherme