Re: [R] R-help Digest, Vol 54, Issue 30

2007-08-31 Thread David Duffy
Ron Crump wrote: Hi, I have a dataframe that contains pedigree information; that is individual, sire and dam identities as separate columns. It also has date of birth. These identifiers are not numeric, or not sequential. Obviously, an identifier can appear in one or two columns,

Re: [R] R Help

2007-08-28 Thread Henrique Dallazuanna
You don't have installed the akima pakage. install.packages(akima, dep=T) -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 28/08/07, Ola Asteman [EMAIL PROTECTED] wrote: I got the Warning message below when I tried to load Locfit. What is wrong? Regards Ola

Re: [R] R Help

2007-08-28 Thread Petr PIKAL
[EMAIL PROTECTED] napsal dne 28.08.2007 13:33:13: You don't have installed the akima pakage. install.packages(akima, dep=T) And wait about two months and update your R version to 2.6.0. Or update now to 2.5.1 Regards Petr -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40

Re: [R] R help

2007-06-07 Thread Jim Lemon
scott flemming wrote: Hi, I wonder whether R can finish the following project: I want to make a chart to represent 10 genes. Each gene has orientation and length. Therefore, a gene can be represented by arrows. Can R be used to draw 10 arrows in one line ? Hi Scott, Maybe the

Re: [R] R help

2007-06-06 Thread Sarah Goslee
On 6/5/07, scott flemming [EMAIL PROTECTED] wrote: Can R be used to draw 10 arrows in one line ? Um, sure. Assuming you actually also want to know how to do it, why don't you take a look at the help for arrows(). Sarah -- Sarah Goslee http://www.functionaldiversity.org

Re: [R] R help

2007-06-06 Thread michael watson \(IAH-C\)
Yes, but you need to be a bit more specific... When it comes to graphs and drawing lines, there isn't much R can't do... -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of scott flemming Sent: 06 June 2007 04:49 To: r-help@stat.math.ethz.ch Subject: [R] R

Re: [R] R help

2007-06-06 Thread John Kane
--- scott flemming [EMAIL PROTECTED] wrote: Hi, I wonder whether R can finish the following project: I want to make a chart to represent 10 genes. Each gene has orientation and length. Therefore, a gene can be represented by arrows. Can R be used to draw 10 arrows in one line ?

Re: [R] R-help with apply and ccf

2007-05-22 Thread ONKELINX, Thierry
You could combine them with cbind, and then split the rows again inside the function you're calling with apply. Mat - cbind(mat1, mat2) apply(Mat, 1, function(x){ row.mat1 - x[seq_len(length(x)/2)] row.mat2 - x[length(x)/2 + seq_len(length(x)/2)] cor(row.mat1, row.mat2) })

Re: [R] R-help with apply and ccf

2007-05-22 Thread Bert Gunter
I understand you to want correlations of corresponding rows (** not ccf, which returns a vector ccf for each pair of rows). If that is so, 1) ... in theory, diag(cor(t(A), t(B)) would work without apply, except 196,000 rows is probably too large, and it is probably too inefficient to compute and

Re: [R] R-help Digest, Vol 46, Issue 27

2006-12-27 Thread Grant Izmirlian
On Wednesday 27 December 2006 06:00, [EMAIL PROTECTED] wrote: jingjiangyan I agree, you can use 'assign'. To be more explicit, you could use the following function. jingjiangyan - function(formula, data) { m - match.call() %,% - function(x,y)paste(x,y,sep=) d.nm - as.character(m$data)

Re: [R] R-Help

2006-12-06 Thread Uwe Ligges
You might want to contact the nortest maintainer, Juergen Gross (CCing), who is not listening to the traffic on this list. Uwe Ligges amna khan wrote: Respected Sir I am a very new user of R. I want to ask a question about the nortest package. In this package how we can write the code of

Re: [R] R Help

2006-11-11 Thread David Barron
d - read.table(lahore.txt, header=TRUE) On 11/11/06, amna khan [EMAIL PROTECTED] wrote: Respected Sir I request you to please fill the following read.table function and read.csvfor my understanding by assuming my data attached with this maiL, because I am fail to run these functions using

Re: [R] R Help

2006-11-11 Thread Benilton Carvalho
data = read.delim(lahore.txt) is enough for what you want to do. b On Nov 11, 2006, at 2:11 PM, amna khan wrote: Respected Sir I request you to please fill the following read.table function and read.csvfor my understanding by assuming my data attached with this maiL, because I am fail to

Re: [R] R Help

2006-11-11 Thread Nicolas Mazziotta
Not sure of what you ask... Does it help? $ read.table(Desktop/lahore.txt, header=T) Years lrmax n1 n2 n3 n4 arranged 1 1980 207.6 25 24 23 22 29.4 2 1981 92.7 24 23 22 21 49.4 3 1982 67.5 23 22 21 20 55.1 4 1983 93.8 22 21 20 19 58.0 5 1984 60.6 21 20 19 18

Re: [R] R Help

2006-11-11 Thread Jeffrey Robert Spies
If you were trying to use the sep=\t argument, you might have encountered an error, as there are three tabs one of the two blank lines at the end of your data file. The default for read.table and read.delim (as has been suggested by David and Benilton) is whitespace, which consumes the

Re: [R] R-help in Newsreader?

2006-11-06 Thread Jens Scheidtmann
Matthias Voigt [EMAIL PROTECTED] writes: Dear list, I am new to this all and therefore have following Newbie question: How can I receive and read R-help mailings in a newsreader like thunderbird? Use gmane.comp.lang.r.general on gmane.org HTH, Jens

Re: [R] R-help in Newsreader?

2006-11-03 Thread Marc Schwartz
On Fri, 2006-11-03 at 17:36 +0100, Matthias Voigt wrote: Dear list, I am new to this all and therefore have following Newbie question: How can I receive and read R-help mailings in a newsreader like thunderbird? Thanks Matthias If you go to the main R web site

Re: [R] R-help in a newsgroup

2006-07-18 Thread Marc Schwartz (via MN)
On Tue, 2006-07-18 at 11:30 -0700, Darren Weber wrote: Hi, I find a lot of the R-help email traffic overloads my inbox. My IT managers are not really happy for me to be subscribed to several high-traffic email lists. I don't want to lose my contact with the R-help emails, so I'm having to

Re: [R] R-help in a newsgroup

2006-07-18 Thread Duncan Murdoch
On 7/18/2006 2:30 PM, Darren Weber wrote: Hi, I find a lot of the R-help email traffic overloads my inbox. My IT managers are not really happy for me to be subscribed to several high-traffic email lists. I don't want to lose my contact with the R-help emails, so I'm having to consider

Re: [R] r-help at stat.math.ethz.ch

2006-05-16 Thread Ben Bolker
Xin jasonshi510 at hotmail.com writes: Dear All: Then error messga there: initial value in 'vmmin' is not finite In addition: There were 38 warnings (use warnings() to see them). Could you give some advice please? Thanks a lot! Xin Shi [[alternative HTML version

Re: [R] R-help Digest, Vol 39, Issue 13

2006-05-13 Thread Alan Cobo-Lewis
r-help@stat.math.ethz.ch on Saturday, May 13, 2006 at 6:00 AM -0500 wrote: lme(biomass~age, random=~woods/age)? Jörn Consult Pinheiro and Bates (2000, Mixed-effects models in S and S-Plus, Springer, ISBN 0-387-98957-0 ref 7 at http://www.r-project.org/doc/bib/R-books.html ) for how to fit

Re: [R] R-help Digest, Vol 38, Issue 30

2006-04-30 Thread isaac . martin
Mi nueva dirección de correo es: [EMAIL PROTECTED] New e-mail address: [EMAIL PROTECTED] __ 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] R help

2006-04-25 Thread Doran, Harold
Ezra I don't know what the elements of your matrix are, but if there are a large proportion of 0s you can work with sparse matrices in the Matrix package. Harold -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Erez Sent: Tuesday, April 25, 2006 8:39

Re: [R] R help

2006-04-25 Thread Gabor Grothendieck
A similar question was just asked. See: http://tolstoy.newcastle.edu.au/R/help/06/04/25898.html On 4/25/06, Erez [EMAIL PROTECTED] wrote: Hello, I'm working with large matrix data and i would like to know if there is any way to reduce the size of it because even that I'm increasing the

Re: [R] R help

2006-04-24 Thread Petr Pikal
Hi more memory new comp new OS think about possibility to reformulate the problem with help of database and loading/processing data in chunks. Couple of similar questions were answered not long ago so check archives. HTH Petr On 24 Apr 2006 at 13:06, Erez wrote: Date sent:

Re: [R] R-Help

2006-04-20 Thread Gabor Grothendieck
Read them in as zoo objects (you can replace textConnection(Lines1) with the filename) and then merge them using all = FALSE to retain only common time points. Note that in my English locale I had to modify your Apl to Apr. Lines1 - Date x y Jan-1,2005120 230

Re: [R] R-Help

2006-04-20 Thread john seers \(IFR\)
I think this does what you require. #Read your data in whatever way you wish: d1-data.frame(Date=c(2005/1/1,2005/2/1,2005/1/3,2005/1/4,2005/ 1/7,2005/3/5), x=c(119,123,-110,114,11,200), y=c(230,-125,300,-21,299,311))

Re: [R] R-help Digest, Vol 38, Issue 19

2006-04-19 Thread isaac . martin
Mi nueva dirección de correo es: [EMAIL PROTECTED] New e-mail address: [EMAIL PROTECTED] __ 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] R-help Digest, Vol 38, Issue 9

2006-04-09 Thread isaac . martin
Mi nueva dirección de correo es: [EMAIL PROTECTED] New e-mail address: [EMAIL PROTECTED] __ 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] R-help Digest, Vol 37, Issue 26

2006-03-26 Thread isaac . martin
Mi nueva dirección de correo es: [EMAIL PROTECTED] New e-mail address: [EMAIL PROTECTED] __ 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] R-help Digest, Vol 37, Issue 15

2006-03-15 Thread isaac . martin
Mi nueva dirección de correo es: [EMAIL PROTECTED] New e-mail address: [EMAIL PROTECTED] __ 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] R-help Digest, Vol 37, Issue 12

2006-03-12 Thread Ferran Carrascosa
Hi r-users, I would like to know if R have any solution to the Address standardization. The problem is to classify a database of addresses with the real addresses of a streets of Spain. Ideally, I would like to assign Postal code, census data and other geographic information. If this is not

Re: [R] R-help Digest, Vol 37, Issue 12

2006-03-12 Thread Liaw, Andy
From: Ferran Carrascosa Hi r-users, I would like to know if R have any solution to the Address standardization. The problem is to classify a database of addresses with the real addresses of a streets of Spain. Ideally, I would like to assign Postal code, census data and other

Re: [R] R-help list: temporary problem in local archives

2006-03-06 Thread Martin Maechler
Dimitri == Dimitri Szerman [EMAIL PROTECTED] on Mon, 6 Mar 2006 14:17:33 -0300 writes: Dimitri Hi, Dimitri It seems the list faced some problems during the Dimitri weekend, so I am re-sending this message. To be specific: The only problems it saw was that the local

Re: [R] R-help list: temporary problem in local archives

2006-03-06 Thread Martin Maechler
Martin == Martin Maechler [EMAIL PROTECTED] on Mon, 6 Mar 2006 18:37:06 +0100 writes: Dimitri == Dimitri Szerman [EMAIL PROTECTED] on Mon, 6 Mar 2006 14:17:33 -0300 writes: Dimitri Hi, Dimitri It seems the list faced some problems during the Dimitri weekend, so I am

Re: [R] R-help Digest, Vol 37, Issue 1

2006-03-01 Thread isaac . martin
Mi nueva dirección de correo es: [EMAIL PROTECTED] New e-mail address: [EMAIL PROTECTED] __ 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] R-help Digest, Vol 36, Issue 21

2006-02-21 Thread Evgeniy Kachalin
Hello, dear R users. I've already sent a question here, but I'm not sure that it had been read. I need to visualize classification of my numerical data based on 2-3 factors. As I suppose, the best way is a tree. With an orbitrary function at the ends (leaves), or at least with means of my data

Re: [R] R-help Digest, Vol 36, Issue 21

2006-02-21 Thread Evgeniy Kachalin
Evgeniy Kachalin wrote: Hello, dear R users. I've already sent a question here, but I'm not sure that it had been read. I need to visualize classification of my numerical data based on 2-3 factors. As I suppose, the best way is a tree. With an orbitrary function at the ends (leaves), or

Re: [R] R-help, specifying the places to decimal

2006-02-13 Thread Dimitris Rizopoulos
have a look at ?round() Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015 Web: http://www.med.kuleuven.be/biostat/

Re: [R] R-help, specifying the places to decimal

2006-02-13 Thread Vivek Satsangi
In addition to round() mentioned earlier, if you are merely looking to *display* your results differently, you may want to check out the digits option, e.g. in summary(): (This is the method signature for data.frame 's): summary(object, maxsum = 7, digits = max(3,

Re: [R] R-help Digest, Vol 35, Issue 24

2006-01-25 Thread Ted Harding
I've been reluctant to step into this topic, but now feel that it may be helpful to make a certain point. On the internet, for the most part, the person behind the email is invisible and intangible. It is therefore possible, when someone puts their foot down, to stamp inadvertently on someone

Re: [R] R-help Digest, Vol 35, Issue 24

2006-01-25 Thread François Pinard
[Gabor Grothendieck] [...] this list is inhabited by some rather rude participants but everyone puts up with them in the hope that they do have some useful remarks. I've been witnessing this list for about one year, and also read *lots* of archived messages. While it is true that a few members

Re: [R] R-help Digest, Vol 35, Issue 24

2006-01-24 Thread Keith . Chamberlain
Dear Prof Ripley, First of all, unless you are an english professor, then I do not think you have any business policing language. I'm still very much a student, both in R, and regarding signal analysis. My competence on the subject as compared too your own level of expertise, or my spelling for

Re: [R] R-help Digest, Vol 35, Issue 24

2006-01-24 Thread Spencer Graves
Dear Mr. Chamberlain: You asked for free consulting, and as near as I can tell, you got pretty good advice. Now you complain that you don't like the packaging. If you can't stand the heat, get out of the kitchen. Professor Brian Ripley has an international reputation

Re: [R] R-help Digest, Vol 35, Issue 24

2006-01-24 Thread François Pinard
[EMAIL PROTECTED], addressing to Brian Ripley] First of all, unless you are an english professor, then I do not think you have any business policing language. We all do mistakes (English or otherwise). I'm very grateful that people forgive my own errors, and I try to be tolerant to others.

Re: [R] R-help Digest, Vol 35, Issue 24

2006-01-24 Thread Gabor Grothendieck
Its not really you. Its a fact of life that this list is inhabited by some rather rude participants but everyone puts up with them in the hope that they do have some useful remarks. This has been discussed repeatedly on the list and there is even a group of thought that feels it is a justifiable

Re: [R] R-help Digest, Vol 35, Issue 23

2006-01-23 Thread Dr. Herwig Meschke
summary.aov(aovRes, split=list(interval = list(i1 vs i2 = 1, i2 vs i3 = 2, i3 vs i4 = 3, i4 vs i5 = 4, i5 vs i6 = 5))) try class(aovRes) #- aovlist ! summary.aovlist(aovRes, spit=...) or simply summary(aovRes, spit=...) Hoping this helps, Herwig -- Dr. Herwig Meschke Wissenschaftliche

Re: [R] r-help, how can i use my own distance matrix without usin g dist()

2006-01-18 Thread Liaw, Andy
Use something like hclust(as.dist(mydist), ...) ought to work. Andy -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Wednesday, January 18, 2006 4:47 PM To: r-help@stat.math.ethz.ch Subject: [R] r-help, how can i use my own

Re: [R] r-help, how can i use my own distance matrix without using dist()

2006-01-18 Thread Marco Geraci
see ?dist there's an example x - matrix(rnorm(100), nrow=5) m - as.matrix(dist(x)) d - as.dist(m) 'as.dist' is what you're probably looking for regards, Marco [EMAIL PROTECTED] wrote: Dear R-helpers, i am a beginner of R and i am using cluster

Re: [R] R-help Digest, Vol 35, Issue 14

2006-01-16 Thread Achim Zeileis
On Sun, 15 Jan 2006, Werner Wernersen wrote: Dear all, Is anybody aware of a tutorial, introduction, overview or alike for cluster analysis with R? I have been searching for something like that but it seems there are only a few rather specialized articles around. As an overview (rather

Re: [R] R-help Digest, Vol 35, Issue 14

2006-01-14 Thread Werner Wernersen
Dear all, Is anybody aware of a tutorial, introduction, overview or alike for cluster analysis with R? I have been searching for something like that but it seems there are only a few rather specialized articles around. I would very much appreciate any hint. Thanks a million, Werner

Re: [R] R-help Digest, Vol 35, Issue 14

2006-01-14 Thread Prof Brian Ripley
On Sun, 15 Jan 2006, Werner Wernersen wrote: Is anybody aware of a tutorial, introduction, overview or alike for cluster analysis with R? I have been searching for something like that but it seems there are only a few rather specialized articles around. Chapter 11 of MASS (the book

Re: [R] R-help Digest, Vol 35, Issue 7

2006-01-08 Thread Evgeniy Kachalin
Uwe Ligges пишет: Evgeniy Kachalin wrote: Hello, dear participants! Could you tip me, is there any simple and nice way to build scatter-plot for three different types of data (, and o and * - signs, for example) with legend. Now i can guess only that way:

Re: [R] R-help Digest, Vol 35, Issue 7

2006-01-08 Thread Kyosti H Kurikka
Hi! Just use your factors for indexing c(15,16,17) and c(red,green,blue). So, with the iris data: with(iris, plot(Sepal.Length, Sepal.Width, pch=c(15,16,17)[as.integer(Species)], col=c(red,green,blue)[as.integer(Species)] )) Best regards, Kyosti Kurikka Evgeniy Kachalin wrote:

Re: [R] R-help Digest, Vol 35, Issue 7

2006-01-08 Thread Uwe Ligges
Evgeniy Kachalin wrote: Uwe Ligges пишет: Evgeniy Kachalin wrote: Hello, dear participants! Could you tip me, is there any simple and nice way to build scatter-plot for three different types of data (, and o and * - signs, for example) with legend. Now i can guess only that way:

Re: [R] R-help Digest, Vol 35, Issue 7

2006-01-07 Thread Evgeniy Kachalin
Hello, dear participants! Could you tip me, is there any simple and nice way to build scatter-plot for three different types of data (, and o and * - signs, for example) with legend. Now i can guess only that way: plot(x~y,data=subset(mydata,factor1=='1'), pch='.',col='blue')

Re: [R] R-help Digest, Vol 35, Issue 7

2006-01-07 Thread Uwe Ligges
Evgeniy Kachalin wrote: Hello, dear participants! Could you tip me, is there any simple and nice way to build scatter-plot for three different types of data (, and o and * - signs, for example) with legend. Now i can guess only that way: plot(x~y,data=subset(mydata,factor1=='1'),

Re: [R] R-help: gls with correlation=corARMA

2005-12-12 Thread gaffigan
Thank you for your extra efforts in pinpointing the source of my problem. That is a smart workaround, to reduce the parameters until invertibility conditions are satisfied. The concern in my case is what effect the constant will have on the covariance matrix and the resulting estimates for the

Re: [R] R-help Digest, Vol 34, Issue 14

2005-12-11 Thread Dominik Schaub
Guten Tag, Ich bin vom 12. bis 23. Dezember 2005 im Militär-WK. Ich werde die Mails somit nur verzögert beantworten können. Für dringende Fälle: Während diesen zwei Wochen bin ich via Natel (am besten per SMS) erreichbar unter der Nummer 079 438 27 68. Mit freundlichem Gruss Dominik Schaub

Re: [R] R-help: gls with correlation=corARMA

2005-12-11 Thread Spencer Graves
The error message is misleading. It should say something like, Error in corARMA(q = 25, value = -ma.coefs, fixed = T) : The moving average process specified is not invertible, having roots outside the unit circle. Instead it says, Error in corARMA(q = 25, value = -ma.coefs, fixed =

Re: [R] R-help Digest, Vol 33, Issue 27

2005-11-27 Thread A.J. Rossini
From: Duncan Murdoch [EMAIL PROTECTED] I'd recommend using the RWinEdt package instead for a different way to integrate winedit with R. winedit and winedt are two different editors, last I checked. best, -tony [EMAIL PROTECTED] Muttenz, Switzerland. Commit early,commit often, and commit in

Re: [R] R-help: conversion of long decimal numbers into hexadeci

2005-11-10 Thread Ted Harding
On 10-Nov-05 Antje Döring wrote: Hi there, could somebody help me to convert a decimal number into a hexadecimal number? I know that there is the function sprintf, but the numbers I want to convert consist of 20 or more numbers. Spintf is not able to convert these big numbers. If I

Re: [R] R-help Digest, Vol 32, Issue 26

2005-10-26 Thread Alan Cobo-Lewis
r-help@stat.math.ethz.ch on Wednesday, October 26, 2005 at 6:00 AM -0500 wrote: Ronaldo, Try Harold's suggestion. The df still won't agree, because lmer (at least in its current version) just puts an upper bound on the df. But that should be OK, because all those t tests are approximations

Re: [R] R-help Digest, Vol 32, Issue 26

2005-10-26 Thread Doran, Harold
PROTECTED] On Behalf Of Alan Cobo-Lewis Sent: Wednesday, October 26, 2005 8:53 AM To: r-help@stat.math.ethz.ch Subject: Re: [R] R-help Digest, Vol 32, Issue 26 r-help@stat.math.ethz.ch on Wednesday, October 26, 2005 at 6:00 AM -0500 wrote: Ronaldo, Try Harold's suggestion. The df still won't agree

Re: [R] R-help Digest, Vol 31, Issue 30

2005-09-30 Thread John Maindonald
With lme4, use of mcmcsamp can be insightful. (Douglas Bates drew my attention to this function in a private exchange of emails.) The distributions of random effects are simulated on a log scale, where the distributions are much closer to symmetry than on the scale of the random effects

Re: [R] R-help Digest, Vol 31, Issue 9

2005-09-10 Thread Wuming Gong
?summary.lm and check the Value section. Wuming On 9/10/05, Ping Yao [EMAIL PROTECTED] wrote: Hi: I use lm (linear model) to analyze 47 variables , 8 responses So I use loop to finish it . I want the program to show the results that P-value is less than 0.05. How can I cite the P-valus from

Re: [R] R-help Digest, Vol 31, Issue 9

2005-09-10 Thread Wuming Gong
Hi Ping, You can use zz$coefficients[,4] to get the p values for each estimated coefficients in your context. Wuming On 9/11/05, Ping Yao [EMAIL PROTECTED] wrote: Wuming: Thanks for your help. I use the fuction: call(fstatistic,zz) call(p-value,zz) I can

Re: [R] R-help Digest, Vol 31, Issue 9

2005-09-09 Thread Ping Yao
Hi: I use lm (linear model) to analyze 47 variables , 8 responses So I use loop to finish it . I want the program to show the results that P-value is less than 0.05. How can I cite the P-valus from lm result ? Ping The code: #using LM to model general fati for (j in 48:52) { for (i in 3:46){

Re: [R] R-help Digest, Vol 30, Issue 26

2005-08-26 Thread Jean-Marc Ottorini
Dear R helpers, For me ( i.e. R 2.1.1 on Mac OS X), using trellis.device (postscript, onefile = F, etc ... with the lattice library within a R function works fine to obtain the desired graph as an EPS file , provided that : 1) the command dev.off() is not included in this function

Re: [R] R-help

2005-08-22 Thread Mike Waters
Diego, Have you checked out the home site for nls2? Specifically the system requirements page? http://www.inra.fr/miaj/public/AB/nls2/available.html That says that nls2 requires a Unix-like operating system. Basically, the script for building the library is for such systems only, it also

Re: [R] R-help Digest, Vol 30, Issue 6

2005-08-08 Thread David Duffy
On Fri, 5 Aug 2005 Julia Reid wrote: Subject: [R] GAP pointer I am trying to do a simple segregation analysis using the GAP package. I have the documentation for pointer but I desperately need an example so that I can see how to format the datfile and the jobfile. For each individual, I

Re: [R] R-help Digest, Vol 28, Issue 28

2005-06-28 Thread A. Mani
On Tuesday 28 June 2005 15:30, [EMAIL PROTECTED] wrote: Re : 37. Re: A. Mani : colours in Silhouette (Mulholland, Tom) Message: 37 Date: Tue, 28 Jun 2005 09:08:24 +0800 From: Mulholland, Tom [EMAIL PROTECTED] Subject: Re: [R] A. Mani : colours in Silhouette To: [EMAIL PROTECTED],

Re: [R] r-help

2005-06-25 Thread Uwe Ligges
Sotdikov Mansor wrote: There is a function 'simMD()' in 'popgen' library which package, not library. simulates a sample of genotype data as follows: library(popgen) x - simMD(20, 2, 2, p = NULL, c(0.09, 0.05), ac = 2, beta = 1) x , , 1 [,1] [,2] [1,]11 [2,]1

Re: [R] r-help

2005-06-25 Thread Peter Dalgaard
Sotdikov Mansor [EMAIL PROTECTED] writes: There is a function 'simMD()' in 'popgen' library which simulates a sample of genotype data as follows: library(popgen) x - simMD(20, 2, 2, p = NULL, c(0.09, 0.05), ac = 2, beta = 1) x ... How can I repeat this function, for example, 1000 times

Re: [R] R-help

2005-06-21 Thread Gabor Grothendieck
On 6/21/05, McClatchie, Sam (PIRSA-SARDI) [EMAIL PROTECTED] wrote: Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 - Colleagues Is there a function in R that is an equivalent of zoom in matlab? This is

Re: [R] R-help

2005-06-21 Thread Francisco J. Zagmutt
Try RSiteSearch(zoom plot). There are some good suggestions there. Cheers Francisco From: McClatchie, Sam (PIRSA-SARDI) [EMAIL PROTECTED] To: R-Help-Request (E-mail) r-help@stat.math.ethz.ch Subject: [R] R-help Date: Wed, 22 Jun 2005 09:05:05 +0930 Background: OS: Linux Mandrake 10.1

Re: [R] R-help Digest, Vol 28, Issue 11

2005-06-11 Thread dwfu
Dear all, I'm new using R and in (geo)statistics. I have a problem with solving my homework questions. We are working with variograms and trying to write down basic equations for different models (spherical, exponential, Gaussian). I tried to use the 'gstat' and 'geoR' packages to solve the

Re: [R] R-help Digest, Vol 28, Issue 11

2005-06-11 Thread Uwe Ligges
dwfu wrote: Dear all, I'm new using R and in (geo)statistics. I have a problem with solving my homework questions. We are working with variograms and trying to write down basic equations for different models (spherical, exponential, Gaussian). I tried to use the 'gstat' and 'geoR' packages

Re: [R] R-help,

2005-05-19 Thread Sean Davis
On May 19, 2005, at 7:01 AM, Jean Eid wrote: I do not fully understand your example but if you need to act on the columns of a dataframe why don't you just call its columns the title you want. something like X-matrix(rnorm(1000), ncol=4) colnames(X)-c(foo, foo1, foo2, foo3) X-as.data.frame(X)

Re: [R] R-help Digest

2005-05-06 Thread Prof Brian Ripley
On Fri, 6 May 2005, Sebastian Schoenherr wrote: Hi folks, I have to create my own time series, Is it possible to generate ARIMA time series, where i can define the range of the values in the y axis. (e.g: Values only between 0 and 1) No. Take a look at the definition of an ARIMA process. Suppose

Re: [R] [R-help] install problem

2005-03-22 Thread Brian D Ripley
On Mon, 21 Mar 2005, Tae-Young Goo wrote: Hello. I've tried to install R to IBM AIX(v.5.1) machine. I've used compile options indicated by R-admin. There are several sets there, so which exactly? Then, I met following error messages. /home/local/R_2.0.1/lib/R/bin/exec/R is unchanged

Re: [R] R-help Time Series

2005-02-26 Thread Prof Brian Ripley
I think 1) You have the units wrong: these appear to be the figures quoted for KB of compressed files, and the compression is nothing like 1024:1. 2) This is not `a series' unless you add a time base, e.g. via a call to ts(). Surely subscribers are aware that they do not get many MB/day and

Re: [R] R-help Time Series

2005-02-26 Thread Ted Harding
On 26-Feb-05 Prof Brian Ripley wrote: I think 1) You have the units wrong: these appear to be the figures quoted for KB of compressed files, and the compression is nothing like 1024:1. Sorry, yes, you are correct: it is KB and not MB (a slip of the eye on my part). 2) This is not `a

Re: [R] R-help Time Series

2005-02-26 Thread Gabor Grothendieck
Ted.Harding at nessie.mcc.ac.uk writes: : : Hi Folks, : : While I was browsing in the R-help archives yesterday, : I got curious about the time series of the sizes of : the monthly archives in MB. : : This turned out to have an unexpected feature or two, : which I leave to readers to explore

Re: [R] R-help Time Series

2005-02-26 Thread Uwe Ligges
Prof Brian Ripley wrote: I think 1) You have the units wrong: these appear to be the figures quoted for KB of compressed files, and the compression is nothing like 1024:1. 2) This is not `a series' unless you add a time base, e.g. via a call to ts(). Surely subscribers are aware that they do

Re: [R] R-help Time Series

2005-02-26 Thread Ted Harding
On 26-Feb-05 Uwe Ligges wrote: So let's be immensely unfair and do some speculation ... Assuming 1000MB/month means a compressed archive file of (very) *roughly* 250MB. Looking at the data with linear models, lm(sqrt(MB) ~ monthindex) seems not to be the worst model (removing the

RE: [R] R-help

2005-02-22 Thread abunn
This is what you are trying to get at I think: dev.off() ## start with a new graphics device # X11() or postscript() #par(mar = c(3,3,1,3), oma = c(0,0,0,0), mgp = c(2, 1, 0), bg = white) par(mar = rep(5,4)) plot(x-rnorm(100),y-rnorm(100), ylab = This is y, xlab = This is x) z-rnorm(100)*250

RE: [R] R-Help : running MIX package

2005-01-28 Thread Ted Harding
On 28-Jan-05 Jeanhee Hong wrote: Hello all. I am inexperienced with R and am clumsily trying to work through it for specific multiple imputations Id like to run for my thesis.In running the MIX package, I keep getting an error message regarding the use of the prelim.mix command. Error in

RE: [R] R-Help : running MIX package

2005-01-28 Thread Ted Harding
On 28-Jan-05 Ted Harding wrote: On 28-Jan-05 Jeanhee Hong wrote: Hello all. I am inexperienced with R and am clumsily trying to work through it for specific multiple imputations Id like to run for my thesis.In running the MIX package, I keep getting an error message regarding the use of the

Re: [R] R help search and Java(tm)?

2005-01-04 Thread Peter Dalgaard
Dan Bolser [EMAIL PROTECTED] writes: I found this great search source(http://www.statslab.cam.ac.uk/~djw1005/Stats/Interests/search.R;) helpHTML() Has (or will) this become the standard search method? I think it got superseded by Jon Baron's RSiteSearch() function which is finding its

Re: [R] R help

2004-11-16 Thread Spencer Graves
I don't understand your question. PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html;. In particular, have you read the R Data Import / Export documentation that some with R [available, from www.r-project.org - Manuals or R - help.start()]? Have you tried read.spss

Re: [R] R-help

2004-09-13 Thread Spencer Graves
Did you do a search at www.r-project.org - R site search? Searching for .dbf there just now exposed a read.shape function in the maptools package. hope this helps. p.s. Did you read the posting guide! http://www.R-project.org/posting-guide.html;? Tips provided there may help

Re: [R] R-help

2004-09-13 Thread Prof Brian Ripley
I believe .dbf files are more commonly DBase files, in which case see package RODBC. On Mon, 13 Sep 2004, Spencer Graves wrote: Did you do a search at www.r-project.org - R site search? Searching for .dbf there just now exposed a read.shape function in the maptools package.

Re: [R] R-help

2004-08-25 Thread Dimitris Rizopoulos
Hi Jim, you could try this: mat - matrix(sample(1:25), 5, 5) mat row(mat)[mat==max(mat)] col(mat)[mat==max(mat)] I hope it helps. Best, Dimitris Dimitris Rizopoulos Doctoral Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35,

Re: [R] R help in Firefox on Windows XP

2004-06-18 Thread Erich Neuwirth
Thank you Brian, Installing 0.8 first and then upgrading solved the problem. I noticed that installing 0.9 from scratch creates a registry key [HKEY_LOCAL_MACHINE\SOFTWARE\mozilla.org\Mozilla] CurrentVersion=1.7 Installing 0.8 creates [HKEY_LOCAL_MACHINE\SOFTWARE\mozilla.org\Mozilla]

Re: [R] R help in Firefox on Windows XP

2004-06-17 Thread Marc Schwartz
On Thu, 2004-06-17 at 12:06, Erich Neuwirth wrote: I had to reinstall my machine, so I installed Firefox 0.9 as browser I am using WinXP and R 1.9.1 beta. Now search in R html help does not work. I checked that the Java VM is working correctlt, Sun's test site says my installation is OK.

Re: [R] R help in Firefox on Windows XP

2004-06-17 Thread Damon Wischik
Erich Neuwirth wrote: I had to reinstall my machine, so I installed Firefox 0.9 as browser I am using WinXP and R 1.9.1 beta. Now search in R html help does not work. A workaround (which is slow, provisional, and as yet untested on your configuration) is to use a different help engine,

Re: [R] R help in Firefox on Windows XP

2004-06-17 Thread James MacDonald
Works for me with Firefox 0.8. James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623 Prof Brian Ripley [EMAIL PROTECTED] 06/17/04 02:10PM Works for me. Did it work with Firefox 0.8?

Re: [R] R-help

2004-05-14 Thread Gabor Grothendieck
If you want pairwise correlations for the cols of matrix x and N is some parameter to your routine that does not vary throughout the calculation then: apply(x,2,function(a)apply(x,2,function(b)Transmer.cor(a,b,N))) Replace both occurrences of 2 with 1 if you want rows. Ruedi Epple ruedi.epple

Re: [R] R-help

2004-02-12 Thread Rolf Turner
Jim Gustafsson [EMAIL PROTECTED] wrote: I have a problem. Boy, do you ever. I would like to put my SAS-code into R. Could I do that, if yes, how? Right. I have a couple of questions for you too. (1) How long is a piece of string? (2) How do I bring peace,

  1   2   >