Re: [R] Calling LAPACK functions directly from R

2007-04-12 Thread Prof Brian Ripley
This was much more appropriate for R-devel, so please move any followup there. On Wed, 11 Apr 2007, William Constantine wrote: I am interested in tapping into LAPACK functionality directly from R. Using R-2.4.1 for Windows, I was able to do so ala: dyn.load(bin/Rlapack.dll)

Re: [R] Looping through series of names

2007-04-12 Thread Petr PIKAL
Hallo [EMAIL PROTECTED] napsal dne 11.04.2007 22:45:19: Hi I am very new to R and have not able to find the answer to my question in the manual or any other post. I have a dataset that has various different sites names with data relating to each site. The data is in one dataset. I want

[R] taskPR

2007-04-12 Thread Benno Willemsen
Hello I tried to load the taskPR-package in R, but it doesn't work at all. Everytime I try /library(taskPR)/, I get the error-message: Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library '/home/bw135690/R-2.4.1/library/taskPR/libs/taskPR.so':

[R] FFT output and phase

2007-04-12 Thread #HU GUANGAN#
I'm doing the time-series data analysis using FFT. After fft(), I used spectrum() to calculate the power spectrum, which is used to define the main frequency. Now I want to get the phase information according to this frequency. For example, input 12 time-point data,

[R] Assignment from list

2007-04-12 Thread Eric Blanc
I have a list of groups of xy positions I want to set to 0 in an array full of 1s. When the assignments are done directly from the list, they are incorrect, while if I use a temporary array derived from the list the assignments are correct. The following example will hopefully make my

Re: [R] taskPR

2007-04-12 Thread Prof Brian Ripley
It would have been helpful to have given the details the R posting guide asks for, such as your OS. That message probably does not mean that /home/bw135690/R-2.4.1/library/taskPR/libs/taskPR.so does not exist. It says 'cannot open shared object file', and this can have many causes. One is that

[R] Tweaking my plot of matrix as image

2007-04-12 Thread Wee-Jin Goh
Greetings list, I have a rectangular 20 x 20 matrix with values in the range of [0 , 1]. I'd like to plot it as an image. To that end, I have used the image() function that seems to do what I want. Now, I just want to tweak it to look perfect. So here is my question: At the moment, the

Re: [R] Assignment from list

2007-04-12 Thread Dimitris Rizopoulos
try the following: z - matrix(1, ncol = 249, nrow = 240) zz - matrix(1, ncol = 249, nrow = 240) for (k in seq_along(x)) { z[rbind(x[[k]])] - 0 } zz[xx] - 0 all.equal(z, zz) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public

[R] LME: incompatible formulas for groups

2007-04-12 Thread Jan Thiele
Dear R-Users, I am currently working with LME to analyse repeated measures data. I encounter a problem when including both a random effect and a correlation structure with different grouping levels into the LME model. The error message is: Error in lme.formula(diameter ~ flowers*timef +

Re: [R] Tweaking my plot of matrix as image

2007-04-12 Thread Petr PIKAL
Hi Did you by chance look at the help page of image? If you did, you could read x,y locations of grid lines at which the values in z are measured. These must be finite, non-missing and in (strictly) ascending order. By default, equally spaced values from 0 to 1 are used. If x is a list, its

Re: [R] Reasons to Use R

2007-04-12 Thread Jim Lemon
Charilaos Skiadas wrote: A new fortune candidate perhaps? On Apr 10, 2007, at 6:27 PM, Greg Snow wrote: Remember, everything is better than everything else given the right comparison. Only if we remove the grammatical blip that turns it into an infinite regress, i.e. Remember, anything

Re: [R] Tweaking my plot of matrix as image

2007-04-12 Thread Jim Lemon
Wee-Jin Goh wrote: Greetings list, I have a rectangular 20 x 20 matrix with values in the range of [0 , 1]. I'd like to plot it as an image. To that end, I have used the image() function that seems to do what I want. Now, I just want to tweak it to look perfect. So here is my

[R] reporting multinomial logistic regression results

2007-04-12 Thread Bob Green
I am hoping to obtain information regarding the reporting of multinomial logistic regression results - either recommendations of published studies employing R or suggestions about what output should be reported. I am used to SPSS which routinely produces a broader range of output than R (e.g

[R] using install.packages via proxy

2007-04-12 Thread Erbe, Nicole
Hello everybody, I'm trying to install rJava to use JRI in a Java program as JRI is a part of rJava. There should be the option to install this package via install.packages(rJava). Unfortunately I'm connected to the internet via a proxy. I have no idea how to tell R to connect via the proxy.

Re: [R] Regression on two time series

2007-04-12 Thread rdporto1
Ron, you're right. It's not legitimate at all. I suggest you to take a look at the HUGE bibliography on cointegration, as a start up. Rogerio Dear all R user, Please forgive me if my question is too simple. My question is related to Statistics rather directly to R. Suppose I have

Re: [R] Positioning in xyplot

2007-04-12 Thread Michael Kubovy
Wow! Thanks to Sundar and to Deepayan for this selfless effort.(I rather suspect they had a good time figuring this out.) The outcome is that I will be able to publish just the right figure for the job. I will be sure to aknowledge your contribtution. Michael On Apr 12, 2007, at 12:48 AM,

Re: [R] using install.packages via proxy

2007-04-12 Thread Vladimir Eremeev
What operating system do you use? If you use Windows, than open the Control Panel, double click on the System icon, go to Advanced tab, press Environment variables button, then press New to create one more new variable, enter http_proxy as the Name of variable, and

[R] data file import - numbers and letters in a matrix(!)

2007-04-12 Thread Felix Wave
Hello, I have a problem with the import of a date file. I seems verry tricky. I have a text file (end of the mail). Every file has a different number of measurments witch start with START OF HEIGHT DATA and ende with END OF HEIGHT DATA. I imported the file in a matrix but the letters before the

Re: [R] using install.packages via proxy

2007-04-12 Thread Richard Gott
I had that problme and ended up downloading the package and installing locally. there is probably a better way Richard Erbe, Nicole wrote: Hello everybody, I'm trying to install rJava to use JRI in a Java program as JRI is a part of rJava. There should be the option to install this package

Re: [R] using install.packages via proxy

2007-04-12 Thread Erbe, Nicole
Hi, I'm (currently forced to) use Win XP. I actually found an easy way to solve the problem of installing the rJava package. To avoid the problems with the proxy I just downloaded the zip file and installed the package in R from a local source. Now I can use the JRI libraries in eclipse by

Re: [R] using install.packages via proxy

2007-04-12 Thread Prof Brian Ripley
On Thu, 12 Apr 2007, Vladimir Eremeev wrote: What operating system do you use? If you use Windows, than open the Control Panel, double click on the System icon, go to Advanced tab, press Environment variables button, then press New to create one more new variable, enter http_proxy as the

Re: [R] Reasons to Use R

2007-04-12 Thread Joel J. Adamson
Douglas Bates writes: One can do data analysis by using the computer as a blunt instrument with which to bludgeon the problem to death but one can't do elegant data analysis like that. One nice thing about a blunt instrument like Stata is the ability to hold an entire dataset in memory

Re: [R] Reasons to Use R

2007-04-12 Thread Lucke, Joseph F
A re-interpretation of Zorn's lemma? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon Sent: Thursday, April 12, 2007 5:14 AM To: [EMAIL PROTECTED] Subject: Re: [R] Reasons to Use R Charilaos Skiadas wrote: A new fortune candidate perhaps?

[R] pattern

2007-04-12 Thread Schmitt, Corinna
Dear R-experts, I have the following function: userInput - function() { ANSWER - readline(saving place of the data (example Z:/Software/test.mat)? ) x = c(.mat) endingTest = x %in% ANSWER print (endingTest) if (endingTest == ??) cat (saving place not

Re: [R] graphs superimposed on pictures?

2007-04-12 Thread Michael Dewey
At 14:39 11/04/2007, Robert Biddle wrote: Hi: I am doing some work that involves plotting points of interest superimposed on photographs and maps. I can produce the plots fine in R, but so far I have had to do the superimposition externally, which makes it tedious to do exploratory work. I have

[R] Method dispatch for print() in package its

2007-04-12 Thread Albrecht, Dr. Stefan (AZ Private Equity Partner)
Dear all, in the package its the print() method does not seem to correctly work in all circumstances: selectMethod(print, its) Method Definition: function (x, ...) { print([EMAIL PROTECTED] mailto:[EMAIL PROTECTED] , ...) } environment: namespace:its Signatures: x

Re: [R] data file import - numbers and letters in a matrix(!)

2007-04-12 Thread Gabor Grothendieck
Try pasting this into an R session: Lines.raw - FILEDATE:02.02.2007 ... START OF HEIGHT DATA S= 0 y=0.0 x=0. S= 0 y=0.1 x=0.00055643 ... S= 9 y=4.9 x=1.67278117 S= 9 y=5.0 x=1.74873257 S=10 y=0.0 x=0. S=10 y=0.1 x=0.00075557 ... S=99 y=5.3 x=1.94719490 END OF HEIGHT DATA ...

Re: [R] Method dispatch for print() in package its

2007-04-12 Thread Gabor Grothendieck
Without a reproducible example I can't tell what the problem is but as a workaround you could print it as a zoo object: library(its) library(zoo) # create an its object example(its) ii - its(mat, times) # print it without using its' print print(as.zoo(ii)) # or from the command line just

Re: [R] Method dispatch for print() in package its

2007-04-12 Thread Seth Falcon
Albrecht, Dr. Stefan (AZ Private Equity Partner) [EMAIL PROTECTED] writes: Dear all, in the package its the print() method does not seem to correctly work in all circumstances: My understanding is that one should define only a method for the show generic for S4 classes and leave print

Re: [R] Reasons to Use R

2007-04-12 Thread Joel J. Adamson
Lucke, Joseph F writes: A re-interpretation of Zorn's lemma? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon Sent: Thursday, April 12, 2007 5:14 AM To: [EMAIL PROTECTED] Subject: Re: [R] Reasons to Use R Charilaos Skiadas

Re: [R] Memory management

2007-04-12 Thread yoooooo
Okay thanks, I'm going through the docs now.. and I came through this.. The named field is set and accessed by the SET_NAMED and NAMED macros, and take values 0, 1 and 2. R has a `call by value' illusion, so an assignment like b - a appears to make a copy of a and refer to it as b.

Re: [R] Method dispatch for print() in package its

2007-04-12 Thread Prof Brian Ripley
Setting S4 methods for print() works correctly, as in your examples. Your first example does not call print(), so it is not suprising that no method dispatch takes place. You seem to suppose that such a method should be used for top-level auto-printing, and it will not be as your example

[R] How to manipulate the pointer of a file?

2007-04-12 Thread Yuchen Luo
Dear friends. With file( ) to obtain a pointer of a file, every time we use scan ( ) to read one row of it, the pointer will point to the next row of the file. In the following example, d1 and d2 are obtained the same way but they correspond to different rows of the same file because the pointer

Re: [R] data file import - numbers and letters in a matrix(!)

2007-04-12 Thread Adaikalavan Ramasamy
Here is the contents of my testdata.txt : - START OF HEIGHT DATA S= 0y=0.0 x=0. S= 0 y=0.1 x=0.00055643 S= 9 y=4.9 x=1.67278117 S= 9 y=5.0 x=1.74873257 S=10 y=0.0 x=0. S=10y=0.1 x=0.00075557 S=99

Re: [R] how to reverse a list

2007-04-12 Thread Weiwei Shi
Thanks all of you! But my next question is, how to improve R programming skills? I never have time in improving it but I feel I need to. Regards, W On 4/11/07, Seth Falcon [EMAIL PROTECTED] wrote: Weiwei Shi [EMAIL PROTECTED] writes: I forgot to add my bad solution here: reverseList -

Re: [R] pattern

2007-04-12 Thread Uwe Ligges
Schmitt, Corinna wrote: Dear R-experts, I have the following function: userInput - function() { ANSWER - readline(saving place of the data (example Z:/Software/test.mat)? ) x = c(.mat) endingTest = x %in% ANSWER print (endingTest) if (endingTest

[R] zoo merge() method

2007-04-12 Thread Matthew.Findley
R users: I'd like to get some insight on an error I encounter when attempting to work with two moderately sized sets of time series data. FYI - I'm using the following versions of R and supporting packages on a Windows 2000 OS: - R version 2.4.1 (2006-12-18) - zoo version 1.2-2 - chron version

Re: [R] How to manipulate the pointer of a file?

2007-04-12 Thread Prof Brian Ripley
On Thu, 12 Apr 2007, Yuchen Luo wrote: With file( ) to obtain a pointer of a file, every time we use scan ( ) to read one row of it, the pointer will point to the next row of the file. In the following example, d1 and d2 are obtained the same way but they correspond to different rows of the

Re: [R] How to manipulate the pointer of a file?

2007-04-12 Thread Greg Snow
?seek -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Yuchen Luo Sent: Thursday, April 12, 2007 9:33 AM To: [EMAIL PROTECTED]

[R] Subsetting list of vectors with list of (boolean) vectors?

2007-04-12 Thread Johannes Graumann
Dear Rologists, I'm stuck with this. How would you do this efficiently: aPGI [[1]] [1] 864 5576 aPGItest [[1]] [1] TRUE FALSE result - [magic box involving subset) result [[1]] [1] 864 Thanks for any hints, Joh __ [EMAIL PROTECTED] mailing

Re: [R] zoo merge() method

2007-04-12 Thread Achim Zeileis
On Thu, 12 Apr 2007 [EMAIL PROTECTED] wrote: attempted to get at this by merging the two time series (by union), interpolating the NAs, and finally, subtracting one vector from the other. The problem is that I can not combine the two zoo time series objects using the merge() or cbind()

Re: [R] Subsetting list of vectors with list of (boolean) vectors?

2007-04-12 Thread Marc Schwartz
On Thu, 2007-04-12 at 18:12 +0200, Johannes Graumann wrote: Dear Rologists, I'm stuck with this. How would you do this efficiently: aPGI [[1]] [1] 864 5576 aPGItest [[1]] [1] TRUE FALSE result - [magic box involving subset) result [[1]] [1] 864 Thanks for any

Re: [R] Subsetting list of vectors with list of (boolean) vectors?

2007-04-12 Thread Seth Falcon
Johannes Graumann [EMAIL PROTECTED] writes: Dear Rologists, I'm stuck with this. How would you do this efficiently: aPGI [[1]] [1] 864 5576 aPGItest [[1]] [1] TRUE FALSE result - [magic box involving subset) result - aPGI[aPGItest[[1]]] result [[1]] [1] 864 + seth -- Seth

Re: [R] Subsetting list of vectors with list of (boolean) vectors?

2007-04-12 Thread Liaw, Andy
From: Marc Schwartz On Thu, 2007-04-12 at 18:12 +0200, Johannes Graumann wrote: Dear Rologists, I'm stuck with this. How would you do this efficiently: aPGI [[1]] [1] 864 5576 aPGItest [[1]] [1] TRUE FALSE result - [magic box involving subset) result

Re: [R] sas.get problem

2007-04-12 Thread John Kane
--- Tim Churches [EMAIL PROTECTED] wrote: John Kane wrote: I have 3 SAS files all in the directory F:/sas, two data files and a format file : form.ea1.sas7bdat form.ea2.sas7bdat sas.fmts.sas7bdat F is a USB. I am trying import them to R using sas.get. I have not

[R] Construct time series objects from raw data stored in csv files

2007-04-12 Thread tom soyer
Hi, I have time series data stored in csv files (see below for an example of the data). I understand that in order to analyze my data in R, I need to first transform it into a ts object. Howeve, I could not find an example on how exactly to do that. Is ts the only function I need? What are the

[R] Putting 2 breaks on Y axis

2007-04-12 Thread Inman, Brant A. M.D.
R plotting experts: I have a bivariate dataset composed of 300 (x,y) continuous datapoints. 297 of these points are located within the y range of [0,10], while 2 are located at 20 and one at 55. No coding errors, real outliers. When plotting these data with a scatterplot, I obviously have a

Re: [R] sas.get problem

2007-04-12 Thread John Kane
Cancel this! . I see that Hmisc is already changed. When I updated today it was changed. My thanks to Mr Dupont the maintainer. I still have a problem with sas.get but that's another post. --- John Kane [EMAIL PROTECTED] wrote: --- Tim Churches [EMAIL PROTECTED] wrote: John Kane

Re: [R] Construct time series objects from raw data stored in csv files

2007-04-12 Thread Achim Zeileis
On Thu, 12 Apr 2007, tom soyer wrote: Hi, I have time series data stored in csv files (see below for an example of the data). I understand that in order to analyze my data in R, I need to first transform it into a ts object. Howeve, I could not find an example on how exactly to do that. Is

Re: [R] Construct time series objects from raw data stored in csv files

2007-04-12 Thread Gabor Grothendieck
On 4/12/07, tom soyer [EMAIL PROTECTED] wrote: Hi, I have time series data stored in csv files (see below for an example of the data). I understand that in order to analyze my data in R, I need to first transform it into a ts object. Howeve, I could not find an example on how exactly to do

Re: [R] Construct time series objects from raw data stored in csv files

2007-04-12 Thread tom soyer
Thanks Gabor! I think your example works, but check this out: as.ts(z) Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 1921 19.0 18.4 18.3 18.1 17.7 17.6 17.7 17.7 17.5 17.5 17.4 17.3 1922 16.9 16.9 16.7 16.7 16.7 16.7 16.8 16.6 16.6 16.7 16.8 16.9 is.ts(z) [1] FALSE How

Re: [R] Construct time series objects from raw data stored in csv files

2007-04-12 Thread Gabor Grothendieck
On 4/12/07, tom soyer [EMAIL PROTECTED] wrote: Thanks Gabor! I think your example works, but check this out: as.ts(z) Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 1921 19.0 18.4 18.3 18.1 17.7 17.6 17.7 17.7 17.5 17.5 17.4 17.3 1922 16.9 16.9 16.7 16.7 16.7 16.7 16.8

Re: [R] Construct time series objects from raw data stored in csv files

2007-04-12 Thread tom soyer
Oh yes, I forgot. Thanks!! On 4/12/07, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 4/12/07, tom soyer [EMAIL PROTECTED] wrote: Thanks Gabor! I think your example works, but check this out: as.ts(z) Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 1921 19.0 18.4

Re: [R] Construct time series objects from raw data stored in csv files

2007-04-12 Thread tom soyer
Gabor, What should I do for weekly and daily data series? Are there functions similar to yearmon() for other time intervals? I see that there is a built-in yearqtr() function for quarterly data, but that's it. Thanks! On 4/12/07, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 4/12/07, tom

Re: [R] Putting 2 breaks on Y axis

2007-04-12 Thread Marc Schwartz
On Thu, 2007-04-12 at 13:41 -0500, Inman, Brant A. M.D. wrote: R plotting experts: I have a bivariate dataset composed of 300 (x,y) continuous datapoints. 297 of these points are located within the y range of [0,10], while 2 are located at 20 and one at 55. No coding errors, real outliers.

[R] making a counter of consecitive positive cases in time series

2007-04-12 Thread Jose Bustos Melo
Hi all..RCounters! I´m working with standarized time series, and i need make a counter of consecutives positves numbers to make a cumulative experimental funtion. I have x: the time series (0,1) and y: my counter, i have this for step. What is wrong?.. any can help me please!

Re: [R] Putting 2 breaks on Y axis

2007-04-12 Thread Stephen Weigand
On 4/12/07, Marc Schwartz [EMAIL PROTECTED] wrote: On Thu, 2007-04-12 at 13:41 -0500, Inman, Brant A. M.D. wrote: R plotting experts: I have a bivariate dataset composed of 300 (x,y) continuous datapoints. 297 of these points are located within the y range of [0,10], while 2 are located

[R] LME: internal workings of QR factorization --repost

2007-04-12 Thread Izmirlian, Grant \(NIH/NCI\) [E]
Hi: I've been reading Computational Methods for Multilevel Modeling by Pinheiro and Bates, the idea of embedding the technique in my own c-level code. The basic idea is to rewrite the joint density in a form to mimic a single least squares problem conditional upon the variance parameters.

Re: [R] Putting 2 breaks on Y axis

2007-04-12 Thread Greg Snow
Try something like this (modify to how you like it): x - runif(100) y - rnorm(100, 5, 2) y[1:3] - c(19, 21, 50) layout(matrix( 3:1, ncol=1 ), heights=c(2,3,4)) par(mar=c(5,4,0,2)+0.1) plot(x,y, ylim=c(0,10), ylab='') par(mar=c(0.5,4,0,2)+0.1) plot(x,y, ylim=c(18,22), xlab='', xaxt='n' )

Re: [R] Putting 2 breaks on Y axis

2007-04-12 Thread Jim Lemon
Inman, Brant A. M.D. wrote: R plotting experts: I have a bivariate dataset composed of 300 (x,y) continuous datapoints. 297 of these points are located within the y range of [0,10], while 2 are located at 20 and one at 55. No coding errors, real outliers. When plotting these data with a

[R] LME: internal workings of QR factorization

2007-04-12 Thread Izmirlian, Grant \(NIH/NCI\) [E]
Hi: I've been reading Computational Methods for Multilevel Modeling by Pinheiro and Bates, the idea of embedding the technique in my own c-level code. The basic idea is to rewrite the joint density in a form to mimic a single least squares problem conditional upon the variance parameters. The

[R] Random Sequence

2007-04-12 Thread Anup Nandialath
Dear Friends, I'm trying to generate a sequence of 100 observations with either a 1 or -1. In other words the sequence should look something like this. y = 1 1 -1 1 -1 -1 -1 1 1 .. Can somebody please give me some direction on how I can do this in R. Thanks Anup

Re: [R] Random Sequence

2007-04-12 Thread Matthew Keller
Hi Anup, (runif(100).5)*1 #would give you 0's and 1's. sample(rep(c(-1,1),50),100) #A bit slower I think, gives you -1's and 1's Best, Matt On 4/12/07, Anup Nandialath [EMAIL PROTECTED] wrote: Dear Friends, I'm trying to generate a sequence of 100 observations with either a 1 or -1.

Re: [R] Random Sequence

2007-04-12 Thread Peter Dalgaard
Anup Nandialath wrote: Dear Friends, I'm trying to generate a sequence of 100 observations with either a 1 or -1. In other words the sequence should look something like this. y = 1 1 -1 1 -1 -1 -1 1 1 .. Can somebody please give me some direction on how I can do this in R.

Re: [R] Random Sequence

2007-04-12 Thread Francisco J. Zagmutt
There are many ways to do this. The first that comes to my mind is sample(c(1,-1),100,TRUE). Notice that sample also has a prob argument that may be useful for you. Francisco Anup Nandialath wrote: Dear Friends, I'm trying to generate a sequence of 100 observations with either a 1 or

Re: [R] Putting 2 breaks on Y axis

2007-04-12 Thread Greg Snow
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Stephen Weigand Sent: Thursday, April 12, 2007 2:32 PM To: Inman, Brant A. M.D. Cc: [EMAIL PROTECTED] Subject: Re: [R] Putting 2 breaks on Y axis [snip] I was thinking plot the data without the

Re: [R] Random Sequence

2007-04-12 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Matthew Keller Sent: Thursday, April 12, 2007 2:32 PM To: Anup Nandialath Cc: [EMAIL PROTECTED] Subject: Re: [R] Random Sequence Hi Anup, (runif(100).5)*1 #would give you 0's and 1's.

[R] GLM with random effects

2007-04-12 Thread Brendan Connors
Hi R-Users, I have 3 replicates ('Replicate) of counts of parasites ('nor.tot.lep') before and after an experiment ('In.Out'). I am trying to treat the three replicates as a random effect in order to determine if the main effect (In.Out) significantly influences my dependent variable

[R] problems in loading MASS

2007-04-12 Thread Weiwei Shi
Hi, there: After I upgraded my R to 2.4.1, it is my first time of trying to use MASS and found the following error message: install.packages(MASS) --- Please select a CRAN mirror for use in this session --- trying URL 'http://cran.cnr.Berkeley.edu/bin/macosx/universal/contrib/2.4/VR_7.2-33.tgz'

[R] Question on ridge regression with R

2007-04-12 Thread Benny Wong
Hi, I am working on a project about hospital efficiency. Due to the high multicolinearlity of the data, I want to fit the model using ridge regression. However, I believe that the data from large hospital(indicated by the number of patients they treat a year) is more accurate than from small

Re: [R] Reasons to Use R

2007-04-12 Thread Ted Harding
On 12-Apr-07 10:14:21, Jim Lemon wrote: Charilaos Skiadas wrote: A new fortune candidate perhaps? On Apr 10, 2007, at 6:27 PM, Greg Snow wrote: Remember, everything is better than everything else given the right comparison. Only if we remove the grammatical blip that turns it into an

Re: [R] GLM with random effects

2007-04-12 Thread Andrew Robinson
Hello Brendan, do you have any particular objection to lmer(nor.tot.lep ~ In.Out + (1|Replicate), data=coho, family=poisson) ? Andrew On Thu, Apr 12, 2007 at 02:58:07PM -0700, Brendan Connors wrote: Hi R-Users, I have 3 replicates ('Replicate) of counts of parasites ('nor.tot.lep')

Re: [R] LME: internal workings of QR factorization

2007-04-12 Thread Douglas Bates
On 4/12/07, Izmirlian, Grant (NIH/NCI) [E] [EMAIL PROTECTED] wrote: Hi: I've been reading Computational Methods for Multilevel Modeling by Pinheiro and Bates, the idea of embedding the technique in my own c-level code. The basic idea is to rewrite the joint density in a form to mimic a

[R] class and method

2007-04-12 Thread jiho.han
Hi, R-experts I defined a new class of object, called alpha, which is basically a data.frame. And (I think) I know how to create a method for alpha, such as summary.alpha, plot.alpha. The problem is, when I try to access alpha object by usual data.frame method, it won't. For example, suppose X

Re: [R] problems in loading MASS

2007-04-12 Thread Charilaos Skiadas
Weiwei, I have never had much success installing packages from within R.app on MacOSX, because the location that it is supposed to save things, / Library/Frameworks/, needs elevated priviledges, which the app doesn't seem to try to get. So it at best ends up saving it in some

Re: [R] problems in loading MASS

2007-04-12 Thread Weiwei Shi
Hi Charilaos: I think the installation has no problem as indicated by install.packages(). Anyway, I tried your way but unfortunately it does not help. Thanks though, -W On 4/12/07, Charilaos Skiadas [EMAIL PROTECTED] wrote: Weiwei, I have never had much success installing packages

Re: [R] Construct time series objects from raw data stored in csv files

2007-04-12 Thread Gabor Grothendieck
On 4/12/07, tom soyer [EMAIL PROTECTED] wrote: What should I do for weekly and daily data series? Are there functions similar to yearmon() for other time intervals? I see that there is a built-in yearqtr() function for quarterly data, but that's it. ts series cannot directly represent daily

[R] NYC Area R / S / S-Plus Users Group Forming, Thursday, April 26th

2007-04-12 Thread Talbot Katz
Hi gang. So I'm actually carrying through on the threat I made a couple of months back, and helping to start a local users group for the NYC Metropolitan area. Here's a copy of the announcement of our introductory / planning meeting, I also can send a flyer version of this to anyone who

[R] trouble getting bw plots using trellis graphics

2007-04-12 Thread Brian Riordan
I'm having trouble getting a black and white plot using the trellis device. I have color graphics, but I'd like to use bw for a paper. For example, trellis.device(color=F) d - read.table(textConnection( A,B 0,1 1,2), header=T, sep=,) xyplot(A ~ B, data=d, type=o, lty=c(1:2), col=2, pch=19:20)

Re: [R] trouble getting bw plots using trellis graphic

2007-04-12 Thread Deepayan Sarkar
On 4/12/07, Brian Riordan [EMAIL PROTECTED] wrote: I'm having trouble getting a black and white plot using the trellis device. I have color graphics, but I'd like to use bw for a paper. For example, trellis.device(color=F) d - read.table(textConnection( A,B 0,1 1,2), header=T, sep=,)

Re: [R] class and method

2007-04-12 Thread Prof Brian Ripley
On Thu, 12 Apr 2007, jiho.han wrote: Hi, R-experts I defined a new class of object, called alpha, which is basically a data.frame. And (I think) I know how to create a method for alpha, such as summary.alpha, plot.alpha. The problem is, when I try to access alpha object by usual data.frame

Re: [R] making a counter of consecitive positive cases in time series

2007-04-12 Thread Marc Schwartz
On Thu, 2007-04-12 at 22:30 +0200, Jose Bustos Melo wrote: Hi all..RCounters! Im working with standarized time series, and i need make a counter of consecutives positves numbers to make a cumulative experimental funtion. I have x: the time series (0,1) and y: my counter, i have this

Re: [R] problems in loading MASS

2007-04-12 Thread Prof Brian Ripley
MASS is part of the standard R installation, as a recommended package. I suggest you try re-installing R and not then trying to download a different version. I would also suggest trying 2.5.0 beta at this stage in preference to 2.4.1 if you do need to re-install R. Beyond that, there is a