答复: [R] fail in adding library in new version.

2005-07-08 Thread Ivy_Li
Dear All, I imitated the web age Description example to make my description. I wrote: Package: example Version: 1.0 Date: 2005-07-08 Title: example library Author: Ivy [EMAIL PROTECTED] Maintainer: Ivy [EMAIL PROTECTED] Description: simple function License: GPL Version 2 or later.

Re: [R] [R-SIG-Mac] r version 2.1.0 and graphics in mac os 10.3.9

2005-07-08 Thread stefano iacus
Hi Luis, please remove /usr/local/lib/libxml* files and in particular /usr/local/lib/libxml2.2.dylib from your system sudo rm /usr/local/lib/libxml2.2.dylib stefano On 04/lug/05, at 20:42, Luis Borda de Agua wrote: I use mac os 10.3.9 and I've installed in my computer R 2.1.0 (I

Re: [R] Problem compiling R 2.1.* on SUSE 9.2

2005-07-08 Thread Sean O'Riordain
Hi! what point exactly in the command sequence does it fail? during make? or make install? You might try saving the output and the error from the build into a file? Something along the lines of ./configure 21 filename or even ./configure 21 filename | tee configure_output.txt etc. cheers!

Re: [R] r: LOOPING

2005-07-08 Thread Vehbi Sinan Tunalioglu
Uwe Ligges wrote: Clark Allan wrote: i know that one should try and limit the amount of looping in R programs. i have supplied some code below. i am interested in seeing how the code cold be rewritten if we dont use the loops. It is not always a good thing to remove loops (without having

[R] [OT] Dispersion in French

2005-07-08 Thread Ted Harding
Greetings, I'm posting this OT query here because of out very international membership! In the French sentence Les taux de tirage sont calculés de manière à ce que la dispersion soit inférieure à 5 % dans chaque strate. it would seem intended that the dispersion is to be calculated in a

Re: [R] xmat[1, 2:3] - NULL

2005-07-08 Thread Mikkel Grum
Thanks for this! This is even simpler than using !is.null(): xmat - as.data.frame(matrix(NA, 2, 3)) try(xmat[i, 2:3] - dbGetQuery(...), silent = TRUE) Best wishes, Mikkel --- Uwe Ligges [EMAIL PROTECTED] wrote: Mikkel Grum wrote: I have a situation where I'm filling out a dataframe

[R] Update on SuSE-9.2 fails

2005-07-08 Thread Thomas Schönhoff
Hi, yesterday, I switched to newest version of R-2.1.1 by using YAst to install necessary SuSE-rpms from a separate local directory. So far , so good. Afterwards I tried to find out if there are any updates for add-on packages, and of cause there were some of them requiring a refreshment.

Re: [R] [OT] Dispersion in French

2005-07-08 Thread Philippe Grosjean
Hello Ted, I would interpret this sentence the same way as you do. for me, in French, dispersion is a general term and there is no clue of which dispersion measurement was used... perhaps elsewhere in your text? Best, Philippe ..°})) ) ) )

[R] about image() function in R and colors

2005-07-08 Thread Mikkel Grum
Javier, on the image and colors, did you try using the breaks option? You will need 21 breaks to display 20 colors (see ?image). Janek, the option asp = 1 should generally solve your scaling problem (point 2 below). I don't remember where I got that from, but it works unless you use functions

Re: [R] Problems with nlme: Quinidine example

2005-07-08 Thread Dieter Menne
rich at mi.fu-berlin.de writes: This concerns the Clinical Study of Quinidine example on page 380 of the book Mixed-Effects Models in S and S-PLUS by Pinheiro and Bates (2000). I have tried to reproduce the example, but get an error: .. system is computationally singular:

[R] Orthogonal regression

2005-07-08 Thread Vito Ricci
Dear R-Users, is there any statement to fit a orthogonal regression in R environment? Many thanks in advance. Best regards, Vito Diventare costruttori di soluzioni Became solutions' constructors The business of the statistician is to catalyze the scientific learning process. George E. P.

Re: [R] Making Package, Chm error, Html Help Workshop

2005-07-08 Thread Lars
Hi Matthias, I had the same problem a couple of weeks ago. It wasn't easy to come around it but i made it. Did you use the skeloton function in R? If yes: Did you fill out both the discripton-form and the help-pages? Did you follow all the further steps recommended? Did you point in the

Re: [R] fail in adding library in new version.

2005-07-08 Thread Duncan Murdoch
Ivy_Li wrote: Dear all, I really appreciate your help. I think I have a little advancement. ^_^ When I enter the Dos environment, at first, into the D:\, I type the following code: cd Program Files\R\rw2011\ bin\R CMD install /example example is in the d:\, which

Re: [R] r: LOOPING

2005-07-08 Thread Duncan Murdoch
Vehbi Sinan Tunalioglu wrote: Uwe Ligges wrote: Clark Allan wrote: i know that one should try and limit the amount of looping in R programs. i have supplied some code below. i am interested in seeing how the code cold be rewritten if we dont use the loops. It is not always a good thing to

[R] Possible Solution to Tempfile error (for documentation)

2005-07-08 Thread Doran, Harold
Dear List: I'm posting this to provide a possible solution and to document to what appears to be an R limitation. The solution is more of a cheap hack that works for now. To provide a little background, I am looping through a dataframe and creating Sweave documents using data from each row in the

Re: [R] pairs() uses col argument for axes coloring

2005-07-08 Thread Duncan Murdoch
Olaf Mersmann wrote: Hi list, not sure if this is the wanted behavior, but running the following code: I'd say it's a bug. version platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major2 minor1.1 year 2005 month06 day

Re: [R] Problems with nlme: Quinidine example

2005-07-08 Thread Douglas Bates
On 7/8/05, Dieter Menne [EMAIL PROTECTED] wrote: rich at mi.fu-berlin.de writes: This concerns the Clinical Study of Quinidine example on page 380 of the book Mixed-Effects Models in S and S-PLUS by Pinheiro and Bates (2000). I have tried to reproduce the example, but get an error:

Re: [R] Possible Solution to Tempfile error (for documentation)

2005-07-08 Thread Duncan Murdoch
Doran, Harold wrote: Dear List: I'm posting this to provide a possible solution and to document to what appears to be an R limitation. The solution is more of a cheap hack that works for now. To provide a little background, I am looping through a dataframe and creating Sweave documents

Re: [R] Making Package, Chm error, Html Help Workshop [solved]

2005-07-08 Thread TEMPL Matthias
Dear Uwe Ligges, There was really a problem with an irregular name of a help topic (%cin% - the % was the problem) as you said. With the new R Version and with the corrected Rd file all works fine. Thanks a lot for your help! Matthias TEMPL Matthias wrote: Hello, When building my

Re: [R] [OT] Dispersion in French

2005-07-08 Thread Renaud Lancelot
(Ted Harding) a écrit : Greetings, I'm posting this OT query here because of out very international membership! In the French sentence Les taux de tirage sont calculés de manière à ce que la dispersion soit inférieure à 5 % dans chaque strate. it would seem intended that the

Re: [R] [OT] Dispersion in French

2005-07-08 Thread Peter Dalgaard
Renaud Lancelot [EMAIL PROTECTED] writes: (Ted Harding) a écrit : Greetings, I'm posting this OT query here because of out very international membership! In the French sentence Les taux de tirage sont calculés de manière à ce que la dispersion soit inférieure à 5 % dans

Re: [R] Orthogonal regression

2005-07-08 Thread davidr
This topic has come up a few times recently, so it must be 'in the wind' these days. Depending on what approach you take and what area you are coming from, it goes under the names Orthogonal [Distance] Regression, Total Least Squares, Errors in Variables, Deming Regression. The middle two

[R] removing factor level represented by less than x rows

2005-07-08 Thread Mikkel Grum
In a number of different situations I'm trying to remove factor levels that are represented by less than a certain number of rows, e.g. if I had the dataset aa below and wanted to remove the species that are represented in less than 2 rows: data(iris) aa - iris[1:101,] In this case, since I can

[R] Garch in a model with explanatory variables

2005-07-08 Thread Carlo Fezzi
Dear helpers, does anyone know a function to fit a model with: - y mean that is regressed on a set of explanatory variables - y variace behaving as a garch or as a garch in mean Thank you so much for your help, Carlo __ R-help@stat.math.ethz.ch

[R] segfault on q() using RODBC ... was q() == Segmentation fault

2005-07-08 Thread Omar Lakkis
-- Forwarded message -- From: Omar Lakkis [EMAIL PROTECTED] Date: Jul 7, 2005 6:25 PM Subject: q() == Segmentation fault To: r-help@stat.math.ethz.ch I created the simple library, attached. When I terminate an R session where the library has been loaded with q() a segmentation

Re: [R] removing factor level represented by less than x rows

2005-07-08 Thread Sundar Dorai-Raj
Mikkel Grum wrote: In a number of different situations I'm trying to remove factor levels that are represented by less than a certain number of rows, e.g. if I had the dataset aa below and wanted to remove the species that are represented in less than 2 rows: data(iris) aa -

[R] exact conditional mantelhaen.test estimate is 0 ?!

2005-07-08 Thread Florian Hahne
Dear listers, I am trying to compute the exact conditional test given strata margins of a 2 by 2 by K array using the mantelhaen.test function to get a common odds ratio estimate. The estimate for the test on the following data is 0, which in my opinion dosen't make any sense. x - array(c(53,

Re: [R] removing factor level represented by less than x rows

2005-07-08 Thread Frank E Harrell Jr
Mikkel Grum wrote: In a number of different situations I'm trying to remove factor levels that are represented by less than a certain number of rows, e.g. if I had the dataset aa below and wanted to remove the species that are represented in less than 2 rows: data(iris) aa - iris[1:101,]

Re: [R] pairs() uses col argument for axes coloring

2005-07-08 Thread Anne York
A work-around for this problem was posted on April 7 by Bill Venables and Deepayan Sarkar using a panel function: Anne On Thu, 7 Apr 2005, Deepayan Sarkar wrote: DS On Thursday 07 April 2005 17:51, Anne York wrote: DS The following command produces red axis line in a pairs DS plot: DS

[R] ANN: BioC2005 Conference, online registration now available

2005-07-08 Thread Seth Falcon
BioC2005: Where Biology and Software Connect August 16 and 17 in Seattle, WA, USA The online registration form is now available. http://www.bioconductor.org/meeting05/ (You will be redirected to our secure server: https://cobra.fhcrc.org/BioC2005/) About BioC2005: This conference will

[R] Plotting a simple subset

2005-07-08 Thread Michael Hopkins
Hi all Just converting from Stata to R and struggling a little to come to terms with the new philosophy/command line. E.g. I want to plot x against y if x 5 In Stata: graph x y, if( x 5 ) How do I do this in R? Have tried most of the obvious options without success. Can I have multiple

Re: [R] Plotting a simple subset

2005-07-08 Thread Berton Gunter
Please first read An Introduction to R (one of the pdf manuals that ships with R) before posting these sorts of questions, as it is written specifically to help you get started (I think fairly clearly). Other (links to) learning resources may be found on the CRAN website. Please take advantage of

Re: [R] Error message NA/NaN/Inf in foreign function call (arg 6) when using knn()

2005-07-08 Thread Uwe Ligges
Kerri-Ann Norton wrote: I am trying to use knn to do a nearest neighbor classification. I tried using my dataset and got an error message so I used a simple example to try and understand what I was doing wrong and got the same message. Here is what I typed into R: try [,1] [,2] [,3]

Re: [R] how to do something like symptoms==c('a', 'e', 'z')

2005-07-08 Thread Uwe Ligges
RenE J.V. Bertin wrote: I find myself doing lots of tests like subset( data, symptoms=='a' | symptoms=='e' | symptoms=='z' ) with symptoms one of the factors contained in the data frame. and I wonder if there is not an existing operator or function which implements this sort

Re: [R] how to do something like symptoms==c('a', 'e', 'z')

2005-07-08 Thread Chuck Cleland
RenE J.V. Bertin wrote: I find myself doing lots of tests like subset( data, symptoms=='a' | symptoms=='e' | symptoms=='z' ) with symptoms one of the factors contained in the data frame. and I wonder if there is not an existing operator or function which implements this sort of

Re: [R] Plotting a simple subset

2005-07-08 Thread Peter Wolf
Michael Hopkins wrote: Hi all Just converting from Stata to R and struggling a little to come to terms with the new philosophy/command line. E.g. I want to plot x against y if x 5 In Stata: graph x y, if( x 5 ) How do I do this in R? Have tried most of the obvious options without success.

Re: [R] Plotting a simple subset

2005-07-08 Thread Gabor Grothendieck
To be fair none of Introduction to R, ?plot nor the reference card really cover this without substantial digging. # test data x - 1:10 y - x*x plot(x[x 5], y[x 5]) # or plot(y ~ x, subset = x 5) # We can have combine conditions like this: plot(y ~ x, subset = x 5 y 50) # also if

Re: [R] Plotting a simple subset

2005-07-08 Thread Michael Hopkins
Thanks to all of the respondents for helpful (except one) and super-fast replies! I think the two packages differ more than I thought in terms of philosophy, but I have been seduced by the potential power and graphics facilities available in R (not to mention cost and cross-platform

Re: [R] comparing strength of association instead of strength of evidence?

2005-07-08 Thread Weiwei Shi
Dear all: I still need some further help since I think the question itself might be very interesting (i hope so:) : the question is on chisq.test, my concern is which criteria should be used here to evaluate the independence. The reason i use this old subject of the email is, b/c I think the

Re: [R] comparing strength of association instead of strength ofevidence?

2005-07-08 Thread Berton Gunter
But isn't this the old chestnut that any effect will be found significant given enough data; and with too few data, not even a large one can be distinguished from noise? If so, it's a good question that has more to do with the philosophy of science than statistics. Bayesians, of course, would

[R] Problem with filter() ?

2005-07-08 Thread Matthew Dowle
Dear list, This is ok : filter(1:5, 1, recursive) Time Series: Start = 1 End = 5 Frequency = 1 [1] 1 3 6 10 15 But this? : filter(c(rep(NA,5),1:5), 1, recursive) Time Series: Start = 1 End = 10 Frequency = 1 [1] NA 0 NA 0 NA 0 2 5 9 14 version _

Re: [R] exact conditional mantelhaen.test estimate is 0 ?!

2005-07-08 Thread Thomas Lumley
On Fri, 8 Jul 2005, Florian Hahne wrote: Dear listers, I am trying to compute the exact conditional test given strata margins of a 2 by 2 by K array using the mantelhaen.test function to get a common odds ratio estimate. The estimate for the test on the following data is 0, which in my

Re: [R] comparing strength of association instead of strength ofevidence?

2005-07-08 Thread Ravi Varadhan
Hi, The following article by William DuMouchel, which takes an empirical Bayes approach, might be helpful to you: Bayesian Data Mining in Large Frequency Tables, with an Application to the FDA Spontaneous Reporting System, The American Statistician, Vol. 53, No. 3 (Aug., 1999), pp. 177-190.

[R] missing data imputation

2005-07-08 Thread Anders Schwartz Corr
Dear R-help, I am trying to impute missing data for the first time using R. The norm package seems to work for me, but the missing values that it returns seem odd at times -- for example it returns negative values for a variable that should only be positive. Does this matter in data analysis,

[R] Help with Mahalanobis

2005-07-08 Thread Jose Claudio Faria
Dear R list, I'm trying to calculate Mahalanobis distances for 'Species' of 'iris' data as obtained below: Squared Distance to Species From Species: Setosa Versicolor Virginica Setosa 0 89.86419 179.38471 Versicolor 89.86419 0 17.20107 Virginica

[R] explained deviance in multinom

2005-07-08 Thread alexbri
Hi: I'm working with multinomial models with library nnet, and I'm trying to get the explained deviance (pseudo R^2) of my models. I am assuming that: pseudo R^2= 1 - dev(model) / dev (null) where dev(model) is the deviance for the fitted model and dev(null) is the deviance for the null

Re: [R] Help with Mahalanobis

2005-07-08 Thread Christian Hennig
Dear Jose, normal mixture clustering (mclust) operates on points times variables data and not on a distance matrix. Therefore it doesn't make sense to compute Mahalanobis distances before using mclust. Furthermore, cluster analysis based on distance matrices (hclust or pam, say) operates on a

[R] extract prop. of. var in pca

2005-07-08 Thread K. Steinmann
Dear R-helpers, Using the package Lattice, I performed a PCA. For example pca.summary - summary(pc.cr - princomp(USArrests, cor = TRUE)) The Output of pca.summary looks as follows: Importance of components: Comp.1Comp.2Comp.3 Comp.4 Standard deviation

Re: [R] extract prop. of. var in pca

2005-07-08 Thread Chuck Cleland
K. Steinmann wrote: Dear R-helpers, Using the package Lattice, I performed a PCA. For example pca.summary - summary(pc.cr - princomp(USArrests, cor = TRUE)) The Output of pca.summary looks as follows: Importance of components: Comp.1Comp.2Comp.3

[R] Overlying a Normal Dist in a Barplot

2005-07-08 Thread Bret Collier
R-Users, Hopefully someone can shed some light on these questions as I had little luck searching the archives (although I probably missed something in my search due to the search phrase). I estimated multinomial probabilities for some count data (number successful offspring) ranging from 0 to 8

Re: [R] extract prop. of. var in pca

2005-07-08 Thread Achim Zeileis
On Fri, 8 Jul 2005, K. Steinmann wrote: Dear R-helpers, Using the package Lattice, I performed a PCA. In my R installation, the function princomp() is contained in the stats package, not lattice (sic!). For example pca.summary - summary(pc.cr - princomp(USArrests, cor = TRUE)) The Output

[R] Multiple assignments in one statement

2005-07-08 Thread Jeffrey Horner
Is this possible? For instance, I have a function that returns a vector length 3. In one statement I'd like to assign each element of the vector to different variables. Syntactically, I hoped this would work: c(x,y,z) - myfun(); Thanks, -- Jeffrey Horner Computer Systems Analyst

Re: [R] Multiple assignments in one statement

2005-07-08 Thread Huntsinger, Reid
Yes and no. Your function doesn't work because c() is a function and R is call-by-value, so c(x,y,z) is just a value, like 7. You could use lists like this: l - list(x=NULL,y=NULL,z=NULL) l[c(x,y,z)] - 1:3 l $x [1] 1 $y [1] 2 $z [1] 3 That, together with attaching lists (help(attach)),

Re: [R] Finding indices of NA values in a data frame

2005-07-08 Thread Sundar Dorai-Raj
Ken Termiso wrote: Hi all, I've got a data frame with NA values scattered throughout. I would like to find the rows which contain an NA. Trying to find the indices of the NAs with grep() doesn't work (apparently you can't convert NA to NA character value, even with as.character()) so I

[R] Finite Mixture Models with logistic regression

2005-07-08 Thread Feng Gao
Do we have any R package that can do analysis on finite mixture model with logistic regression? Thanks Faith Feng Gao Dept. of Statistics Virginia Tech. Email: [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list

[R] time series regression

2005-07-08 Thread yyan liu
Hi: I have two time series y(t) and x(t). I want to regress Y on X. Because Y is a time series and may have autocorrelation such as AR(p), so it is not efficient to use OLS directly. The model I am trying to fit is like Y(t)=beta0+beta1*X(t)+rho*Y(t-1)+e(t) e(t) is iid normal random error.

Re: [R] missing data imputation

2005-07-08 Thread Frank E Harrell Jr
Anders Schwartz Corr wrote: Dear R-help, I am trying to impute missing data for the first time using R. The norm package seems to work for me, but the missing values that it returns seem odd at times -- for example it returns negative values for a variable that should only be positive. Does

[R] help with ARIMA and predict

2005-07-08 Thread Brian Scholl
I'm trying to do the following out of sample regression with autoregressive terms and additional x variables: y(t+1)=const+B(L)*y(t)+C(1)*x_1(t)...+C(K)*x_K(t) where: B(L) = lag polynom. for AR terms C(1..K) = are the coeffs. on K exogenous variables that have only 1 lag Question 1:

Re: [R] time series regression

2005-07-08 Thread Gabor Grothendieck
On 7/8/05, yyan liu [EMAIL PROTECTED] wrote: Hi: I have two time series y(t) and x(t). I want to regress Y on X. Because Y is a time series and may have autocorrelation such as AR(p), so it is not efficient to use OLS directly. The model I am trying to fit is like

[R] merge

2005-07-08 Thread Ling Jin
Hi all, I have two data frames to merge by a column containing the site names (as characters). However, somehow, one of the site names of one data frame have fixed length, say 8, so the names sometimes have spaces at the end. For example, the site name is ST, but in one data frame, it is ST

Re: [R] more and tab functionalities in R under linux

2005-07-08 Thread Gabor Grothendieck
Here are some possibilities: - head(iris) will show the first few rows of the data frame - edit(iris) will put up a spreadsheet with the data frame in it that you can scroll - In JGR (a GUI front end for R) you can use the object browser (ctrl-B) - If the object is a file rather than a data frame

Re: [R] Multiple assignments in one statement

2005-07-08 Thread Gabor Grothendieck
Check out: http://tolstoy.newcastle.edu.au/R/help/04/06/1430.html http://tolstoy.newcastle.edu.au/R/help/04/06/1406.html On 7/8/05, Jeffrey Horner [EMAIL PROTECTED] wrote: Is this possible? For instance, I have a function that returns a vector length 3. In one statement I'd like to assign

Re: [R] merge

2005-07-08 Thread Gabor Grothendieck
trim in package gdata will trim spaces off the beginning and end. On 7/8/05, Ling Jin [EMAIL PROTECTED] wrote: Hi all, I have two data frames to merge by a column containing the site names (as characters). However, somehow, one of the site names of one data frame have fixed length, say 8,

Re: [R] merge

2005-07-08 Thread Marc Schwartz
One other option during the import is to set 'strip.white = TRUE' in read.csv(). See ?read.csv for more information. Bear in mind that this will strip both leading and trailing white space in all columns, which may have unintended consequences. Yet another post-import option, would be to use

[R] It was a sad day for the Statistics profession

2005-07-08 Thread Liaw, Andy
For those of you who have not heard: Prof. Breiman passed away on July 5th. http://www.berkeley.edu/news/media/releases/2005/07/07_breiman.shtml Andy __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Overlying a Normal Dist in a Barplot

2005-07-08 Thread Gabor Grothendieck
On 7/8/05, Bret Collier [EMAIL PROTECTED] wrote: R-Users, Hopefully someone can shed some light on these questions as I had little luck searching the archives (although I probably missed something in my search due to the search phrase). I estimated multinomial probabilities for some count

Re: [R] Overlying a Normal Dist in a Barplot

2005-07-08 Thread Marc Schwartz
On Fri, 2005-07-08 at 21:58 -0400, Gabor Grothendieck wrote: On 7/8/05, Bret Collier [EMAIL PROTECTED] wrote: R-Users, Hopefully someone can shed some light on these questions as I had little luck searching the archives (although I probably missed something in my search due to the search

Re: [R] Overlying a Normal Dist in a Barplot

2005-07-08 Thread Gabor Grothendieck
Offline, Marc pointed out to me that boxplot has an at= argument. This suggests that we could substitute a boxplot command for the rect command since a boxplot of c(0,a) looks like a bar from 0 to a if we use medlty=0 (which omits the median line) and boxwex=1 (which eliminates the space between

Re: [R] help with ARIMA and predict

2005-07-08 Thread Spencer Graves
What class is newx when you get the error message? Is it a vector or a 1xK array? If the former, force it to be an array. (Hint: array(1:4, dim=c(2,2))[1,] is a vector of length 2, while array(1:4, dim=c(2,2))[1,,drop=FALSE] is a 1 x 2 matrix..( If this does not solve