Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-10 Thread Thomas Lumley
On Thu, 5 Apr 2007, John Kane wrote: > Heck. I'd be happy to get an answer to what is > happening here: >> mac <- spss.get("H:/ONTH/Raw.data/Follow.sav") > Warning message: > H:/ONTH/Raw.data/Follow.sav: Unrecognized record type > 7, subtype 16 encountered in system file > It means that your file

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-10 Thread John Kane
--- Thomas Lumley <[EMAIL PROTECTED]> wrote: > On Thu, 5 Apr 2007, John Kane wrote: > > Heck. I'd be happy to get an answer to what is > > happening here: > >> mac <- spss.get("H:/ONTH/Raw.data/Follow.sav") > > Warning message: > > H:/ONTH/Raw.data/Follow.sav: Unrecognized record > type > > 7, su

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-10 Thread John Kane
--- Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > On Thu, 5 Apr 2007, Michael Conklin wrote: > > > Not being the developer I cannot answer > definitively but, as a frequent > > user of SPSS files I can give you my experience. > > > > 1) The unrecognized coding is perhaps due to the > locale of

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-10 Thread Prof Brian Ripley
On Thu, 5 Apr 2007, Michael Conklin wrote: Not being the developer I cannot answer definitively but, as a frequent user of SPSS files I can give you my experience. 1) The unrecognized coding is perhaps due to the locale of the SPSS installation. I have had success reading in files from versio

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-10 Thread Prof Brian Ripley
On Thu, 5 Apr 2007, Michael Conklin wrote: Not being the developer I cannot answer definitively but, as a frequent user of SPSS files I can give you my experience. 1) The unrecognized coding is perhaps due to the locale of the SPSS installation. I have had success reading in files from versio

Re: [R] plot log scale, axis original scale

2007-04-10 Thread Petr PIKAL
Hi Dean Sonneborn <[EMAIL PROTECTED]> napsal dne 10.04.2007 18:28:43: > Petr, > This is great! Thank you so much for responding. Could I get one more > point clarified. My A values range from 1 to 35. I would really like to > use something like > AT=1 to 35 by 5 instead of AT=log(a). at=log(

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-10 Thread Thomas Lumley
On Thu, 5 Apr 2007, John Kane wrote: > Heck. I'd be happy to get an answer to what is > happening here: >> mac <- spss.get("H:/ONTH/Raw.data/Follow.sav") > Warning message: > H:/ONTH/Raw.data/Follow.sav: Unrecognized record type > 7, subtype 16 encountered in system file > It means that your file

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-10 Thread John Kane
--- Thomas Lumley <[EMAIL PROTECTED]> wrote: > On Thu, 5 Apr 2007, John Kane wrote: > > Heck. I'd be happy to get an answer to what is > > happening here: > >> mac <- spss.get("H:/ONTH/Raw.data/Follow.sav") > > Warning message: > > H:/ONTH/Raw.data/Follow.sav: Unrecognized record > type > > 7, su

Re: [R] read.spss (package foreign) and SPSS 15.0 files

2007-04-10 Thread John Kane
--- Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > On Thu, 5 Apr 2007, Michael Conklin wrote: > > > Not being the developer I cannot answer > definitively but, as a frequent > > user of SPSS files I can give you my experience. > > > > 1) The unrecognized coding is perhaps due to the > locale of

Re: [R] Reasons to Use R

2007-04-10 Thread Greg Snow
I think SAS has the database part built into it. I have heard 2nd hand of new statisticians going to work for a company and asking if they have SAS, the reply is "Yes we use SAS for our database, does it do statistics also?" Also I heard something about SAS is no longer considered an acronym, the

Re: [R] Reasons to Use R

2007-04-10 Thread Wensui Liu
Greg, As far as I understand, SAS is more efficient handling large data probably than S+/R. Do you have any idea why? On 4/10/07, Greg Snow <[EMAIL PROTECTED]> wrote: > > -Original Message- > > From: [EMAIL PROTECTED] > > [mailto:[EMAIL PROTECTED] On Behalf Of > > Bi-Info (http://members.h

[R] Using Sampling Weights in R

2007-04-10 Thread Thomas W. Volscho
Dear List, I have a dataset that provides sampling weights (National Survey of Family Growth 2002). I want to produce a cross-tabulation and use the provided sampling weights to obtain representative population estimates. (I believe they are simply frequency weights but codebook is uninformative).

Re: [R] Function knn.dist from knnflex library

2007-04-10 Thread Prof Brian Ripley
On Wed, 11 Apr 2007, [EMAIL PROTECTED] wrote: > Hello, > > I am feeling that this question can have a very simple answer, but I > can't find it. > > I need to use the function knn.dist from knnflex library. > > Whatever I try, I get the error: > Error in as.vector.dist(x, "character") : unused arg

Re: [R] Memory management

2007-04-10 Thread Prof Brian Ripley
On Tue, 10 Apr 2007, yoo wrote: > > Hi all, I'm just curious how memory management works in R... I need to run an > optimization that keeps calling the same function with a large set of > parameters... so then I start to wonder if it's better if I attach the > variables first vs passing them i

Re: [R] Rserve and R to R communication

2007-04-10 Thread AJ Rossini
On Tuesday 10 April 2007 23:17, Ramon Diaz-Uriarte wrote: > Of course, you are right there. I think that might still be the case. > At the time we made our decision, and decided to go for MPI, MPI 2 was > already out, and MPI seemed "more like the current/future standard" > than PVM. That's alwa

Re: [R] crosstabs in R?

2007-04-10 Thread Marc Schwartz
On Wed, 2007-04-11 at 10:18 +0700, A.R. Criswell wrote: > Hello All, > > I looked through the archives and I've seen this question asked > before, but with no answer: Is there an R equivalent to crosstabs() in > S-Plus? > > Andrew Hey Andrew, The closest is probably the CrossTable() function in

[R] crosstabs in R?

2007-04-10 Thread A.R. Criswell
Hello All, I looked through the archives and I've seen this question asked before, but with no answer: Is there an R equivalent to crosstabs() in S-Plus? Andrew __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEAS

[R] Function knn.dist from knnflex library

2007-04-10 Thread jmoreira
Hello, I am feeling that this question can have a very simple answer, but I can't find it. I need to use the function knn.dist from knnflex library. Whatever I try, I get the error: Error in as.vector.dist(x, "character") : unused argument(s) ("character") First example: > a<-NULL > a<-rbind

Re: [R] problems understanding sapply

2007-04-10 Thread Marc Schwartz
On Tue, 2007-04-10 at 21:19 -0400, John Sorkin wrote: > R 2.4.1 > windows XP > > I am trying to fill in a matrix with binomial probabilities without > using a for loop. I am trying to obtain a value for pbinom using the > value stored in column one of the matrix delete. Clearly I am doing > somet

Re: [R] problems understanding sapply

2007-04-10 Thread Floris Van Ogtrop
Dear John, Try the following? delete<-matrix(nrow=31,ncol=2) delete[1:31,1]<-1:31 f<- function(X) { pbinom(X, 30, 0) } delete[,2]<-sapply(delete[,1], f) delete HTH Floris University of Sydney -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of

[R] problems understanding sapply

2007-04-10 Thread John Sorkin
R 2.4.1 windows XP I am trying to fill in a matrix with binomial probabilities without using a for loop. I am trying to obtain a value for pbinom using the value stored in column one of the matrix delete. Clearly I am doing something wrong. Please help me understand my error. > delete<-matri

Re: [R] Positioning in xyplot

2007-04-10 Thread Michael Kubovy
On Apr 10, 2007, at 7:07 PM, Deepayan Sarkar wrote: > On 4/10/07, Michael Kubovy <[EMAIL PROTECTED]> wrote: >> Dear R-helpers, >> >> I have an xyplot >> >> of the following kind: >> xYplot(y ~ x | p, groups = factor(gg, levels = c('1', '5', '2', '4', >> '3')), >> strip = strip.custom(var.

Re: [R] Reasons to Use R

2007-04-10 Thread Frank E Harrell Jr
Taylor, Z Todd wrote: > On Monday, April 09, 2007 3:23 PM, someone named Wilfred wrote: > >> So what's the big deal about S using files instead of memory >> like R. I don't get the point. Isn't there enough swap space >> for S? (Who cares anyway: it works, isn't it?) Or are there >> any problems w

[R] p-value and significance code

2007-04-10 Thread Paul Suckling
Dear R list members. Is there anyone out there that might be able to explain the difference between p-values and significance codes for terms in a linear regression model (as given by summary.lm for example) for me please? Do they both give the probability that the observed variation in the depen

[R] tapply and storage

2007-04-10 Thread Floris Van Ogtrop
Dear R-Users, I have the following problem of which I have provided a simple example. Using the tapply command I can efficiently run the function genflo for all months and years. I am new to R and I do not understand how I can store the results of f such that as the function loops through the mont

Re: [R] Reasons to Use R

2007-04-10 Thread Andrew Robinson
Hi Todd, I guess I don't see the difference between that strategy and using make to look after scripts, raw data, Sweave files, and (if necessary) images. I find that I can get pretty fine-grained control over what parts of a project need to be rerun by breaking the analysis into chapters. I sup

Re: [R] Reasons to Use R

2007-04-10 Thread Taylor, Z Todd
On Monday, April 09, 2007 3:23 PM, someone named Wilfred wrote: > So what's the big deal about S using files instead of memory > like R. I don't get the point. Isn't there enough swap space > for S? (Who cares anyway: it works, isn't it?) Or are there > any problems with S and large datasets? I do

Re: [R] Positioning in xyplot

2007-04-10 Thread Deepayan Sarkar
On 4/10/07, Michael Kubovy <[EMAIL PROTECTED]> wrote: > Dear R-helpers, > > I have an xyplot > > of the following kind: > xYplot(y ~ x | p, groups = factor(gg, levels = c('1', '5', '2', '4', > '3')), > strip = strip.custom(var.name = 'g', factor.levels = c(1', '5', '2', > '4', '3'), >

Re: [R] Reasons to Use R

2007-04-10 Thread Gabor Grothendieck
I think SAS was developed at a time when computer memory was much smaller than it is now and the legacy of that is its better usage of computer resources. On 4/10/07, Wensui Liu <[EMAIL PROTECTED]> wrote: > Greg, > As far as I understand, SAS is more efficient handling large data > probably than S

Re: [R] sas.get problem

2007-04-10 Thread Tim Churches
John Kane wrote: > I have 3 SAS files all in the directory F:/sas, two > data files > and a format file : > form.ea1.sas7bdat > form.ea2.sas7bdat > sas.fmts.sas7bdat > > F is a USB. > > I am trying import them to R using "sas.get". > > I have not used SAS since I was downloading data from > ma

Re: [R] Positioning in xyplot

2007-04-10 Thread Sundar Dorai-Raj
Michael Kubovy said the following on 4/10/2007 3:21 PM: > Dear R-helpers, > > I have an xyplot > > of the following kind: > xYplot(y ~ x | p, groups = factor(gg, levels = c('1', '5', '2', '4', > '3')), > strip = strip.custom(var.name = 'g', factor.levels = c(1', '5', '2', > '4', '3'),

[R] Positioning in xyplot

2007-04-10 Thread Michael Kubovy
Dear R-helpers, I have an xyplot of the following kind: xYplot(y ~ x | p, groups = factor(gg, levels = c('1', '5', '2', '4', '3')), strip = strip.custom(var.name = 'g', factor.levels = c(1', '5', '2', '4', '3'), strip.levels = T, strip.names = T, sep = ' = ', shingle.

Re: [R] sas.get problem

2007-04-10 Thread Cody_Hamilton
John, According to the sas.get documentation, the format library option must specify the directory containing the file formats.sct (I don't think it can be a sas dataset). Also, do you need a forward slash after C: in the sasprog option below? Finally, sas may not like the period in the file na

[R] Memory management

2007-04-10 Thread yoooooo
Hi all, I'm just curious how memory management works in R... I need to run an optimization that keeps calling the same function with a large set of parameters... so then I start to wonder if it's better if I attach the variables first vs passing them in (coz that involves a lot of copying.. ) Thu

[R] sas.get problem

2007-04-10 Thread John Kane
I have 3 SAS files all in the directory F:/sas, two data files and a format file : form.ea1.sas7bdat form.ea2.sas7bdat sas.fmts.sas7bdat F is a USB. I am trying import them to R using "sas.get". I have not used SAS since I was downloading data from mainframe and having to write JCL. I had for

Re: [R] Rserve and R to R communication

2007-04-10 Thread Ramon Diaz-Uriarte
On 4/10/07, AJ Rossini <[EMAIL PROTECTED]> wrote: > On Monday 09 April 2007 23:02, Ramon Diaz-Uriarte wrote: > > > (Yes, maybe I should check snowFT, but it uses PVM, and I recall a > > while back there was a reason why we decided to go with MPI instead of > > PVM). > > There is no reason that you

[R] Computing fundamental harmonics from a periodogram

2007-04-10 Thread hadley wickham
Dear all, I'm trying to finding the fundamental harmonics (ie. peaks in a periodogram) from a time series (extracted from an mp3). For example, if I look at spectrum(fdeaths, spans = c(3,3)) I'd say the fundamental harmonics are about 1, 2, 3.5 and 4.5 - but how can I extract these automaticall

[R] Kmeans cluster analysis

2007-04-10 Thread nathaniel Grey
Hello, I have a data-set containing 22 variables, after appropriate transformations etc I ran a kmeans cluster analysis for 4 clusters , I ran it 20 times to find a result with the lowest within sum of squares. My question is how best do I go about finding out what the characteristics are

Re: [R] Barplot legend position

2007-04-10 Thread Peter Danenberg
> However, the legend does not reproduce the color/shading used in the > original barplot, are those available somehow? Actually, Ingmar, there's a more elegant way to recre- ate the original colors; to expand on your example: data <- 1:10 rows <- 2 cols <- 5 labels <-

Re: [R] sem vs. LISREL: sem fails

2007-04-10 Thread John Fox
Dear Dimitri, Using the tools provided in the sem package [e.g., the debug argument to sem()], I took a closer look at this convergence problem, which arises because of poor start values selected for the variances of the latent variables. It turns out that you can get a solution for your original

Re: [R] Reasons to Use R

2007-04-10 Thread Greg Snow
> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of > Bi-Info (http://members.home.nl/bi-info) > Sent: Monday, April 09, 2007 4:23 PM > To: Gabor Grothendieck > Cc: Lorenzo Isella; r-help@stat.math.ethz.ch > Subject: Re: [R] Reasons to Use R [snip] >

Re: [R] Problem installing gnomeGUI in Ubuntu: "HAVE_ORBIT" was never defined

2007-04-10 Thread fsando
Thanks! It worked after a some trial and error: I Installed the pmg package through R. When loading pmg it complained about missing gWidgets which I installed through R. Then it complained about missing guitoolkit I guessed that would be gWidgetsRGtk2 (turned out to be correct). Installation fai

Re: [R] Rserve and R to R communication

2007-04-10 Thread AJ Rossini
On Monday 09 April 2007 23:02, Ramon Diaz-Uriarte wrote: > (Yes, maybe I should check snowFT, but it uses PVM, and I recall a > while back there was a reason why we decided to go with MPI instead of > PVM). There is no reason that you can't run both MPI and PVM on the same cluster. There is a pa

Re: [R] How do I back transforme ordinary log-krigged prediction values?

2007-04-10 Thread Zia Uddin Ahmed
Farrar, I did ordinary kriging (OK) after log10 transformation. I need to backtransfer my predicted data to original unit. You are right! Anti-log of the predictions for bcaktransformation might give biased values. For minimize the biasness, it is necessary know the lagrange-multiplier (y) and pre

Re: [R] How do I back transforme ordinary log-krigged prediction values?

2007-04-10 Thread Zia Uddin Ahmed
Farrar, I did ordinary kriging (OK) after log10 transformation. I need to backtransfer my predicted data to original unit. You are right! Anti-log of the predictions for bcaktransformation might give biased values. For minimize the biasness, it is necessary know the lagrange-multiplier (y) and pre

[R] Clustered vectors

2007-04-10 Thread RK
Dear R users, Just wondering how to extract clustered data in a file. More precisely after running any clustering algorithms, you get information about the clusterd vectors. I want to save that data after clustering is done in the same order as the clusters are produced. How can I do that After run

Re: [R] How do I back transforme ordinary log-krigged prdiction values?

2007-04-10 Thread David Farrar
Zia, I'm not an expert on Kriging, or on the particular package. It sounds like you just want to know how to anti-log the predictions,. which shouldn't be hard. However, depending on what you actually need to predict, it seems your predictions might be biased, right? regards, Farrar

Re: [R] nlm() and optim()

2007-04-10 Thread Bert Gunter
Numerical optimization is sensitive to (at least) the method chosen, control/convergence specifications, and the parameterization of the function being optimized (all of this is well known). Defining what you mean by "reproduce" in a precise, operational way is therefore essential. You have not don

[R] Lagrange multiplier for ordinary kriging in gstat in R??

2007-04-10 Thread Edzer J. Pebesma
Zia, this is not trivial, as gstat uses the kriging solution where Lagrange parameters are not involved. You can get the separate components of it and work from there. Another angle would be the difference between the ordinary and simple kriging variance, but I'm not 100% sure that equals the La

[R] nlm() and optim()

2007-04-10 Thread Silvia Lucato
Dear R-users, I have just joint the list and much appreciate any thoughts on 2 issues. Firstly, I want to reproduce some minimization results conducted in MATLAB. I have suceeded with nlm and optim-method CG. I have been told that I should get also with other optim methods. Actually, I found th

[R] Estimate new values from GLM (familiy=Gamma)

2007-04-10 Thread Antonio Punzon
Hi, we have this model glm(CPUE~Year+Month+Length+Year:Month, family=Gamma (link="inverse")) where, CPUE is a numerical variable >0, We use this model for standardization of CPUE and we need the new estimates for each year. How can we obtain it? Thank you ___

Re: [R] Reasons to Use R

2007-04-10 Thread Greg Snow
For a previous version of SAS we had parts installed on each computer where it was used, but there were key pieces located on a network drive (not internet, but local network) such that if you tried to start SAS while someone else was using it you would get an error message. We had troubles with t

Re: [R] Matlab import

2007-04-10 Thread Spencer Graves
Have you considered exporting from Matlab using the '-ascii' option, then reading it with 'read.table', after perhaps checking the first few rows using, e.g., 'readLines', and confirming that all records have the same length using 'count.fields'? Hope this helps. Spencer Gra

[R] How to do spatial-temporal regression

2007-04-10 Thread 董莉娟
hello all, If I want to do regression of spatial-temporal data, what package or function should I use? Thanks. ¶­Àò¾ê 2007-04-10 [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mai

Re: [R] Problem installing gnomeGUI in Ubuntu: "HAVE_ORBIT" was never defined

2007-04-10 Thread Michael Lawrence
If you want a GTK2-based R GUI, try the pmg package. On 4/10/07, fsando <[EMAIL PROTECTED]> wrote: > > Thanks a lot - Yes I have gnome 2.18, that must be it. > So no gnomeGUI for me :( > > Michael Lawrence wrote: > > My guess is that you have Gnome >=2.0, while the gnomeGUI only works > > with the

Re: [R] Computational speed question

2007-04-10 Thread Prof Brian Ripley
On Tue, 10 Apr 2007, David Vonka wrote: > Hello, > > I'm doing bootstrap in R 2.4.1 in order to compute standard errors of an > estimator. I'm running the same program on 4 computers. The core of the > program is a nlm minimization of a function, which is (I believe) nicely > coded using (t)apply

[R] Computational speed question

2007-04-10 Thread David Vonka
Hello, I'm doing bootstrap in R 2.4.1 in order to compute standard errors of an estimator. I'm running the same program on 4 computers. The core of the program is a nlm minimization of a function, which is (I believe) nicely coded using (t)apply and all the vectorized stuff. I'm slightly puzzl

Re: [R] Matlab import

2007-04-10 Thread Schmitt, Corinna
Hallo, >With readMat, don't use the "-ascii" option (which you didn't have in your >first posting). I've never tried reading matlab's ascii format. In any >case, >readMat reads matlab's binary format. >- Tom I did the saving again without 'ascii' option but the import also did not work. I get

Re: [R] sem vs. LISREL: sem fails

2007-04-10 Thread John Fox
Dear Dimitri, > -Original Message- > From: John Smith [mailto:[EMAIL PROTECTED] > Sent: Tuesday, April 10, 2007 8:11 AM > To: John Fox > Cc: r-help@stat.math.ethz.ch > Subject: Re: [R] sem vs. LISREL: sem fails > > Dear Scott, (John, actually.) > thanks a lot - all your comments are ex

Re: [R] Matlab import

2007-04-10 Thread tshort
With readMat, don't use the "-ascii" option (which you didn't have in your first posting). I've never tried reading matlab's ascii format. In any case, readMat reads matlab's binary format. - Tom Tom Short EPRI Schmitt, Corinna wrote: > > Hallo, > >>I've used Henrik Bengtsson's R.matlab pac

Re: [R] Reasons to Use R

2007-04-10 Thread Jeffrey J. Hallman
"halldor bjornsson" <[EMAIL PROTECTED]> writes: > ... > Now, R does not have everything we want. One thing missing is a decent > R-DB2 connection, for windows the excellent RODBC works fine, but ODBC > support on Linux is a hassle. > A hassle? I use RODBC on Linux to read data from a mainf

Re: [R] Matlab import

2007-04-10 Thread Schmitt, Corinna
Hallo, >I've used Henrik Bengtsson's R.matlab package several times to successfully >read in matlab data files. It's normally as easy as: >library(R.matlab) >mats <- readMat("matrixM.txt") >- Tom I have imported this package, too. And tried your commands with the new txt-file as mentioned in m

[R] WG: Matlab import

2007-04-10 Thread Schmitt, Corinna
Dear R-Experts, >here again a question concerning matlab. With the command "matrixM=[1 2 >3;4 5 6]" a matrix under Matlab was constructed. It was than stored with >the command "save('matrixM.txt','matrixM')". now I found my error in the matlab command for storing the data in a txt file. You need

Re: [R] sem vs. LISREL: sem fails

2007-04-10 Thread John Smith
Dear Scott, thanks a lot - all your comments are extremely valuable! One quick follow-up question - the thing is, I was really trying to compare sem and LISREL solutions, which is why I was trying to specify exactly the same model. Hence, I'd prefer to stick to the Var-Covar matrix, rather than

Re: [R] Matlab import

2007-04-10 Thread tshort
I've used Henrik Bengtsson's R.matlab package several times to successfully read in matlab data files. It's normally as easy as: library(R.matlab) mats <- readMat("matrixM.txt") - Tom Tom Short EPRI Schmitt, Corinna wrote: > > Dear R-Experts, > > here again a question concerning matlab. Wi

[R] string edit distance

2007-04-10 Thread David Meyer
It's in package cba (sdists()). David -- I have a column of words, for example "DOG" "DOOG" "GOD" "GOOD" "DOOR" ... and I am interested in creating a matrix that contains the string edit distances between each pair of words. I am this close -> ' ' <- to writing the algorit

[R] Matlab import

2007-04-10 Thread Schmitt, Corinna
Dear R-Experts, here again a question concerning matlab. With the command "matrixM=[1 2 3;4 5 6]" a matrix under Matlab was constructed. It was than stored with the command "save('matrixM.txt','matrixM')". Now I tried to import the data in R with the help of the command "Z=matrix(scan("Z:/Softwar

Re: [R] Problem installing gnomeGUI in Ubuntu: "HAVE_ORBIT" was never defined

2007-04-10 Thread fsando
Thanks a lot - Yes I have gnome 2.18, that must be it. So no gnomeGUI for me :( Michael Lawrence wrote: > My guess is that you have Gnome >=2.0, while the gnomeGUI only works > with the ancient Gnome 1.4, although I am not that familiar with the > package. > > Michael > > On 4/9/07, * fsando*

Re: [R] Modified Sims test

2007-04-10 Thread Millo Giovanni
Dear Chris, I do not have the references here, but AFAIR: if x and y are two time series, we say that "x does not Granger-cause y" (x ngc y) if the models (1) y~y(-1)+y(-2)+...+x(-1)+x(-2)+... and (2) y~y(-1)+y(-2)+... are "equivalent", i.e. if past values of x "do not help explaining" y. Th

[R] Alerte de l'antivirus de l'IPC: file85.eml\raglans.png Malformed attachment blocked

2007-04-10 Thread SV-IPC-8
L'anti-virus de l'Institut Paoli Calmettes a détecté un virus dans ce message. Time: 10 Apr 2007 12:27:17 Scan result: Mail modified to remove malicious content Protocol: SMTP in File Name\Mail Subject: R-help Digest, Vol 50, Issue 10 Source: 192.168.4.13 Destination: 192.168.3.5 Mail Sender: r-

[R] Multiple output from cor.test

2007-04-10 Thread Joseph Wakeling
Hello all, I'm trying to use cor.test across a range of values via the apply command, c <- apply(X,1,cor.test,y) where X is an array, giving me multiple output from cor.test. I can access individual elements of this by e.g. c[[1]] or c[[1]]$estimate or c[[n]] or whatever, but is there a way to

Re: [R] &500 and above ignored

2007-04-10 Thread Jim Lemon
Constant Depièreux wrote: > Hello, > > I have a data crunching program which causes me a problem. > > Some of the data are identified with an & sign followed by a integer. > No problem is experienced when analysing data below &500, but value > above are ignored. > > Is there some limitation

Re: [R] &500 and above ignored

2007-04-10 Thread jim holtman
Read the data in as a string and then remove the "&" with 'gsub'. What does your data file look like and how big is it? It would be useful to provide a sample of your data so we could answer the question easier. On 4/10/07, Constant Depièreux <[EMAIL PROTECTED]> wrote: > Hello, > > I have a data

[R] TXTSETUP.OEM and CD-ROM?

2007-04-10 Thread Keisha
Content-Type: multipart/related; boundary="DB3A3FF1.D25D383D" --DB3A3FF1.D25D383D Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable We awaken in others the same attitude of mind we hold toward them=2E An unattempted lady cou

[R] Question regarding the mvtnorm R package

2007-04-10 Thread shlomi lifshits
Hi, I am trying to compare two methods for the calculation of a trivariate normal cdf: 1) "direct" - using the original covariance matrix 2) "indirect" - transforming the integration interval using the Mahalanobis transformation and then integrating as a standard random variable. The followin

Re: [R] When to use quasipoisson instead of poisson family

2007-04-10 Thread Prof Brian Ripley
On Tue, 10 Apr 2007, ronggui wrote: > On 4/10/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: >> On Tue, 10 Apr 2007, Achim Zeileis wrote: >> >> > On Tue, 10 Apr 2007, ronggui wrote: >> > >> >> It seems that MASS suggest to judge on the basis of >> >> sum(residuals(mode,type="pearson"))/df.resid

Re: [R] &500 and above ignored

2007-04-10 Thread Ted Harding
On 10-Apr-07 07:48:38, Constant Depièreux wrote: > Hello, > > I have a data crunching program which causes me a problem. > > Some of the data are identified with an & sign followed by a integer. > No problem is experienced when analysing data below &500, but value > above are ignored. > > Is

Re: [R] When to use quasipoisson instead of poisson family

2007-04-10 Thread ronggui
On 4/10/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > On Tue, 10 Apr 2007, Achim Zeileis wrote: > > > On Tue, 10 Apr 2007, ronggui wrote: > > > >> It seems that MASS suggest to judge on the basis of > >> sum(residuals(mode,type="pearson"))/df.residual(mode). > > Not really; that is the convent

Re: [R] plot log scale, axis original scale

2007-04-10 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 09.04.2007 22:10:22: > I want to produce some boxplots and plot the logged values but have the > axis scale in the original, not-logged scale. It seeming like I have the > first few steps but I'm having trouble with the last. Here's what I'm > doing (which I go

Re: [R] tapply

2007-04-10 Thread Petr PIKAL
Hallo Seems to me that you can make a summary table using aggregate(RESPONSE, list(TREATMENT, MEASUREMENT, BLOCK, STUDY), mean) and then if you want you can use reshape function or melt/cast function from reshape package to get wide form of your table. Regards Petr Pikal [EMAIL PROTECTED] [E

Re: [R] When to use quasipoisson instead of poisson family

2007-04-10 Thread Prof Brian Ripley
On Tue, 10 Apr 2007, Achim Zeileis wrote: > On Tue, 10 Apr 2007, ronggui wrote: > >> It seems that MASS suggest to judge on the basis of >> sum(residuals(mode,type="pearson"))/df.residual(mode). Not really; that is the conventional moment estimator of overdispersion and it does not suffer from th

[R] &500 and above ignored

2007-04-10 Thread Constant Depièreux
Hello, I have a data crunching program which causes me a problem. Some of the data are identified with an & sign followed by a integer. No problem is experienced when analysing data below &500, but value above are ignored. Is there some limitation in R on this topic? Many thanks for your he

Re: [R] When to use quasipoisson instead of poisson family

2007-04-10 Thread Achim Zeileis
On Tue, 10 Apr 2007, ronggui wrote: > It seems that MASS suggest to judge on the basis of > sum(residuals(mode,type="pearson"))/df.residual(mode). My question: Is > there any rule of thumb of the cutpoiont value? > > The paper "On the Use of Corrections for Overdispersion" suggests > overdispersi

Re: [R] median polishing

2007-04-10 Thread Martin Maechler
> "Sorn" == Sorn Norng <[EMAIL PROTECTED]> > on Tue, 10 Apr 2007 16:11:38 +1000 writes: Sorn> Hi, In SPlus there is a function called "twoway" for Sorn> median polishing gridded data. Is there an equivalent Sorn> function in R? yes Sorn> I have been searching for it in