Re: [R] random uniform sample of points on an ellipsoid (e.g. WG

2007-03-02 Thread R Heberto Ghezzo, Dr
OK, I stand corrected, my previous post generated a point on an ellipsoid as a proyection of an uniform random point in a sphere. now, to generate a random point in a sphere I generate an uniform on one axis and an uniform angle in the second, since the sphere has the property that the perifery

Re: [R] random uniform sample of points on an ellipsoid (e.g. WG

2007-03-01 Thread R Heberto Ghezzo, Dr
Hello, I do not know if I am completely out of it but . . . if x,y,z is a point in a sphere and [u,v,w]'A[u,v,w] = 1 is the equation of an ellipsoid and A = T'T (cholesky) then T.[x,y,z] should be a point in the ellipsoid ? isn't it? Heberto Ghezzo McGill University Montreal Canada _

[R] lmer formula

2006-12-18 Thread R Heberto Ghezzo, Dr
Hello, I have some problems trying to write up the formula for lmer. I have 43 subjects ( random factor) which were seen twice ( Visit : repeated measure - fixed). on each visit the patient performed a graded effort exercise ( effort : repeated measures, ordered, fixed 4 levels). So Subject is a

[R] Compiling a package

2006-09-01 Thread R Heberto Ghezzo, Dr
Hello, I am in Win-XP R:2.3.0 latest rtools and Perl - of today I got Rcmdr.HH source code and tried to compile it myself copy all directory to R/R-2.3.0/src/library/Rcmdr.HH from R/R-2.3.0/src/library I typed: ..\..\bin\R CMD build --force --binary --auto-zip Rcmdr.HH * checking for file 'Rcmdr.HH

Re: [R] Correlation adjusted Bonferroni? (was: Multiple tests on repeated measurements)

2006-08-02 Thread R Heberto Ghezzo, Dr
HI, just my 2 cents. Bonferroni et al, assume independent tests, thus p ~ p*k with k the number of tests, in repeated measures, each measure is correlated with the previous, so k is not 2 but 2-q. I do not know q but it should be a function of the correlation between measures, ie the Sigma in ge

[R] install.packages

2004-04-15 Thread R. Heberto Ghezzo
Hello, I just downloaded RW1090. No problems. My thanks to everybody involved in the project. I work in Win98 I updated my library, some problems with some files that were in the PACKAGES list but not in 1.9/ site, now all are. I tried to install "Zelig" from Harvard install.packages("Zelig",CRAN="

[R] problems with library in 1.7.1

2003-06-26 Thread R. Heberto Ghezzo
Hello. I am using R 1.7.1 just downloaded on Win98. With the old 1.6.2 I had modified the etc/Rprofile file as etc/Rprofile # Things you might want to change # options(width=80) # options(papersize="a4") # options(editor="notepad") # options(pager="internal") # to prefer Compiled HTML help options

[R] RODBC and Excel in Widows

2003-03-13 Thread R. Heberto Ghezzo
Hello, I have some problems with RODBC and Excel in Win98 I am using R 1.6.2 and just upgraded RODBC to the last version on CRAN. I have an Excel file with columns Number, Name, Sex, Age, FEV1 on Sheet 1 and Number, Age, FEV1, Name, Sex on Sheet 2. Now I open the channel to the file > chan1 <- odbc