[R] Urgent Help needed

2007-08-16 Thread AbouEl-Makarim Aboueissa
Dear All:

Urgent help is needed.


I have a data set in matrix format  of three columns: X, Y and index of four 
groups (1,2,3,4). What I need to do is the following;

1- How I can subtract the sample mean of each group indexed 1,2,3,4 from the 
 corresponding data values of this group and create new columns say 
X-sample mean 
  and Y-sample mean? I tried to use the tapply but I have some 
difficulties to restore the new data


2- How I can use the “tapply” if possible or any other R-function to find the 
correlation 
 coefficient between the X and Y columns for each group indexed 1,2,3,4.? 
Could not use the tapply.


I attached part of the data as txt file.


Thank you so much for your attention to this matter, and I look forward to hear 
from you soon.

Regards,

Abou


Data:

x   y   index
15807.2412.54
15752.5133.54
12893.7601.53
8426.88 22.23
5706.24 333 3
3982.08 560 2
3642.62 670 2
295.68  124 1
215.40  104 1
195.40  204 1
4240.21 22.42
1222.72 45.92
1142.26 23.62
63.00   90.11
1216.00 82.42
2769.60 111 2
1790.46 34.72
26.10   26.10   1
19676.830.994
10920.60203 3
6144.00 46  3
4534.48 4534.48 3
4.0065  4
29500.0056  4
17100.0077  4
9000.00 435 3
6300.00 84  3
3962.88 334 2
5690.00 653 3
3736.00 233 2
2750.00 22  2
1316.00 345 2
4595.00 4595.00 3
5928.00 45  3
2645.70 0.002
2580.24 454 2
6547.34 6547.34 3
1615.68 5   2
194.06  55  1
184.80  6   1
82.94   44  1
16649.0056  4
4500.00 74  3
1600.00 744 2

=



==
AbouEl-Makarim Aboueissa, Ph.D.
Assistant Professor of Statistics
Department of Mathematics  Statistics
University of Southern Maine
96 Falmouth Street
P.O. Box 9300
Portland, ME 04104-9300

Tel: (207) 228-8389
Email: [EMAIL PROTECTED]
  [EMAIL PROTECTED]
Office: 301C Payson Smith

x   y   index
15807.2412.54
15752.5133.54
12893.7601.53
8426.88 22.23
5706.24 333 3
3982.08 560 2
3642.62 670 2
295.68  124 1
215.40  104 1
195.40  204 1
4240.21 22.42
1222.72 45.92
1142.26 23.62
63.00   90.11
1216.00 82.42
2769.60 111 2
1790.46 34.72
26.10   26.10   1
19676.830.994
10920.60203 3
6144.00 46  3
4534.48 4534.48 3
4.0065  4
29500.0056  4
17100.0077  4
9000.00 435 3
6300.00 84  3
3962.88 334 2
5690.00 653 3
3736.00 233 2
2750.00 22  2
1316.00 345 2
4595.00 4595.00 3
5928.00 45  3
2645.70 0.002
2580.24 454 2
6547.34 6547.34 3
1615.68 5   2
194.06  55  1
184.80  6   1
82.94   44  1
16649.0056  4
4500.00 74  3
1600.00 744 2
__
R-help@stat.math.ethz.ch 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.


Re: [R] Urgent Help needed

2007-08-16 Thread Weiwei Shi
try this:

t0 = read.table(datatest.txt, header=T)
X.mean = ave(t0[,1], as.factor(t0[,3]))

you do the rest of Y.mean and make them into a data.fame or whatever.

HTH,

Weiwei

On 8/16/07, AbouEl-Makarim Aboueissa [EMAIL PROTECTED] wrote:
 Dear All:

 Urgent help is needed.


 I have a data set in matrix format  of three columns: X, Y and index of four 
 groups (1,2,3,4). What I need to do is the following;

 1- How I can subtract the sample mean of each group indexed 1,2,3,4 from the
  corresponding data values of this group and create new columns say 
 X-sample mean
   and Y-sample mean? I tried to use the tapply but I have some 
 difficulties to restore the new data


 2- How I can use the tapply if possible or any other R-function to find the 
 correlation
  coefficient between the X and Y columns for each group indexed 1,2,3,4.? 
 Could not use the tapply.


 I attached part of the data as txt file.


 Thank you so much for your attention to this matter, and I look forward to 
 hear from you soon.

 Regards,

 Abou


 Data:
 
 x   y   index
 15807.2412.54
 15752.5133.54
 12893.7601.53
 8426.88 22.23
 5706.24 333 3
 3982.08 560 2
 3642.62 670 2
 295.68  124 1
 215.40  104 1
 195.40  204 1
 4240.21 22.42
 1222.72 45.92
 1142.26 23.62
 63.00   90.11
 1216.00 82.42
 2769.60 111 2
 1790.46 34.72
 26.10   26.10   1
 19676.830.994
 10920.60203 3
 6144.00 46  3
 4534.48 4534.48 3
 4.0065  4
 29500.0056  4
 17100.0077  4
 9000.00 435 3
 6300.00 84  3
 3962.88 334 2
 5690.00 653 3
 3736.00 233 2
 2750.00 22  2
 1316.00 345 2
 4595.00 4595.00 3
 5928.00 45  3
 2645.70 0.002
 2580.24 454 2
 6547.34 6547.34 3
 1615.68 5   2
 194.06  55  1
 184.80  6   1
 82.94   44  1
 16649.0056  4
 4500.00 74  3
 1600.00 744 2

 =



 ==
 AbouEl-Makarim Aboueissa, Ph.D.
 Assistant Professor of Statistics
 Department of Mathematics  Statistics
 University of Southern Maine
 96 Falmouth Street
 P.O. Box 9300
 Portland, ME 04104-9300

 Tel: (207) 228-8389
 Email: [EMAIL PROTECTED]
   [EMAIL PROTECTED]
 Office: 301C Payson Smith


 __
 R-help@stat.math.ethz.ch 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.





-- 
Weiwei Shi, Ph.D
Research Scientist
GeneGO, Inc.

Did you always know?
No, I did not. But I believed...
---Matrix III

__
R-help@stat.math.ethz.ch 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.


Re: [R] Urgent Help needed

2007-08-16 Thread Henrique Dallazuanna
For the 2nd item, perhaps:

by(df[,1:2], df$index, FUN=cor)

where df is your data.frame.

-- 
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O

On 16/08/07, AbouEl-Makarim Aboueissa [EMAIL PROTECTED] wrote:

 Dear All:

 Urgent help is needed.


 I have a data set in matrix format  of three columns: X, Y and index of
 four groups (1,2,3,4). What I need to do is the following;

 1- How I can subtract the sample mean of each group indexed 1,2,3,4 from
 the
  corresponding data values of this group and create new columns say
 X-sample mean
   and Y-sample mean? I tried to use the tapply but I have some
 difficulties to restore the new data


 2- How I can use the tapply if possible or any other R-function to find
 the correlation
  coefficient between the X and Y columns for each group indexed
 1,2,3,4.? Could not use the tapply.


 I attached part of the data as txt file.


 Thank you so much for your attention to this matter, and I look forward to
 hear from you soon.

 Regards,

 Abou


 Data:
 
 x   y   index
 15807.2412.54
 15752.5133.54
 12893.7601.53
 8426.88 22.23
 5706.24 333 3
 3982.08 560 2
 3642.62 670 2
 295.68  124 1
 215.40  104 1
 195.40  204 1
 4240.21 22.42
 1222.72 45.92
 1142.26 23.62
 63.00   90.11
 1216.00 82.42
 2769.60 111 2
 1790.46 34.72
 26.10   26.10   1
 19676.830.994
 10920.60203 3
 6144.00 46  3
 4534.48 4534.48 3
 4.0065  4
 29500.0056  4
 17100.0077  4
 9000.00 435 3
 6300.00 84  3
 3962.88 334 2
 5690.00 653 3
 3736.00 233 2
 2750.00 22  2
 1316.00 345 2
 4595.00 4595.00 3
 5928.00 45  3
 2645.70 0.002
 2580.24 454 2
 6547.34 6547.34 3
 1615.68 5   2
 194.06  55  1
 184.80  6   1
 82.94   44  1
 16649.0056  4
 4500.00 74  3
 1600.00 744 2

 =



 ==
 AbouEl-Makarim Aboueissa, Ph.D.
 Assistant Professor of Statistics
 Department of Mathematics  Statistics
 University of Southern Maine
 96 Falmouth Street
 P.O. Box 9300
 Portland, ME 04104-9300

 Tel: (207) 228-8389
 Email: [EMAIL PROTECTED]
   [EMAIL PROTECTED]
 Office: 301C Payson Smith


 __
 R-help@stat.math.ethz.ch 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.




[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch 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.


Re: [R] Urgent Help needed

2007-08-16 Thread Marc Schwartz
On Thu, 2007-08-16 at 12:33 -0400, AbouEl-Makarim Aboueissa wrote:
 Dear All:
 
 Urgent help is needed.
 
 
 I have a data set in matrix format  of three columns: X, Y and index
 of four groups (1,2,3,4). What I need to do is the following;
 
 1- How I can subtract the sample mean of each group indexed 1,2,3,4
 from the 
  corresponding data values of this group and create new columns
 say X-sample mean 
   and Y-sample mean? I tried to use the tapply but I have some
 difficulties to restore the new data
 
 
 2- How I can use the “tapply” if possible or any other R-function to
 find the correlation 
  coefficient between the X and Y columns for each group indexed
 1,2,3,4.? Could not use the tapply.
 
 
 I attached part of the data as txt file.
 
 
 Thank you so much for your attention to this matter, and I look
 forward to hear from you soon.
 
 Regards,
 
 Abou
 
 
 Data:
 
 x y   index
 15807.24  12.54
 15752.51  33.54
 12893.76  01.53
 8426.88   22.23
 5706.24   333 3
 3982.08   560 2
 3642.62   670 2
 295.68124 1
 215.40104 1
 195.40204 1
 4240.21   22.42
 1222.72   45.92
 1142.26   23.62
 63.00 90.11
 1216.00   82.42
 2769.60   111 2
 1790.46   34.72
 26.10 26.10   1
 19676.83  0.994
 10920.60  203 3
 6144.00   46  3
 4534.48   4534.48 3
 4.00  65  4
 29500.00  56  4
 17100.00  77  4
 9000.00   435 3
 6300.00   84  3
 3962.88   334 2
 5690.00   653 3
 3736.00   233 2
 2750.00   22  2
 1316.00   345 2
 4595.00   4595.00 3
 5928.00   45  3
 2645.70   0.002
 2580.24   454 2
 6547.34   6547.34 3
 1615.68   5   2
 194.0655  1
 184.806   1
 82.94 44  1
 16649.00  56  4
 4500.00   74  3
 1600.00   744 2
 
 =


I might be tempted to take the following approach:

If your data is a matrix, coerce it to a data frame first. Let's call
that 'DF'.

 str(DF)
'data.frame':   44 obs. of  3 variables:
 $ x: num  15807 15753 12894  8427  5706 ...
 $ y: num  12.5 33.5 1.5 22.2 333 560 670 124 104 204 ...
 $ index: int  4 4 3 3 3 2 2 1 1 1 ...


Now use split() to break up the data frame into a list of 4
sub-dataframes, based upon the index value.  We can use scale() within a
lapply() loop to center the 'x' and 'y' columns for each sub-dataframe:


DF.ctr - lapply(split(DF[, -3], DF$index), scale, scale = FALSE)


 str(DF.ctr)
List of 4
 $ 1: num [1:8, 1:2]  138.5   58.2   38.2  -94.2 -131.1 ...
  ..- attr(*, dimnames)=List of 2
  .. ..$ : chr [1:8] 8 9 10 14 ...
  .. ..$ : chr [1:2] x y
  ..- attr(*, scaled:center)= Named num [1:2] 157.2  81.7
  .. ..- attr(*, names)= chr [1:2] x y
 $ 2: num [1:16, 1:2]  1469  1129  1727 -1291 -1371 ...
  ..- attr(*, dimnames)=List of 2
  .. ..$ : chr [1:16] 6 7 11 12 ...
  .. ..$ : chr [1:2] x y
  ..- attr(*, scaled:center)= Named num [1:2] 2513  230
  .. ..- attr(*, names)= chr [1:2] x y
 $ 3: num [1:13, 1:2]  5879  1413 -1308  3906  -870 ...
  ..- attr(*, dimnames)=List of 2
  .. ..$ : chr [1:13] 3 4 5 20 ...
  .. ..$ : chr [1:2] x y
  ..- attr(*, scaled:center)= Named num [1:2] 7014 1352
  .. ..- attr(*, names)= chr [1:2] x y
 $ 4: num [1:7, 1:2] -6262 -6317 -2393 17931  7431 ...
  ..- attr(*, dimnames)=List of 2
  .. ..$ : chr [1:7] 1 2 19 23 ...
  .. ..$ : chr [1:2] x y
  ..- attr(*, scaled:center)= Named num [1:2] 2206943
  .. ..- attr(*, names)= chr [1:2] x y


Now, create a new single DF comprised of the sub-dataframes from DF.ctr:

DF.new - do.call(rbind, DF.ctr)


Define colnames:

colnames(DF.new) - c(x-mean, y-mean)


 str(DF.new)
 num [1:44, 1:2]  138.5   58.2   38.2  -94.2 -131.1 ...
 - attr(*, dimnames)=List of 2
  ..$ : chr [1:44] 8 9 10 14 ...
  ..$ : chr [1:2] x-mean y-mean


Now, use merge() to join DF and DF.new by the rownames:

DF.final - merge(DF, DF.new, by = row.names)

 DF.final
   Row.namesx   y index  x-mean   y-mean
1  1 15807.24   12.50 4 -6262.12857   -30.498571
2 10   195.40  204.00 138.22750   122.35
3 11  4240.21   22.40 2  1726.93188  -208.037500
4 12  1222.72   45.90 2 -1290.55812  -184.537500
5 13  1142.26   23.60 2 -1371.01812  -206.837500
6 1463.00   90.10 1   -94.17250 8.45
7 15  1216.00   82.40 2 -1297.27812  -148.037500
8 16  2769.60  111.00 2   256.32188  -119.437500
9 17  1790.46   34.70 2  -722.81812  -195.737500
101826.10   26.10 1  -131.07250   -55.55
1119 19676.830.99 4 -2392.53857   -42.008571
12 2 15752.51   33.50 4 -6316.85857-9.498571
1320 10920.60  

Re: [R] Urgent Help needed

2007-08-16 Thread AbouEl-Makarim Aboueissa
Thanks all.

it works.

Just one more thing: if you look to this out put,

 by(data1[,2:3], data1[,4], cor)[1]
$`1`
XY
X   1.0000.4400451
Y   0.4400451  1.000

Q. How I can just pick the value of the correlation 0.4400451 from this output 
and call it sat corxy.


Once again thank you all so much for your helps.


Abou


==
AbouEl-Makarim Aboueissa, Ph.D.
Assistant Professor of Statistics
Department of Mathematics  Statistics
University of Southern Maine
96 Falmouth Street
P.O. Box 9300
Portland, ME 04104-9300

Tel: (207) 228-8389
Email: [EMAIL PROTECTED]
  [EMAIL PROTECTED]
Office: 301C Payson Smith

 Henrique Dallazuanna [EMAIL PROTECTED] 8/16/2007 2:05 PM 
For the 2nd item, perhaps:

by(df[,1:2], df$index, FUN=cor)

where df is your data.frame.

-- 
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O

On 16/08/07, AbouEl-Makarim Aboueissa [EMAIL PROTECTED] wrote:

 Dear All:

 Urgent help is needed.


 I have a data set in matrix format  of three columns: X, Y and index of
 four groups (1,2,3,4). What I need to do is the following;

 1- How I can subtract the sample mean of each group indexed 1,2,3,4 from
 the
  corresponding data values of this group and create new columns say
 X-sample mean
   and Y-sample mean? I tried to use the tapply but I have some
 difficulties to restore the new data


 2- How I can use the tapply if possible or any other R-function to find
 the correlation
  coefficient between the X and Y columns for each group indexed
 1,2,3,4.? Could not use the tapply.


 I attached part of the data as txt file.


 Thank you so much for your attention to this matter, and I look forward to
 hear from you soon.

 Regards,

 Abou


 Data:
 
 x   y   index
 15807.2412.54
 15752.5133.54
 12893.7601.53
 8426.88 22.23
 5706.24 333 3
 3982.08 560 2
 3642.62 670 2
 295.68  124 1
 215.40  104 1
 195.40  204 1
 4240.21 22.42
 1222.72 45.92
 1142.26 23.62
 63.00   90.11
 1216.00 82.42
 2769.60 111 2
 1790.46 34.72
 26.10   26.10   1
 19676.830.994
 10920.60203 3
 6144.00 46  3
 4534.48 4534.48 3
 4.0065  4
 29500.0056  4
 17100.0077  4
 9000.00 435 3
 6300.00 84  3
 3962.88 334 2
 5690.00 653 3
 3736.00 233 2
 2750.00 22  2
 1316.00 345 2
 4595.00 4595.00 3
 5928.00 45  3
 2645.70 0.002
 2580.24 454 2
 6547.34 6547.34 3
 1615.68 5   2
 194.06  55  1
 184.80  6   1
 82.94   44  1
 16649.0056  4
 4500.00 74  3
 1600.00 744 2

 =



 ==
 AbouEl-Makarim Aboueissa, Ph.D.
 Assistant Professor of Statistics
 Department of Mathematics  Statistics
 University of Southern Maine
 96 Falmouth Street
 P.O. Box 9300
 Portland, ME 04104-9300

 Tel: (207) 228-8389
 Email: [EMAIL PROTECTED] 
   [EMAIL PROTECTED] 
 Office: 301C Payson Smith


 __
 R-help@stat.math.ethz.ch 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.




__
R-help@stat.math.ethz.ch 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.


Re: [R] Urgent Help needed

2007-08-16 Thread Henrique Dallazuanna
Hi, try this:

by(df[,1:2], df$index, FUN=function(x)cor(x[1],x[2]))

-- 
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O

On 16/08/07, AbouEl-Makarim Aboueissa [EMAIL PROTECTED] wrote:

 Thanks all.

 it works.

 Just one more thing: if you look to this out put,

  by(data1[,2:3], data1[,4], cor)[1]
 $`1`
 XY
 X   1.0000.4400451
 Y   0.4400451  1.000

 Q. How I can just pick the value of the correlation 0.4400451 from this
 output and call it sat corxy.


 Once again thank you all so much for your helps.


 Abou


 ==
 AbouEl-Makarim Aboueissa, Ph.D.
 Assistant Professor of Statistics
 Department of Mathematics  Statistics
 University of Southern Maine
 96 Falmouth Street
 P.O. Box 9300
 Portland, ME 04104-9300

 Tel: (207) 228-8389
 Email: [EMAIL PROTECTED]
   [EMAIL PROTECTED]
 Office: 301C Payson Smith

  Henrique Dallazuanna [EMAIL PROTECTED] 8/16/2007 2:05 PM 
 For the 2nd item, perhaps:

 by(df[,1:2], df$index, FUN=cor)

 where df is your data.frame.

 --
 Henrique Dallazuanna
 Curitiba-Paraná-Brasil
 25° 25' 40 S 49° 16' 22 O

 On 16/08/07, AbouEl-Makarim Aboueissa [EMAIL PROTECTED] wrote:
 
  Dear All:
 
  Urgent help is needed.
 
 
  I have a data set in matrix format  of three columns: X, Y and index of
  four groups (1,2,3,4). What I need to do is the following;
 
  1- How I can subtract the sample mean of each group indexed 1,2,3,4 from
  the
   corresponding data values of this group and create new columns say
  X-sample mean
and Y-sample mean? I tried to use the tapply but I have some
  difficulties to restore the new data
 
 
  2- How I can use the tapply if possible or any other R-function to
 find
  the correlation
   coefficient between the X and Y columns for each group indexed
  1,2,3,4.? Could not use the tapply.
 
 
  I attached part of the data as txt file.
 
 
  Thank you so much for your attention to this matter, and I look forward
 to
  hear from you soon.
 
  Regards,
 
  Abou
 
 
  Data:
  
  x   y   index
  15807.2412.54
  15752.5133.54
  12893.7601.53
  8426.88 22.23
  5706.24 333 3
  3982.08 560 2
  3642.62 670 2
  295.68  124 1
  215.40  104 1
  195.40  204 1
  4240.21 22.42
  1222.72 45.92
  1142.26 23.62
  63.00   90.11
  1216.00 82.42
  2769.60 111 2
  1790.46 34.72
  26.10   26.10   1
  19676.830.994
  10920.60203 3
  6144.00 46  3
  4534.48 4534.48 3
  4.0065  4
  29500.0056  4
  17100.0077  4
  9000.00 435 3
  6300.00 84  3
  3962.88 334 2
  5690.00 653 3
  3736.00 233 2
  2750.00 22  2
  1316.00 345 2
  4595.00 4595.00 3
  5928.00 45  3
  2645.70 0.002
  2580.24 454 2
  6547.34 6547.34 3
  1615.68 5   2
  194.06  55  1
  184.80  6   1
  82.94   44  1
  16649.0056  4
  4500.00 74  3
  1600.00 744 2
 
  =
 
 
 
  ==
  AbouEl-Makarim Aboueissa, Ph.D.
  Assistant Professor of Statistics
  Department of Mathematics  Statistics
  University of Southern Maine
  96 Falmouth Street
  P.O. Box 9300
  Portland, ME 04104-9300
 
  Tel: (207) 228-8389
  Email: [EMAIL PROTECTED]
[EMAIL PROTECTED]
  Office: 301C Payson Smith
 
 
  __
  R-help@stat.math.ethz.ch 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.
 
 
 



[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch 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.


[R] Urgent: How to obtain the Consistent Standard Errors after apply 2SLS through tsls() from sem or systemfit(2SLS) without this error message !!!!!!!!!!!!!

2007-02-19 Thread Guillermo Julián San Martín
Hi,



I am trying to obtain the heteroskedasticity consitent standard errors
(HCSE) after apply 2SLS. I obtain 2SLS through tsls from package sem or
systemfit:



 tsls 

library (sem)

Reg2SLS -tsls(LnP~Sc+Ag+Ag2+Var+R+D,~I2+Ag+Ag2+Var+R+D)

summary (Reg2SLS)



 systemfit 

library (systemfit)

RS - LnP~Sc+Ag+Ag2+Var+R+D

Inst - ~I2+Ag+Ag2+Var+R+D

labels -list(RS)

system -list(RS)

Reg2SLS - systemfit(2SLS, system, labels, Inst, saveMemory=TRUE)

summary (Reg2SLS)



If I try to obtain the HCSE with robcov I obtain the following error message
in both cases:



Error in rep.default(1, p) : rep() incorrect type for second argument



If I tried to apply vcovHAC or vcovHC in the systemfit output of 2SLS I
receive the following messages of errors respectively:



Error in estfun(x) : no applicable method for estfun (for vcovHAC)



Error in terms.default(object) : no terms component (for vcovHC)



If somebody can help me to solve this problem I will thank you a lot.



Best regards



Guillermo

[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch 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.


[R] Urgent: How to obtain the Consistent Standard Errors after apply 2SLS through tsls() from sem or systemfit(2SLS) without this error message !!!!!!!!!!!!!

2007-02-18 Thread Guillermo Julián San Martín
Hi,

I am trying to obtain the heteroskedasticity consitent standard errors
(HCSE) after apply 2SLS. I obtain 2SLS through tsls from package sem or
systemfit:


 tsls 

library (sem)

Reg2SLS -tsls(LnP~Sc+Ag+Ag2+Var+R+D,~I2+Ag+Ag2+Var+R+D)

summary (Reg2SLS)



 systemfit 

library (systemfit)

RS - LnP~Sc+Ag+Ag2+Var+R+D

Inst - ~I2+Ag+Ag2+Var+R+D

labels -list(RS)

system -list(RS)

Reg2SLS - systemfit(2SLS, system, labels, Inst, saveMemory=TRUE)

summary (Reg2SLS)

If I try to obtain the HCSE with robcov I obtain the following error message
in both cases:


Error in rep.default(1, p) : rep() incorrect type for second argument



If I tried to apply vcovHAC or vcovHC in the systemfit output of 2SLS I
receive the following messages of errors respectively:



Error in estfun(x) : no applicable method for estfun (for vcovHAC)



Error in terms.default(object) : no terms component (for vcovHC)



If somebody can help me to solve this problem I will thank you a lot.



Best regards



Guillermo

[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch 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.


[R] urgent

2005-12-06 Thread Subhabrata

Hello R Users,

I have two sets of values

x - c(7, 7 , 8, 9, 15, 17, 18)

y - c(7, 8, 9, 15, 17, 19, 20, 20, 25, 23, 22)

I am able to create multi histogram using
multhist(). But not able to control the 'xlim'.
ie the xaxis is showing 7.5, 13, 18, 23

1st on what basis it is calculated

2nd I want it to be like 7 8 9 15 17 and so on


Can any one help me


With Regards
Subhabrata Pal
[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] urgent

2005-12-06 Thread Adaikalavan Ramasamy
1) R-help mailing list is run entirely by volunteers, so requests such
as urgent may sound rude

2) Use an informative subject line please !

3) Please state which package multhist comes from.

4) Please show your call to multhist.

5) multhist does _histograms_ by aggregating points within certain
intervals. In your case, you simply want a plot of your raw data. You
can use barplot directly via


 multi.barplot - function( mylist, ... ){ 
   u   -  unique( unlist( mylist ) )
   tb  -  t(sapply( mylist, function(v) table(factor(v, levels=u)) ) ) 
   barplot( tb, beside=TRUE, ... )
   return(tb)
 }
 

 x - c(7, 7 , 8, 9, 15, 17, 18)
 y - c(7, 8, 9, 15, 17, 19, 20, 20, 25, 23, 22)
 z - c(8, 9, 9, 9, 31)
 multi.barplot( list(x, y, z), col=1:3 )
 legend( topright, legend=c(one, two, three), fill=1:3 )


Regards, Adai



On Tue, 2005-12-06 at 15:32 +0530, Subhabrata wrote:
 Hello R Users,
 
 I have two sets of values
 
 x - c(7, 7 , 8, 9, 15, 17, 18)
 
 y - c(7, 8, 9, 15, 17, 19, 20, 20, 25, 23, 22)
 
 I am able to create multi histogram using
 multhist(). But not able to control the 'xlim'.
 ie the xaxis is showing 7.5, 13, 18, 23
 
 1st on what basis it is calculated
 
 2nd I want it to be like 7 8 9 15 17 and so on
 
 
 Can any one help me
 
 
 With Regards
 Subhabrata Pal
   [[alternative HTML version deleted]]
 
 __
 R-help@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Re: [R] urgent

2005-12-06 Thread P Ehlers
I don't have an answer to your query, but I do have
three suggestions:

1. Use a sensible subject line. This may be urgent to
you, but I doubt that it is to anyone else.

2. Do indicate what package contains multhist(). I have
no idea (nor do I know what a 'multi histogram' is).

3. Don't send HTML mail.

People are very willing to help, but you do have to
make it easy to do so.

Peter Ehlers


Subhabrata wrote:

 Hello R Users,
 
 I have two sets of values
 
 x - c(7, 7 , 8, 9, 15, 17, 18)
 
 y - c(7, 8, 9, 15, 17, 19, 20, 20, 25, 23, 22)
 
 I am able to create multi histogram using
 multhist(). But not able to control the 'xlim'.
 ie the xaxis is showing 7.5, 13, 18, 23
 
 1st on what basis it is calculated
 
 2nd I want it to be like 7 8 9 15 17 and so on
 
 
 Can any one help me
 
 
 With Regards
 Subhabrata Pal
   [[alternative HTML version deleted]]
 
 __
 R-help@stat.math.ethz.ch mailing list
 https://stat.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] urgent request

2005-03-09 Thread h0394018
Hi guys,
I want to do weighted conditional logistic regression, but clogit do not accept
weights. Would you tell me whether there are any ways to do this?
Thanks a lot in advance,
Spring

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] urgent

2004-11-25 Thread Jennifer Wilson
Hello,

My name is Mrs. Jennifer Wilson i am a dying woman who have decided to donate 
what i have to you/ church.I am 59 years old and i was diagnosed for cancer for 
about 2 years ago,immediately after the death of my husband, who has left me 
everything he worked for.

I have been touched by God to donate from what i have inherited from my late 
husband to the you for the god work of God,rather than allow my relatives to 
use my husband hard earned funds ungodly.Please pray,that the good Lord forgive 
me my sins.I have asked God to forgive me and i beleive he has because He is a 
merciful God. I will be going in for an operation in less than one hour.

I decided to WILL/donate the sum of $1,500,000 (One million five hundred 
thousand dollars) to you for the good work of the lord, and also to help the 
motherless and less privilege and also for the assistance of the widows 
according to (JAMES 1:27).

At the moment i cannot take any telephone calls right now due to the fact that 
my relatives are around me and my health status.I have adjusted my WILL and my 
lawyer is aware i have changed my will you and he will arrange the transfer of 
the funds from my account to you.

I wish you all the best and may the good Lord bless you abundantly, and please 
use the funds well and always extend the good work to others. Contact my lawyer 
with this specified email   [EMAIL PROTECTED]   and tell him that i have WILLED 
($1,500,000.00) to you and i have also notified him that i am WILLING that 
amount to you for a specific and good work.I know i dont know you but i have 
been directed to do this.Thanks and God bless.

NB: I will appreciate your utmost confidentiality in this matter until the task 
is accomplished as I don't want anything that will Jeopardize my last wish. And 
Also I will be contacting with you by email as I don't want my relation or 
anybody to know because they are always around me. 

Regards,
Jennifer Wilson

__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Spam {Re: [R] urgent}

2004-11-25 Thread Martin Maechler
I'm sorry for this spam that astonishingly came through to the
mailing list.

It seems these guys have been exercising against known spam
filters and achieved more than in the past.
Also, recent versions of our spamfilter have been tuned such as
to rather produce a few false negatives {spam not detected}
with hardly ever any false positive {non-spam not delivered to recipient},
which makes a lot of sense.

We (I and local e-mail administrators) do keep an eye on this, 
me spending a little time for manual tuning, but we do not
want to allocate too much time for this. 

PLEASE do not reply to this e-mail (at least not to R-help!). 
It is *not* relevant to R and not worth the time (also since too
many people think they know what they are talking about :-).

If you want, reply to me privately.

Martin Maechler, ETH Zurich

__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


[R] Urgent Help on Cumulant computation on two R.V's

2003-02-12 Thread Feng Zhang
Hey

Now I am going to check the independence of random variables using cumulant
function.

So if R has such package or functions to calculate
the sample cumulant of a random vector?

Thanks a lot.

Fred

__
[EMAIL PROTECTED] mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help