[R] RPgSQL W2K

2003-02-10 Thread Christian Schulz
Hi, i found tis message in the archive and have got the same problems ? John, perhaps you have found now a way install RPgSQL for windows2000, or anybody other ? many thanks for advance christian [SNIP] I wonder whether anyone has succeeded in building either the RPgSQL or the Rdbi.PGSQL

[R] Covariance matrix for GMM

2003-02-10 Thread Feng Zhang
Hey, All Now I generate a data vector X (d-dimension column vector) from a Gaussian Mixture Model (GMM). That is, the pdf of vector X is f(X) = a1*N(u1, Cov1) + a2*(u2, Cov2) where a1+a2 = 1, N is multidimensional normal distribution, ui is the mean vecotr, Covi is the covariance matrix, i=1, 2.

Re: [R] installation on FreeBSD

2003-02-10 Thread Peng Zhang
Dr Mccarter and Dr Gentleman, Thank you for your reply. I reinstalled the system and found out that you can not install R on FreeBSD 5.0 (I mean in the default setting when installing R and also I tried to open the --enable-R-shlib. I chose installing full binary, kernel src, & X during the instal

[R] Parallel Processing Interface for lapply()

2003-02-10 Thread John Maindonald
In co-operation with Markus Hegland and myself, Zhongwen Ding has written a package, based partly on Markus Hegland's code, that provides a parallel processing interface to a remote multi-processsor system. Pyro (Python Remote Objects) and R must both be installed, both on the client machine and o

Re: [R] geoR question from new R user

2003-02-10 Thread Roger Peng
Is `D' the original matrix or a `geodata' object? You need to assign the output of `as.geodata to another variable and pass that to `variog'. For example, B <- as.geodata(D, coords.col=1:2, data.col=3) variog(B, option = "cloud", max.dist = 1) It seems like the problem is that you're passing th

[R] geoR question from new R user

2003-02-10 Thread Sharon Lambert
Hi, I'm a new R user. My goal is to do a variogram using geoR. I started by trying to do the example in the geoR Illustrative Session using my own data. I am able to read in my Ascii data using: D <- matrix(scan("file.dat", n=530*3), 530,3, byrow=TRUE). Then I use: as.geodata(D, coords.col=1:2

Re: [R] multilm for simseg/acm

2003-02-10 Thread Prof Brian D Ripley
On Mon, 10 Feb 2003, Christian Schulz wrote: > hi, > for working with the simseg/acm approach i need multilm and it seems > that the last version is 0.1-4.tar ? > What's wrong in the description file, if i attempt install the package ? > > c:\DataMining\rw1062\bin>Rcmd Install multilm > Malforme

[R] multilm for simseg/acm

2003-02-10 Thread Christian Schulz
hi, for working with the simseg/acm approach i need multilm and it seems that the last version is 0.1-4.tar ? What's wrong in the description file, if i attempt install the package ? c:\DataMining\rw1062\bin>Rcmd Install multilm Malformed DCF file (file multilm/DESCRIPTION, line 6) [R1.6.2 /W2

Re: [R] Type of multi-valued variable

2003-02-10 Thread Frank E Harrell Jr
On Mon, 10 Feb 2003 21:51:50 +0100 Fan <[EMAIL PROTECTED]> wrote: > Hi, > > I've read in the past a thead in the R discussion list > about the multi-valued type variable (what was called checklist). > At the moment Gregory had intention to add some general code > in his gregmisc package. > > I'm

RE: [R] non-SQL sqlQuery error

2003-02-10 Thread Pijus Virketis
Dear All, Reverting to RODBC 0.9.1 that I had archived resolved the problem. Thanks to Professor Ripley for noting that I had a non-release version of RODBC (how I might have obtained it is still beyond me ;)). Pijus __ [EMAIL PROTECTED] mailing list

[R] Type of multi-valued variable

2003-02-10 Thread Fan
Hi, I've read in the past a thead in the R discussion list about the multi-valued type variable (what was called checklist). At the moment Gregory had intention to add some general code in his gregmisc package. I'm wondering if there's some general code / packages available ? A general class for

RE: [R] non-SQL sqlQuery error

2003-02-10 Thread Pijus Virketis
Professor Ripley, > There is no statement 3 in the released version of sqlGetResults ... In other words, I have an outdated/wrong version of RODBC? > You can. R are RODBC and Open Source, so why not open it and read it? Because I don't know C/C++. I realise I should learn it, but I hope that t

Re: [R] non-SQL sqlQuery error

2003-02-10 Thread ripley
On Mon, 10 Feb 2003, Pijus Virketis wrote: > I've encountered a curious problem. I am trying to run an SQL query > using sqlQuery() function in RODBC. > The query works fine when run in a stand-alone SQL browser (Microsoft > Query Analyzer, in particular). > However, when I use the exact same th

[R] non-SQL sqlQuery error

2003-02-10 Thread Pijus Virketis
Dear all, I've encountered a curious problem. I am trying to run an SQL query using sqlQuery() function in RODBC. The query works fine when run in a stand-alone SQL browser (Microsoft Query Analyzer, in particular). However, when I use the exact same thing from sqlQuery() function, I get the fo

[R] hierarchical partitioning in multiple regression

2003-02-10 Thread David Nogués
Hello to everybody: Im starting to use R and this is my first message to the list. My first question is: Can R-Program calculate hierarchical partitioning analysis to measure the "independent effect" of each independent variable in a multiple regression (lm, glm) Thanks in advance -- David Nogué

Re: [R] problems using lqs()

2003-02-10 Thread Luca Scrucca
Many thanks to Kjetil Halvorsen for explaining me the reason I got different results. Regarding the second point, Both Kjetil and Brian Ripley suggest to go through debugging and see what the reason for the message. I'm not so expert to fully understand why, but I have found a patch that makes the

Re: [R] Memory allocation, IBM-AIX and R-1.6.2 - addendum

2003-02-10 Thread Peter Dalgaard BSA
Laurent Gautier <[EMAIL PROTECTED]> writes: > ...sorry for the spam, but answers I get to my previous > question suggest that I should specify that the > machine has a *lot* of memory and should be able > to the instanciation... > > Anybody with an IBM mainframe, R-1.6.2 and large > matrices ? F

Re: [R] Zero rows/cols in the hessian matrix

2003-02-10 Thread Timur Elzhov
On Mon, Feb 10, 2003 at 05:25:33PM +0100, Ott Toomet wrote: > ind <- c(1, 3) # indices you need > varcov <- matrix(0, 4, 4) > varcovar[ind,ind] <- solve(hessian[ind,ind]) > > I.e. invert only parameter-depending part of the hessian and put it > into the corresponding elements in the full varcovar

Re: [R] lda

2003-02-10 Thread ripley
R per se does not have lda() function. Package MASS does, and MASS (the book) describes it in detail. If you use a package supporting a book (there are several) do expect to read the book for the fine details On Mon, 10 Feb 2003, Luis Silva wrote: > There are some versions of lda. I would

Re: [R] Zero rows/cols in the hessian matrix

2003-02-10 Thread Ott Toomet
Hi, | From: Timur Elzhov <[EMAIL PROTECTED]> | Date: Mon, 10 Feb 2003 19:06:18 +0300 | | Dear R experts! | | I try to minimize a function with external C fitting function. | I get the hessian matrix. Here it is: | | [,1] [,2] [,3] [,4] | [1,] 1.88166310 0.885980

[R] (no subject)

2003-02-10 Thread Tomer Maymon
my name is Tomer Maimon, and i'm a student in the israeli Institue for science ( "Technion") .i have some problem with the R program. i want to write formula with two factors to predict(in rpart algorithm). i'll be more then thankfull if you could send me an answer & internet link, where can i find

Re: [R] problems using lqs()

2003-02-10 Thread kjetil brinchmann halvorsen
On 10 Feb 2003 at 15:40, Luca Scrucca wrote: I am not sure about this, but lqs have an extra argument `adjust'. from the help page: adjust: should the intercept be optimized for each sample? Defaults to TRUE When you use formula, the estimation algorithm knows the columns of 1's is the intercep

Re: [R] Factor level comparisons in lme

2003-02-10 Thread Douglas Bates
Martin Biuw <[EMAIL PROTECTED]> writes: > Hello, > I'm trying to fit a linear mixed effects model of the form: > > lme(y ~ x * Sex * Year, random=x|subject) Did you mean lme(y ~ x * Sex * Year, random= ~ x|subject) The random argument should be a formula or a list. > where Sex and Year are f

[R] lda

2003-02-10 Thread Luis Silva
Dear helpers There are some versions of lda. I would like to know which one is behind R's lda function. Luis -- Ofereça a si uma prenda no Dia dos Namorados: SAPO.ADSL.PT, e receba ainda uma oferta especial para a sua cara-metade. Tudo só por 50€ http://adsl.sapo.pt/destaque_namorados.html __

[R] Zero rows/cols in the hessian matrix

2003-02-10 Thread Timur Elzhov
Dear R experts! I try to minimize a function with external C fitting function. I get the hessian matrix. Here it is: [,1] [,2] [,3] [,4] [1,] 1.88166310 0.88598030 [2,] 0.0000 0.0000 [3,] 0.88598030 0.48599830 [4,] 0.0000 0.000

Re: [R] problems using lqs()

2003-02-10 Thread ripley
On Mon, 10 Feb 2003, Luca Scrucca wrote: > I found a strange behaviour in the lqs function. Why is this strange? It *is* covered in the Note on the help page: please read it and the reference there. > Suppose I have the following data: [...] > Now, instead of using the formula, I want to prov

Re: [R] shapiro.test

2003-02-10 Thread Richard A. Bilonick
Ernesto Jardim wrote: Hi It's a bootstrap empirical distribution and has 1000 replicates. EJ On Mon, 2003-02-10 at 15:33, Richard A. Bilonick wrote: It depends on how non-Normal the distribution and the size of the sample. A t-distribution with df = 30 isn't Normal but it is close to bein

Re: [R] shapiro.test

2003-02-10 Thread Richard A. Bilonick
Ernesto Jardim wrote: Ok, let me put it the other way around. On another test I have W = 0.9907, p-value = 6.024e-06. The same question stands, with such huge W should it be expected to be normal ? EJ You have it backwards. The null hypothesis is that the distribution is Normal. You reject

Re: [R] shapiro.test

2003-02-10 Thread Ernesto Jardim
Ok, let me put it the other way around. On another test I have W = 0.9907, p-value = 6.024e-06. The same question stands, with such huge W should it be expected to be normal ? EJ On Mon, 2003-02-10 at 14:47, Richard A. Bilonick wrote: > Ernesto Jardim wrote: > > >Hi > > > >The shapiro.test func

[R] Memory allocation, IBM-AIX and R-1.6.2 - addendum

2003-02-10 Thread Laurent Gautier
...sorry for the spam, but answers I get to my previous question suggest that I should specify that the machine has a *lot* of memory and should be able to the instanciation... Anybody with an IBM mainframe, R-1.6.2 and large matrices ? L. __ [EMAIL PR

[R] Importing e-prime data

2003-02-10 Thread Hedderik van Rijn
Before I'm reinventing the wheel: did anyone already work on an E-Prime (experiment presentation and data collection tool used in Psychology, http://www.pstnet.com/e-prime/default.htm) data import function? - Hedderik. (Googling with "e-prime site:r-project.org" didn't result in any hits.)

[R] problems using lqs()

2003-02-10 Thread Luca Scrucca
Dear List-members, I found a strange behaviour in the lqs function. Suppose I have the following data: y <- c(7.6, 7.7, 4.3, 5.9, 5.0, 6.5, 8.3, 8.2, 13.2, 12.6, 10.4, 10.8, 13.1, 12.3, 10.4, 10.5, 7.7, 9.5, 12.0, 12.6, 13.6, 14.1, 13.5, 11.5, 12.0, 13.0, 14.1, 15.1) x1 <- c(8.2, 7.6,, 4.6, 4.3,

Re: [R] shapiro.test

2003-02-10 Thread Richard A. Bilonick
Ernesto Jardim wrote: Hi The shapiro.test function outputs a value of the W statistic, which should be 1 if the distribution is normal, and a p-value for the test (as the documentation states). I'm a bit confused with some results. I'm getting a W=0.9977 and a p-value=0.1889. I was expecting t

[R] shapiro.test

2003-02-10 Thread Ernesto Jardim
Hi The shapiro.test function outputs a value of the W statistic, which should be 1 if the distribution is normal, and a p-value for the test (as the documentation states). I'm a bit confused with some results. I'm getting a W=0.9977 and a p-value=0.1889. I was expecting that a W of 0.9977 would

[R] DSC 2003: Workshop on Distributed Statistical Computing

2003-02-10 Thread Achim Zeileis
Dear R users, the preliminary program for the 3rd international workshop on 'Distributed Statistical Computing' (DSC 2003) is now available at http://www.ci.tuwien.ac.at/Conferences/DSC-2003/ The workshop will take place at the Technische Universität Wien in Vienna, Austria from March 20-22.

Re: [R] nls

2003-02-10 Thread John Fox
Dear Gregory, As Brian Ripley pointed out (and as the help page for nls mentions), you shouldn't test nls on zero-residual data. But, in any event, you can fit this model with linear least squares, as lm(y ~ I(1/x) - 1). I hope that this helps, John At 10:11 AM 2/10/2003 +0100, Gregory BENMEN

[R] Factor level comparisons in lme

2003-02-10 Thread Martin Biuw
Hello, I'm trying to fit a linear mixed effects model of the form: lme(y ~ x * Sex * Year, random=x|subject) where Sex and Year are factors with two and three levels respectively. I want to compare the fixed effects for each level to the overall mean, but the default in R is to compare to the f

[R] hashmaps/hashtables and R

2003-02-10 Thread Steffen Durinck
Dear, Do hashmaps or hashtables exist in R? If so how can I implement them? Kind Regards, Steffen __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] Memory allocation, IBM-AIX and R-1.6.2

2003-02-10 Thread Laurent Gautier
Dear all, I compiled R-1.6.2 for IBM-AIX (using the native compilers) and I am facing problems to instanciate (rather) large matrices. I have: > m <- matrix(0, 640*640, 102) Error: cannot allocate vector of size 326400 Kb I am not truly familiar with AIX, but this does not seem t

Re: [R] nls

2003-02-10 Thread Prof Brian D Ripley
Don't use nls for exact fits, as that's not what it is intended for, and you are trying to achieve a fit to machine accuracy with rounding errors. You will find several discussions of this in the archives. On Mon, 10 Feb 2003, Gregory BENMENZER wrote: > I want to estimate parameters of the model

[R] nls

2003-02-10 Thread Gregory BENMENZER
Hello, I want to estimate parameters of the model y=a/x with x=1:10 and y=3/x. I tested the NLS function. Could you tell me why it doesn't converge ? Regards, Grégory Benmenzer > x=1:10 > y=3/x > nls(y~a/x, start=list(a=2), control=nls.control(maxiter=100)) Error in nls(y ~ a/x, start = list(a