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 column

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

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? > > Re

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 th

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 l

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 he

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 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 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 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(

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 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 u

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 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 fa

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

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 (http://www.r-proj

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 conside

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 havi

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

2006-05-16 Thread Ben Bolker
Xin 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 dele

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 mo

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 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 th

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 AM

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: Mon

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)) d2<-data.frame(Date=c("2005/1/3"

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 Jan-

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 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 > ot

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 pos

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

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

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
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

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 18

2006-02-18 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 17

2006-02-17 Thread Jeffrey Barrett
I will be out of the office from Friday (Feb 17) through Wednesday (Feb 22). __ 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 - Problem in drawing braplot with a huge value of data

2006-02-15 Thread Marc Schwartz
On Thu, 2006-02-16 at 09:20 +0530, Subhabrata wrote: > Hello R-experts, > > I am facing a strange problem while creating a barplot. I have serise of > data of which the first > on is around 162589 while the remaining data are around 0-150. so when I am > ploting the barplot > with all the data I a

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, getOption("digi

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/ http://w

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 membe

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 else

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 w

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 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 base

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 th

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 Wissenscha

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, 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 distance

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 (

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 d

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 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 ca

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 K

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: >> >> plot(x~y,

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==

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') points(x~y,data=s

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: 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 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 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 commi

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

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

2005-10-26 Thread Doran, Harold
ED] [mailto:[EMAIL 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. T

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 anywa

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 themselves

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) >

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-valu

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){ g

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 depends

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

2005-08-22 Thread A Mani
Re: A. Mani : Avoiding loops (Petr Pikal) > Message: 9 > Date: Mon, 22 Aug 2005 06:40:45 +0200 > From: "Petr Pikal" <[EMAIL PROTECTED]> > Subject: Re: [R] A. Mani : Avoiding loops > To: "A. Mani" <[EMAIL PROTECTED]>, r-help > > > On 20 Aug 2005 at 3:26, A. Mani wrote: > > > On Friday 19 August

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

2005-08-07 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,

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 PROTECTE

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 exam

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,]1

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)" Subject: [R] R-help Date: Wed, 22 Jun 2005 09:05:05 +0930 Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor

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? T

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' packag

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

2005-06-11 Thread Ruben Roa
> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] Behalf Of dwfu > Sent: 11 June 2005 13:59 > To: r-help@stat.math.ethz.ch > Subject: Re: [R] R-help Digest, Vol 28, Issue 11 > > > Dear all, > I'm new using R and in (geo)statis

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 questi

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.fr

Re: [R] R-help,

2005-05-19 Thread Jean Eid
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) par(mfrow=c(2,2)) lapply(colnames(X), f

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 e

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 unchange

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

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 no

Re: [R] R-help Time Series

2005-02-26 Thread Gabor Grothendieck
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 for themselves

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 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 tha

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

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

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. > > Err

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 fin

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 "re

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 packag

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 ma

RE: [R] R-help

2004-08-25 Thread Wayne Jones
Or you could use: mat <- matrix(sample(1:25), 5, 5) mat which(mat==max(mat),arr.ind = T) Wayne -Original Message- From: Jim Gustafsson [mailto:[EMAIL PROTECTED] Sent: 25 August 2004 10:07 To: [EMAIL PROTECTED] Subject: [R] R-help Dear R users, I have just start working with R

  1   2   >