Hello Sachin,
a sequential testing procedure is described in the useR! book:
@Book{,
title = {Analysis of Integrated and Cointegrated Time Series with R},
author = {B. Pfaff},
publisher = {Springer},
edition = {First},
address = {New York},
year = {2006},
note = {ISB
Dear list,
In some places (for example,
http://en.wikipedia.org/wiki/Support_vector_machine) , the polynomail
kernel in SVM is written as (u'*v + 1)^d, while in the document of
svm() in e1071 package, the polynomial kernel is written as
(gamma*u'*v + coef0)^d. I am a little confused here:
When d
> "JohnT" == Thaden, John J <[EMAIL PROTECTED]>
> on Thu, 6 Jul 2006 12:29:42 -0500 writes:
JohnT> Martin Maechler replied to my query "Warning while subsetting...":
MartinM> > "JohnT" == Thaden, John J <[EMAIL PROTECTED]>
MartinM> > on Thu, 6 Jul 2006 00:02:10
Dear friends,
The "s" in the following argument don't have a variable name, how should i
give it a name?
> s<-data.frame(seq(1,6,by=2))
> s
seq.1..6..by...2.
1 1
2 3
3 5
thanks very much!
--
Kind Regards,
Zhi Jie,Zhang ,PHD
Department of Epidem
Hi,
try this:
formatEng <- function(x) {
s<-as.numeric(strsplit(format(x, scientific=T),"e")[[1]])
return(paste(s[1]*10^(s[2]%%3),as.integer(s[2]-(s[2]%%3)),sep="e"))
}
>>
>> Some examples:
>>
>> 1635 000 000 => 1.635E9
>> 163 500 000 => 163.5E6
>> 0.000 000 000 135 != 135E-9
>>
s<-data.frame(var.name=seq(1,6,by=2))
s
var.name
11
23
35
2006/7/7, zhijie zhang <[EMAIL PROTECTED]>:
Dear friends,
The "s" in the following argument don't have a variable name, how should i
give it a name?
> s<-data.frame(seq(1,6,by=2))
> s
seq.1..6..by...2.
1
I have not seen a reply, so I will offer a few suggestions, even
though I've never used the 'polspline' package. I scanned several of
the help pages and looked in "~\library\polspline' where R is installed
on my hard drive and found no further documentation, and I found nothing
new f
flash johny gmail.com> writes:
>
> Hi,
> I am working on project to build a model for analysing Mortgage-backed
> securities.
> I am doing this as a part of my summer project.
Try
http://cran.r-project.org/src/contrib/Views/
Dieter
__
R-help@stat.m
Bi-Info (http://members.home.nl/bi-info) wrote:
> Dear Users,
>
> I'm a complete novice to R.
>
> I need to do a crosstabs in R, but my data is almost completely
> alphanumeric (with some variables scaled). The Table routine does not
> seem to accept alphanumeric data. What should I do? Do I ne
Wuming Gong wrote:
> Dear list,
>
> In some places (for example,
> http://en.wikipedia.org/wiki/Support_vector_machine) , the polynomail
> kernel in SVM is written as (u'*v + 1)^d, while in the document of
> svm() in e1071 package, the polynomial kernel is written as
> (gamma*u'*v + coef0)^d. I a
Dear List,
I apologize in advance if this is silly. I tried to replicate an analysis I
did previously in SPSS using R, and was surprised to find different results.
So my question is: shouldn't the following SPSS syntax
REGRESSION
DEPENDENT INC89
/METHOD=ENTER hiedyrs experien SE93rec.
Yeld
"Celso Barros" <[EMAIL PROTECTED]> writes:
> Dear List,
>
> I apologize in advance if this is silly. I tried to replicate an analysis I
> did previously in SPSS using R, and was surprised to find different results.
>
>
> So my question is: shouldn't the following SPSS syntax
>
>
> REGRESSION
If results are *very* different, I would suspect that R treated "experien" as a
factor (when you use "read.spss" the default is to use value labels, which
makes R treat such variables as factors --> set "use.value.labels = F").
Another explanation could be missing data in SPSS (not sysmis, but co
Le 23.06.2006 10:42, Martin Maechler a écrit :
> {BTW: Did you know that to *search* mailing list archives of
> such R-SIG-foo mailing lists, you can use google very
> efficiently by prepending the mailing list name and 'site:stat.ethz.ch'?
> e.g., use google search on
> "R-
How to convert the following ms() in Splus to Optim in R? The "Calc" function
is also attached.
ms(~ Calc(a.init, B, v, off, d, P.a, lambda.a, P.y, lambda.y,
10^(-8), FALSE, 20, TRUE)$Bic,
start = list(lambda.a = 0.5, lambda.y = 240),
control = list(maxiter = 10, tol = 0.1))
Ca
Madam,
I've read your answer :
///
Here are other alternatives that are easy to implement and that you
should consider depending on what you want to do:
1- one is to use R in batch mode ie you create a file in your java code
with all the R commands a
Clever...
Good observation on the last example, mistake on my part.
Thanks for all the suggestions.
Sam
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Hans-Joerg Bibiko
Sent: Friday, July 07, 2006 4:27 AM
To: r-help@stat.math.ethz.ch
Subject: Re: [R] e
Dear all,
I am trying to find a suitable R-function for
parametric proportional hazard regressions. The
package survival contains the coxph() function which
performs a Cox regression which leaves the base hazard
unspecified, i.e. it is a semi-parametric method. The
package Design contains the func
Hi all,
I have a three columned list that I have imported into R. The first column
is a plot (ex. Plot1), the second is a species name (ex ACERRUB) and the
third a numeric value. I want to replace some of the second column names
with other names (for example replace ACERRUB with ACERDRU). The o
Hi,
The following produces almost exactly what I needed. The problems are
that the 'panel.dotplot' call (commented) generates the error 'Error in
NextMethod("[") : argument "subscripts" is missing, with no default'. The
other problem is that the colors alternate between the levels of the 'site'
On Fri, 2006-07-07 at 11:20 -0400, Wade Wall wrote:
> Hi all,
>
> I have a three columned list that I have imported into R. The first column
> is a plot (ex. Plot1), the second is a species name (ex ACERRUB) and the
> third a numeric value. I want to replace some of the second column names
> wit
Hi,
I am trying to plot odds ratios and the corresponding confidence
intervals in horizontal segments. It would be ideal if the confidence
interval segment can be drawn with little vertical bars at both ends. I
have tried very hard to change the type of ends by using 'lend'
arguments, but cannot m
Naras presented something with similar functionalities using SVG a couple of
years ago, and since then the svgdevice package had been made available. I
do not know SVG myself, but according to Naras, it's all possible.
Best,
Andy
From: Terry Therneau
>
> I have an application where the Splus
Dear Peter and Florian,
I think you have spotted the problem, thank you very much. I checked
the frequencies, and R is not understanding that the (-1) value in a binary
variable was defined in SPSS as "sysmis". This was precisely the variable
that showed the larger difference (it's the only
On Fri, 7 Jul 2006, Valentin Dimitrov wrote:
> I am trying to find a suitable R-function for
> parametric proportional hazard regressions. The
> package survival contains the coxph() function which
> performs a Cox regression which leaves the base hazard
> unspecified, i.e. it is a semi-parametric
Valentin,
Have you tried survreg() in the Design library?
Regards,
-Cody
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Valentin Dimitrov
Sent: Friday, July 07, 2006 10:1 AM
To: r-help@stat.math.ethz.ch
Subject: [R] parametric proportional hazard re
Hello,
I have a question about Rpad Server Installation on Windows XP.
Here is a guide
http://www.rpad.org/Rpad/ServerNotes.html
and I have been reading this more than ten times to figure this out, but
still having trouble executing a Rpad website properly.
Problem: for example, if I go to
On Fri, 2006-07-07 at 11:47 -0400, Lu, Jiang Jane wrote:
> Hi,
>
> I am trying to plot odds ratios and the corresponding confidence
> intervals in horizontal segments. It would be ideal if the confidence
> interval segment can be drawn with little vertical bars at both ends. I
> have tried very ha
You need angle=90. From the help page:
angle: angle from the shaft of the arrow to the edge of the arrow
head.
On Fri, 7 Jul 2006, Lu, Jiang Jane wrote:
> Hi,
>
> I am trying to plot odds ratios and the corresponding confidence
> intervals in horizontal segments. It would be idea
Lu, Jiang Jane wrote:
> Hi,
>
> I am trying to plot odds ratios and the corresponding confidence
> intervals in horizontal segments. It would be ideal if the confidence
> interval segment can be drawn with little vertical bars at both ends. I
> have tried very hard to change the type of ends by us
It seems to work if you handle the col= argument explicitly.
Suggest you double check this since I have not carefully done so:
barley$yield2 <- with(barley, yield + 5)
dotplot(site ~ yield | variety, data=barley, groups=year,
yield2=barley$yield2, col=c("gray", "black"),
panel=funct
On Fri, 2006-07-07 at 11:47 -0400, Lu, Jiang Jane wrote:
> Hi,
>
> I am trying to plot odds ratios and the corresponding confidence
> intervals in horizontal segments. It would be ideal if the confidence
> interval segment can be drawn with little vertical bars at both ends. I
> have tried very ha
Thank you all for such a prompt response. I misunderstood the meaning of
'lend' in par(). Now I have got the ideal confidence segments. I
appreciate all the specific explanations you have offered.
Sincerely,
Jane
-Original Message-
From: Gavin Simpson [mailto:[EMAIL PROTECTED]
Sent: Fri
The format is like this.
Plot species Value
P1 ACERRUB 3
P2 MAGNVIR2
P3 ARONARB 2
etc.
imported using x<-read.table(file="filename.txt")
I want to replace a list of values in the 2nd column with another list. For
example, I want to replace ARONARB with PHOTPYR.
list<-read.table(fi
Dear list,
I know this is really basic question but I just couldn't get anything
to work. (I did a R site search with keywords "zoo" and "data frame" but
the server timed out on me.)
I have a time series which has the following (typical) format,
DATE Open High Low
Thanks for re-posting onlist. My offlist reply:
The 'list' argument, as per ?replace, needs to be a vector of one or
more indices into another vector, generally the source vector argument
'x', not the actual values that you want to replace. 'list' can either
be explicit integers as the indices, o
Check out read.zoo in the zoo package.
On 7/7/06, Horace Tso <[EMAIL PROTECTED]> wrote:
> Dear list,
>
> I know this is really basic question but I just couldn't get anything
> to work. (I did a R site search with keywords "zoo" and "data frame" but
> the server timed out on me.)
>
> I have a time
Wade,
Given that you appear to have multiple search and replace items to deal
with, here is a possible loop based "Global Search and Replace"
solution:
gsr <- function(Source, Search, Replace)
{
if (length(Search) != length(Replace))
stop("Search and Replace Must Have Equal Number of Items\
Claire,
I know you emailed me a question earlier and I just haven't gotten around to
answering it. For both perl and rgui.exe (and rterm.exe) your path needs to
state the location of these files. (I am not sure about this, but after you
change the path you may even need to reboot for the change
I have a large R program that I am constantly running over and over again. At
the beginning of this program, I create a hige matrix and a huge dataframe but
these are constant. What I mean by constant is that, if I run the program over
later, I really should just use the old matrix and datafram
Thanks Gabor. I figured out what went wrong. The culprit turns out to be
the headers in my data. read.zoo doesn't recognize column headers and
complains
Error in read.zoo("C:\\...\\table.csv", :
index contains NAs
Or is there an option as in read.table(x, header=...) ?
After the
read.zoo supports all the arguments that read.table supports including
header=. From the read.zoo help file:
...: further arguments passed to 'read.table'.
Regards.
On 7/7/06, Horace Tso <[EMAIL PROTECTED]> wrote:
> Thanks Gabor. I figured out what went wrong. The culprit turns out to be
> the
There is an example of saving just two objects in the example
section of ?save
To check whether an object of a given name exists see ?exists
On 7/7/06, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
>
> I have a large R program that I am constantly running over and over again. At
> the beginning o
[EMAIL PROTECTED] wrote:
> I have a large R program that I am constantly running over and over again. At
> the beginning of this program, I create a hige matrix and a huge dataframe
> but these are constant. What I mean by constant is that, if I run the program
> over later, I really should just
I am using the as.factor command to use with glm. When I use the command
>maj <- as.factor(data.logistic$Majors)
>maj
I receive the following output:
[1] M M N M M M M N N M M M N M M M M M M M M M M M N M N N M M N M
M N M M M M M
[40] N M N M M N M M M N M N M N M N N N M N M M M M M M N M
I think the ggplot package is extremely promising.
Parts of the dokumentation are very good alread, e.g I recently managed
to write my first own grob function.
One thing I am missing in the dokumentation is a little more detail on
how to modify the colours of a plot.
Specifically, if I have:
x
On Fri, 2006-07-07 at 12:57 -0500, [EMAIL PROTECTED] wrote:
>
> 1) i can't find an example somewhere of
> just saving two objects rather than the whole session.
> i've looked and looked and i can't find it.
?save
> 2) if i am able to save these two objects, i was hoping that it
> would be possi
I'm not sufficiently familiar with 'trellis' / 'lattice' to provide
an easy, complete answer to your question, but I can explain the current
behavior and provide a hack to get you what you want. With luck,
someone else will suggest an improvement.
First, let's decompose you
On 7/7/06, Sebastian Luque <[EMAIL PROTECTED]> wrote:
> Hi,
>
> The following produces almost exactly what I needed. The problems are
> that the 'panel.dotplot' call (commented) generates the error 'Error in
> NextMethod("[") : argument "subscripts" is missing, with no default'.
It's just as it s
Could you explain what panel.groups= does and what the difference
is between panel.groups= and panel= ? In ?xyplot it just says:
panel.groups: useful mostly for 'xyplot' and 'densityplot'. Applies
when 'panel' is 'panel.superpose' (which happens by default
in these cases if 'g
Dear WizaRds, dear Thomas,
First of all, I want to tell you how grateful I am for all your
support. I wish I will be able to help others along one day the same way
you do. Thank you so much. I am struggling with a multistage sampling
design:
library(survey)
multi3 <- data.frame(cluster=c(
>From: Mark Hempelmann <[EMAIL PROTECTED]>
>Date: Fri Jul 07 14:05:29 CDT 2006
>To: r-help@stat.math.ethz.ch
>Subject: [R] Multistage Sampling
i also find it an truly amazing group also. the general kindness
and generosity of everyone is beyond belief. it will be a long time coming but
i hope i
On 7/7/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> Could you explain what panel.groups= does and what the difference
> is between panel.groups= and panel= ? In ?xyplot it just says:
>
> panel.groups: useful mostly for 'xyplot' and 'densityplot'. Applies
> when 'panel' is 'panel.s
jdrapp,
By default, R fits full rank models. If you are coming from SAS, you're
probably used to less than full rank model parameterizations.
>From Section 11.1.1 of "An Introduction to R" at
http://cran.r-project.org/doc/manuals/R-intro.html#Contrasts
there is this:
"What about a k-level f
One correction... since you are fitting a logistic model, it is
technically correct to say the "mean value of the linear predictor,"
instead of "mean response".
20 lashes for me.
Max
-Original Message-
From: Kuhn, Max
Sent: Friday, July 07, 2006 4:11 PM
To: 'r-help@stat.math.ethz.ch
"Kuhn, Max" <[EMAIL PROTECTED]> writes:
>
> One correction... since you are fitting a logistic model, it is
> technically correct to say the "mean value of the linear predictor,"
> instead of "mean response".
>
> 20 lashes for me.
...plus five more for forgetting that the link is nonlinear a
On 7/7/06, Deepayan Sarkar <[EMAIL PROTECTED]> wrote:
> On 7/7/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> > Could you explain what panel.groups= does and what the difference
> > is between panel.groups= and panel= ? In ?xyplot it just says:
> >
> > panel.groups: useful mostly for 'xyplot'
On Fri, 7 Jul 2006, Mark Hempelmann wrote:
> library(survey)
> multi3 <- data.frame(cluster=c(1,1,1,1 ,2,2,2, 3,3), id=c(1,2,3,4,
> 1,2,3, 1,2),
> nl=c(4,4,4,4, 3,3,3, 2,2), Nl=c(100,100,100,100, 50,50,50, 75,75),
> M=rep(23,9),
> y=c(23,33,77,25, 35,74,27, 37,72) )
>
> dmulti3 <- svydesign(id=~cl
Since I haven't seen an answer to this, I'll offer a couple of
comments. I don't recall having heard the term 'harmonic regression'
prior to seeing your email, but it sounded interesting, so I did some
searching. First, RSiteSearch("harminic regression") produced 7 hits,
one of whic
I have not seen it myself (I have the first edition which uses FORTRAN)
but I believe the second edition of Peter Bloomfield's book on Fourier
Analysis contains harmonic regression code in S-Plus and that may
work in R.
On 7/7/06, Spencer Graves <[EMAIL PROTECTED]> wrote:
Since I haven'
Hello. I am a beginer in R and I can not implement the KhmaladzeTest in the
following command. Please help me!!!
PD: I attach thw results and the messages of the R program
R : Copyright 2006, The R Foundation for Statistical Computing
Version 2.3.1 (2006-06-01)
ISBN 3-900051-07-0
>
> Those are not proportional hazards families of
> distributions. That is, if
> the distribution is gaussian for one value of the
> hazard ratio parameters
> it will not be gaussian for any other value.
>
> You can get accelerated failure models with these
> distributions using
> survreg()
Cody,
I have tried the survreg() in the Design library, it
is analogous to survreg() in the survival library and
it seems to me it is designed only for accelerated
time models like the accelerated failure model (or
accelerated lifetime model) and not for proportional
hazard models. Correct me if I
What problem are you trying to solve? The mixture proportion you are
trying to estimate violates the assumptions that provide a normal
approximation to the distribution of maximum likelihood estimators.
These situations even violate the assumptions for 2*log(likelihood
ratio) being a
Valentin Dimitrov wrote:
>
>> Those are not proportional hazards families of
>> distributions. That is, if
>> the distribution is gaussian for one value of the
>> hazard ratio parameters
>> it will not be gaussian for any other value.
>>
>> You can get accelerated failure models with these
>> di
Valentin Dimitrov wrote:
> Cody,
>
> I have tried the survreg() in the Design library, it
> is analogous to survreg() in the survival library and
> it seems to me it is designed only for accelerated
> time models like the accelerated failure model (or
> accelerated lifetime model) and not for prop
Dear all,
is there existing code to fit a ordinal regression model with probit link where
the variance can be modeled? For example, the polr function from the MASS
package fits:
probitP(Y <= k | x) = zeta_k - eta
Instead of,
probitP(Y <= k | x) = (zeta_k - eta)/scale, where scale = exp(b*z
67 matches
Mail list logo