However, note that I DID get it WRONG:
The method invoked by the NextMethod() call in the boxcox.formula
method is boxcox.default, NOT boxcox.lm.
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374
Data is not information. Information is not knowledge. And knowledge
On Feb 13, 2014, at 11:36 AM, Ross Boylan wrote:
If R changes from 3.0.1 to 3.0.2, or more generally from m.n.p to m.n.q,
is it necessary to refresh libraries to match the version, e.g., with
update.packages(checkBuilt=TRUE, ask=FALSE)? The R Windows FAQ 2.8 says
For those with a personal
Dear all,
I have a RNA-seq cancer dataset (2 classes) with all raw integer counts,
does NOT contain any NaN or missing value.
When I run samseq and Poission-seq, I got the following error:
Error in quantile.default(prop, c(0.25, 0.75)) :
missing values and NaN's not allowed if 'na.rm' is FALSE
Yes I read the help on NextMethod. In fact, since people frequently
respond with did you read the help I mentioned that I had read the help
in my original post. I'm very grateful for the time and effort that people
put into answering questions, so I always try to answer the question myself
first
Dear R-help,
I try to estimate alpha for a factor that loads onto two items, using
the psych-package (newest version).
The data is from the boston housing data set.
The problem I face can be reproduced by the following approximation of
the correlation matrix:
library(psych)
tt -
In searching for NextMethod on http://www.rseek.org/ I found some helpful
tutorials on S3 and S4 methods.
Between your answer and the tutorials, I think I'm starting to understand.
The NextMethod is just a dispatcher type of thing that doesn't do anything
directly.
I think you're saying that
Hi!
In 'An Introduction to R', section 11.7 on nonlinear least squares fitting,
the following example is given for obtaining the standard errors of the
estimated parameters:
To obtain the approximate standard errors (SE) of the estimates we do:
sqrt(diag(2*out$minimum/(length(y) - 2) *
I found some good explanations of what the Box Cox transform for anyone who
happens upon this thread while searching for boxcox details:
This is the best explanation I found of what how the Box Cox transformation
is calculated. It explains why QR decomposition is used, has nice code
examples
On Fri, Feb 14, 2014 at 7:38 PM, Bert Gunter gunter.ber...@gene.com wrote:
However, note that I DID get it WRONG:
The method invoked by the NextMethod() call in the boxcox.formula
method is boxcox.default, NOT boxcox.lm.
The method chain is slightly unusual here, because of the liberal way
Hi,
Try ?multhist() from library(plotrix)
set.seed(49)
vec1 - rnorm(1e4,0,0.5)
set.seed(9453)
vec2 - rnorm(1e4,0,0.5)
library(plotrix)
multhist(list(vec1,vec2))
A.K.
I've met a problem recently. Have you ever been met a plot like the belowing?
In
this picture, as you see, the solid
On Feb 13, 2014, at 2:07 PM, Gil Gamesh wrote:
Hi,
I'm printing a bunch of summary tables to a latex file using latexVerbatim
from the Hmisc package.
An example looks like this...
x
Visit N Min. 1st Qu. Median Mean 3rd Qu. Max. NAs
[1,] 1 92 25 27.28 29.05 29.47
beer - read.csv(C:/Users/Administrator/Desktop/beer.csv)
View(beer) attach(beer)The following objects are masked from chem_1:
C_SALE, CFO, Company, DISEX, PPE, REV, ROA, SALE, TA, TAC, Year
library(plm) y-cbind(TAC) x-cbind(TA,REV,PPE,ROA)
pdata-plm.data(beer,index=c(Company,Year))
On 02/15/2014 03:19 PM, Ramesh Chandra Das wrote:
beer- read.csv(C:/Users/Administrator/Desktop/beer.csv)
View(beer) attach(beer)The following objects are masked from chem_1:
C_SALE, CFO, Company, DISEX, PPE, REV, ROA, SALE, TA, TAC, Year
library(plm) y-cbind(TAC)
Hello
I have to convert character strings into POSIXt format. And would like to
combine two of them. The following code does not what I expect. The single
conversions of the character strings, gives a the date and time with time zone
GMT as I expect. However if I combine two date time with c()
Dear all,
a simple question about boxplot() function,
on how not to display a set of values (outliers), but keep the same median
of the original dataset on the display ;
more precisely :
given a dataset : 1,2,3,4,5,6, 100, 200,300
the median is : 5
after removing the extreme values : 100,
Hi Raphael
Bug or not that for others to say. Here is an explanation and a work around.
This is on my Windows 8 laptop and R-3.0.2:
## strptime gives you POSIXlt objects (see ?DateTimeClasses or ?striptime)
t1 - strptime(paste(01.01.2013, 20:00:00),format=%d.%m.%Y %H:%M:%S,
tz=GMT)
t2 -
On 02/14/2014 07:52 PM, Bogdan Tanasa wrote:
Dear all,
a simple question about boxplot() function,
on how not to display a set of values (outliers), but keep the same median
of the original dataset on the display ;
more precisely :
given a dataset : 1,2,3,4,5,6, 100, 200,300
the median is :
Thanks Dr. I did read the help page but maybe I just don't understand enough
about contrasts to take what I need from it. It says:
If value supplies more than how.many contrasts, the first how.many are
used. If too few are supplied, a suitable contrast matrix is created by
extending value after
Hi all,
I have a very large number of vectors that I want first to look fast which
distribution might be considered candidate for fitting.
I made a simple loop that checks for all vector (the code below is for one
vector and being called for each vector separately). If a good fit is found
this
that is a great suggestion ! thanks a lot Jim !
On Fri, Feb 14, 2014 at 1:14 AM, Jim Lemon j...@bitwrit.com.au wrote:
On 02/14/2014 08:08 PM, Bogdan Tanasa wrote:
thank you Jim ...the idea behind would be to subset a set of values ?
I would need to automate it the display for a list of
Thanks a lot for the fast answer. Sys.setenv(TZ=GMT) is a good solution for
me.
Best regards
Raphael
-Ursprüngliche Nachricht-
Von: Frede Aakmann Tøgersen [mailto:fr...@vestas.com]
Gesendet: Freitag, 14. Februar 2014 10:02
An: Felber Raphael Agroscope; r-help@r-project.org
Betreff:
Hi
See
?try
?tryCatch
So you still want to do these distribution tests even though you have been
warned not to?
Yours sincerely / Med venlig hilsen
Frede Aakmann Tøgersen
Specialist, M.Sc., Ph.D.
Plant Performance Modeling
Technology Service Solutions
T +45 9730 5135
M +45 2547 6050
Hello,
I do data analysis on json data (Twitter). An example of the data:
**
\id\: 433662713886429200,
\id_str\: \433662713886429184\,
\text\: \Hond vast in water in Bargerveen bij Zwartemeer -
http://t.co/FqbkOMzYd1 #Zwartemeer #bargerveen #hond #innood\,
Hi thanks for the answer.
I have so vast majority of vectors that some distributions will not work for
some cases but they work for others pretty well. This is the reason I store
only those that worked to helpe get a coarse understanding how to proceed. If
you have a better idea let me know
Thank you Jim !
I made a little change in the script when assigning the median values - and
it works ! thank you very much !
x_median-median(x,na.rm=T)
y_median-median(y,na.rm=T)
x-x[x10]
y-y[y10]
xy_box-boxplot(list(x,y),plot=FALSE)
xy_box$stats[3,1]-x_median
xy_box$stats[3,2]-y_median
Hi everyone!I am trying to make some synthetic data using two AR(1) models,
but I am having some troubles.I want to make data from:x_t = \alpha x_{t-1}
+ a_{1t}y_t = \beta y_{t-1} + \gamma a_{1t} + \sqrt{1-\gamma^2} a_{2t}But I
don't know how to set a fixed error term in the arima.sim()
HI ,
I am using R on a very huge datasets which contains lot of text. I prepared
a all word vector of words by using strsplit function. Now I want to compute
frequency of unique words from all word vector. For doing
so, I used two ways
1) as.data.frame(table(x))
2) sapply(x,x,length)
x contains
That sounds like the operating system is terminating R due to overloading the
system. You have not supplied the information requested in the Posting Guide,
so you may not get very specific responses on how to solve this. I would
suggest running partial data sets of larger and larger size to
http://cran.r-project.org/web/packages/MASS/citation.html
---
Jeff NewmillerThe . . Go Live...
DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...
You need to use the JSON library or equivalent to solve this problem. I don't
understand why you think that having the data in the clipboard prevents you
from doing this since that is just another file (but I usually avoid using the
clipboard for reproducible analysis anyway).
Hello,
On Feb 13, 2014, at 2:16 PM, Ross Boylan r...@biostat.ucsf.edu wrote:
Redefining methods of a refClass does not change the methods for existing
instances. I found this somewhat surprising, and also pretty inconvenient
for debugging.
The documentation for reference classes does
Sir,
can you solve my problem. I did not get why this type of problem is there.
duplicate couples (time-id)Error in pdim.default(index[[1]], index[[2]]) :
Thanks and regards
Ramesh Chandra Das
Doctoral Student
Vinod Gupta School of
Hi,
In my function, I want to allow input to be a vector or a data.frame.
Certain operations need to be done if the length or nrows exceeds one, but
since nrow doesn't work for vectors, I cannot simply use
if( nrow(input)1 | length(input)1 ) ...
So is there a more elegant way to do this then
Sorry, psychic powers too weak.
Suggest you read the Posting Guide for this mailing list. You may also benefit
from reading
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example.
---
Jeff
Hi Berry,
What about using
NROW(input)
?
Best,
Jorge.-
On Sat, Feb 15, 2014 at 2:26 AM, Berry Boessenkool
berryboessenk...@hotmail.com wrote:
Hi,
In my function, I want to allow input to be a vector or a data.frame.
Certain operations need to be done if the length or nrows exceeds one,
Yes.
And there's also a possibility of keeping the vector as a one column
data frame (using a drop=FALSE argument in indexing to avoid
conversion to vectors perhaps) so that nrow() can be used. This has
the disadvantage of the inefficiency of using a data frame when a
vector will do, and the
Hi,
Try:
x - 1985:1990
set.seed(45)
y - sample(70:90,length(x),replace=TRUE)
plot(x,y,xlab=Year,ylab=expression(Temperature~degree*C),main=April Average
Temperature 1985-1990,type=p)
A.K.
I want to provide a label for my x,y plot as such: ylab=Temperature
(^oC), where the o is meant to be
URL - http://r.789695.n4.nabble.com/file/n4685033/cylinder.dat;
cyl - read.table(URL, header = TRUE)
plot(vol ~ ht, data = cyl)
Since importing data is one of the more error-prone operations in any
computer system, I like to recommend that people always take a look
at what they read before
Have you tried setting all the rownames to be the empty string? Its a
bit of a kludge, but a simple one.
On Thu, Feb 13, 2014 at 3:07 PM, Gil Gamesh games...@gmail.com wrote:
Hi,
I'm printing a bunch of summary tables to a latex file using latexVerbatim
from the Hmisc package.
An example
Hi,
Check these links:
http://r.789695.n4.nabble.com/quot-duplicate-couples-time-id-quot-Problem-td4636316.html
http://stackoverflow.com/questions/15960129/r-plm-year-fixed-effects-year-and-quarter-data
A.K.
On Friday, February 14, 2014 7:57 AM, Ramesh Chandra Das
rameshchandrada...@gmail.com
The given solution is the right one for the degree symbol, but does not
answer the overall question of how to superscript, in particular how to
superscript with nothing actually preceding the superscript.
In general ^ is used to obtain superscripts in plotmath (RTFM!) and if
nothing comes
Well, since this is really a question about understanding how S3
methods work, and this is not the place for a tutorial, I think what
you need to do is search out a tutorial that you understand.
But very briefly, it does what it says. The object argument is
supplied to the boxcox generic; lm()
Inline.
Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374
Data is not information. Information is not knowledge. And knowledge
is certainly not wisdom.
H. Gilbert Welch
On Fri, Feb 14, 2014 at 4:19 PM, Gene Leynes gleyne...@gmail.com wrote:
In searching for NextMethod on
43 matches
Mail list logo
