On Sun, 20 Jan 2013, Geoffrey Smith wrote:
Can anyone please tell me how to get the r-squared output from a
piecewise (segmented) regression using the strucchange package? Here is
the R code I have tried thus far.
library(lmtest)
library(strucchange)
data <- ts(c(rnorm(30), runif(30)), freq
HI,
Your example data had only 4 columns
mydata<-read.table(text="
3.014505 62.96425 3.014505 138.0673
2.817503 56.03400 2.817503 133.3411
2.976227 47.12192 2.976227 139.2438
2.825495 75.05284 2.825495 129.2959
2.793500 52.75190 2.793500 130.9874
3.006333 54.56210 3.006333 136.2982
Hello All,
I have a data file in a text format and there are two data sets. The data
set are continuous.
For each data set there is a header which has the number of data rows and
the name of data series.
For example first data set has "6240 Terry Cove-Model". Then the data for
that series follows
On 01/21/2013 04:39 PM, Roslina Zakaria wrote:
Dear r-users,
I have these data below:
I would like to compare each column with a certain value and count how many in
each column less than that specified value.
odd column (1,3,5) will compare with 1.61 and even column (2,4,6) will compare
with
Hi Henrik,
that's perfect, thanks!
Greetings
Christof
Am 18-01-2013 19:26, schrieb Henrik Bengtsson:
Christof,
I've added support for this to the R.filesets package. In your case,
then all you need to do is:
library("R.filesets")
dlf <- readDataFrame(filename, skip="^year")
No need to spec
Dear r-users,
I have these data below:
I would like to compare each column with a certain value and count how many in
each column less than that specified value.
odd column (1,3,5) will compare with 1.61 and even column (2,4,6) will compare
with 75 and I would like to count how many for e
On Jan 20, 2013, at 3:05 PM, Geoffrey Smith wrote:
Can anyone please tell me how to get the r-squared output from a
piecewise
(segmented) regression using the strucchange package? Here is the R
code I
have tried thus far.
library(lmtest)
library(strucchange)
data <- ts(c(rnorm(30), runif
On Jan 20, 2013, at 3:52 PM, Duncan Murdoch wrote:
On 13-01-20 6:23 PM, David Winsemius wrote:
On Jan 20, 2013, at 11:46 AM, Duncan Murdoch wrote:
On 13-01-20 2:28 PM, Steve Taylor wrote:
From: Duncan Murdoch
Maybe we just need a manual on how to use the existing help
system. But
I suspec
require(lattice)
require(latticeExtra)
data(mtcars)
mt2 <- mtcars
mt2$vs <- factor(mtcars$vs)
mt2$cyl <- factor(mtcars$cyl)
mt2$gear <- factor(mtcars$gear)
lp2 <- levelplot(mpg ~ cyl*gear | vs, data=mt2,
strip=strip.custom(strip.names=c(TRUE,TRUE)),
scales=list(x=
Probably good enough, thanks!
Nick
On 1/20/13 2:03 PM, Duncan Murdoch wrote:
On 13-01-20 4:51 PM, Nick Matzke wrote:
Hi, here's an obscure question someone can hopefully help
with.
I have some R code that uses stuff from parallel (now a part
of the R core in 2.15 I believe), especially cluste
On 01/21/2013 08:59 AM, Murat Tasan wrote:
hi all - i'm having some difficulty figuring out how to convert
between "user units" (which i can't find a definition for in the
plotrix package) and either (a) device units (e.g. inches with PDFs)
or (b) user coordinates along any particular axis.
as a
On 13-01-20 6:23 PM, David Winsemius wrote:
On Jan 20, 2013, at 11:46 AM, Duncan Murdoch wrote:
On 13-01-20 2:28 PM, Steve Taylor wrote:
From: Duncan Murdoch
Maybe we just need a manual on how to use the existing help
system. But
I suspect people who won't read the existing manuals won't rea
You can't.
You shouldn't.
Segmented regression with unknown breakpoints in nonlinear regression,
and R-squared doesn't make sense for nonlinear regression.
See:
https://stat.ethz.ch/pipermail/r-help/2002-July/023461.html
-- Bert
On Sun, Jan 20, 2013 at 3:05 PM, Geoffrey Smith wrote:
> Can anyo
On Jan 20, 2013, at 11:46 AM, Duncan Murdoch wrote:
On 13-01-20 2:28 PM, Steve Taylor wrote:
From: Duncan Murdoch
Maybe we just need a manual on how to use the existing help
system. But
I suspect people who won't read the existing manuals won't read that
one, either.
Duncan, I assume you
On Sun, Jan 20, 2013 at 5:54 PM, John Sorkin
wrote:
> Why are the help pages not right? The ... construct is a fundamental part of
> the language syntax. Information about this fundamental construct should be
> easily available!
This strikes me as a bit harsh. The information is easily availabl
Can anyone please tell me how to get the r-squared output from a piecewise
(segmented) regression using the strucchange package? Here is the R code I
have tried thus far.
library(lmtest)
library(strucchange)
data <- ts(c(rnorm(30), runif(30)), frequency = 12, start = c(2005, 01))
bpts <- breakp
On 13-01-20 5:54 PM, John Sorkin wrote:
Why are the help pages not right? The ... construct is a fundamental
part of the language syntax. Information about this fundamental
construct should be easily available!
Yes, I agree it should be easily available. Why does that mean it has
to be on the
Why are the help pages not right? The ... construct is a fundamental part of
the language syntax. Information about this fundamental construct should be
easily available!
John
John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division
I am trying to create a LaTex table based on a multinom (nnet) object using
the stargazer command.
I have created a small data frame to demonstration the problem:
data <- data.frame(age=1:21, hight=20:40,
ed=c(1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,2,3,1,2,3))
data$ed <- as.factor(data$ed)
I then make a
Perfect! That did it, and I learned a new function :)
On Sun, Jan 20, 2013 at 2:49 PM, John Fox wrote:
> Dear Dustin and Steve,
>
> For some reason, I didn't see Dustin's original message.
>
> As documented in ?sem, the warn argument to sem() covers only warnings
> directly produced by the optim
On 13-01-20 2:56 PM, Steve Taylor wrote:
On 13-01-20 2:28 PM, Steve Taylor wrote:
From: Duncan Murdoch
Maybe we just need a manual on how to use the existing help system. But
I suspect people who won't read the existing manuals won't read that
one, either.
Duncan, I assume you're being faceti
On 13-01-20 4:51 PM, Nick Matzke wrote:
Hi, here's an obscure question someone can hopefully help with.
I have some R code that uses stuff from parallel (now a part
of the R core in 2.15 I believe), especially clusterApply.
However, this seems to cause problems in R.app, and I've
seen advice to
hi all - i'm having some difficulty figuring out how to convert
between "user units" (which i can't find a definition for in the
plotrix package) and either (a) device units (e.g. inches with PDFs)
or (b) user coordinates along any particular axis.
as an example, suppose i set up a PDF device with
Hi, here's an obscure question someone can hopefully help with.
I have some R code that uses stuff from parallel (now a part
of the R core in 2.15 I believe), especially clusterApply.
However, this seems to cause problems in R.app, and I've
seen advice to not use these multicore functions, e.g
Dear Dustin and Steve,
For some reason, I didn't see Dustin's original message.
As documented in ?sem, the warn argument to sem() covers only warnings
directly produced by the optimizer,
"warn: if TRUE, warnings produced by the optimization function will
be printed. This should generally
Have you tried
suppressWarnings(sem(mod, S = as.matrix(dataset), N = 1000, maxiter = 1,
warn = FALSE))
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Dustin Fife
> Sent: Saturday, 19 January 2013 5:32a
> To: r-help
> Su
> On 13-01-20 2:28 PM, Steve Taylor wrote:
> >> From: Duncan Murdoch
> >> Maybe we just need a manual on how to use the existing help system. But
> >> I suspect people who won't read the existing manuals won't read that
> >> one, either.
> >
> > Duncan, I assume you're being facetious. Every R us
On 13-01-20 2:28 PM, Steve Taylor wrote:
From: Duncan Murdoch
Maybe we just need a manual on how to use the existing help system. But
I suspect people who won't read the existing manuals won't read that
one, either.
Duncan, I assume you're being facetious. Every R user soon learns to use
hel
> From: Duncan Murdoch
> Maybe we just need a manual on how to use the existing help system. But
> I suspect people who won't read the existing manuals won't read that
> one, either.
Duncan, I assume you're being facetious. Every R user soon learns to use
help() and help.search() or their equiv
Dear Ilya,
Please check these links.
http://stackoverflow.com/questions/4556524/whats-the-way-to-learn-r
http://www.r-bloggers.com/learn-to-use-r-for-free-with-coursera/
http://stackoverflow.com/questions/192369/books-for-learning-the-r-language
You may also benefit from "An introduction to
On Jan 20, 2013, at 9:27 AM, David Winsemius wrote:
On Jan 20, 2013, at 8:26 AM, Sam Steingold wrote:
* Bert Gunter [2013-01-19 22:26:46 -0800]:
But David W. and Bill Dunlap gave you solutions that also work and
are
much faster, no?!
Yes, indeed, and I am now using David's solution as
On Jan 20, 2013, at 8:26 AM, Sam Steingold wrote:
* Bert Gunter [2013-01-19 22:26:46 -0800]:
But David W. and Bill Dunlap gave you solutions that also work and
are
much faster, no?!
Yes, indeed, and I am now using David's solution as it is fast
(enough), simple and concise.
I am a bit
> * Bert Gunter [2013-01-19 22:26:46 -0800]:
>
> But David W. and Bill Dunlap gave you solutions that also work and are
> much faster, no?!
Yes, indeed, and I am now using David's solution as it is fast
(enough), simple and concise.
Thanks a lot to David, Bill, Rui, and arun for their answers (t
On Sun, Jan 20, 2013 at 1:59 AM, ronny wrote:
> Hi,
>
> I am using levelplot, and would like remove from each panel (condition) its
> unused x levels.
>
Uneven scales on categorical axes lead to distortion and a miss
representation (as in your example - the area for levels 4,6 in vs1 will be
big
Chris,
I've fixed Surv in rms. The fix will be in the next release. For now you
can do source('http://biostat.mc.vanderbilt.edu/tmp/Surv.s') after issuing
require(rms).
Frank
Frank Harrell wrote
> Chris,
>
> Thanks for sending the specifics. It appears that I've let Surv in rms
> fall behind
Deat R users,
I am looking for a way to get standard errors of marginal effects of
multinomial logit model.
So far I have estimated coeficients of multinomial logit and covariance matrix
of coefficients.
I also managed to obtain marginal effects. However, I can not find a way to
find standard
Hi,
thank you very much for your kind answer.
>If you look a bit further down the manual page you will see
>### using a model formula to specify the same model
>rma(yi, vi, mods=~factor(alloc)+year+ablat, data=dat, method="REML",
>btt=c(2,3))
>which is much easier.
I have seen the possibili
At 17:14 19/01/2013, Alma Wilflinger wrote:
Hi,
I am quite new to R and in need of some advice. I am trying to
conduct a meta regression over a some studies with about 7 mod
variables which I have to dummy encode.
Alma, although you can generate your own dummy variables by hand you
do not h
Hi,
I am using levelplot, and would like remove from each panel (condition) its
unused x levels.
e.g.
Remove from panel vs=1 the cyl level=8.
data(mtcars)
levelplot(mpg~factor(cyl)*factor(gear)|factor(vs))
Thanks for your help,
Ronny
--
View this message in context:
http://r.789695.n4.nabbl
39 matches
Mail list logo