Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Jon Olav Skoien
On 12/17/2010 6:22 PM, Duncan Murdoch wrote: On 17/12/2010 11:13 AM, Jon Olav Skoien wrote: Dear list, (R 2.12.0, Windows 7, 64bit) I recently tried to install a new package (spacetime), that depends on sp among others. I already had the last one installed, but there was probably a newer

Re: [R] Alternative to extended recode sintax? Bug?

2010-12-20 Thread Luca Meyer
All right, I get it now: lubridate's week() define weeks from Thursday till the following Wednesday. You'd probably agree with me that it's a bit strange what it is going to do over the turn of the year: y -

[R] R is not running well

2010-12-20 Thread Roslina Zakaria
Hi, I'm working on windows 7.  Recently I install the latest R and also use together with Tinn-R but it is not working well.  I got the following message.  I have tried saving the Rprofile as suggested by the R forum, but still it is not running.  What should I do? source(.trPaths[5],

[R] incremental learning for LOESS time series model

2010-12-20 Thread ying zhang
Hi All, I am currently working on some time series data, I know I can use LOESS/ARIMA model. The data is written to a vector whose length is 1000, which is a queue, updating every 15 minutes, Thus the old data will pop out while the new data push in the vector. I can rerun the

Re: [R] Alternative to extended recode sintax? Bug?

2010-12-20 Thread Uwe Ligges
On 20.12.2010 06:54, Luca Meyer wrote: All right, I get it now: lubridate's week() define weeks from Thursday till the following Wednesday. You'd probably agree with me that it's a bit strange what it is going to do over the turn of the year: y-

Re: [R] Layout of mulitpage conditioned lattice plots

2010-12-20 Thread Dieter Menne
David Winsemius wrote: Here's my latest guess at what you may want: pdf(file=multpage.pdf) xyplot(val~time|subj + comp, data=dt,type=l, layout=c(3,5, 3), skip=rep(c(rep(FALSE,13), TRUE, TRUE), 3) ) dev.off() Not really, but skip was the right idea. I added

Re: [R] barplot: width of label

2010-12-20 Thread Jim Lemon
On 12/20/2010 02:13 AM, fransiepansiekevertje wrote: Hello, I try to make barplots with rather wide labels. A simplified example of this: x- c(12, 33, 56, 67, 15, 66) names(x)- c('Richard with a long surname','Minnie with a long name,'Albert','Helen','Joe','Kingston') barplot(x, las = 2) Now

[R] contourplot help

2010-12-20 Thread Andre Nathan
Hello I'm using the following call to create a contourplot: library(lattice) m - as.matrix(read.table(data.txt)) contourplot(m[,3] ~ m[,2] * -m[,1], at = c(1e-6, 1e-5, 1e-4, 1e-3, 1e-2, 1e-1), scales = list(x = list(log = 10, labels =

[R] For-loop

2010-12-20 Thread Anne-Christine Mupepele
Hi, I have the following problem: I have a data.frame with 36 sample sites (colums) for which I have covariates in 3 categories: Area, Month and River. Each Area consists of 3 rivers, which were sampled over 3 month. Now I want to fuse River 1-3 for one area in one month. To get a data.frame

[R] transposing panel data

2010-12-20 Thread taby gathoni
I am currently trying to transpose some large panel data set ie transposing multiple rows in into a single column. instead the transpose functionality transposes all rows into columns. my sample data set looks like below: ACCT_NUM ACCOUNT_NAME TRAN_AMT DATE EMPLOYER

[R] Modifying effect display (effects package)

2010-12-20 Thread Anna Berthinussen
Hi, I am using the effects package and plot() to produce an effect plot for a higher order term from a GLM. However, I need to modify the graph it produces by removing the tick marks and axes from the right and top of the graph, and moving the remaining tick marks inside the axes. Does

[R] Metafor package

2010-12-20 Thread petretta
I have some question about metafor package. I'm interest to perform a random effect meta-analysis of proportion (single group summary of prevalence of disease in a population as reported by different study) It ask: 1. PFT: The Freeman-Tukey double arcsine transformed proportion is

Re: [R] R is not running well

2010-12-20 Thread Uwe Ligges
On 20.12.2010 10:29, Roslina Zakaria wrote: Hi, I'm working on windows 7. Recently I install the latest R and also use together with Tinn-R but it is not working well. I got the following message. I have tried saving the Rprofile as suggested by the R forum, but still it is not running.

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Uwe Ligges
On 20.12.2010 09:41, Jon Olav Skoien wrote: On 12/17/2010 6:22 PM, Duncan Murdoch wrote: On 17/12/2010 11:13 AM, Jon Olav Skoien wrote: Dear list, (R 2.12.0, Windows 7, 64bit) I recently tried to install a new package (spacetime), that depends on sp among others. I already had the last one

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Duncan Murdoch
Jon Olav Skoien wrote: On 12/17/2010 6:22 PM, Duncan Murdoch wrote: On 17/12/2010 11:13 AM, Jon Olav Skoien wrote: Dear list, (R 2.12.0, Windows 7, 64bit) I recently tried to install a new package (spacetime), that depends on sp among others. I already had the last one installed, but there

[R] How to read japanese characters from a database/file

2010-12-20 Thread vikrant
Hi, I am using rodbc package to connect with Sql Server 2005. I have a table which contains some japanese characters. Now I want to use these japanese characters as labels in simple x y plot? Here is the data in the table サービスデスク 35355 データネットワーク 8417 I want to read this data from

Re: [R] Alternative to extended recode sintax? Bug?

2010-12-20 Thread David Winsemius
On Dec 20, 2010, at 12:54 AM, Luca Meyer wrote: All right, I get it now: lubridate's week() define weeks from Thursday till the following Wednesday. You'd probably agree with me that it's a bit strange what it is going to do over the turn of the year: y -

[R] Odp: For-loop

2010-12-20 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 20.12.2010 11:48:51: Hi, I have the following problem: I have a data.frame with 36 sample sites (colums) for which I have covariates in 3 categories: Area, Month and River. Each Area consists of 3 rivers, which were sampled over 3 month. Now

Re: [R] For-loop

2010-12-20 Thread jim holtman
try this: # create indexing vector to select 3 adjacent columns indx - sapply(seq(1, 36, 3), seq, length = 3) # process each row of the table ans - t(apply(data, 1, function(.row){ + # use indx to sum up the columns + apply(indx, 2, function(.indx){ + sum(.row[.indx]) + })

Re: [R] transposing panel data

2010-12-20 Thread Ben Bolker
taby gathoni tabieg at yahoo.com writes: I am currently trying to transpose some large panel data set ie transposing multiple rows in into a single column. instead the transpose functionality transposes all rows into columns. my sample data set looks like below: Try the melt and cast

Re: [R] Turning a Variable into String

2010-12-20 Thread Paolo Rossi
Thank you very much to the both of you Paolo On 20 December 2010 00:35, Duncan Murdoch murdoch.dun...@gmail.com wrote: On 19/12/2010 7:21 PM, Paolo Rossi wrote: I would like to know how to turn a variable into a string. I have tried as.symbol and as.name but it doesnt work for what I'd

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Jon Olav Skoien
On 12/20/2010 1:43 PM, Duncan Murdoch wrote: Jon Olav Skoien wrote: On 12/17/2010 6:22 PM, Duncan Murdoch wrote: On 17/12/2010 11:13 AM, Jon Olav Skoien wrote: Dear list, (R 2.12.0, Windows 7, 64bit) I recently tried to install a new package (spacetime), that depends on sp among others. I

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Jon Olav Skoien
On 12/20/2010 1:30 PM, Uwe Ligges wrote: On 20.12.2010 09:41, Jon Olav Skoien wrote: On 12/17/2010 6:22 PM, Duncan Murdoch wrote: On 17/12/2010 11:13 AM, Jon Olav Skoien wrote: Dear list, (R 2.12.0, Windows 7, 64bit) I recently tried to install a new package (spacetime), that depends on

Re: [R] Defining Two-dimensional Separable Variance-Covariance Structure in nlme Package

2010-12-20 Thread Muhammad Yaseen
*Hello All,* * * *I'd highly appreciate if someone can explain the way of defining Two-dimensional Separable Variance-Covariance Structure in nlme Package. Thanks* * * * * -- * Muhammad Yaseen * [[alternative HTML version deleted]] __

Re: [R] box-and-whisker plots based on summary not data

2010-12-20 Thread Matthew Vernon
Matthew Vernon m.c.ver...@warwick.ac.uk writes: Is it possible to produce box-and-whisker plots given that I have the median, interquartile and 5/95th centile values, but not the data from which they come? The answer, which came from Steve Ellison (thanks!) is to note that stats is the only

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Duncan Murdoch
On 20/12/2010 9:03 AM, Jon Olav Skoien wrote: On 12/20/2010 1:43 PM, Duncan Murdoch wrote: Jon Olav Skoien wrote: On 12/17/2010 6:22 PM, Duncan Murdoch wrote: On 17/12/2010 11:13 AM, Jon Olav Skoien wrote: Dear list, (R 2.12.0, Windows 7, 64bit) I recently tried to install a new

Re: [R] box-and-whisker plots based on summary not data

2010-12-20 Thread David Winsemius
On Dec 20, 2010, at 8:33 AM, Matthew Vernon wrote: Matthew Vernon m.c.ver...@warwick.ac.uk writes: Is it possible to produce box-and-whisker plots given that I have the median, interquartile and 5/95th centile values, but not the data from which they come? Please note that the whiskers

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Uwe Ligges
On 20.12.2010 15:19, Duncan Murdoch wrote: On 20/12/2010 9:03 AM, Jon Olav Skoien wrote: On 12/20/2010 1:43 PM, Duncan Murdoch wrote: Jon Olav Skoien wrote: On 12/17/2010 6:22 PM, Duncan Murdoch wrote: On 17/12/2010 11:13 AM, Jon Olav Skoien wrote: Dear list, (R 2.12.0, Windows 7,

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Duncan Murdoch
On 20/12/2010 9:26 AM, Uwe Ligges wrote: On 20.12.2010 15:19, Duncan Murdoch wrote: On 20/12/2010 9:03 AM, Jon Olav Skoien wrote: On 12/20/2010 1:43 PM, Duncan Murdoch wrote: Jon Olav Skoien wrote: On 12/17/2010 6:22 PM, Duncan Murdoch wrote: On 17/12/2010 11:13 AM, Jon Olav

Re: [R] moving average with gaps in time series

2010-12-20 Thread Josef . Kardos
Thank you for the suggestion Jannis. I ended up using the function 'running' in package gregmisc. 'rollyapply' was not able to process all the data on my computer. Happy holidays to you and the R community Jannis bt_jan...@yahoo.de 12/16/2010 06:56 PM To josef.kar...@phila.gov cc

Re: [R] lower/upper case question

2010-12-20 Thread Bos, Roger
toupper() and tolower() -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Erin Hodgess Sent: Friday, December 17, 2010 4:10 PM To: R help Subject: [R] lower/upper case question Dear R People: Is there a function to convert a

Re: [R] dotchart for matrix data

2010-12-20 Thread e-letter
On 18/12/2010, e-letter inp...@gmail.com wrote: On 18/12/2010, Peter Ehlers ehl...@ucalgary.ca wrote: On 2010-12-18 07:50, e-letter wrote: Ben Bolker Sat, 18 Dec 2010 07:07:24 -0800 [... snip ...] I am trying to create a chart like this

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Duncan Murdoch
On 20/12/2010 9:29 AM, Duncan Murdoch wrote: On 20/12/2010 9:26 AM, Uwe Ligges wrote: On 20.12.2010 15:19, Duncan Murdoch wrote: On 20/12/2010 9:03 AM, Jon Olav Skoien wrote: [ lots deleted ] Yes, I had it open. In this case it was intentional to give a reproducible example

Re: [R] contourplot help

2010-12-20 Thread David Winsemius
On Dec 20, 2010, at 5:36 AM, Andre Nathan wrote: Hello I'm using the following call to create a contourplot: library(lattice) m - as.matrix(read.table(data.txt)) contourplot(m[,3] ~ m[,2] * -m[,1], at = c(1e-6, 1e-5, 1e-4, 1e-3, 1e-2, 1e-1), scales = list(x = list(log

Re: [R] R.matlab memory use

2010-12-20 Thread Stefano Ghirlanda
Hi Ben, Thanks for your reply. My data structure is about 2 x 2000 so one order of magnitude the one you tried. I have no problem saving and reading smaller data structures (even large ones, just not his large) between octave and R using octave's save -7 (which saves MATLAB v5 files) and

Re: [R] help on timeseries

2010-12-20 Thread cameron
Thanks Frederic. that solved the problem. Thanks again. Cameron -- View this message in context: http://r.789695.n4.nabble.com/help-on-timeseries-tp3076866p3095633.html Sent from the R help mailing list archive at Nabble.com. __

[R] After heteroskedasticity correction, how can I get new confidential interval?

2010-12-20 Thread JoonGi
I just corrected std.error of my 'model'(Multi Regression). Then how can I get new t and p-values? Isn't there any R command which shows new t and p values? -- View this message in context:

[R] How to optimize function parameters?

2010-12-20 Thread zbynek.jano...@gmail.com
Hi, I have a dataset and I want to fit a function to it. The function is variogram model (http://en.wikipedia.org/wiki/Variogram) The variogram model is defined by three parameters and I want them to be automatically optimized for real time data. I tried to use gafit {gafit} for this, but there

Re: [R] How to optimize function parameters?

2010-12-20 Thread Jonathan P Daily
Without an example, I'm not sure how much help I can be. Does your Variogram calculation return a single value that could be minimized? Or could it be coerced to be so? I would suggest, if you haven't done so, a look at ?nlm and ?nlminb. Alternatively, the CRAN package 'optimx' serves as a

Re: [R] Alternative to extended recode sintax? Bug?

2010-12-20 Thread Luca Meyer
Right, I appreciate the first day of the year start date. I am just wondering why then the cut off day is not the same for the rest of the year...but it's all right to use other packages. Thanks, Luca Il giorno 20/dic/2010, alle ore 14.16, David Winsemius ha scritto: On Dec 20, 2010, at

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Jon Olav Skoien
This sounds like a good solution for the case I described in my first email. Thanks a lot! Jon On 12/20/2010 4:05 PM, Duncan Murdoch wrote: On 20/12/2010 9:29 AM, Duncan Murdoch wrote: On 20/12/2010 9:26 AM, Uwe Ligges wrote: On 20.12.2010 15:19, Duncan Murdoch wrote: On 20/12/2010

Re: [R] Alternative to extended recode sintax? Bug?

2010-12-20 Thread David Winsemius
On Dec 20, 2010, at 10:58 AM, Luca Meyer wrote: Right, I appreciate the first day of the year start date. I am just wondering why then the cut off day is not the same for the rest of the year...but it's all right to use other packages. Are you saying it shifts within the year? I am not

Re: [R] After heteroskedasticity correction, how can I get new confidential interval?

2010-12-20 Thread Andrew Miles
It is hard to say without knowing more about the type of model you are running. A good place to look would be at the coeftest() function in the package lmtest. Andrew Miles On Dec 20, 2010, at 10:04 AM, JoonGi wrote: I just corrected std.error of my 'model'(Multi Regression). Then

[R] ideas, modeling highly discrete time-series data

2010-12-20 Thread Mike Williamson
Hello all, First of all, thanks so those of you who helped me a week or so ago managing a time series with varying gaps between the data series in 'R'. (My final preferred solution was to use its function then forecast(Arima( ) ). ) My next question is a general statistical question

[R] Passing parameter to a function

2010-12-20 Thread Luca Meyer
I am trying to pass a couple of variable names to a xtabs formula: tab - function(x,y){ xtabs(time~x+y, data=D) } But when I run: tab(A,B) I get: Error in eval(expr, envir, enclos) : object A not found I am quite sure that there is some easy way out, but I have tried with different

Re: [R] Alternative to extended recode sintax? Bug?

2010-12-20 Thread Luca Meyer
Yes, I am seeing that at the end of 2010-beginning 2011. Try: weekdays(as.POSIXct(2010-12-25)+(0:20)*24*60*60) week(as.POSIXct(2010-12-25)+(0:20)*24*60*60) Week 1 (2011) is made up of 6 days Luca Il giorno 20/dic/2010, alle ore 17.54, David Winsemius ha scritto: On Dec 20, 2010, at

Re: [R] Complicated nls formula giving singular gradient message

2010-12-20 Thread Jared Blashka
Though my topic is slightly old already, I feel that it is necessary to post an update on my situation. I ended being able to estimate the parameters for this problem without having to worry as much about initial parameter estimates using AD Model Builder. It calculates the exact gradient using

Re: [R] Passing parameter to a function

2010-12-20 Thread Duncan Murdoch
On 20/12/2010 1:13 PM, Luca Meyer wrote: I am trying to pass a couple of variable names to a xtabs formula: tab- function(x,y){ xtabs(time~x+y, data=D) } But when I run: tab(A,B) I get: Error in eval(expr, envir, enclos) : object A not found I am quite sure that there is some easy

Re: [R] package survey

2010-12-20 Thread Martin Maechler
JS == Joel Schwartz j...@joelschwartz.com on Sat, 18 Dec 2010 19:09:24 -0800 writes: -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Saturday, December 18, 2010 5:54 PM To: Joel Schwartz Cc: r-help@r-project.org Subject: Re: [R]

Re: [R] sum with times

2010-12-20 Thread Martin Maechler
RRJ == Ronaldo Reis Junior chrys...@gmail.com on Sun, 19 Dec 2010 12:10:54 -0200 writes: RRJ Hi, Forget, the chron package work with this No, don't forget, rather solve the problem: Date / dateTime classes are native R entities; extension package 'chron' objects are not.

Re: [R] R.matlab memory use

2010-12-20 Thread Ben Bolker
Stefano Ghirlanda dr.ghirlanda at gmail.com writes: Hi Ben, Thanks for your reply. My data structure is about 2 x 2000 so one order of magnitude the one you tried. I have no problem saving and reading smaller data structures (even large ones, just not his large) between octave and R

Re: [R] Passing parameter to a function

2010-12-20 Thread Luca Meyer
Hi Duncan, Yes, A and B are columns in D. Having said that I and trying to avoid tab(D$A,D$B) and I would prefer: tab(A,B) Unfortunately the syntax you suggest is giving me the same error: Error in eval(expr, envir, enclos) : object A not found I have tried to add some deparse() but I have

[R] survexp - unable to reproduce example

2010-12-20 Thread Heinz Tuechler
Dear All, when I try to reproduce an example of survexp, taken from the help page of survdiff, I receive the error message Error in floor(temp) : Non-numeric argument to mathematical function . It seems to come from match.ratetable. I think, it has to do with character variables in a

[R] sample() issue

2010-12-20 Thread cory n
length(sample(25000, 25000*(1-.55))) [1] 11249 25000*(1-.55) [1] 11250 length(sample(25000, 11250)) [1] 11250 length(sample(25000, 25000*.45)) [1] 11250 So the question is, why do I get 11249 out of the first command and not 11250? I can't figure this one out. Thanks Cory

Re: [R] After heteroskedasticity correction, how can I get new confidential interval?

2010-12-20 Thread JoonGi
First of all, thanks for your guide! Let me be more specific, here. Using data(Housing) from Ecdat library I ran a regression raw.model-lprice~llot+lbed+lbath+lsto+factor(driveway)+factor(recroom)+factor(fullbase)+factor(gashw)+factor(airco)+factor(prefarea)+factor(garagepl)

Re: [R] How to optimize function parameters?

2010-12-20 Thread zbynek.jano...@gmail.com
Thanks for your suggestions, I will certainly look at that To answer your question... I can calculate root square error between empirical data a those predicted by model (I used this to optimize the parameters using gafit()). -- View this message in context:

[R] performance issue with merge

2010-12-20 Thread Maxim
Hi, I'm trying to merge a couple of matrices together. The matrices are created from a couple of input files that have very similar (but not in all cases) values in first and second column. 1st and 2nd column are used to generate rownames, the actually interesting value is in column 3. I'd like

Re: [R] sample() issue

2010-12-20 Thread David Winsemius
On Dec 20, 2010, at 2:04 PM, cory n wrote: length(sample(25000, 25000*(1-.55))) [1] 11249 25000*(1-.55) [1] 11250 length(sample(25000, 11250)) [1] 11250 length(sample(25000, 25000*.45)) [1] 11250 So the question is, why do I get 11249 out of the first command and not 11250? I

Re: [R] Passing parameter to a function

2010-12-20 Thread Bert Gunter
Duncan's solution works for me. I strongly suspect that you have not read the Help file for xtabs carefully enough. Does the time column in D consist of a vector of counts? A - factor(c(A,A,B,B)) B - factor(c(1,2,1,2)) D - data.frame(A,B) D$time - c(10,12,14,16) tab(A,B) B A1 2

Re: [R] After heteroskedasticity correction, how can I get new confidential interval?

2010-12-20 Thread Achim Zeileis
On Mon, 20 Dec 2010, JoonGi wrote: First of all, thanks for your guide! Let me be more specific, here. Using data(Housing) from Ecdat library I ran a regression

Re: [R] sample() issue

2010-12-20 Thread Viechtbauer Wolfgang (STAT)
See FAQ 7.31: http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f And try this: 25000*(1-.55) - 11250 25000*.45 - 11250 Notice a difference? Best, -- Wolfgang Viechtbauer Department of Psychiatry and Neuropsychology School for Mental Health and

Re: [R] sample() issue

2010-12-20 Thread Peter Langfelder
On Mon, Dec 20, 2010 at 11:04 AM, cory n corynis...@gmail.com wrote: length(sample(25000, 25000*(1-.55))) [1] 11249 25000*(1-.55) [1] 11250 length(sample(25000, 11250)) [1] 11250 length(sample(25000, 25000*.45)) [1] 11250 So the question is, why do I get 11249 out of the first command

Re: [R] sample() issue

2010-12-20 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On Behalf Of cory n Sent: Monday, December 20, 2010 11:04 AM To: r-help@r-project.org Subject: [R] sample() issue length(sample(25000, 25000*(1-.55))) [1] 11249 25000*(1-.55) [1]

Re: [R] sample() issue

2010-12-20 Thread Jorge Ivan Velez
Hi Cory, Check out http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f HTH, Jorge On Mon, Dec 20, 2010 at 2:04 PM, cory n wrote: length(sample(25000, 25000*(1-.55))) [1] 11249 25000*(1-.55) [1] 11250 length(sample(25000, 11250)) [1]

Re: [R] How to optimize function parameters?

2010-12-20 Thread Jonathan P Daily
If you can somehow add the sign of the output to the root square error, nlm, nlminb, and many options within optimx accept bounds. -- Jonathan P. Daily Technician - USGS Leetown Science Center 11649 Leetown Road Kearneysville WV, 25430 (304) 724-4480 Is the

Re: [R] Complicated nls formula giving singular gradient message

2010-12-20 Thread Bert Gunter
Jared: You realize, of course, that just because you get estimates of the parameters from the software is no guarantee that the estimates mean anything? Nor does it mean that they mean nothing, I hasten to add. If, as one might suspect, the model is overparameterized, the estimates may be so

Re: [R] Passing parameter to a function

2010-12-20 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Luca Meyer Sent: Monday, December 20, 2010 11:03 AM To: Duncan Murdoch Cc: R-help@r-project.org Subject: Re: [R] Passing parameter to a function Hi Duncan, Yes, A and B

Re: [R] sample() issue

2010-12-20 Thread Viechtbauer Wolfgang (STAT)
Trying to make the fastest referal to FAQ 7.31 is becoming quite a sport around here. David beat us all to it on this round. Chances are there will be many more rounds to come. Best, Wolfgang Original Message From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On

[R] Sine Regression in R

2010-12-20 Thread Paolo Rossi
Hi everyone, I am trying to fit a sine function on one year of wind data. I have two questions below. Looking around on the net I managed to get the following: Sine Equation: y = a + b * sin( c + d*x ) b is the amplitude, c is the phase shift, d is something deal with periodicty of data*.*

Re: [R] Sine Regression in R

2010-12-20 Thread Dennis Murphy
Hi: Try using package sos to address your question re which packages perform sine/trig regression: library(sos) findFn('sine regression')# 12 matches on my system findFn('trigonometric regression') # 9 matches The results point largely to the mFilter, oce, circStats and circular packages.

Re: [R] breaking an inner loop in R

2010-12-20 Thread Berend Hasselman
kayj wrote: Hi All, is there a way to break an inner loop in R? for example for ( i in 1: 100){ if ( condition){ break } } I have search and I have seen exampels where they add a semicolon after break, do I need to do that? Also when I type ?break in R in returns nothing

Re: [R] Metafor package

2010-12-20 Thread Viechtbauer Wolfgang (STAT)
Dear Mario, 1) Whether to use 1/2* or not is arbitary. It's just a multiplicative constant. 2) vi - 1/(4*ni + 2) If one leave off the 1/2* part, then it would be vi - 1/(ni + 1/2). 3) Yes, but it takes a bit more work. Fortunately, the same issue was raised not very long ago, so the following

Re: [R] sample() issue

2010-12-20 Thread Petr Savicky
On Mon, Dec 20, 2010 at 01:04:18PM -0600, cory n wrote: length(sample(25000, 25000*(1-.55))) [1] 11249 25000*(1-.55) [1] 11250 Besides FAQ 7.31, look also at http://rwiki.sciviews.org/doku.php?id=misc:r_accuracy for further examples and hints concerning rounding errors in simple

[R] RExcel doesn't get active dataset

2010-12-20 Thread jryan . daniels
Hey everyone When I try to 'get active dataframe' from the context menu in excel using Rcmdr menus in excel, only the header and first row of data gets copied to the excel spread sheet. No clue why this is happening. I followed the instructions from 'R through Excel' pg 25-26. I'm using excel

[R] finding and opening C source called from R functions.

2010-12-20 Thread Chiquoine, Ben
Hi, As will soon be very clear I'm an R novice. I'm trying to better understand the ks.test function in the stats package. When I look at the source code there are several calls to C functions (for example .C(pkstwo, is.integer(length(x[IND])), p=as.double(x[IND]), as.double(tol), PACKAGE =

Re: [R] finding and opening C source called from R functions.

2010-12-20 Thread Peter Langfelder
On Mon, Dec 20, 2010 at 12:48 PM, Chiquoine, Ben bchiquo...@tiff.org wrote: Hi, As will soon be very clear I'm an R novice.  I'm trying to better understand the ks.test function in the stats package.  When I look at the source code there are several calls to C functions (for example

Re: [R] Sine Regression in R

2010-12-20 Thread David Winsemius
On Dec 20, 2010, at 3:01 PM, Paolo Rossi wrote: Hi everyone, I am trying to fit a sine function on one year of wind data. I have two questions below. Looking around on the net I managed to get the following: Sine Equation: y = a + b * sin( c + d*x ) b is the amplitude, c is the phase

Re: [R] finding and opening C source called from R functions.

2010-12-20 Thread Duncan Murdoch
On 20/12/2010 3:56 PM, Peter Langfelder wrote: On Mon, Dec 20, 2010 at 12:48 PM, Chiquoine, Benbchiquo...@tiff.org wrote: Hi, As will soon be very clear I'm an R novice. I'm trying to better understand the ks.test function in the stats package. When I look at the source code there are

Re: [R] R Code and the Pygments Python SyntaxHighlighter

2010-12-20 Thread Liviu Andronic
(reviving from the grave) On Thu, Mar 27, 2008 at 5:44 PM, Hans W. Borchers hwborch...@gmail.com wrote: is someone going to write a R/S language lexer for the Pygments Python syntax highlighter http://pygments.org/? As it is used now by Trac, Django, or the Python documentation tool Sphinx,

[R] .Rd file for S4-method warning

2010-12-20 Thread Mark Heckmann
Dear R users, I want to create a proper .Rd file for the show method for an S4 class. I am encountering problems in the \usage{} line, I guess. An example: setClass(testClass, representation(a=character)) setMethod(show, testClass, function(object){ }) The .Rd file:

Re: [R] .Rd file for S4-method warning

2010-12-20 Thread Duncan Murdoch
On 20/12/2010 5:18 PM, Mark Heckmann wrote: Dear R users, I want to create a proper .Rd file for the show method for an S4 class. I am encountering problems in the \usage{} line, I guess. An example: setClass(testClass, representation(a=character)) setMethod(show, testClass,

Re: [R] R.matlab memory use

2010-12-20 Thread Stefano Ghirlanda
Hi, The data set is indeed highly redundant, mostly zeros (which means I could use sparse matrices, but that's another issue). I have posted an octave ascii file here: http://dl.dropbox.com/u/4000198/freq.txt and a .mat file here: http://dl.dropbox.com/u/4000198/freq.mat    A few more

Re: [R] Plot a matrix recursively

2010-12-20 Thread AlexZolot
for data.frame: for(j in grep('Laser_', names(m)) lines(m[,j]) for matrix: for(j in grep('Laser_', colnames(m)) lines(m[,j]) or for(j in 2:4) lines(m[,j]) Shorter could be worse if you insert additional columns later. -- View this message in context:

Re: [R] plot more plots from one matrix

2010-12-20 Thread AlexZolot
See library(reshape) -- View this message in context: http://r.789695.n4.nabble.com/plot-more-plots-from-one-matrix-tp3069545p3107761.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] source file for rnorm

2010-12-20 Thread Yue Zhang
Dear Sir/Madam: I am sorry if this is a simple question. I have downloaded the source code for R but could not find the source file that has the rnorm function definition. Could you please let me know which file has the function? Also, in general, how do I find the source file of a given

Re: [R] Plot a matrix recursively

2010-12-20 Thread AlexZolot
for data.frame: for(j in grep('Laser_', names(m)) lines(m[,j]) for matrix: for(j in grep('Laser_', colnames(m)) lines(m[,j]) or: for(j in 2:4) lines(m[,j]) Shorter could be worse if you insert additional columns later. alcesgabbo wrote: Hi, I have the following matrix (named

Re: [R] Filter data

2010-12-20 Thread AlexZolot
You missed the first, the most important line of str() output. What type is your data? -- View this message in context: http://r.789695.n4.nabble.com/Filter-data-tp3070069p3108011.html Sent from the R help mailing list archive at Nabble.com. __

[R] logistic regression or not?

2010-12-20 Thread array chip
Hi, I have a dataset where the response for each person on one of the 2 treatments was a proportion (percentage of certain number of markers being positive), I also have the number of positive negative markers available for each person. what is the best way to analyze this kind of data? I can

Re: [R] package arules - 'transpose' of the transactions

2010-12-20 Thread Kohleth Chia
Just found 'ngCMatrix' class, it must be the way to go? - Message: 47Date: Mon, 20 Dec 2010 17:41:20 +1100 From: Kohleth Chia kohl...@gmail.com To: r-help@r-project.org Subject: [R] package arules - 'transpose' of the transactions Message-ID:

Re: [R] install.packages() - old version deleted, new version did not install

2010-12-20 Thread Duncan Murdoch
I've just tentatively put code to fix this in place in R-devel. By default it is not enabled, because if it goes wrong it could really mess things up. To enable it, set options(install.lock=TRUE) before installing or updating binary packages in Windows. It will slow down all installs, so

Re: [R] source file for rnorm

2010-12-20 Thread David Winsemius
On Dec 20, 2010, at 5:40 PM, Yue Zhang wrote: Dear Sir/Madam: I am sorry if this is a simple question. I have downloaded the source code for R but could not find the source file that has the rnorm function definition. Could you please let me know which file has the function? Also, in