Re: [R] Forcing best-fit lines to go through the origin

2005-04-18 Thread Henric Nilsson
Jim Milks said the following on 2005-04-18 03:08: Dear All, I have a rather unusual problem. I have a set of data for a class in subsurface processes. From that dataset, I must calculate the slope of the best-fit line (which is the parameter of interest). The problem I have is twofold: 1)

[R] lmer question

2005-04-18 Thread Markus Bean
Hi -- I'm using lmer for binomial data. I am trying to replicate estimates provided by Agresti (2002, Categorical data analysis, Wiley) using abortion data in table 10.13 (estimates provided in table 12.3 p. 505). I fit the same model using these three commands: a1 - lmer(resp ~ sex + option1

[R] how to use sink.number()

2005-04-18 Thread Michael S
Dear All R-helper, my question is how to use sink.number. if I want to label the output was sinked by sink function how should I do ? example: zz - file(c:\\sinktest1.txt,open=wt) sink(zz) g- 1:10 cat(g,\n) f-list(1:100) h-capture.output(f) cat(h) sink() I want there are marks in sinktest1.txt

[R] Mac GUI Slow?

2005-04-18 Thread Samuel E. Kemp
Hi, I have R version 2.0.1 on my PowerBook G4 1.5Ghz, 512RAM, running OS X 10.3.9. I have noticed that the official R GUI runs considerably slower than R in the terminal and also JGR. Whilst the differences in speed are unimportant for simple parametric stats, it becomes a bit of a pig when

[R] [R-pkgs] RWinEdt 1.7-0 released

2005-04-18 Thread Uwe Ligges
Dear useRs, RWinEdt_1.7-0 has been released. Windows users of R-2.1.0 (which is scheduled to be released today) who are going to use a translated (i.e. non-english) version of RGui together with WinEdt will have to upgrade. It is not advisable to upgrade RWinEdt for R versions 2.1.0. RWinEdt

[R] the graph gallery strikes back

2005-04-18 Thread Romain Francois
Hello useRs and helpRs, Some time ago, in a gallaxy far away (here is the thread : http://finzi.psych.upenn.edu/R/Rhelp02a/archive/46532.html ) we discussed about a graph gallery showing the power of R in that domain. I did some work around that, and there is a (pretty advanced) draft here :

[R] Rd.sty problems.

2005-04-18 Thread Federico Calboli
Hi All, I am trying to build a new R package to submit, but it's failing to create a tex manual: R CMD check Biodem * checking for working latex ... OK * using log directory '/home/greatsage/Fede/R-packages/temp/Biodem.Rcheck' * checking for file 'Biodem/DESCRIPTION' ... OK * checking if this is

[R] R2.0.1 for Mac OS X 10.3 problem

2005-04-18 Thread Greg Arnold
This combination was operating satisfactorily until I tried updating lme4 and Matrix. My attempts to do this ultimately broke R. The R console appears briefly, then collapses. I have tried downloading and reinstalling R without success. Typing 'R' into Terminal gives the error message

[R] Invalid *.Rd in package {was Rd.sty problems}

2005-04-18 Thread Martin Maechler
Federico == Federico Calboli [EMAIL PROTECTED] on Mon, 18 Apr 2005 10:33:22 +0100 writes: Federico I am trying to build a new R package to submit, Federico but it's failing to create a tex manual: Federico R CMD check Biodem Federico * checking for working latex ... OK

[R] Rcmdr's problem in East Asian Languages Version

2005-04-18 Thread ronggui
maybe this is a problem specific to the East Asian Languages Version.when i reinstall R without checking the Version for East Asian Languages option ,it works.but if i install R with the Version for East Asian Languages option, it show the following error msg. library(Rcmdr) Loading

Re: [R] how to use sink.number()

2005-04-18 Thread Prof Brian Ripley
On Mon, 18 Apr 2005, Michael S wrote: Dear All R-helper, my question is how to use sink.number. if I want to label the output was sinked by sink function how should I do ? example: zz - file(c:\\sinktest1.txt,open=wt) sink(zz) g- 1:10 cat(g,\n) f-list(1:100) h-capture.output(f) cat(h) sink() I

Re: [R] R2.0.1 for Mac OS X 10.3 problem

2005-04-18 Thread Prof Brian Ripley
It did not `break R': just rename your .RData which appears to have become corrupted. Please always start R with --vanilla before concluding it is broken. On Mon, 18 Apr 2005, Greg Arnold wrote: This combination was operating satisfactorily until I tried updating lme4 and Matrix. My attempts

Re: [R] the graph gallery strikes back

2005-04-18 Thread Gabor Grothendieck
On 4/18/05, Romain Francois [EMAIL PROTECTED] wrote: Hello useRs and helpRs, Some time ago, in a gallaxy far away (here is the thread : http://finzi.psych.upenn.edu/R/Rhelp02a/archive/46532.html ) we discussed about a graph gallery showing the power of R in that domain. I did some work

Re: [R] Rcmdr's problem in East Asian Languages Version

2005-04-18 Thread Prof Brian Ripley
The issue here is not the `East Asian Languages Version' but the locale you are using, which you have not told us. Try using the C locale or an English locale: Rcmdr works there for me. Also be sure to tell us the version of the package (here Rcmdr) you are using. The posting guide does ask

Re: [R] the graph gallery strikes back

2005-04-18 Thread Romain Francois
Oops :) That should be allright now. Romain Le 18.04.2005 13:52, Naji a crit : Hi Romain, Great job.. Thanks a lot. When pointing at the graph or the link (ie Kernel density estimator in R, Perspective plot and contour plot, etc) the browser is sending back a SQL error (Safari, MacOSX3.8) Erreur

Re: [R] the graph gallery strikes back

2005-04-18 Thread Romain Francois
Hello Hans, For the moment, you can't put yourself the graph in it. I'll do it for you. Just send me the source code and i'll put your graph when i found the time (probably tonight). You may also want to add comments on Andrew curves. Romain Le 18.04.2005 13:37, Hans Werner Borchers a écrit :

[R] R-2.1.0 is released

2005-04-18 Thread Peter Dalgaard
I've rolled up R-2.1.0.tar.gz a short while ago. This version contains several substantial changes and additions, perhaps most notably the UTF8 and Internationalization changes. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.1.0.tar.gz or wait

Re: [R] the graph gallery strikes back

2005-04-18 Thread Ales Ziberna
Great work! Nice collection! For me everything worked, although I did't click on all the graphs! Thans! - Original Message - From: Romain Francois [EMAIL PROTECTED] To: RHELP R-help@stat.math.ethz.ch; [EMAIL PROTECTED] Sent: Monday, April 18, 2005 11:29 AM Subject: [R] the graph gallery

Re: [R] the graph gallery strikes back

2005-04-18 Thread Romain Francois
Hello, I'm not sure i got you right, i added a link to index on the top right corner of each page. Is that what you meant ? Romain Le 18.04.2005 14:33, BXC (Bendix Carstensen) a écrit : On my windows system, I cannot get back to the overview page once I clicked on one of the items. Is this yours

Re: [R] lmer question

2005-04-18 Thread Douglas Bates
Markus Bean wrote: Hi -- I'm using lmer for binomial data. I am trying to replicate estimates provided by Agresti (2002, Categorical data analysis, Wiley) using abortion data in table 10.13 (estimates provided in table 12.3 p. 505). I fit the same model using these three commands: a1 - lmer(resp

Re: [R] the graph gallery strikes back

2005-04-18 Thread Gabor Grothendieck
Yes, it works for me now. Thanks. Very impressive I might add. On 4/18/05, Romain Francois [EMAIL PROTECTED] wrote: Oops :) That should be allright now. Romain Le 18.04.2005 13:52, Naji a écrit : Hi Romain, Great job.. Thanks a lot. When pointing at the graph or the link (ie

Re: [R] R2.0.1 for Mac OS X 10.3 problem

2005-04-18 Thread Horacio Montenegro
I have the same problem, in Windows, and I think the .Rdata is not corrupted. Load R without loading the.Rdata file - move or rename it. Then library(lme4) and load the .Rdata - it should work. cheers, Horacio -- Horacio Montenegro PhD Student in Genetics and

Re: [R] generalized linear mixed models - how to compare?

2005-04-18 Thread Peter Spreeuwenberg
Liset Dat moet lukken en hoe sneller je het anlevert hoe eerder je resultaten terug ziet. Dus, Ik wacht af. groet Peter S Date sent: Sun, 17 Apr 2005 18:07:28 +0100 (BST) From: Prof Brian Ripley [EMAIL PROTECTED] To: Deepayan Sarkar

[R] nls error in formula

2005-04-18 Thread Anaid Diaz
Hi, I'm a new R user, with a lot of questions. At the moment I'm stoped on an error traying to fit a model: x - sandeel ## numeric data (2500-6) y - Noss ## numeric data (0-1.2) A - 0.8 B - 0.6 C - 1/4 nls( y ~ A-B*exp(-C*x)) Error in match.call(definition, call, expand.dots) :

Re: [R] Generating a binomial random variable correlated with a

2005-04-18 Thread Giovanni Petris
You may generate a single standard normal random variable Z and set X = (Zx). According to back-of-the-envelope calculations, the two have a correlation of exp(-x^2/2)/sqrt(2*pi*pnorm(x)*(1-pnorm(x))) which goes from a maximum of about 0.79 at x=0 to 0 for x going to infinity. If you aim at

Re: [R] the graph gallery strikes back

2005-04-18 Thread Peter Dalgaard
Romain Francois [EMAIL PROTECTED] writes: Oops :) That should be allright now. Romain Yep. You might want to antialias the small images though. And, btw, just for fun you might want to include the real Maunga Whau among the volcano pictures:

Re: [R] nls error in formula

2005-04-18 Thread Douglas Bates
Anaid Diaz wrote: Hi, I'm a new R user, with a lot of questions. At the moment I'm stoped on an error traying to fit a model: x - sandeel ## numeric data (2500-6) y - Noss ## numeric data (0-1.2) A - 0.8 B - 0.6 C - 1/4 nls( y ~ A-B*exp(-C*x)) Error in match.call(definition, call,

Re: [R] nls error in formula

2005-04-18 Thread Anaid Diaz
thank you, it worked. Sylvia --- Douglas Bates [EMAIL PROTECTED] wrote: Anaid Diaz wrote: Hi, I'm a new R user, with a lot of questions. At the moment I'm stoped on an error traying to fit a model: x - sandeel ## numeric data (2500-6) y - Noss ## numeric data (0-1.2) A -

Re: [R] the graph gallery strikes back

2005-04-18 Thread Gabor Grothendieck
Just one suggestion. Maybe you could provide on a link on that page to http://www.r-project.org . On 4/18/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: Yes, it works for me now. Thanks. Very impressive I might add. On 4/18/05, Romain Francois [EMAIL PROTECTED] wrote: Oops :) That

[R] Barplot and colors for legend

2005-04-18 Thread Werner Wernersen
Hi all! One quick question: How do I get the standard colors used by barplot? I have some stacked bars and would like to add a horizontal legend via legend() but I don't know how to find the colors for the fills of the legend points. Thanks! Werner

Re: [R] Barplot and colors for legend

2005-04-18 Thread Uwe Ligges
Werner Wernersen wrote: Hi all! One quick question: How do I get the standard colors used by barplot? I have some stacked bars and would like to add a horizontal legend via legend() but I don't know how to find the colors for the fills of the legend points. Thanks! Werner Type

Re: [R] the graph gallery strikes back

2005-04-18 Thread Romain Francois
Hi, I don't know about that web page, what is it ? Romain Le 18.04.2005 16:33, Gabor Grothendieck a écrit : Just one suggestion. Maybe you could provide on a link on that page to http://www.r-project.org . On 4/18/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: Yes, it works for me now.

Re: [R] the graph gallery strikes back

2005-04-18 Thread Gabor Grothendieck
Its the R home page. I think the graph gallery might be used to direct people to R before they have used it since its quite attractive. At that point they will want to know how to get to the rest of R. On 4/18/05, Romain Francois [EMAIL PROTECTED] wrote: Hi, I don't know about that web page,

Re: [R] Barplot and colors for legend

2005-04-18 Thread Achim Zeileis
On Mon, 18 Apr 2005 16:48:42 +0200 Uwe Ligges wrote: Werner Wernersen wrote: Hi all! One quick question: How do I get the standard colors used by barplot? I have some stacked bars and would like to add a horizontal legend via legend() but I don't know how to find the colors

[R] colClasses = Date in read.delim, how to pass date-format?

2005-04-18 Thread Christoph Lehmann
Hi I have a huge data-set with one column being of type date. Of course I can import the data using this column as factor and then convert it later to dates, using: sws.bezuege$FaktDat - dates(as.character(sws.bezuege$FaktDat), format = c(dates = d.m.y)) But the

Re: [R] generalized linear mixed models - how to compare?

2005-04-18 Thread Spencer Graves
No puedo entender. Nicht versteh. Je ne comprend pas. Peter Spreeuwenberg wrote: Liset Dat moet lukken en hoe sneller je het anlevert hoe eerder je resultaten terug ziet. Dus, Ik wacht af. groet Peter S Date sent: Sun, 17 Apr 2005 18:07:28 +0100 (BST) From:

Re: [R] the graph gallery strikes back

2005-04-18 Thread Romain Francois
Gabor, I was obviously kidding. Do you really think I would build a gallery for R without knowing its website. I spend a lot of time around http://www.r-project.org/ The link has been added on the minute after i did read your mail. Romain. Le 18.04.2005 16:54, Gabor Grothendieck a écrit : Its

Re: [R] colClasses = Date in read.delim, how to pass date-format?

2005-04-18 Thread Gabor Grothendieck
On 4/18/05, Christoph Lehmann [EMAIL PROTECTED] wrote: Hi I have a huge data-set with one column being of type date. Of course I can import the data using this column as factor and then convert it later to dates, using: sws.bezuege$FaktDat - dates(as.character(sws.bezuege$FaktDat),

[R] refitting lm() with same x, different y

2005-04-18 Thread William Valdar
Dear All, Is there is a fast way of refitting lm() when the design matrix stays constant but the response is different? For example, y1 ~ X y2 ~ X y3 ~ X ...etc. where y1 is the 1st instance of the response vector. Calling lm() every time seems rather wasteful since the QR-decomposition of X

Re: [R] R2.0.1 for Mac OS X 10.3 problem

2005-04-18 Thread Jari Oksanen
On Mon, 2005-04-18 at 06:39 -0700, Horacio Montenegro wrote: I have the same problem, in Windows, and I think the .Rdata is not corrupted. Load R without loading the.Rdata file - move or rename it. Then library(lme4) and load the .Rdata - it should work. This also is my experience (in

[R] a small problem with data

2005-04-18 Thread Yoko Nakajima
Hello, I have a strange problem. I may have specified and scanned the data set (I receive a message: Read 384134 items). I can check the dimension of the data set by dim(xyz). But when I try to load the data set or to list the data set by data(xyz), then I receive an error message saying that

[R] Install problem on Solaris 9

2005-04-18 Thread Shields, Rusty (IMS)
I'm trying to install R-2.0.1 on Solaris 9 and I'm receiving the following error messages during make. begin installing recommended package foreign make[2]: *** [foreign.ts] Error 1 make[2]: Leaving directory `/opt/net/source/R-2.0.1/src/library/Recommended' make[1]: ***

RE: [R] Creating packages with windows (accessing data)

2005-04-18 Thread Alan T. Arnholt
Thank you, The DESCRIPTION is what I needed...all works as I want now. Alan- Alan T. Arnholt Associate Professor Department of Mathematical Sciences Appalachian State University Boone, NC 28608 (828) 262 2863 www1.appstate.edu/~arnholta -Original Message- From: Prof Brian Ripley

[R] Latent Class Analyisis in R

2005-04-18 Thread Dennis, Helen
Hi, I am new to R. I would like to know if there is a package available that will estimate a 2-class Latent Class Analysis, with trichotomous observed variables, some of which are free to covary. I have read the PDF files regarding the e1071 and gllm packages, but it appears that these only

[R] Retrieving column descriptions into Stata

2005-04-18 Thread AJ Bostian
Dear All, I am working on a project with Stata 7 datafiles. I'm trying to do some of my analysis in R since I don't have Stata handy. My problem is with accessing the Stata column descriptions. (Part of the project is to clean up these column descriptions.) I imported the Stata files ok

Re: [R] Latent Class Analyisis in R

2005-04-18 Thread Tim F Liao
Helen, You may consider using LCA 1.1, written by Niels Waller at Vanderbilt University. You may need to get it from his website at Vanderbilt instead of CRAN, at least it was the case in the past. Tim Liao Original message Date: Mon, 18 Apr 2005 12:29:41 -0400 From: Dennis, Helen

Re: [R] a small problem with data

2005-04-18 Thread Douglas Bates
Yoko Nakajima wrote: Hello, I have a strange problem. I may have specified and scanned the data set (I receive a message: Read 384134 items). I can check the dimension of the data set by dim(xyz). But when I try to load the data set or to list the data set by data(xyz), then I receive an error

Re: [R] colClasses = Date in read.delim, how to pass date-format?

2005-04-18 Thread Prof Brian Ripley
You are confusing class Date (part of R) with class dates (part of package chron). There is no as() method for class dates, so you can't do this. You can read the column as character (not factor) and convert later, but it sounds like the `huge amount of memory (and time)' is in fact taken by

Re: [R] Very Slow Gower Similarity Function

2005-04-18 Thread Jari Oksanen
On 18 Apr 2005, at 19:10, Tyler Smith wrote: Hello, I am a relatively new user of R. I have written a basic function to calculate the Gower similarity function. I was motivated to do so partly as an excercise in learning R, and partly because the existing option (vegdist in the vegan package)

Re: [R] refitting lm() with same x, different y

2005-04-18 Thread Prof Brian Ripley
William, As a first shot, use lm with a matrix response. That fits them all at once with one QR-decomposition. No analogue for glm or lmer, though, since for those the iterative fits run do depend on the response. Brian On Mon, 18 Apr 2005, William Valdar wrote: Dear All, Is there is a fast

[R] Storing vectors as vectors and iterating through them

2005-04-18 Thread Ken Termiso
Hi all, I have a bunch of int vectors. Each vector holds a bunch of ints that correspond to row numbers of an existing matrix. I use the int vectors to pull out rows of data from a matrix, i.e. data - my_matrix[int_vector,] I would like to store these int vectors in some sort of data structure

RE: [R] Storing vectors as vectors and iterating through them

2005-04-18 Thread Wiener, Matthew
Ken -- try using lists: vec.list - list(1:5, 2:4, 3:8) vec.list - c(vec.list, list(7:9)) vec.list [[1]] [1] 1 2 3 4 5 [[2]] [1] 2 3 4 [[3]] [1] 3 4 5 6 7 8 [[4]] [1] 7 8 9 Then you can use lapply or sapply, or just a for loop, to iterate over the list, applying your function to each

Re: [R] Very Slow Gower Similarity Function

2005-04-18 Thread Anon.
Jari Oksanen wrote: On 18 Apr 2005, at 19:10, Tyler Smith wrote: Hello, I am a relatively new user of R. I have written a basic function to calculate the Gower similarity function. I was motivated to do so partly as an excercise in learning R, and partly because the existing option (vegdist in

RE: [R] Pearson corelation and p-value for matrix

2005-04-18 Thread Dren Scott
Hi all, Thanks Andy, Mark and John for all the help. I really appreciate it. I'm new to both R and statistics, so please excuse any gaffes on my part. Essentially what I'm trying to do, is to evaluate for each row, how many other rows would have a p-value 0.05. So, after I get my N x N

Re: [R] colClasses = Date in read.delim, how to pass date-format?

2005-04-18 Thread Christoph Lehmann
so what do you recommend: I just need to be able to sort a data.frame according to the date entry, and e.g. compute differences between subsequent dates. Shall I stay with dates (thanks for the hint about confusion of Date and dates) or is there a better way for this kind of task? thanks a lot

Re: [R] colClasses = Date in read.delim, how to pass date-format?

2005-04-18 Thread Prof Brian Ripley
On Mon, 18 Apr 2005, Christoph Lehmann wrote: so what do you recommend: I just need to be able to sort a data.frame according to the date entry, and e.g. compute differences between subsequent dates. Shall I stay with dates (thanks for the hint about confusion of Date and dates) or is there a

[R] dynamic function loading

2005-04-18 Thread Omar Lakkis
How can I do dynamic function loading in R? I have a list of R functions in the database and want to, dynamicly, execute teh function that corresponds to my in comming data. Is there a way to do that without a big if/else? __ R-help@stat.math.ethz.ch

RE: [R] Generating a binomial random variable correlated with a

2005-04-18 Thread Mohammad A. Chaudhary
Dear: Thank you very much for your 'generous' help. Overall I have learnt the limits of the solution to my problem. I was not feeling compfortable how I was going around with the problem. Bill and Ted provided a wonderful insight into what I was accomplishing and provided alternative solutions.

[R] R 2.1.0 GUI language

2005-04-18 Thread Stephan Tolksdorf
Hi How do I set the language of the GUI? With R 2.1.0 it suddenly changed to German and I want it to stay English. I couldn't find any setting to change. Regards, Stephan __ R-help@stat.math.ethz.ch mailing list

[R] when can we expect Prof Tierney's compiled R?

2005-04-18 Thread Jason Liao
I am excited to learn that Prof. Tierney is bringing to us compiled R. I would like to learn when it will be available. This information will be useful in scheduling some of my projects. Thanks. Jason Jason Liao, http://www.geocities.com/jg_liao Dept. of Biostatistics,

Re: [R] Very Slow Gower Similarity Function

2005-04-18 Thread Jari Oksanen
On 18 Apr 2005, at 20:36, Anon. wrote: Jari Oksanen wrote: On 18 Apr 2005, at 19:10, Tyler Smith wrote: Hello, I am a relatively new user of R. I have written a basic function to calculate the Gower similarity function. I was motivated to do so partly as an excercise in learning R, and partly

RE: [R] Pearson corelation and p-value for matrix

2005-04-18 Thread John Fox
Dear Dren, Since cor(), on which Andy's solution is based, can compute pairwise-present correlations, you could adapt his function -- you'll have to adjust the df for each pair. Alternatively, you could probably save some time (programming time + computer time) by just using my solution: R -

Re: [R] R 2.1.0 GUI language

2005-04-18 Thread Prof Brian Ripley
On Mon, 18 Apr 2005, Stephan Tolksdorf wrote: How do I set the language of the GUI? With R 2.1.0 it suddenly changed to German and I want it to stay English. I couldn't find any setting to change. Which OS is this? PLEASE do read the posting guide. This *is* documented in the `R Installation and

Re: [R] Very Slow Gower Similarity Function

2005-04-18 Thread Martin Maechler
Tyler == Tyler Smith [EMAIL PROTECTED] on Mon, 18 Apr 2005 12:10:34 -0400 writes: Tyler Hello, I am a relatively new user of R. I have Tyler written a basic function to calculate the Gower Tyler similarity function. I was motivated to do so partly Tyler as an excercise in

RE: [R] when can we expect Prof Tierney's compiled R?

2005-04-18 Thread Liaw, Andy
Do you mean the byte code compiler? You can find it at: http://www.stat.uiowa.edu/~luke/R/compiler/ Andy From: Jason Liao I am excited to learn that Prof. Tierney is bringing to us compiled R. I would like to learn when it will be available. This information will be useful in scheduling

RE: [R] Pearson corelation and p-value for matrix

2005-04-18 Thread Liaw, Andy
I believe this will do: cor.pval2 - function(x, alternative=two-sided) { corMat - cor(x, use=if (any(is.na(x))) pairwise.complete.obs else all) df - crossprod(!is.na(x)) - 2 STATISTIC - sqrt(df) * corMat / sqrt(1 - corMat^2) p - pt(STATISTIC, df) p - if (alternative == less)

Re: [R] R 2.1.0 GUI language

2005-04-18 Thread Stephan Tolksdorf
Thanks for your prompt reply. Which OS is this? PLEASE do read the posting guide. Win XP, SP2 (German), which I forgot to mention despite heaving read the guide. This *is* documented in the `R Installation and Adminstration' Manual, section `Internationalization'. Looked in the wrong place

Re: [R] Very Slow Gower Similarity Function

2005-04-18 Thread Tyler Smith
Quoting Martin Maechler [EMAIL PROTECTED]: I don't know what exactly you want. The Gower coefficient I am referring to comes from his 1971 article in Biometrics (27(4):857-871). It differs from most commonly used measures (but not, apparently, daisy!) by allowing the incorporation of

[R] Construction of a large sparse matrix

2005-04-18 Thread Doran, Harold
Dear List: I'm working to construct a very large sparse matrix and have found relief using the SparseM package. I have encountered an issue that is confusing to me and wonder if anyone may be able to suggest a smarter solution. The matrix I'm creating is a covariance matrix for a larger research

Re: [R] R 2.1.0 GUI language

2005-04-18 Thread Prof Brian Ripley
On Mon, 18 Apr 2005, Stephan Tolksdorf wrote: [...] LC_ALL doesn't seem to change the message language. LC_MESSAGES, which is mentioned in the installation and administration guide, doesn't seem to exist in R 2.1.0 and isn't listed in the locales help. Right: it cannot be set inside R. For as

RE: [R] Construction of a large sparse matrix

2005-04-18 Thread Huntsinger, Reid
Your statements for sample.size=10,000 try to construct a matrix with 3 dense 10,000 x 10,000 blocks. That's approximately 10*10*8 MB each and very likely explains your error message. Since you have a simple formula for the matrix, why not define a function to implement multiplication by this

Re: [R] Construction of a large sparse matrix

2005-04-18 Thread roger koenker
The dense blocks are too big as Reid has already written -- for smaller instances of this sort of thing I would suggest that the the kronecker product %x% operator in SparseM, would be more convenient. url:www.econ.uiuc.edu/~rogerRoger Koenker email [EMAIL PROTECTED]

[R] R-squared in summary(lm...)

2005-04-18 Thread Benjamin M. Osborne
What is the difference between the two R-squareds returned for a linear regression by summary(lm...)? When might one report multiple vs. adjusted R-squared? Thank you, Ben Osborne -- Botany Department University of Vermont 109 Carrigan Drive Burlington, VT 05405 [EMAIL PROTECTED] phone:

RE: [R] R-squared in summary(lm...)

2005-04-18 Thread Berton Gunter
Please consult any basic book on linear regression/linear models; for example, APPLIED REGRESSION ANALYSIS by Draper and Smith. Or google on adjusted R-Squared. In brief, Adjusted R-squared arrempts to compensate for the property that adding extra regressors -- even random ones -- to a linear

[R] Discrepancy between gam from gam package and gam in S-PLUS

2005-04-18 Thread John Fox
Dear Trevor, I've noticed a discrepancy in the degrees of freedom reported by gam() from the gam package in R vs. gam() in S-PLUS. The nonparametric df differ by 1; otherwise (except for things that depend upon the df), the output is the same: - snip *** From R (gam

RE: [R] Pearson corelation and p-value for matrix

2005-04-18 Thread John Fox
Dear Andy, Very nice! (My point was that if this is a one-time thing, for Dren to puzzle over it is probably more time-consuming than simply doing it inefficiently.) Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S

[R] longer object length, is not a multiple of shorter object length in: kappa * gcounts

2005-04-18 Thread Hui Han
Hi, I was using a density estimation function as follows: est - KernSmooth::bkde(x3, bandwidth=10) When setting bandwidth less than 5, I got the error longer object length, is not a multiple of shorter object length in: kappa * gcounts . I wonder if there is anybody who can explain the error

[R] Re: longer object length, is not a multiple of shorter object length in: kappa * gcounts

2005-04-18 Thread Hui Han
Hi, Another phenomenon is that if I don't sort the data vector x3, I get the same error regardless of the bandwidth value: longer object length, is not a multiple of shorter object length in: kappa * gcounts . Thanks in advance for any help you can give me! Hui Hui Han wrote: Hi, I was using

[R] Re: longer object length, is not a multiple of shorter object length in: kappa * gcounts

2005-04-18 Thread Hui Han
Hui Han wrote: Hi, Another phenomenon is that if I don't sort the data vector x3, I get the same error regardless of the bandwidth value: Actually for really large bandwidth,e.g. 120, it deosn't report errors. Sorry for the confusing description. longer object length, is not a multiple

[R] Re-release of R 2.1.0

2005-04-18 Thread Peter Dalgaard
Due to a mishap in the build procedure, the originally released sources erroneously had 2.2.0 in several places within the configure file. A re-released version has been put up in http://biostat.ku.dk/~pd/R-release/ and should find its way to CRAN and its mirrors some time tomorrow. The new

[R] standard normal gaussian

2005-04-18 Thread simone gabbriellini
Hello List, I need to calculate the p-value (in a standard normal gaussian) for a Z value I have. I guess if I had to calculate by myself the value of the area in my function or refer to some already made function in R. Any hints? Thank you very much, Simone

RE: [R] Pearson corelation and p-value for matrix

2005-04-18 Thread Liaw, Andy
From: John Fox Dear Andy, Very nice! (My point was that if this is a one-time thing, for Dren to puzzle over it is probably more time-consuming than simply doing it inefficiently.) There's a memory/speed tradeoff. The loop avoids creating multiple 6000x6000 matrices that my version

RE: [R] longer object length, is not a multiple of shorter object length in: kappa * gcounts

2005-04-18 Thread Liaw, Andy
I suspect you set the bandwidth too small. Consider: x - sort(runif(100)) * 100 min(diff(x)) [1] 0.001700642 library(KernSmooth) KernSmooth 2.22 installed Copyright M. P. Wand 1997 dx - bkde(x, bandwidth=.001) Warning message: longer object length is not a multiple of shorter object

Re: [R] standard normal gaussian

2005-04-18 Thread Spencer Graves
?pnorm For more information, try help.start() - An Introduction to R - Probability distributions. spencer graves simone gabbriellini wrote: Hello List, I need to calculate the p-value (in a standard normal gaussian) for a Z value I have. I guess if I had to calculate by myself

Re: [R] longer object length, is not a multiple of shorter object length in: kappa * gcounts

2005-04-18 Thread Hui Han
I see. Thanks, Andy! Hui Liaw, Andy wrote: I suspect you set the bandwidth too small. Consider: x - sort(runif(100)) * 100 min(diff(x)) [1] 0.001700642 library(KernSmooth) KernSmooth 2.22 installed Copyright M. P. Wand 1997 dx - bkde(x, bandwidth=.001) Warning message:

RE: [R] Barplot and colors for legend

2005-04-18 Thread Mulholland, Tom
Well I don't know how I can live with myself. I guess I can't wait for the site to mirror itself in case someone thinks I'm yesterday's man. ;-) -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Achim Zeileis Sent: Monday, 18 April 2005 10:56 PM ...

[R] Compatability with Tiger OS X 10.4

2005-04-18 Thread Hemant Ishwaran
Does anyone know if R will install and run properly on the soon to be released Tiger operating system (Mac OS X 10.4)? There must be several Mac users out there (like myself) who will be considering the upgrade. -- Hemant Ishwaran Staff, Dept of Quantitative Health Sciences, Cleveland Clinic

[R] Listas de E-mails para Mala Direta , Marketing Direto

2005-04-18 Thread marcosjbd_g5a
Cadastros para mala direta livre de spam programas grátis para e-mail marketing Cadastros para mala direta livre de spam e-mails para mala direta por classe social: http://www.gueb.de/divulgamail e-mails para mala direta por classe social Listas de E-mails para Mala Direta , Marketing Direto

Re: [R] Compatability with Tiger OS X 10.4

2005-04-18 Thread David L. Van Brunt, Ph.D.
don't know about the binaries, but I'd be surprised if you couldn't compile it. You've given me something to think about, though, as to whether I want to upgrade on my analytic box. On 4/18/05, Hemant Ishwaran [EMAIL PROTECTED] wrote: Does anyone know if R will install and run properly on

[R] using imported tables

2005-04-18 Thread Owen Buchner
I've recently found out using read.table i can insert a table into R from a .dat file. The problem i'm having is now that i've got R to read the package i want to use the data from each column. I've tried adding titles to each column using col.name, but it assigns a name to the entire table and

[R] Odd diagnostic plots in mixed-effects models

2005-04-18 Thread Andrew Robinson
Dear R community, In the excellent nlme package the default diagnostic plot graphs the innermost residuals against innermost fitted values. I recently fit a mixed-effects model in which there was a very clear positive linear trend in this plot. I inferred that this trend occurred because my

Re: [R] Compatability with Tiger OS X 10.4

2005-04-18 Thread Prof Brian Ripley
No one _can_ know. Mac OS 10.4 is unreleased, and things can easily break before release. It is working on pre-releases. One problem is that Mac OS 10.4 (don't think you need two tens) will use gcc 4.0.0, which is also unreleased and its release candidates have problems compiling R (see the