[R] Odp: xtabs - missing combination

2009-12-13 Thread Petr PIKAL
Hi I am not sure if it solves your problem tab1<-xtabs(value ~ r07 + r08, df1) tab1[xtabs( ~ r07 + r08, df1)==0]<-NA tab1 r08 r07 A B C D A 1 2 4 B 5 6 7 8 C 9 10 11 12 D 13 0 15 16 > however this is slightly different from tapply output. As I understand help page if

[R] R square in NLS-urgent help

2009-12-13 Thread ruchita gupta
Hello I need one urgent help I am trying to fit the Sigmod curve of logistic growth model using NLS estimation. But i do not get the R square value in that even after getting the "Summary" In that case how to compare the fit for 3 models and find which one is better fit?? How to get this R Square

[R] R square in NLS-urgent help

2009-12-13 Thread ruchita gupta
Hello I need one urgent help I am trying to fit the Sigmod curve of logistic growth model using NLS estimation. But i do not get the R square value in that even after getting the "Summary" In that case how to compare the fit for 3 models and find which one is better fit?? How to get this R Square

[R] RPART - printing full splitting rule number on tree plot

2009-12-13 Thread Rebecca O'Leary
Dear R-users I am using RPART package to get regression trees. However having trouble getting the text function to put the full splitting rule number on the plot, instead to puts it in scientific notation. When a covariate has 1e4 or greater number of digits then the splitting rule number

[R] hierarchical model in R

2009-12-13 Thread albina
Hi! I have the following situation of mixture distributions: The number of eggs, N, laid by a female of certain species of butterfly follows a Poisson(F) distribution, where the fertility F follows an exponential distribution with expectation 1/lambda. Each egg hatches with probability p, and ap

[R] Error with hgu133a2.db

2009-12-13 Thread mkna005 mkna005
Hi, I'm trying to get this library working under R2.10.0. I downloaded it and installed it from zip file. i get the following output. > library(hgu133a2.db) Error in library(hgu133a2.db) : there is no package called 'hgu133a2.db' > utils:::menuInstallLocal() package 'hgu133a2.db' successfully unpa

[R] how to use Rengine instance to parse R script String

2009-12-13 Thread hobartlul
Hello everyone I am currently developing a web based R script editor. My idea is to pass the R script command as a string into the backend , then use the Rengine instance to parse the R script command and get the resutls. Do anyone know how to use Rengine instance to parse a R script String? if s

Re: [R] help: forest plots

2009-12-13 Thread Kim Jung Hwa
Thanks for this link. It helps, but I have to make lots of forest plots and these R scripts are not generic. Are you aware of similar functions as forest(), which can take input in the form of estimates, lower_limit, upper_limit? Thanks a lot! ~Kim On Sun, Dec 13, 2009 at 8:12 PM, C.H. wrote:

Re: [R] Need help to complete missing value (Date and Time) in Sp500 Data

2009-12-13 Thread ivan popivanov
Alright, one more idea - try the zoo merge function. According to the documentation, it allows intersection based on the index. Later you can remove the columns that you don't need. Something like: C = merge(A, B, all=FALSE) --- the last parameter request intersect instead of union now

[R] Confused on using expand.grid(), array(), image() and npudens(np) in my case

2009-12-13 Thread rusers.sh
Hi all, I want to use the npudens() function in the np package (multivariate kernel density estimation), but was confused by the several functions in the following codes,expand.grid(),array(),image() and npudensbw(). This confusion will only be generated in >=3 dimensions. I marked the four pla

Re: [R] help with graphing -- Points in my graph are not apparent, always displayed in steps

2009-12-13 Thread philip robinson
Thank you very much, -- Philip David Winsemius wrote: > > > -- View this message in context: http://n4.nabble.com/help-with-graphing-Points-in-my-graph-are-not-apparent-always-displayed-in-steps-tp961629p963227.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Duncan Murdoch
On 13/12/2009 5:22 PM, Carl Witthoft wrote: Duncan Murdoch wrote: "This is by design. Carriage returns have syntactic meaning in R, so the editor shouldn't display phantom ones. In Linux there's no choice because R doesn't control how text is displayed, but the Windows GUI gets it right." I

Re: [R] Fwd: Error with hgu133a2.db

2009-12-13 Thread Martin Morgan
Hi Christoph -- ask questions about Bioconductor packages on the Bioconductor mailing list (provide the output of sessionInfo(), as some of the error message below suggests a seriously out-of-date installation). http://bioconductor.org/docs/mailList.html Here the likely problem is that some of y

Re: [R] help: forest plots

2009-12-13 Thread C.H.
This one does required the metafor package. http://tables2graphs.com/doku.php?id=04_regression_coefficients On Mon, Dec 14, 2009 at 8:13 AM, Kim Jung Hwa wrote: > Hi All, > > I'm fitting a Poisson regression. And I want to plot 95% Confidence Interval > of Regression Estimates. > > After comin

Re: [R] plot merging/compression with R

2009-12-13 Thread Felix Andrews
Hi, If your PDFs are huge, I suggest you use a different type of plot that does not explicitly plot millions of data points, but rather produces some sort of summary. Examples are box plots, or hexagonal binning (hexbin package). -Felix 2009/12/14 Walther, Alexander : > Dear list, > > i just en

[R] Fwd: Error with hgu133a2.db

2009-12-13 Thread mkna005 mkna005
Hi, I'm trying to get this library working under R2.10.0. I downloaded it and installed it from zip file. i get the following output. > library(hgu133a2.db) Error in library(hgu133a2.db) : there is no package called 'hgu133a2.db' > utils:::menuInstallLocal() package 'hgu133a2.db' successfully unpa

Re: [R] how to replace a single backward slash with a doublebackward slash?

2009-12-13 Thread William Dunlap
Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com > -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Sean Zhang > Sent: Sunday, December 13, 2009 4:01 PM > To: David Winsemius > Cc: r-help@r-project.org > Subject: Re: [R

[R] help: forest plots

2009-12-13 Thread Kim Jung Hwa
Hi All, I'm fitting a Poisson regression. And I want to plot 95% Confidence Interval of Regression Estimates. After coming back to original scale (using following formula): exponential(estimate +/- 1.96*SE), at best I can get the output in the form of estimates, lower_limit, upper_limit values.

Re: [R] how to replace a single backward slash with a double backward slash?

2009-12-13 Thread Sean Zhang
David and William, Thanks for your reply which make me know the concept of escape symbols. As David guessed, I was trying to write a function which will accept a path cut from windows explorer. and as you know windows explorer uses "\". e.g., c:\temp\function.r I originally would like that the f

[R] CRAN (and crantastic) updates this week

2009-12-13 Thread Crantastic
CRAN (and crantastic) updates this week New packages * Bergm (1.0) Alberto Caimo http://crantastic.org/packages/Bergm Functions implementing Bayesian estimation for exponential random graph models via exchange algorithm Updated packages lmtest (0.9-26), l

Re: [R] Need help to complete missing value (Date and Time) in Sp500 Data

2009-12-13 Thread Hamid
Hi Ivan, thanks for the code. I think the other way would be to compare the both data series and remove (shorted) the longest. For that I need a function to compare this both Matrices and then cut the rows. I did not find any proper packege in R to compare vectors or matrices( in respect to one

Re: [R] help with graphing -- Points in my graph are not apparent, always displayed in steps

2009-12-13 Thread David Winsemius
On Dec 13, 2009, at 6:06 PM, philip robinson wrote: Sorry to have not provided enough information in my prior post. Here is a summary of the calculation set that I am doing. The histograms have changed since my last post, but I am not so concerned about them. The scatter plots do not sho

Re: [R] forest() in metafor package.

2009-12-13 Thread Kim Jung Hwa
Please ignore this email. Got this link: http://cran.r-project.org/web/packages/metafor/metafor.pdf Thanks, Kim On Sun, Dec 13, 2009 at 6:14 PM, Kim Jung Hwa wrote: > Hi All, > > I'm using forest() from metafor package to plot forest plots. Here is the > code > > e<-1:6 > v<-seq(2,3,.20) > fore

[R] forest() in metafor package.

2009-12-13 Thread Kim Jung Hwa
Hi All, I'm using forest() from metafor package to plot forest plots. Here is the code e<-1:6 v<-seq(2,3,.20) forest(e,v) I want to edit default labels "Study 1", "Study 2" and so on... how can I do it? I tried ?forest(), but couldn't find any details. Any suggestions? Thanks, Kim [[al

Re: [R] Need help to complete missing value (Date and Time) in Sp500 Data

2009-12-13 Thread ivan popivanov
Hi Hamid, Not sure why would you need missing data, but you can definitely generate it:). A few random thoughts: First, you can generate all possible timestamps from the beginning to the end of the interval (skipping weekends and skipping the stamps when the markets are closed). Then for

Re: [R] help with graphing -- Points in my graph are not apparent, always displayed in steps

2009-12-13 Thread philip robinson
Sorry to have not provided enough information in my prior post. Here is a summary of the calculation set that I am doing. The histograms have changed since my last post, but I am not so concerned about them. The scatter plots do not show all of the points, they seem to round off to the most common

Re: [R] Export R output to Word/RTF?

2009-12-13 Thread Wenjie Lee
Oh okay :), thanks a lot! On Sun, Dec 13, 2009 at 5:53 PM, Remko Duursma wrote: > http://en.wikipedia.org/wiki/Software_release_life_cycle#Beta > > > > r > > > - > Remko Duursma > Post-Doctoral Fellow > > Centre for Plants and the Environment > Univ

Re: [R] Export R output to Word/RTF?

2009-12-13 Thread Remko Duursma
http://en.wikipedia.org/wiki/Software_release_life_cycle#Beta r - Remko Duursma Post-Doctoral Fellow Centre for Plants and the Environment University of Western Sydney Hawkesbury Campus Richmond NSW 2753 Dept of Biological Science Macquarie Uni

Re: [R] Export R output to Word/RTF?

2009-12-13 Thread Wenjie Lee
Thanks, can anyone explain a little more on what "beta" version means? Thanks! On Sat, Dec 12, 2009 at 10:59 AM, Frank Bloos wrote: > I am using SWord from statconn: http://rcom.univie.ac.at/download.html > It allows you to put r-commands into Word (similar to odfWeave). Output is > directed to

Re: [R] how to replace a single backward slash with a double backwardslash?

2009-12-13 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Sean Zhang > Sent: Sunday, December 13, 2009 2:11 PM > To: r-help@r-project.org > Subject: [R] how to replace a single backward slash with a > double backwardslash? > > Dear R-h

Re: [R] how to replace a single backward slash with a double backward slash?

2009-12-13 Thread David Winsemius
On Dec 13, 2009, at 5:11 PM, Sean Zhang wrote: Dear R-helpers: Hours ago, I asked how to replace a single forward slash with a double backward slash and recieved great help. Thanks again for all the repliers. In the meantime, I wonder how to replace a single backward slash with a double

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Carl Witthoft
Duncan Murdoch wrote: "This is by design. Carriage returns have syntactic meaning in R, so the editor shouldn't display phantom ones. In Linux there's no choice because R doesn't control how text is displayed, but the Windows GUI gets it right." I have to disagree here. Soft-wrapping is the

[R] plot merging/compression with R

2009-12-13 Thread Walther, Alexander
Dear list, i just encountered a problem concerning the export of multiple plots. is it possible to merge several PNGs into one PDF document? i know that this could be easily done by pdf(), but the outcome of this is /way/ too huge (> 15 MB, four plots) and to my knowledge there's no way to compres

[R] how to replace a single backward slash with a double backward slash?

2009-12-13 Thread Sean Zhang
Dear R-helpers: Hours ago, I asked how to replace a single forward slash with a double backward slash and recieved great help. Thanks again for all the repliers. In the meantime, I wonder how to replace a single backward slash with a double backward slash? e.g., I want change "c:\test" into "c:\

[R] Drawing an axis with a subset of the data

2009-12-13 Thread Wells Oliver
I have a vector of dates: > data$date [1] "2009-04-06 CDT" "2009-04-07 CDT" "2009-04-08 CDT" "2009-04-09 CDT" "2009-04-10 CDT" "2009-04-11 CDT" "2009-04-12 CDT" "2009-04-14 CDT" [9] "2009-04-15 CDT" "2009-04-16 CDT" "2009-04-17 CDT" "2009-04-18 CDT" "2009-04-19 CDT" "2009-04-21 CDT" "2009-04-22

Re: [R] Need help to complete missing value (Date and Time) in Sp500 Data

2009-12-13 Thread Hamid
Hi uwe, here the example: 1) Dataset sp500 has values at 01.05.2008: 30.04.2008 21:56;1.385,63;1.385,63;1.385,14;1.385,14;0 30.04.2008 21:57;1.385,67;1.385,67;1.384,48;1.384,48;0 30.04.2008 21:58;1.384,74;1.384,91;1.384,27;1.384,91;0 30.04.2008 21:59;1.385,05;1.385,81;1.385,05;1.385,08;0 01.05.20

Re: [R] How to resolve include Rcpp.h problems?

2009-12-13 Thread Amine Jallouli
Thank you Romain and Dirk. The problem is solved. It works well. I created the file src/Makevars, then I added these lines: PKG_CPPFLAGS=`Rscript -e 'Rcpp:::CxxFlags()'` PKG_LIBS=`Rscript -e 'Rcpp:::LdFlags()'` Good night. Dirk Eddelbuettel wrote: > > > On 13 December 2009 at 19:37, Romai

[R] Reshape a data set

2009-12-13 Thread dadrivr
I am trying to reshape a data set. Could someone please help me with the reshape, cast, and melt functions? I am new to R and I have tried reading up on how to use the reshape package, but I am very confused. Here is an example of what I am trying to do: subject coder score time [1,]

Re: [R] different randomForest performance for same data

2009-12-13 Thread Uwe Ligges
Häring wrote: Hello, I came across a problem when building a randomForest model. Maybe someone can help me. I have a training- and a testdataset with a discrete response and ten predictors (numeric and factor variables). The two datasets are similar in terms of number of predictor, name of

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Duncan Murdoch
On 13/12/2009 1:30 PM, Viju Moses wrote: Thanks again to all for the persuasive expert statements. @ Duncan, I tried Rterm.exe as advised by you in the other post, hoping for my desired (soft wrapping) behavior. I didn't find it there either. Right, I rarely use it and didn't notice that. R

Re: [R] How to resolve include Rcpp.h problems?

2009-12-13 Thread Dirk Eddelbuettel
On 13 December 2009 at 19:37, Romain Francois wrote: | Hi, | | The file lives in the lib directory of the installed package Rcpp. | | Rscript -e "cat(system.file( 'lib', 'Rcpp.h', package = 'Rcpp' ))" Yes, and we typically use similar tricks to build with Rcpp -- suppose your file is called myf

Re: [R] option for multi thread or multi core machines

2009-12-13 Thread Uwe Ligges
Standard answer: It is possible to make use of more than 1 CPU cores by linking against some optimized BLAS (for matrix computations) or by parallelizing your code. Some packages might be of help, they are listed in the "high performance computing" task view on CRAN: http://cran.r-project.or

Re: [R] Need help creating multiple datasets using loops or vectorization

2009-12-13 Thread Uwe Ligges
Forget the loop and use vocallGrp <- split(vocallsub, vocallsub$Group) Uwe Ligges Anderson, Chris wrote: I am trying to create multiple dataset by group like the following using either a loop or vectorization: vocallGp1<-subset(vocallsub, Group==1) vocallGp2<-subset(vocallsub, Group==2) vo

[R] odfWeave produces output file that OO can't open

2009-12-13 Thread Michael Kubovy
Dear R-helpers, I'm trying to learn how to use odfSweave. Here is my source file (in /Users/mk/myTeach/2010-1-7720/odfWeave): Analysis of the iris Data Created on \Sexpr{date()} <>= # I usually load the libraries first so that any output produced by load

Re: [R] how to replace a single forward slash with a double backward slash in a string?

2009-12-13 Thread David Winsemius
On Dec 13, 2009, at 1:58 PM, Peter Dalgaard wrote: David Winsemius wrote: On Dec 13, 2009, at 1:00 PM, Sean Zhang wrote: Dear R-helpers. Can someone kindly tell me how to replace a single forward slash with double backward slash in a string? i.e., from "a/b" to "a\\b" > sub("/", "

Re: [R] Need help to complete missing value (Date and Time) in Sp500 Data

2009-12-13 Thread Uwe Ligges
Hamid wrote: Hi Uwe, yes you are right, but I just wanted to show the date format I use. The whole data is very huge (about 2 rows). Nevertheless I like first to find the missing dates/Times. Then put 0. I know they are some different approach to impute missing data, but for me is importa

Re: [R] how to replace a single forward slash with a double backward slash in a string?

2009-12-13 Thread Peter Dalgaard
Peter Dalgaard wrote: David Winsemius wrote: On Dec 13, 2009, at 1:00 PM, Sean Zhang wrote: Dear R-helpers. Can someone kindly tell me how to replace a single forward slash with double backward slash in a string? i.e., from "a/b" to "a\\b" > sub("/", "","a/b") #the backslashes n

Re: [R] Need help to complete missing value (Date and Time) in Sp500 Data

2009-12-13 Thread Hamid
Hi Uwe, yes you are right, but I just wanted to show the date format I use. The whole data is very huge (about 2 rows). Nevertheless I like first to find the missing dates/Times. Then put 0. I know they are some different approach to impute missing data, but for me is important to find a way

Re: [R] how to replace a single forward slash with a double backward slash in a string?

2009-12-13 Thread Peter Dalgaard
David Winsemius wrote: On Dec 13, 2009, at 1:00 PM, Sean Zhang wrote: Dear R-helpers. Can someone kindly tell me how to replace a single forward slash with double backward slash in a string? i.e., from "a/b" to "a\\b" > sub("/", "","a/b") #the backslashes need to be doubled becau

Re: [R] do while loop

2009-12-13 Thread annie Zhang
Thanks, Uwe. It's very useful. Annie 2009/12/13 Uwe Ligges > > > annie Zhang wrote: > >> Hi, All, >> >> Is there a 'do while' loop in R for which I can check conditions after? I >> checked and it seems there is only the 'while' function. >> > > > Right, but you can use repeat{} and check condit

Re: [R] How to resolve include Rcpp.h problems?

2009-12-13 Thread Romain Francois
Hi, The file lives in the lib directory of the installed package Rcpp. Rscript -e "cat(system.file( 'lib', 'Rcpp.h', package = 'Rcpp' ))" The maintainer of Rcpp (cc'ed, just in case) will tell you more tricks. Note that Rcpp has its own mailing list. https://lists.r-forge.r-project.org/cgi-bi

Re: [R] do while loop

2009-12-13 Thread Uwe Ligges
annie Zhang wrote: Hi, All, Is there a 'do while' loop in R for which I can check conditions after? I checked and it seems there is only the 'while' function. Right, but you can use repeat{} and check conditions whereever using if() and exit the loop with the "break" control word. Uwe L

Re: [R] Array of legend text with math symbols from predefined variables

2009-12-13 Thread Uwe Ligges
See Ligges, U. (2002): R Help Desk: Automation of Mathematical Annotation in Plots. R News 2 (3), 32-34. http://www.r-project.org/doc/Rnews/Rnews_2002-3.pdf Uwe Ligges Koen wrote: Hello, I am trying to include legend text with math symbols from a predefined character variable that is read i

[R] do while loop

2009-12-13 Thread annie Zhang
Hi, All, Is there a 'do while' loop in R for which I can check conditions after? I checked and it seems there is only the 'while' function. Thanks, Annie [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.

Re: [R] code for [[.data.frame

2009-12-13 Thread Romain Francois
On 12/13/2009 02:42 PM, Guillaume Yziquel wrote: baptiste auguie a écrit : `[[.data.frame` Note that this will only be the source code if you have options "keep.source" and "keep.source.pkgs" set to TRUE, and the environment variable "R_KEEP_PKG_SOURCE" equal to yes. Otherwise you see the

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Viju Moses
Thanks again to all for the persuasive expert statements. @ Duncan, I tried Rterm.exe as advised by you in the other post, hoping for my desired (soft wrapping) behavior. I didn't find it there either. Shouldn't be difficult to get used to breaking lines manually, I guess. Regards, Viju

[R] cross validation/GAM/package Daim

2009-12-13 Thread Kim Vanselow
Dear r-helpers, I estimated a generalized additive model (GAM) using Hastie's package GAM. Example: gam1 <- gam(vegetation ~ s(slope), family = binomial, data=aufnahmen_0708, trace=TRUE) pred <- predict(gam1, type = "response") vegetation is a categorial, slope a numerical variable. Now I want

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Duncan Murdoch
On 13/12/2009 12:43 PM, Viju Moses wrote: Thanks for responding. Problem: When I am typing a command, it does not wrap like it does in this email while typing. As I reach the right border, a horizontal scroll bar appears, and as I keep typing further, the initial part of the line starts disappe

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Ista Zahn
Many applications have what is known as "soft wrapping", where the text wraps visually without inserting a new line. I don't use Windows, so I don't know if the R gui can do this or not, but I still maintain that the proper solution is to break your long input lines manually at places that both mak

Re: [R] how to replace a single forward slash with a double backward slash in a string?

2009-12-13 Thread David Winsemius
On Dec 13, 2009, at 1:00 PM, Sean Zhang wrote: Dear R-helpers. Can someone kindly tell me how to replace a single forward slash with double backward slash in a string? i.e., from "a/b" to "a\\b" > sub("/", "","a/b") #the backslashes need to be doubled because they are the escape

[R] how to replace a single forward slash with a double backward slash in a string?

2009-12-13 Thread Sean Zhang
Dear R-helpers. Can someone kindly tell me how to replace a single forward slash with double backward slash in a string? i.e., from "a/b" to "a\\b" Many thanks in advance. -Sean [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Viju Moses
Thanks for answering. But this is the default behavior when using R (or anything else) in any linux terminal. -Original Message- From: Ista Zahn [mailto:istaz...@gmail.com] Sent: Sunday, December 13, 2009 19:01 To: Viju Moses Cc: Jorge Ivan Velez; r-help@r-project.org Subject: Re: [R] li

Re: [R] How to control the skewness of a heteroscedastic variable?- The solution

2009-12-13 Thread Karl-Oskar Lindgren
Just to make sure that you don't spend your valuable time on the problem that I posted earlier today, I want to let you know that I have found my mistake. I simply forgot to take account of the higher order dependence bw my independent variable and the squared error term. Regards Karl-Oskar _

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Viju Moses
Thanks for responding. Problem: When I am typing a command, it does not wrap like it does in this email while typing. As I reach the right border, a horizontal scroll bar appears, and as I keep typing further, the initial part of the line starts disappearing below the left border of the console.

Re: [R] Repeated Measures Analysis - GLM

2009-12-13 Thread John Fox
Dear Ingo, One approach would be to use the Anova() function in the car package. See ?Anova and in particular the O'Brien and Kaiser example, which is for a more complicated repeated-measures design. If you want to get "type-III" tests (as opposed to the default "type-II" tests), be careful with t

Re: [R] Need help to complete missing value (Date and Time) in Sp500 Data

2009-12-13 Thread Uwe Ligges
Hamid wrote: Dear all, here my problem my be somone can help to solve this. I have tow timeseries from different stock market with different length (diff ca. 4000 ) Now I would to add the missing times of the one of this series with proper time (they are minute data) and set the value to 0 sin

Re: [R] simple ts.plot question

2009-12-13 Thread Uwe Ligges
arka sinha wrote: *Respected Sir, I have a simple question regarding plots of time series in R. I have to plot "conc" against "time" **for each individual and display in the same panel for the in-built dataset "Indometh" in R. * *I have six time series, say subject1.ts, subject2.ts, .

Re: [R] too large dimension problem

2009-12-13 Thread David Winsemius
On Dec 12, 2009, at 11:35 PM, David Winsemius wrote: On Dec 12, 2009, at 8:31 PM, Moohwan Kim wrote: Dear R family When I run the command below, the error message came up. It seems like the problem is about computer capacity. It would be appreciated if anyone could give me a solution.

[R] Repeated Measures Analysis - GLM

2009-12-13 Thread protonringe
Hello to the R world... I have some problems regarding a GLM - repeated measures analysis. I want to test overall differences between AgeClass and Treatment (between subject) with OpenR1+OpenR2+OpenR3 (repeated measures, within subject). The table looks kind like this: AgeClassTreatm

Re: [R] random effects in mixed model not that 'random'

2009-12-13 Thread Thomas Mang
HI, Thanks for your input; see below On 12/13/2009 4:41 PM, Robert A LaBudde wrote: I think what you are finding is that calling a grouping variable a "random effect" is not the same thing as it actually being a random effect. An effect is really only random when it is chosen randomly. Just be

Re: [R] Subtitle & conova post hoc

2009-12-13 Thread Simon . Gingins
Hi, first of all thanks a lot for answering (that was very quick)! Sorry but I actually made a spelling mistake, the post hoc test I want to do (post hoc to a friedman test) is "Conover" and not "conova". My bad! Simon - Original Message - Expéditeur: David Winsemius à: simon.ging..

Re: [R] Subtitle & conova post hoc

2009-12-13 Thread David Winsemius
On Dec 13, 2009, at 9:53 AM, David Winsemius wrote: On Dec 13, 2009, at 9:33 AM, Gabor Grothendieck wrote: Google for r-help title problem and read the entire thread for a number of approaches. This is a minor modification of Gabor's code in that thread that shows it is possible to have

Re: [R] random effects in mixed model not that 'random'

2009-12-13 Thread Robert A LaBudde
I think what you are finding is that calling a grouping variable a "random effect" is not the same thing as it actually being a random effect. An effect is really only random when it is chosen randomly. Just because you don't want to deal with it as a fixed effect (e.g., too many levels) doesn

Re: [R] How to control the skewness of a heteroscedastic variable? - A Correction

2009-12-13 Thread Karl-Oskar Lindgren
When going through my earlier post I find a mistake in the example that I provided. The correct version is provided below. I also start to suspect that my problem is that although the cumulant of a sum of independent variable is the sum of the cumulants, the moments of a sum is not the sum of t

Re: [R] Subtitle & conova post hoc

2009-12-13 Thread David Winsemius
On Dec 13, 2009, at 9:33 AM, Gabor Grothendieck wrote: Google for r-help title problem and read the entire thread for a number of approaches. This is a minor modification of Gabor's code in that thread that shows it is possible to have varying font sizes within the main arguments. Unfort

Re: [R] Subtitle & conova post hoc

2009-12-13 Thread Gabor Grothendieck
Google for r-help title problem and read the entire thread for a number of approaches. On Sun, Dec 13, 2009 at 8:54 AM, wrote: > Hello, > > I just have two questions I cannot figure out. > > 1) I would like to have a subtitle just below the main title > on my graphs and I cannot figure out how

Re: [R] Subtitle & conova post hoc

2009-12-13 Thread David Winsemius
On Dec 13, 2009, at 8:54 AM, simon.ging...@unil.ch wrote: Hello, I just have two questions I cannot figure out. 1) I would like to have a subtitle just below the main title on my graphs and I cannot figure out how to do this, is it possible? or is it only possible to have it at the bottom?

[R] Subtitle & conova post hoc

2009-12-13 Thread Simon . Gingins
Hello, I just have two questions I cannot figure out. 1) I would like to have a subtitle just below the main title on my graphs and I cannot figure out how to do this, is it possible? or is it only possible to have it at the bottom? 2) I need to run Conova post hoc test, do you know if it is cod

Re: [R] code for [[.data.frame

2009-12-13 Thread Guillaume Yziquel
baptiste auguie a écrit : `[[.data.frame` and more generally see Rnews Volume 6/4, October 2006 "Accessing the Sources". HTH, baptiste Thank you so much. Indeed, I've been able to look at the source code of the function from the source code of R. But I was quite keen on knowing how to do

Re: [R] random effects in mixed model not that 'random'

2009-12-13 Thread Thomas Mang
HI, Thanks for your response; yes you are right it's not fully on topic, but I chose this list not only because I am using R for all my stats and so read it anyway, but also because here many statisticians read too. Do you know another list where my question is more appropriate ? For what it's

Re: [R] help with graphing -- Points in my graph are not apparent, always displayed in steps

2009-12-13 Thread David Winsemius
On Dec 12, 2009, at 7:32 PM, philip robinson wrote: http://n4.nabble.com/file/n962784/egraph_rules_list_2.png The red X's that are listed above represent data points that have been completely ignored in production of my graph. I know that in a specified sentence length there are variati

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Ista Zahn
H Viju, Automatically breaking _input_ lines is dangerous, because line breaks have meaning in R. It sounds to me like you need to use the Return key to insert line breaks at the appropriate places. -Ista On Sat, Dec 12, 2009 at 11:29 PM, Viju Moses wrote: > Dear Jorge, Thanks. > > But. > > (My

Re: [R] debug an error that incapacitates R?

2009-12-13 Thread Liviu Andronic
Hello And thank you for the quick answer. On 12/13/09, Duncan Murdoch wrote: > I think you need to go back to old-fashioned debugging methods. Identify > the line the triggers the error, by using debug() and single stepping > through sosInit(), or by adding print() or cat() statements to it (or

Re: [R] lattice - 2 graphs per panel with 2 differents colours

2009-12-13 Thread S Devriese
On 12/11/2009 02:03 PM, PtitBleu wrote: > Hello, > > I would like to get a lattice plot of 8 panels (unique(df$fact)=8) with 2 > graphs (df$y1 and df$y2 as a function of df$x) and 1 red point at (500, > ymax) per panel. > > The script below is quite ok but I'm not able to define two different color

Re: [R] Regularized gamma function/ incomplete gamma function

2009-12-13 Thread Stefan Evert
Indeed, it seems that the author of zipfR has neither been aware that the (scaled / aka regularized) incomplete gamma (and beta, for that matter!) functions have been part of R all along. ... ... well , inspecting his code reveals he did know it. But why then on earth provide all the n

[R] How to control the skewness of a heteroscedastic variable?

2009-12-13 Thread Karl-Oskar Lindgren
Dear listusers, I don't know whether my problem is statistical or computational, but I hope I could recieve some help in either case. I'm currently working on a MC-simulation in which I would like to control the skewness of a heteroscedastic dependent variable defined as: y=d*z+sqrt(.5+.5*x^

Re: [R] debug an error that incapacitates R?

2009-12-13 Thread Duncan Murdoch
On 13/12/2009 5:31 AM, Liviu Andronic wrote: Dear all How should I attempt debugging this error? sosInit() Error: invalid connection When this happens R refuses to run anything: 2+2 Error: invalid connection sessionInfo () Error in stdout() : invalid connection Debugging it seems impossi

Re: [R] Easily switch between R GUI and R term on Windows?

2009-12-13 Thread Duncan Murdoch
On 12/12/2009 11:16 PM, Viju Moses wrote: Is there an easy way to switch between R GUI and R terminal? I'm currently using R GUI 2.10.0 on Windows Vista and would like to see something like the R terminal found in Linux. But I may also want to switch back to GUI. I don't remember seeing the termi

Re: [R] lines don't wrap. must scroll horizontally to see/edit a long line in R GUI

2009-12-13 Thread Duncan Murdoch
On 12/12/2009 10:34 PM, Viju Moses wrote: I'm facing this problem on R GUI version 2.10.0 on Windows Vista. I have not changed Windows settings or R GUI settings much except to change from MDI to SDI. Someone else reported this problem a few months ago: https://stat.ethz.ch/pipermail/r-help/20

[R] debug an error that incapacitates R?

2009-12-13 Thread Liviu Andronic
Dear all How should I attempt debugging this error? > sosInit() Error: invalid connection When this happens R refuses to run anything: > 2+2 Error: invalid connection > sessionInfo () Error in stdout() : invalid connection Debugging it seems impossible. > options(error=recover) > sosInit() Error

Re: [R] random effects in mixed model not that 'random'

2009-12-13 Thread Daniel Malter
Hi, you are unlikely to (or lucky if you) get a response to your question from the list. This is a question that you should ask your local statistician with knowledge in stats and, optimally, your area of inquiry. The list is (mostly) concerned with solving R rather than statistical problems. Best

Re: [R] code for [[.data.frame

2009-12-13 Thread baptiste auguie
`[[.data.frame` and more generally see Rnews Volume 6/4, October 2006 "Accessing the Sources". HTH, baptiste 2009/12/13 Guillaume Yziquel : > Hello. > > I'm currently trying to wrap up data frames into OCaml via OCaml-R, and I'm > having trouble with data frame subsetting: > >> # x#column 1;;

Re: [R] A random number from any distribution??

2009-12-13 Thread Peter Dalgaard
ivan popivanov wrote: Hello, I have some data, and I want to generate random numbers following the distribution of this data (in other words, to generate a synthetic data set sharing the same stats as a given data set). Reading an old thread I found the following text: If you can compute the

Re: [R] confint for glm (general linear model)

2009-12-13 Thread Peter Dalgaard
David Winsemius wrote: On Dec 12, 2009, at 8:19 PM, casperyc wrote: for an example, counts <- c(18,17,15,20,10,20,25,13,12) outcome <- gl(3,1,9); treatment <- gl(3,3) glm.D93 <- glm(counts ~ outcome + treatment, family=poisson()) confint(glm.D93) confint.default(glm.D93) # based on asympto

[R] need a solution to an R-problem: consultant available?

2009-12-13 Thread Andre Zimmermann
I am trying to get confidence bands for a non-linear power function (y=mx^b). I thought I should be able to figure it out, but can't. Are there any R consultant? I would be willing to pay some amount of money, but not sure such consultants exist. I fit power functions to lots of data, and this

[R] A random number from any distribution?‏

2009-12-13 Thread ivan popivanov
Hello, I have some data, and I want to generate random numbers following the distribution of this data (in other words, to generate a synthetic data set sharing the same stats as a given data set). Reading an old thread I found the following text: >If you can compute the quantile function o

Re: [R] help with graphing -- Points in my graph are not apparent, always displayed in steps

2009-12-13 Thread philip robinson
http://n4.nabble.com/file/n962784/egraph_rules_list_2.png The red X's that are listed above represent data points that have been completely ignored in production of my graph. I know that in a specified sentence length there are variations in number of Conjunctions -- View this message in context:

[R] code for [[.data.frame

2009-12-13 Thread Guillaume Yziquel
Hello. I'm currently trying to wrap up data frames into OCaml via OCaml-R, and I'm having trouble with data frame subsetting: # x#column 1;; Erreur dans (function(x, i, exact) if (is.matrix(i)) as.matrix(x)[[i]] else .subset2(x, : l'élément 1 est vide ; la partie de la liste d'argument

[R] A random number from any distribution?

2009-12-13 Thread ivan popivanov
Hello, I have some data, and I want to generate random numbers following the distribution of this data (in other words, to generate a synthetic data sets having the stats of a give data set). Reading an old thread I found the following text: >If you can compute the quantile function of

[R] How to resolve include Rcpp.h problems?

2009-12-13 Thread Amine Jallouli
Hi, I am Linux Ubuntu 9.04 user. I'm using R-2.6. Actually, I am developing R package for reading/writing 3D images. I needed to use the R package Rcpp in order to make profit of available C/C++ codes. My problem is that my C code is not able to include the Rcpp.h In fact, when I run this comma

  1   2   >