Re: [R] How to get R to ignore certain values when analyzing a column in a data table ?

Hi On 9 Jul 2006 at 20:01, Daryl Manning wrote: To: r-help@stat.math.ethz.ch From: Daryl Manning [EMAIL PROTECTED] Date sent: Sun, 9 Jul 2006 20:01:10 -0700 Subject:[R] How to get R to ignore certain values when analyzing a

Re: [R] Query : Modification of graphs

Hi On 12 Jul 2006 at 16:56, priti desai wrote: Date sent: Wed, 12 Jul 2006 16:56:29 +0530 From: priti desai [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject:[R] Query : Modification of graphs How to change the background

Re: [R] Usage of repeat,while or loop. Documetation?

Hi nobody has answered your question yet, so I try. There is no reproducible code in your mail. Besides you use names of bult-in functions for your data which gives it really messy lookup. Such tasks can be solved by do.call or maybe *apply but without more information, especially a

[R] Scalling/Centering the Data by an Index

Dear All: I would like to center the data in 'x' by 'group'. The following code scale the data and I have not been able to figure out how to change it so I get the centered data. x - c(1, 2, 3, 4, 5, 6, 7, 8) group - c(1,1,1,2,2,2,2,2) unsplit(lapply(split(x,group),scale),group) I would

Re: [R] Scalling/Centering the Data by an Index

If you want to centre the data, but not scale as well, turn off scale: unlist(tapply(x, group, scale, scale=FALSE)) HTH, Simon. Ashraf Chaudhary wrote: Dear All: I would like to center the data in 'x' by 'group'. The following code scale the data and I have not been able to figure out how

Re: [R] shapiro.test() output

-Original Message- From: [EMAIL PROTECTED] [mailto:r-help- [EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Wednesday, July 12, 2006 11:14 PM To: r-help@stat.math.ethz.ch Subject: [R] shapiro.test() output R Users: My question is probably more about elementary

[R] simple question about variables....

Dear R users, I have a simple question on variable manipulation. Imagine I have an object OBJ that has toto as one of its variables. I would like to understand why if I do varname - toto OBJ$varname returns no results whereas OBJ[varname]

Re: [R] simple question about variables....

see ?$ 'x$name' is equivalent to 'x[[name]]' so you need use : eval(parse(text = paste(OBJ$, varname))) --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33

Re: [R] simple question about variables....

Stéphane Cruveiller wrote: Dear R users, I have a simple question on variable manipulation. Imagine I have an object OBJ that has toto as one of its variables. I would like to understand why if I do varname - toto OBJ$varname returns no results whereas

Re: [R] simple question about variables....

Thx for the tip Stéphane. Jacques VESLOT a écrit : see ?$ 'x$name' is equivalent to 'x[[name]]' so you need use : eval(parse(text = paste(OBJ$, varname))) --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du

Re: [R] simple question about variables....

OBJ$toto works to... b - as.data.frame(matrix(1:4,2)) b V1 V2 1 1 3 2 2 4 b$V1 [1] 1 2 but varname is not evaluated in OBJ$varname. --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P.

[R] Matrix binary for Mac OS X

Hello all, I can't seem to find package Matrix for Mac OS X. The R package installer, when pointed to Bristol (UK) does not have Matrix amongst its list. Thank you __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Keep value lables with data frame manipulation

At 13:14 12.07.2006 -0500, Marc Schwartz (via MN) wrote: On Wed, 2006-07-12 at 17:41 +0100, Jol, Arne wrote: Dear R, I import data from spss into a R data.frame. On this rawdata I do some data processing (selection of observations, normalization, recoding of variables etc..). The result is

[R] SAR with count data

I would like to undertake a Spatial autoregression with count data, how can I specify a negative binomial distribution? I am using the lagsarlm function in the spdep package. Thanks Julie -- View this message in context: http://www.nabble.com/SAR-with-count-data-tf1936174.html#a5304944 Sent

[R] package GLLMGibbs

Tried to look for package GLLMGibbs-(package for fitting mixed models by Gibbs samppling contributed by Myles' and Clayton) in CRAN but not found it, Does anybody know if it is still available? thanks Oarabile __ R-help@stat.math.ethz.ch mailing

[R] VERY TINY question: missing function to clear the console?

Hi, for presentation purposes I would like to clear to whole console window (like in a UNIX terminal: 'clear'). Is there such a function? If not, I could image that is not too hard to write such a function. Cheers, Hans __

[R] R --gui=GNOME problem

I've installed GNU R statistical comupting language and environment 2.2.1-2 by Synaptic (on Ubunti Linux 6.06 Dapper Drake). I've installed also r-gnome to have a GUI. I can't start my GUI. Why R --gui=GNOME don't start? It cannot find rgnome. R --gui=tk command instead is good. What's the best

Re: [R] VERY TINY question: missing function to clear the console?

Sorry, I forgot to mention that I'm using R 2.3.1 on Mac OSX 10.4.7 via R.app GUI. On 13 Jul 2006, at 12:53, Joris De Wolf wrote: ctrl+l Thanks, but unfortunately this doesn't work on my Mac. Hans Hi, for presentation purposes I would like to clear to whole console window (like in

Re: [R] VERY TINY question: missing function to clear the console?

On Thu, Jul 13, 2006 at 12:47:09PM +0200, Hans-Joerg Bibiko wrote: Hi, for presentation purposes I would like to clear to whole console window (like in a UNIX terminal: 'clear'). Is there such a function? If not, I could image that is not too hard to write such a function. At the

Re: [R] package GLLMGibbs

On Thu, 2006-07-13 at 11:32 +0100, Oarabile Molaodi wrote: Tried to look for package GLLMGibbs-(package for fitting mixed models by Gibbs samppling contributed by Myles' and Clayton) in CRAN but not found it, Does anybody know if it is still available? thanks Oarabile You will find it

Re: [R] package GLLMGibbs

On Thu, 13 Jul 2006, Oarabile Molaodi wrote: Tried to look for package GLLMGibbs-(package for fitting mixed models by Gibbs samppling contributed by Myles' and Clayton) in CRAN but not found it, Does anybody know if it is still available? I presume you mean GLMMGibbs It is where it

[R] Extracting Phi from gls/lme

I am trying to extract into a scalar the value of Phi from the printed output of gls or lme using corAR1 correlation. ie I want the estimate of the autocorrelation. I can't see how to do this and haven't seen it anywhere in str(model.lme). I can get all the other information - fixed and random

Re: [R] RODBC, missing values, and Excel

Hi, Gabor, Thanks for the code. When I tried this I get an error when trying to use a relative path name: read.excel - function(file, sheet, ...) { require(rcom) require(gdata) oxl - comCreateObject('Excel.Application') comSetProperty(oxl, Visible, TRUE) # this line optional owb

Re: [R] RODBC, missing values, and Excel

Perhaps the Excel API requires an absolute path name. Try: chartr(/, \\, file.path(getwd(), myfile.xls)) where we make sure its using backslashes in case that's also required. On 7/13/06, Sundar Dorai-Raj [EMAIL PROTECTED] wrote: Hi, Gabor, Thanks for the code. When I tried this I get an

[R] Access to conditioning values in xyplot

Dear R-help subscribers, Many thanks for all the answers I received either by mail or through the list, and that were most helpful. For the sake of the list records, I wanted to post the solution, I eventually could obtain, to the problem I have submitted. I am summarizing the

Re: [R] Extracting Phi from gls/lme

John Logsdon [EMAIL PROTECTED] writes: I am trying to extract into a scalar the value of Phi from the printed output of gls or lme using corAR1 correlation. ie I want the estimate of the autocorrelation. I can't see how to do this and haven't seen it anywhere in str(model.lme). I can get

Re: [R] Extracting Phi from gls/lme

Peter This looks very promising: x-mod.lme$model$Struct$corStruct Correlation structure of class corAR1 representing Phi -0.1996813 which is the value I want. Yippee (save the bricks) But: coef(x,unconstrained=FALSE) [1] -0.4048011 and any attempt to coerce x into a scalar always

[R] Problem installing ROracle

Hello I have tried to install ROracle (V. 0.5-7 and 0.5-5) the usual way 1 - set/export variables of Oracle environment (ORACLE_HOME, LD_LIBRARY_PATH, PATH) 2 - R CMD INSTALL ROracle_0.5-7.tar.gz (R Version 2.3.1, patched per 11. July, and other versions) but I get following error: checking for

[R] Problems plotting a function defined as a product

In order to define a function f as: f - function(x) (x+1)*(x+2) I want to use the notation: v = c(1,2) g - function(x) prod((v+x)) That apparently works and, for instance, the loop: for (i in 1:100) { print(f(i)-g(i)) } Produces a sequence of zeros. Nevertheless, if I try to plot the

[R] step method in glm()

Hello, I estimaded two logit models via glm(). A null model (called glm00) and full model with all accessible covariates and interactions between them (glm1). Then I tried to get even better model by step procedure. I tried the following code: step(glm00, scope=formula(glm1), method=both)

Re: [R] Problems plotting a function defined as a product

plot(t, sapply(t,g)) --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33 (0)3.20.87.10.44 Fax : 33 (0)3.20.87.10.31 http://www-good.ibl.fr

Re: [R] Keep value lables with data frame manipulation

Heinz Tuechler wrote: At 13:14 12.07.2006 -0500, Marc Schwartz (via MN) wrote: On Wed, 2006-07-12 at 17:41 +0100, Jol, Arne wrote: Dear R, I import data from spss into a R data.frame. On this rawdata I do some data processing (selection of observations, normalization, recoding of variables

Re: [R] step method in glm()

As the posting guide requested of you: 0) Do not send HTML mail. 1) Tell us your version of R. Some problems in this area have been corrected recently: if it is not at least 2.3.1, update as asked. 2) Give a reproducible example. On Thu, 13 Jul 2006, Birycki Konrad DomBank Warszawa wrote:

[R] ts and stl functions - still a problem

Hi I am still having problems with using the stl function, when I read the csv file into R into a file called tkr and use dim(tkr) the result is 132 1 which is fine. When coerce it into a trime series using ts either: tstkr - ts(t(tkr), deltat=1/12) or tstkr -

[R] AICc vs AIC for model selection

Hi, I am using 'best.arima' function from forecast package to obtain point forecast for a time series data set. The documentation says it utilizes AIC value to select best ARIMA model. But in my case the sample size very small - 26 observations (demand data). Is it the right to use AIC

[R] colors on graph

Hy all, I need to draw something in 2 dimension that has 3 dimension, the choice has been made to use colors to display the third dimension into the graph. Has someone done something like that, i can't figure out how to parametize the colors. Thks for all ideas, COMTE Guillaume

Re: [R] Extracting Phi from gls/lme

John Logsdon [EMAIL PROTECTED] writes: Peter This looks very promising: x-mod.lme$model$Struct$corStruct Correlation structure of class corAR1 representing Phi -0.1996813 which is the value I want. Yippee (save the bricks) But: coef(x,unconstrained=FALSE) [1]

Re: [R] colors on graph

?image Cf. See Also --- Jacques VESLOT CNRS UMR 8090 I.B.L (2ème étage) 1 rue du Professeur Calmette B.P. 245 59019 Lille Cedex Tel : 33 (0)3.20.87.10.44 Fax : 33 (0)3.20.87.10.31 http://www-good.ibl.fr

[R] sem question

Dear all, I am trying to estimate simultaneous equation model concerning growth in russian regions. I run the analysis by means of FIML in R sem package. I am not familiar with SEM yet, but I've just got several suitable estimated specifications. Nevertheless, sometimes R gives the following

Re: [R] Keep value lables with data frame manipulation

At 08:11 13.07.2006 -0500, Frank E Harrell Jr wrote: Heinz Tuechler wrote: At 13:14 12.07.2006 -0500, Marc Schwartz (via MN) wrote: On Wed, 2006-07-12 at 17:41 +0100, Jol, Arne wrote: Dear R, I import data from spss into a R data.frame. On this rawdata I do some data processing (selection of

Re: [R] VERY TINY question: missing function to clear the console?

Hi, many thanks for the hints! Hi, for presentation purposes I would like to clear to whole console window (like in a UNIX terminal: 'clear'). Is there such a function? If not, I could image that is not too hard to write such a function. At the risk of this being a stupid answer: An

[R] Question for LM intercept

Hello, I am having the following silly problem with lm. Let X be a dataframe with X[,1] the dependent variable and X[,-1] the independent variables. I want to run the following but without including an intercept. for(i in 1:100 ){ lm( X[,100-i] ) # this works fine but it returns an

Re: [R] Extracting Phi from gls/lme

Peter Yes spelling would help! I transcribed this from one screen to another and still the penny didn't drop. And the transformation works too. getAnywhere() is yet another useful function that I hadn't heard about... Many thanks for saving my house. Best wishes John John Logsdon

Re: [R] Matrix binary for Mac OS X

Just checked the list from Bristol. Matrix 0.995-11 is on there. Regards, Rob On Jul 13, 2006, at 1:47 AM, A.R. Criswell wrote: Hello all, I can't seem to find package Matrix for Mac OS X. The R package installer, when pointed to Bristol (UK) does not have Matrix amongst its list. Thank

Re: [R] colors on graph

What's not mentioned in ?image See Also is levelplot and friends from library(lattice) . I'm using it myself at the moment instead of image for some reason I forgot. Worth a try. best Niels On Thu, 13 Jul 2006, Jacques VESLOT wrote: ?image Cf. See Also

Re: [R] colors on graph

Is filled.contour() what you are looking for? MHP on 7/13/2006 9:26 AM COMTE Guillaume said the following: Hy all, I need to draw something in 2 dimension that has 3 dimension, the choice has been made to use colors to display the third dimension into the graph. Has someone done

[R] set the bahavior that R deal with missing values?

Dear Rusers, The default behavior in R when performing a regression model with missing values is to exclude any case that contains a missing value? How could i set the bahavior that R deal with missing values? e.g.: exclude cases listwise exclude cases pairwise replace with mean Thanks very much!

[R] how to make a contour plot using data in long format

Hi List, There is a dataset with Latitude Longitude information and corresponding value shaped like this: Lati Longi Value 18.0001 -159.6667 123 18.0023 -159.6665 124 19.1001 -158.2988

Re: [R] R --gui=GNOME problem

I had the same problem, and I fixed following Dirk Eddelbuettel's suggestion on R-sig-Debian some time ago: install also the package gnomeGUI to have it working. HTH, --- Fernando Mayer Grupo de Estudos Pesqueiros - GEP Centro de Ciências Tecnológicas, da Terra e do Mar - CTTMar Universidade do

[R] set the bahavior that R deal with missing values?

Dear Rusers, The default behavior in R when performing a regression model with missing values is to exclude any case that contains a missing value? How could i set the bahavior that R deal with missing values? e.g.: exclude cases listwise exclude cases pairwise replace with mean Thanks very much!

[R] Inspecting a dataframe inside a dataframe

I have a rather peculiar dataframe and I have problems at inspecting it (I'm working with R 2.3.0 on Win XP). Let say that example is a dataframe with two variables: first and second; now while first is a quite common numeric variable, second is a variable which contains a dataframe (this of

[R] Clique- Method-Package

Dear R' Users, Is anyone knows, which package in R'Gui is suitable to use Clique technique in a dataset? Thanks in Advance JJ -- Lecturer J. Joshua Thomas KDU College Penang Campus Research Student, University Sains Malaysia [[alternative HTML version deleted]]

[R] Problem installing local additional packages of Rcmdr (WinXP)

Hi R users: I download the .zip files and also PACKAGES and PACKAGES.gz files in a local directory, but when I call library(Rcmdr) It shows me a tcl/tk windows telling me to install the additional packages, the it takes me to another window where I must specify the path where the .zip files are

Re: [R] set the bahavior that R deal with missing values?

You can deal with the missing values before performing any regression model. if you dataset is : D X1 X2 X3 1 NA 2 4 2 1 NA NA 3 2 2 6 4 NA 34 NA 5 3 NA 7 then the following will replace all NA with the mean value f-function(col){ + col[is.na(col)]-mean(col,na.rm=TRUE) + col}

Re: [R] Rgraphviz: Setting the edge width

Søren Højsgaard [EMAIL PROTECTED] writes: I create an undirected graph with Rgraphviz (see code below). I would like to make the edges thicker. Can anyone help on this?? Questions on Rgraphviz are best directed to the bioconductor mail list. I think you already discovered that, but for the

Re: [R] Access values in kpssstat-class

Have you tried 'str(test)'? This tells me it is Formal class 'kpssstat' [package uroot] with 5 slots. I can access the slot 'ltrunc' and the component 'rank' of the 'lmkpss' slot as follows: [EMAIL PROTECTED] [1] 7 [EMAIL PROTECTED] [1] 2 Hope this helps.

[R] Clique-Method-Package-Help

Dear R' Users, Is anyone knows, which package in R'Gui is suitable to use Clique technique in a dataset? Thanks in Advance JJ -- Lecturer J. Joshua Thomas KDU College Penang Campus Research Student, University Sains Malaysia [[alternative HTML version deleted]]

Re: [R] Keep value lables with data frame manipulation

Heinz Tuechler wrote: At 08:11 13.07.2006 -0500, Frank E Harrell Jr wrote: Heinz Tuechler wrote: At 13:14 12.07.2006 -0500, Marc Schwartz (via MN) wrote: On Wed, 2006-07-12 at 17:41 +0100, Jol, Arne wrote: Dear R, I import data from spss into a R data.frame. On this rawdata I do some data

[R] TR: Latent Class Analysis

_ De : Pousset [mailto:[EMAIL PROTECTED] Envoyé : mardi 4 juillet 2006 18:38 À : 'r-help@stat.math.ethz.ch' Objet : Latent Class Analysis Hello everybody, I am working on latent class analysis and have already used the ‘R’ function « lca » (in the e1071 package). I ‘ve got

Re: [R] Access to conditioning values in xyplot

On 7/13/06, Jean-Marc Ottorini [EMAIL PROTECTED] wrote: Dear R-help subscribers, Many thanks for all the answers I received either by mail or through the list, and that were most helpful. For the sake of the list records, I wanted to post the solution, I eventually could obtain, to

Re: [R] which data structure for a set of time series ?

A data.frame is perfect for storing both qualitative and quantitative variables together when they all have more or less the same length and where the observations are connected by something like time. Without substantive common connection, people usually put variables in different

Re: [R] Question for LM intercept

On Thu, 2006-07-13 at 15:02 +0100, Pavlidis, Efthymios wrote: Hello, I am having the following silly problem with lm. Let X be a dataframe with X[,1] the dependent variable and X[,-1] the independent variables. I want to run the following but without including an intercept. for(i in

Re: [R] Keep value lables with data frame manipulation

Further I do not see a simple method to label numerical variables. I often encounter discrete, but still metric data, as e.g. risk scores. Usually it would be nice to use them in their original coding, which may include zero or decimal places and to label them at the same time. ## For this

[R] R-2.3.1 --enable-R-shlib

Hi, I've upgraded our server to the latest R version, everything went ok until I wanted to install RServe for our java clients. RServe needs the shared library which can be compiled using configure --enable-R-shlib. All went ok but after running make I got this error: gcc -shared

Re: [R] convert ms() to optim()

Have you worked through the examples in the help pages for 'optim' and 'nlminb'? Also, 'nls', 'nlm', and 'mle{stats4}' provide related capabilities. If you'd like further help from this groups, please submit another post. To increase your chances of a quick, informative

[R] Invoke operating system command

Hi all, How can I invoke an operating system command in R? I mean something like exclamation mark (!) inside Matlab. Thanks, Gang __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] References verifying accuracy of R for basic statistical calculations and tests

Do you know of any references that verify the accuracy of R for basic statistical calculations and tests. The results of these studies should indicate that R results are the same as the results of other statistical packages to a certain number of decimal places on some benchmark calculations.

Re: [R] Invoke operating system command

?system __ 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] Invoke operating system command

Gang Chen wrote: Hi all, How can I invoke an operating system command in R? I mean something like exclamation mark (!) inside Matlab. See ?system (and additionally ?shell if you are on Windows). Uwe Ligges Thanks, Gang __

Re: [R] References verifying accuracy of R for basic statistical calculations and tests

Actually, you may not want R to agree so precisely with some of the packages since some well-known packages have some not-so-well-known inaccuracies. Clint BowmanINTERNET: [EMAIL PROTECTED] Air Dispersion Modeler INTERNET: [EMAIL PROTECTED] Air Quality

Re: [R] Invoke operating system command

?system Clint BowmanINTERNET: [EMAIL PROTECTED] Air Dispersion Modeler INTERNET: [EMAIL PROTECTED] Air Quality Program VOICE: (360) 407-6815 Department of Ecology FAX:(360) 407-7534 USPS: PO Box

Re: [R] Problem installing local additional packages of Rcmdr (WinXP)

Kenneth Cabrera wrote: Hi R users: I download the .zip files and also PACKAGES and PACKAGES.gz files in a local directory, but when I call library(Rcmdr) It shows me a tcl/tk windows telling me to install the additional packages, the it takes me to another window where I must specify

[R] MLE and QR classes

Hi, I load my data set and separate it as folowing: presu - read.table(C:/_Ricardo/Paty/qtdata_f.txt, header=TRUE, sep=\t, na.strings=NA, dec=., strip.white=TRUE) dep-presu[,3]; exo-presu[,4:92]; Now, I want to use it using the wls and quantreg packages. How I change the data classes for mle

[R] List 'sssitalk' closed to public posts

Sorry, but your posting is being returned to you. The list is temporarily restricted. If your post is substantive, dealing with SSSI issues, it will be forwarded to the list. If it is intended as private email, it will be forwarded accordingly. If you wish to subscribe to SSSITALK, send a

[R] Write a summary or a longer text to a graphical device

How can I write a text, such as a result of a function or an explanation to a graphic device? When I try plot, I'm unable to reduce the axes. Greetings torsten __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Write a summary or a longer text to a graphical device

?text ?plotmath Cheers Francisco Dr. Francisco J. Zagmutt College of Veterinary Medicine and Biomedical Sciences Colorado State University From: Torsten Mathies [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject: [R] Write a summary or a longer text to a graphical device Date: Thu, 13

[R] comparison between parametric and semiparametric models

Dear Listers, I have a question somehow off-topic but think I can find a good answer here. Using the same dataset, I've built a logistic regression and generalized additive model. Might I use likelihood ratio test to see if there is significant improvement in GAM model? Something like: (deviance

[R] collapsing plot lines

All, When plotting a group data object, is there a way to have the plots of all the groups collapsed onto one plot, i.e., instead of separate panels as would be produced from example below? grouped.data.example - groupedData(Y ~ Time | Patient_no, data = example.frm)

Re: [R] set the bahavior that R deal with missing values?

why dont you use Amelia? It's an R package which deal with missing data in R http://gking.harvard.edu/amelia/ Ferdinand Alimadhi wrote: You can deal with the missing values before performing any regression model. if you dataset is : D X1 X2 X3 1 NA 2 4 2 1 NA NA 3 2 2 6 4 NA 34 NA 5

[R] looping using combinatorics

I have a problem where I need to loop over the total combinations of vectors (combined once chosen via combinatorics). Here is a simplification of the problem: STEP 1: Define three vectors a, b, c. STEP 2: Combine all possible pairwise vectors (i.e., 3 choose 2 = 3 possible pairs of

Re: [R] Keep value lables with data frame manipulation

At 12:36 13.07.2006 -0400, Richard M. Heiberger wrote: Further I do not see a simple method to label numerical variables. I often encounter discrete, but still metric data, as e.g. risk scores. Usually it would be nice to use them in their original coding, which may include zero or decimal

[R] ols/gls or systemfit (OLS, WLS, SUR) give identical results

I might be sorry for asking this question :-) I have two equations and I tried to estimate them individually with lm and gls, and then in a system (using systemfit) with OLS, WLS and SUR. Quite surprisingly (for myself at least) the results are identical to the last digit. Could

Re: [R] References verifying accuracy of R for basic statistical calculations and tests

This might be a place to start: http://www.burns-stat.com/pages/Tutor/spreadsheet_addiction.html Among the references listed there are: Assessing the Reliability of Statistical Software: Part I by B. D. McCullough (1998) http://www.amstat.org/publications/tas/mccull-1.pdf Assessing the

Re: [R] looping using combinatorics

I assume your question is given 3 vectors of the same length: a, b and c how do we loop over pairs of them. In the following each iteration displays one pair: library(combinat) DF - data.frame(a = 1:4, b = 5:8, c = 9:12) for(idx in as.data.frame(combn(3,2))) print(DF[,idx]) On 7/13/06,

Re: [R] weights in glmrob

Have you worked through the several examples in the 'glmrob' help page? Those include at least one case that produces sensible standard errors and at least one that return NAs for standard errors. If this does not answer your question, I have two further suggestions:

[R] regarding filtering the data

Hi all a very happy greeting to you all. i have just joined the list, i am newbie in R and this is my first message ( before this i was a member of Linux kernel mailing list) and hence i woulk like to add that i know how to behave in the mailing lists. I have a data something like this which

Re: [R] regarding filtering the data

1. read.csv just calls read.table. Enter this into R to see its source code: read.csv 2. See ?subset . For example using the builtin swiss data frame: subset(swiss, Fertility Agriculture) or swiss[swiss$Fertility swiss$Agriculture, ] On 7/14/06, [EMAIL PROTECTED] [EMAIL