[R] Solving a multinomial gompertz partial differential equation in r

2012-12-02 Thread Brandon Breitling
I haven't used r in quite a while but would like to get back into it. I have a problem that I would like to solve with r. I have some multinomial data that looks to follow an asymmetric sigmoidal growth pattern. Solving a multinomial gompertz partial differential equation in r is what I’m after. Wo

Re: [R] comparing two vectors

2012-12-02 Thread soham chakraborty
Berend, I am working with a time series data.The log-likelihood is to be maximized for a set of parameters. Lets say that x1 and x2 are the output from two successive iterations.I am trying to find out, *if(all(abs(x1-x2)<.0001,say))* -- View this message in context: http://r.789695.n4.n

Re: [R] Warning message: In scan(file, what, nmax...)

2012-12-02 Thread David L Carlson
count.fields() is reporting a single field in each row of your file. Are you sure the separator is a tab ("\t")? In your original query you used a semicolon. Try readLines("/Users/kama/Analysis/GDP10.csv", n=10) to look at the first 10 lines. -- Dav

Re: [R] R beginner

2012-12-02 Thread andrewcd
My $.02: re-do an analysis that you did in another software package in R, making sure that you get the same results. A good way to learn any language. -- View this message in context: http://r.789695.n4.nabble.com/R-beginner-tp4651719p4651759.html Sent from the R help mailing list archive at N

[R] How to make the cell background of a table informative?

2012-12-02 Thread dweeb
I wish to create a table in which the foreground is a number (abundance) but the background is a visual representation of a ratio (shallow:deep). In other words, the background of each cell of the table would be two colours or patterns set in a similar fashion to a stacked bar in a bar graph based

Re: [R] error of installing/building an R package (PortfolioAnalytics) on Win 7

2012-12-02 Thread Jack Bryan
Thanks for your reply. On http://cran.r-project.org/bin/windows/Rtools/ Frozen Rtools216.exe R >2.15.1 to R 2.16.x No Frozen means "no available" ? So, PortfolioAnalytics cannot be used on Linux or Win until Mar. 2012 ? Are there substitutes ? Any help will be appreciated. > Date: S

[R] need help for R's installation

2012-12-02 Thread 蔡良良的gmail邮箱
Hi everyone, I met a erro on the HP-UX system when I compile with the R-2.15.2.tar.gz. the command and erro are shown as follows, I try many times hope anyone can help me. Thanks! $ cd R-2.15.2 $ bash ./configure --with-readline=no --with-x=no checking build system type... ia64-hp-hpux

[R] How to find the lenth of more than ten variables?

2012-12-02 Thread killerkarthick
Hi I have one data set with 20 variables. I want to find the length of each variables at a time. Please help me .. Thanks in advance -- View this message in context: http://r.789695.n4.nabble.com/How-to-find-the-lenth-of-more-than-ten-variables-tp4651751.html Sent from the R he

Re: [R] [mgcv][gam] Manually defining my own knots?

2012-12-02 Thread Andrew Crane-Droesch
Hi Simon, Thanks for your help. I've got another question if you don't mind -- is it possible to "swap out" a set of coefficients of a gamObject in order to change the results when that gamObject is plotted? The (silly) example below illustrates that this is possible with the Vp matrix. But

Re: [R] Fwd: How to calculate different groups of varialbes importance level?

2012-12-02 Thread Bert Gunter
What does this have to do with R?! This is nonsense. -- Bert On Sun, Dec 2, 2012 at 6:34 PM, Jim Lemon wrote: > On 12/02/2012 08:18 PM, Solmaz Filiz KARABAĞ wrote: > >> Dear R user! >> I have a small question! >> I have calculated the relative importance of the variables. >> >> However I woul

Re: [R] Problem with figures

2012-12-02 Thread Yihui Xie
fig=TRUE is irrelevant here, and knitr does not need fig=TRUE at all (plots are automatically recorded). The real problem is the [preview] option; remove it and you are all set. Next time if you have problems with tikzDevice, you can take a look at the LaTeX log to know what exactly is wrong, e.g

Re: [R] Probable spammer Re: Problems with nls

2012-12-02 Thread Jim Lemon
On 12/03/2012 02:48 PM, David Winsemius wrote: It appears to me that rahul143 is a spammer who has automatated extraction of old posts to Nabble as a way of delivering advertisements. Hi David, Very good detective work. I hadn't followed the link and seen the advertisement that has been paste

[R] Probable spammer Re: Problems with nls

2012-12-02 Thread David Winsemius
It appears to me that rahul143 is a spammer who has automatated extraction of old posts to Nabble as a way of delivering advertisements. This is an exact copy of a posting to Nabble on Mar 28, 2011; 2:33pm Worik StantonSullivan, Timmerman and White 1999: TA rules, and R [Apologies if I

Re: [R] model selection with spg and AIC (or, convert list to fitted model object)

2012-12-02 Thread Ben Bolker
Adam Zeilinger umn.edu> writes: > > Dear R Help, > > I have two nested negative log-likelihood functions that I am optimizing > with the spg function [BB package]. I would like to perform model > selection on these two objective functions using AIC (and possibly > anova() too). However, th

Re: [R] function stays in loop

2012-12-02 Thread Jim Lemon
On 12/02/2012 10:57 PM, Mkkl wrote: function<- function(p){ i<- 2 r<- 0 while(i< p) { q<- (p/i) while(q>=1) { q=(q-1) } if (q==0) { r=1 } else i=i+1 } return(r) } This f

Re: [R] Spammer? [wa Re: Fitting binomial lmer-model, high deviance and low logLik]

2012-12-02 Thread Jim Lemon
On 12/03/2012 12:18 PM, Ben Bolker wrote: rahul143 gmail.com> writes: I have a strong suspicion that this person is a spammer. When I saw the almost unbroken list of posts, I wondered what might be going on myself. We are already familiar with those who repeatedly ask ill-formed ques

[R] How to mix three different classification techniques

2012-12-02 Thread Tania
Hello all, do you know how to mix three different classification techniques using R. For example if I have lm regression, svm and neural nets how can I mix them to have something like ensemble learning? Thank you, Tania Sent from my iPod __ R-help@r-

Re: [R] Fwd: How to calculate different groups of varialbes importance level?

2012-12-02 Thread Jim Lemon
On 12/02/2012 08:18 PM, Solmaz Filiz KARABAĞ wrote: Dear R user! I have a small question! I have calculated the relative importance of the variables. However I would like to compare the relative importance of two different groups of variables (i.e Strategy and industry) For example let me say t

[R] Spammer? [wa Re: Fitting binomial lmer-model, high deviance and low logLik]

2012-12-02 Thread Ben Bolker
rahul143 gmail.com> writes: > I have a strong suspicion that this person is a spammer. Evidence: 1. 24 e-mails to the r-help list in a little over an hour, which is excessive even for a real person. 2. The questions range extremely widely, especially for a beginner in R. 3. All of

[R] trying to understand the names of the loadings and factor scores from fa()

2012-12-02 Thread Brent Caldwell
Hi I'm trying to understand the names of the output I get from extracting the loadings and the factor scores from fa() when I export to csv the $loadings I get the following columns in this order: MR4, MR3, MR2, MR6, MR7, MR5, MR1, MR8, MR10, MR9, MR11 When I export to csv the $scores I get th

Re: [R] R beginner

2012-12-02 Thread Bert Gunter
Google "Learn R" (What else!) -- And please do not post here further until you have done some studying. -- Bert On Sun, Dec 2, 2012 at 3:55 PM, Rui Barradas wrote: > Hello, > > An Introduction to R, file R-intro.pdf comes with every installation of R, > folder R/doc. > And there are many boo

Re: [R] Warning message: In scan(file, what, nmax...)

2012-12-02 Thread Pascal Oettli
Hello, Search the option 'fill' in the help of the function 'read.table'. Regards, Pascal Le 03/12/2012 04:19, F86 a écrit : Dear R-users, When i try - Data1<-read.table("/Users/kama/Analysis/GDP10.csv",header=TRUE,sep=";") i am getting this error: Warning message: In scan(file, what, nmax,

[R] error of installing/building an R package (PortfolioAnalytics) on Win 7

2012-12-02 Thread Jack Bryan
Hi, I am trying to install a package (PortfolioAnalytics) of R 2.15.2 on Win 7. I have tried the following instructions on : http://stackoverflow.com/questions/11105131/cannot-install-r-forge-package-using-install-packages I used svn checkout svn://svn.r-forge.r-project.org/svnroot/returnana

Re: [R] R beginner

2012-12-02 Thread Rui Barradas
Hello, An Introduction to R, file R-intro.pdf comes with every installation of R, folder R/doc. And there are many books available at http://www.r-project.org/. (Choose the link Other at the left.) Good luck. Hope this helps, Rui Barradas Em 02-12-2012 22:45, Akrem Zoghlami escreveu: Dear

Re: [R] function stays in loop

2012-12-02 Thread Rui Barradas
Hello, Because p/i == 4/2 == 2 thenyou subtract 1 unil q == 0 and i = i+1 is never executed. (Use print statements to check it.) Hope this helps, Rui Barradas Em 02-12-2012 11:57, Mkkl escreveu: function <- function(p){ i <- 2 r <- 0 while(i < p) { q <- (p/i) while(q

Re: [R] Object Browser

2012-12-02 Thread R. Michael Weylandt
Don't double post. I believe R Studio provides this. MW On Sun, Dec 2, 2012 at 4:14 PM, rahul143 wrote: > What's the best object browser? > > Dear all, > > I have tried all the popular R IDE or editors like Eclipse, Komodo, JGR, > Revolution... > They all have fancy fucntions like auto completi

Re: [R] R beginner

2012-12-02 Thread R. Michael Weylandt
On Sun, Dec 2, 2012 at 10:45 PM, Akrem Zoghlami wrote: > Dear R-help group, > > I'm a R beginner and I find difficulty in manipulating R. Could you send to > me a support that helps me to be familiar with R. Thank you in advance > Try typing help.start() at your prompt. This will bring up the m

Re: [R] Linking C/C++ GUI to R.dll

2012-12-02 Thread R. Michael Weylandt
Look at Dirk's RInside package. MW On Sun, Dec 2, 2012 at 4:38 PM, rahul143 wrote: > Hi, > > I am in the process of linking a C/C++ application to the R.dll > directly. I have obtained the R source code and compiled it > successfully. I have also successfully linked the R.dll directly into > o

Re: [R] Repeated-measures anova with a within-subject covariate (or varying slopes random-effects?)

2012-12-02 Thread Bert Gunter
You should have better response posting to the r-sig-mixed-models list (repeated measures are mixed effects models). -- Bert On Sun, Dec 2, 2012 at 12:18 PM, Giuseppe Pagnoni wrote: > Dear all, > > I am having quite a hard time in trying to figure out how to correctly > spell out a model in R (a

Re: [R] finding index of maximum value in vector

2012-12-02 Thread R. Michael Weylandt
?which.max MW On Sun, Dec 2, 2012 at 4:44 PM, rahul143 wrote: > I found: > max.col(matrix(c(1,3,2),nrow=1)) > Is there a more concise/elegant way? > Thanks, > > > > > - >TO GET MORE DETAILS CLICK HERE > > > > > > -- > View this message in context: > http://r.789695.n4.nabble.com/findi

Re: [R] finding index of maximum value in vector

2012-12-02 Thread Jorge I Velez
?which.max On Mon, Dec 3, 2012 at 3:44 AM, rahul143 wrote: > I found: > max.col(matrix(c(1,3,2),nrow=1)) > Is there a more concise/elegant way? > Thanks, > > > > > - >TO GET MORE DETAILS CLICK HERE > > > > > > -- > View this message in context: > http://r.789695.n4.nabble.com/finding

Re: [R] Problem with figures

2012-12-02 Thread Shige Song
Easiest way: copy and paste the code into Rstudio and hit "compile pdf". >From the command line, I believe you can do "knit2pdf example.Rnw". Shige On Sun, Dec 2, 2012 at 6:12 PM, Duncan Murdoch wrote: > On 12-12-02 5:42 PM, Shige Song wrote: > >> I am having problem making ggplot2, tikzDevice,

Re: [R] RMySQL install on windows

2012-12-02 Thread Gabor Grothendieck
On Sun, Dec 2, 2012 at 11:35 AM, rahul143 wrote: > I have been trying to install RMySQL on Windows 7 following the > procedure at: > http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL > > I think I have properly installed RTools and created a proper > Renviron.site file saying: > MYSQL_HOME="C:/Pro

Re: [R] Can you turn a string into a (working) symbol?

2012-12-02 Thread R. Michael Weylandt
I believe I answered this a few weeks ago at this link: https://stat.ethz.ch/pipermail/r-help/2012-November/328053.html and following. Michael On Sun, Dec 2, 2012 at 2:43 PM, rahul143 wrote: > Dear folks-- > > Suppose I have an expression that evaluates to a string, and that that > string, were

Re: [R] Problem with figures

2012-12-02 Thread Duncan Murdoch
On 12-12-02 5:42 PM, Shige Song wrote: I am having problem making ggplot2, tikzDevice, and knitr working together. I used a very simple example: I don't use knitr so I can't really help, but you didn't tell us how you passed this file to knitr, so maybe nobody can. However, if you were using

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
I finally see what you and David W mean...I was under the assumption that by subsetting Year==2012 it would also extract the max for that year.I guess I have to review 'subset' again :) Thank you much Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wild

Re: [R] Warning message: In scan(file, what, nmax...)

2012-12-02 Thread F86
Also tries this: x<-count.fields("/Users/kama/Analysis/GDP10.csv", sep="\t") > x [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [38] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [75] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

Re: [R] Warning message: In scan(file, what, nmax...)

2012-12-02 Thread F86
In addition, I tried Data1<-read.table("/Users/kama/Analysis/GDP10.csv", sep="\t", fill=TRUE) the result: *** caught segfault *** address 0xd020, cause 'memory not mapped' Possible actions: 1: abort (with core dump, if enabled) 2: normal R exit 3: exit R without saving workspace 4: exit

Re: [R] problem with convergence in mle2/optim function

2012-12-02 Thread Adam Zeilinger
Dear Berend, Thank you so much for your help! I was able to write the gradient function for my NLL function. For you're and others' possible interest, here is my final gradient function: Following from my description of the problem below, gr.p1, gr.p2, gr.mu1, and gr.mu2 are the (very larg

[R] R beginner

2012-12-02 Thread Akrem Zoghlami
Dear R-help group, I'm a R beginner and I find difficulty in manipulating R. Could you send to me a support that helps me to be familiar with R. Thank you in advance Sincerely, [[alternative HTML version deleted]] __ R-help@r-project.org m

Re: [R] model selection with spg and AIC (or, convert list to fitted model object)

2012-12-02 Thread Adam Zeilinger
Dear Ravi, Thank you so much for the help. I switched to using the optimx function but I continue to use the spg method (for the most part) because I found that only spg consistently converges give different datasets. I also decided to use AIC rather that a likelihood ratio test. I have a ne

Re: [R] Labeling a range of bars in barplot?

2012-12-02 Thread rahul143
Marc Schwartz (via MN) wrote: > On Tue, 2005-12-13 at 10:53 +, Dan Bolser wrote: > >>Hi, I am plotting a distribution of (ordered) values as a barplot. I >>would like to label groups of bars together to highlight aspects of the >>distribution. The label for the group should be the range of

Re: [R] superimpose density line over hist

2012-12-02 Thread rahul143
> "Romain" == Romain Francois <[hidden email]> > on Tue, 13 Dec 2005 15:40:59 +0100 writes: Romain> A few comments : Romain> - your code should be reproductible, otherwise it is useless. (that Romain> recommandation is on the posting guide) Ro

[R] concatenating expressions and standard text

2012-12-02 Thread rahul143
Hi all, is it possible to concatenate expressions and basic text when for instance labeling axis of a plot? I would like to see something like the concatenation of expression(C[0]) and "for case 1" on my x axis. Obviously a plot(x, y, xlab=paste(expression(C[0])," in case1")) will not work.

[R] How to simulate correlated data

2012-12-02 Thread rahul143
Hello there, I would like to simulate X --Normal (20, 5) Y-- Normal (40, 10) and the correlation between X and Y is 0.6. How do I do it in R? Thank you very much - TO GET MORE DETAILS CLICK HERE -- View this message in context: http://r.789695.n4.n

[R] Fitting binomial lmer-model, high deviance and low logLik

2012-12-02 Thread rahul143
Hello I have a problem when fitting a mixed generalised linear model with the lmer-function in the Matrix package, version 0.98-7. I have a respons variable (sfox) that is 1 or 0, whether a roe deer fawn is killed or not by red fox. This is expected to be related to e.g. the density of red fo

[R] Kalman Filter Forecast using 'SSPIR'

2012-12-02 Thread rahul143
Dear R Users, I am new to state-space modeling. I am using SSPIR package for Kalman Filter. I have a data set containing one dependent variable and 7 independent variables with 250 data points. I want to use Kalman Filter for forecast the future values of the dependen

[R] Hmisc latex cell background color

2012-12-02 Thread rahul143
Dear latex/R-Sweavers, Using the codel below, I can color text in individual cells for latex output. Is there a similar way to get a background shading? My attempts failed because I did not get the closing brace at the right place with Hmisc/latex. library(Hmisc) x <- as.data.frame(diag(rn

[R] Fitting binomial lmer-model, high deviance and low logLik

2012-12-02 Thread rahul143
Hello I have a problem when fitting a mixed generalised linear model with the lmer-function in the Matrix package, version 0.98-7. I have a respons variable (sfox) that is 1 or 0, whether a roe deer fawn is killed or not by red fox. This is expected to be related to e.g. the density of red fo

[R] finding index of maximum value in vector

2012-12-02 Thread rahul143
I found: max.col(matrix(c(1,3,2),nrow=1)) Is there a more concise/elegant way? Thanks, - TO GET MORE DETAILS CLICK HERE -- View this message in context: http://r.789695.n4.nabble.com/finding-index-of-maximum-value-in-vector-tp4651663.html Sent from the R help mailing list a

[R] residual and null deviance of an lme object with correlation structure

2012-12-02 Thread rahul143
Hello, I am attempting to calculate the residual and null deviance of an lme object that includes a corAR1 correlation structure. I tried deviance(lme.object) and it only returned NULL. Can anyone help? Thank you. - TO GET MORE DETAILS CLICK HERE -- View this message in conte

[R] Sullivan, Timmerman and White 1999: TA rules, and R

2012-12-02 Thread rahul143
Friends I am trying to save myself some tedious work. I am processing a paper from "The Journal Of Finance * Vol. LIV, No. 5 October 1999" by Sullivan, Timmerman and White. "Data-Snooping, Technical Trading Rule Performance, and the Bootstrap" I am aiming to reproduce their results us

[R] Linking C/C++ GUI to R.dll

2012-12-02 Thread rahul143
Hi, I am in the process of linking a C/C++ application to the R.dll directly. I have obtained the R source code and compiled it successfully. I have also successfully linked the R.dll directly into our application and have made calls successfully into the R.dll that are included in the samp

[R] Object Browser

2012-12-02 Thread rahul143
What's the best object browser? Dear all, I have tried all the popular R IDE or editors like Eclipse, Komodo, JGR, Revolution... They all have fancy fucntions like auto completion, syntax highlight BUT, I JUST WANT A OBJECT BROWSER! The easiest way to view objects in R console is fix(),

[R] RMySQL install on windows

2012-12-02 Thread rahul143
I have been trying to install RMySQL on Windows 7 following the procedure at: http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL I think I have properly installed RTools and created a proper Renviron.site file saying: MYSQL_HOME="C:/Program Files/MySQL/MySQL Server 5.5" When I try to install

[R] Problems of metafile plots when converting word to pdf file

2012-12-02 Thread rahul143
Hi all, I copy metafile boxplot from R to Word. Then save as Pdf file. But I found there are some unexpected black lines in some plots within this PDF file. Please give me some advice how to get rid of these black lines in PDF files. Thanks a lot, Jim - TO GET MORE DETAILS CLICK HER

[R] e1071 SVM: Cross-validation error confusion matrix

2012-12-02 Thread rahul143
Hi, I ran two svm models in R e1071 package: the first without cross-validation and the second with 10-fold cross-validation. I used the following syntax: #Model 1: Without cross-validation: > svm.model <- svm(Response ~ ., data=data.df, type="C-classification", > kernel="linear", cost=1) >

[R] Sullivan, Timmerman and White 1999: TA rules, and R

2012-12-02 Thread rahul143
Friends I am trying to save myself some tedious work. I am processing a paper from "The Journal Of Finance * Vol. LIV, No. 5 October 1999" by Sullivan, Timmerman and White. "Data-Snooping, Technical Trading Rule Performance, and the Bootstrap" I am aiming to reproduce their results us

[R] postForm() in RCurl and library RHTMLForms

2012-12-02 Thread rahul143
Hi RUsers, Suppose I want to see the data on the website url <- "http://www.nseindia.com/content/indices/ind_histvalues.htm"; for the index "S&P CNX NIFTY" for dates "FromDate"="01-11-2010","ToDate"="02-11-2010" then read the html table from the page using readHTMLtable() I am using this

[R] Problems with nls

2012-12-02 Thread rahul143
I'm trying to fit the Bass Diffusion Model using the nls function in R but I'm running into a strange problem. The model has either two or three parameters, depending on how it's parameterized, p (coefficient of innovation), q (coefficient of immitation), and sometimes m (maximum market share).

[R] Reading PDF files

2012-12-02 Thread rahul143
I need to do text mining on PDF files. I understand there is a readPDF command in tm that can be used. Have read the 2008 posts on converting PDF files to text by Tony Breyal and others. Wondering if the procedure has been standardized in any tutorial or otherwise? Being new to R, I was ab

[R] postForm() in RCurl and library RHTMLForms

2012-12-02 Thread rahul143
Hi RUsers, Suppose I want to see the data on the website url <- "http://www.nseindia.com/content/indices/ind_histvalues.htm"; for the index "S&P CNX NIFTY" for dates "FromDate"="01-11-2010","ToDate"="02-11-2010" then read the html table from the page using readHTMLtable() I am using this

[R] Stepwise analysis with fixed variables

2012-12-02 Thread rahul143
Hello, How can I run a backward stepwise regression with part of the variables fixed, while the others participate in the backward stepwise analysis? Thank you - TO GET MORE DETAILS CLICK HERE -- View this message in context: http://r.789695.n4.nabble.com/Stepwise-analysis-with-fixed-vari

[R] Conjoint Analysis in R??

2012-12-02 Thread rahul143
Pls advise how I can use R in conjoint analysis?? regds Faisal Afzal Siddiqui Karachi, Pakistan Looking for last minute shopping deals? - TO GET MORE DETAILS CLICK HERE -- View this messag

[R] Can you turn a string into a (working) symbol?

2012-12-02 Thread rahul143
Dear folks-- Suppose I have an expression that evaluates to a string, and that that string, were it not a character vector, would be a symbol. I would like a function, call it doppel(), that will take that expression as an argument and produce something that functions exactly like the symbol wou

[R] Object Browser

2012-12-02 Thread rahul143
Dear all, I have tried all the popular R IDE or editors like Eclipse, Komodo, JGR, Revolution... They all have fancy fucntions like auto completion, syntax highlight BUT, I JUST WANT A OBJECT BROWSER! The easiest way to view objects in R console is fix(), but you have no global view of al

Re: [R] simple subset question

2012-12-02 Thread David L Carlson
As David W. guessed. The maximum is in year 2005 not 2012 so no row from 2012 matches the maximum. > subset(winter,Year==2012 & Total==max(Total)) [1] IDWeek Total FrySmolt FryEq Year <0 rows> (or 0-length row.names) > winter[which(winter$Total==max(winter$Total)),] IDWeek Total

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Steve Lianoglou
Similar to Duncan's example, if you have a script "test.R" which looks like so: start script = a1 <- 1:10 a2 <- 101:122 plot(a1, a1) plot(a1, a2) end script == You can source it one way: R> source('test.R', keep.source=TRUE) Error in xy.coords(x, y, xlabel, ylabel, log) : 'x'

[R] Problem with figures

2012-12-02 Thread Shige Song
I am having problem making ggplot2, tikzDevice, and knitr working together. I used a very simple example: ---example.Rnw- \documentclass[preview]{standalone} \begin{document} \begin{figure} <>= library(ggplot2) qplot(displ, hwy, data = mpg, colo

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
Actually, I dput(winter) and is not that big.so here is the whole thing I am trying to accomplish.  winter <- structure(list(IDWeek = c(27L, 28L, 29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 27L, 28L, 29L, 30L, 31L, 32L

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Duncan Murdoch
On 12-12-02 5:02 PM, John Sorkin wrote: Gentleman, This thread has been of great interest. Perhaps I missed part of it, but do far I have not seen an example of code that has line numbers that demonstrates how one can (in some instances) recover the line number of an error. Can I impose upon the

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread John Sorkin
Gentleman, This thread has been of great interest. Perhaps I missed part of it, but do far I have not seen an example of code that has line numbers that demonstrates how one can (in some instances) recover the line number of an error. Can I impose upon the people who contributed to this thread to p

Re: [R] simple subset question

2012-12-02 Thread arun
Hi, >From the ddply() output, you could get the whole row by:  fish1 <- structure(list(Year = 2002:2012, maxTotal = c(1464311L, 1071051L, 714837L, 2115018L, 850491L, 207537L, 321195L, 935599L, 194429L, 157260L, 303259L)), .Names = c("Year", "maxTotal"), row.names = c(NA, -11L), class = "data.fr

[R] marginal structural modeling

2012-12-02 Thread John Sorkin
Does anyone know an R package that implements marginal structural modeling as described by Jamie Robins? A similar question was asked of the list approximately two-years ago and it does not appear that anyone responded. Thanks, John Sorkin John David Sorkin M.D., Ph.D. Chief, Biostatistics an

[R] Change case of factor in data frame

2012-12-02 Thread Audrey
I am trying to write a function to change the case of all of the text in a data frame to lower case. I do not have foreknowledge of the data frame names or the data types of each column. It seems that if one references the data frame by index, then it returns class "data.frame" but if it is refe

[R] Repeated-measures anova with a within-subject covariate (or varying slopes random-effects?)

2012-12-02 Thread Giuseppe Pagnoni
Dear all, I am having quite a hard time in trying to figure out how to correctly spell out a model in R (a repeated-measures anova with a within-subject covariate, I guess). Even though I have read in the posting guide that statistical advice may or may not get an answer on this list, I decided t

[R] Warning message: In scan(file, what, nmax...)

2012-12-02 Thread F86
Dear R-users, When i try - Data1<-read.table("/Users/kama/Analysis/GDP10.csv",header=TRUE,sep=";") i am getting this error: Warning message: In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : number of items read is not a multiple of the number of columns I wonder what Ia

Re: [R] How to calculate mean of every nth time series data with zoo or xts ?

2012-12-02 Thread arun
Hi, You could use: set.seed(5) dat1<-data.frame(Date=seq(from=as.POSIXct("2012-11-01 00:00:00"),to=as.POSIXct("2012-11-03 23:59:00"),  by="1 min"),col2=rnorm(4320,0,1)) library(xts) dat2<-xts(dat1[,-1],order.by=dat1[,1]) res<-ave(dat2[,1],format(index(dat2),"%H:%M"),FUN=mean) #or res<-tapply(dat2

Re: [R] simple subset question

2012-12-02 Thread arun
Hi, I am getting this: x<-subset(fish,Year==2012 & Total==max(Total))  x #   IDWeek  Total   Fry  Smolt  FryEq Year #21 47 303259 34008 269248 491733 2012 A.K. - Original Message - From: Felipe Carrillo To: R. Michael Weylandt Cc: "r-help@r-project.org" Sent: Sunday, December 2,

Re: [R] CreateThread failure since R 2.15.2 (32-bit)

2012-12-02 Thread Martin Studer
I tried with Java 1.6.0_20 and 1.7.0_09. It actually doesn't matter what command I'm running with system. Any command will fail with the same error message. But for the example I posted, yes, I was using "ls" from Rtools. I found that I can only reproduce this with R 2.15.2+ 32-bit on Windows 32-

[R] help setting up crossed data

2012-12-02 Thread derrick
Hello, and thanks for your time reading this. I'm trying to test interactions of my dataset, in which the all of the factors are within the same column. TypeVol 1 CMass -4.598 2 BBack -4.605 3 BMass -4.602 4 CMass -4.601 5 CBack -4.605 6

Re: [R] Daily Time Series, patterns.

2012-12-02 Thread arun
HI, I guess you can check library(forecast). You could also use ?stl(), though not sure whether this will help you in this case. For e.g. set.seed(5) quantity<-sample(c(120:220,NA),699,replace=TRUE) Date=seq(as.Date("2011-01-01"),len=699,by="1 day") dat3<-data.frame(Date=Date,quantity=quantity) l

[R] Fwd: How to calculate different groups of varialbes importance level?

2012-12-02 Thread Solmaz Filiz KARABAĞ
Dear R user! I have a small question! I have calculated the relative importance of the variables. However I would like to compare the relative importance of two different groups of variables (i.e Strategy and industry) For example let me say that strategy has 2 sub varialbes and industry has four

[R] function stays in loop

2012-12-02 Thread Mkkl
function <- function(p){ i <- 2 r <- 0 while(i < p) { q <- (p/i) while(q>=1) { q=(q-1) } if (q==0) { r=1 } else i=i+1 } return(r) } This function stays in some kind of loo

[R] What is "ts" command for?

2012-12-02 Thread Hard Core
Hello, i'm here to have an explenation about time series I have to analize a time series of a stock prices and i don't understand the meaning of "ts" .. i load, for instance mps <- get.hist.quote(instrument = "bmps.mi" , quote="Close", start="2006-01-01", end="2012-09-30", compression="d") an

Re: [R] cubic spline

2012-12-02 Thread Steve Stephenson
Dear all, many thanks for your answer and support!! Steve -- View this message in context: http://r.789695.n4.nabble.com/cubic-spline-tp4651537p4651698.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing lis

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Milan Bouchet-Valat
Le dimanche 02 décembre 2012 à 14:21 -0500, Duncan Murdoch a écrit : > On 12-12-02 9:52 AM, Milan Bouchet-Valat wrote: > > Le dimanche 02 décembre 2012 à 09:02 -0500, Duncan Murdoch a écrit : > >> On 12-12-02 8:33 AM, Milan Bouchet-Valat wrote: > >>> Le dimanche 02 décembre 2012 à 06:02 -0500, Stev

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
I also thought that I could use it but ggplot didn't like it by itself,,,it wants to whole rowI guess it is hard to understand because I didn't post all the ggplot code of what I am trying to accomplish. The strangest thing is that ddply can extract all my Totals by year but when I try to sub

Re: [R] simple subset question

2012-12-02 Thread David Winsemius
The reason I suggested the alternative that I did was because your code could fail when the max(Total) value was not in the subset where Year==2012. -- David On Dec 2, 2012, at 11:34 AM, Felipe Carrillo wrote: Using my whole dataset I get: library(plyr) ddply(winter,"Year",summarise,maxTot

Re: [R] e1071 SVM: Cross-validation error confusion matrix

2012-12-02 Thread Frank Harrell
What do you mean by accuracy? Proportion classified correctly is not a good index of accuracy if that's the problem. Frank rahul143 wrote > Hi, > > I ran two svm models in R e1071 package: the first without > cross-validation and the second with 10-fold cross-validation. > > I used the follow

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
  Using my whole dataset I get: library(plyr) ddply(winter,"Year",summarise,maxTotal=max(Total))    fish <- structure(list(Year = 2002:2012, maxTotal = c(1464311L, 1071051L, 714837L, 2115018L, 850491L, 207537L, 321195L, 935599L, 194429L, 157260L, 303259L)), .Names = c("Year", "maxTotal"), row.nam

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Duncan Murdoch
On 12-12-02 9:52 AM, Milan Bouchet-Valat wrote: Le dimanche 02 décembre 2012 à 09:02 -0500, Duncan Murdoch a écrit : On 12-12-02 8:33 AM, Milan Bouchet-Valat wrote: Le dimanche 02 décembre 2012 à 06:02 -0500, Steve Lianoglou a écrit : Hi, On Sun, Dec 2, 2012 at 12:31 AM, Worik R wrote: What

Re: [R] simple subset question

2012-12-02 Thread William Dunlap
> I am > still getting an error message > >with : > > x <- subset(fish,Year==2012 & Total==max(Total));x > >I get: > >[1] IDWeek Total FrySmolt FryEq Year > ><0 rows> (or 0-length row.names) The above is not an error message. It says that there are no rows satisfying your criteria. Note

Re: [R] simple subset question

2012-12-02 Thread R. Michael Weylandt
On Sun, Dec 2, 2012 at 6:46 PM, Felipe Carrillo wrote: > Works with the small dataset (2 years) but I get the error message with the > whole dataset (12 years of data). I am going to have > to check what's wrong with it...Thanks Off the cuff guess: there's a NA in Total so max(Total) is NA which

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
Works with the small dataset (2 years) but I get the error message with the whole dataset (12 years of data). I am going to have to check what's wrong with it...Thanks Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA http://www

Re: [R] How to calculate mean of every nth time series data with zoo or xts ?

2012-12-02 Thread Rui Barradas
Hello, You should include a data example, using ?dput. Anyway, making up some data, t1 <- as.POSIXct("2012-1-1 0:0:0") t2 <- as.POSIXct("2012-1-4 0:0:0") d <- seq(t1, t2, by = "1 min") x <- rnorm(length(d)) z <- zoo(x, d) aggregate(z, list(format(index(z), "%H:%M")), FUN = mean) Hope this h

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
Sorry, I was trying it to subset from a bigger dataset called 'winter' and forgot to change the variable names when I asked the question. David W suggestion works but the strange part is that I am still getting an error message with :   x <- subset(fish,Year==2012 & Total==max(Total));x I get: [1

Re: [R] missed values

2012-12-02 Thread Joshua Ulrich
It's simpler, easier, and more correct to use zoo::na.approx. R> library(zoo) R> dat1<-read.table(text= + " date value + 80 2006-10 0.1577647 + 81 2006-11 NaN + 82 2006-12 NaN + 83 2007-01 NaN + 84 2007-02 NaN + 85 2007-03 0.2956429", + header=TRUE,stringsAsFact

[R] How to calculate mean of every nth time series data with zoo or xts ?

2012-12-02 Thread 박상규
Hello, I have 1-minute time series stock data and I'd like to calculate mean of every n-th candle data of m-days. result = c(mean of 1th data, mean of 2nd data, ...) mean of 1th data = (1th data of 2012-1-1 + 1th data of 2012-1-2 + 1th data of 2012-1-3) / 3 mean of 2nd data = (2nd data of

Re: [R] simple subset question

2012-12-02 Thread Gerrit Eichner
Hi, Felipe, two typos? See below! On Sun, 2 Dec 2012, Felipe Carrillo wrote:  Hi, Consider the small dataset below, I want to subset by two variables in one line but it wont work...it works though if I subset separately. I have to be missing something obvious that I did not realize before whil

  1   2   >