On Jan 24, 2023, at 5:00 AM, r-help-requ...@r-project.org wrote:
Send R-help mailing list submissions to r-help@ r-project. org To subscribe or
unsubscribe via the World Wide Web, visit https: //urldefense. proofpoint.
com/v2/url?u=https-3A__stat. ethz.
ch_mailman_listinfo_r-2Dhelp&d=DwIGaQ
s from the deviance
residuals, but I am not sure this is correct. Also I wonder if there is a
better way doing this. Can anyone help?
y <- -(z$observed_target)^(2-p)
y <- y + (fitted(out))^(2-p)
y <- y / (2-p)
x <- (z$observed_target)^(2-p)
x <- x - z$observed_target * ((fitted(ou
Dear mailing list members,
My question is maybe very basic, but I could not find the solution.
I would like to do the following things
1)
colnames(V1)[2] <- par$V2[1]
colnames(V2)[2] <- par$V2[2]
colnames(V3)[2] <- par$V2[3]
...
colnames(V37)[2] <- par$V2[37]
2)
V1 <- V1[,-1]
V2 <- V2[,-1]
V3 <
Dear mailing list members,
Sorry for the repeated posting, but I have to add two more things.
1) I changed the .Rdata into .rdata in the subject to avoid
misunderstanding of readers.
2) To save data in csv file;
> write.csv(WV6_Data_R, "input/ch2/WV6_Data_R_v_2016_01_01.csv")
Best regards,
Yoh
Dear mailing list members,
A member found a problem of lines of codes in my last e-mail.
I submit the correction to this mailing list as follows.
In the final part of my last e-mail message, a csv file is created.
But, I cannot open this file, so please ignore the following part.
Dear mailing list members,
Thank you for all members who gave me advice.
This question is solved.
I tried to create an object using the load function by analogy based on the
read.csv function like this;
> d <- load("~/docdis/input/ch2/WV6_Data_R_v_2016_01_01.rdata")
but, load function cannot be
Dear mailing list members,
I am a beginner of this community.
I would like to analyze data in a .rdata file.
I ran the following code, but the object "d" remained empty.
d <- load("~/docdis/input/ch2/WV6_Data_R_v_2016_01_01.rdata")
I would like to know how to convert this file format into .csv
I'm trying to use the following command.
arima (x, order = c(p,d,q), seasonal =list(order=c(P,D,Q), period=s)
How can I write an estimated seasonal ARIMA model from the outputs. To be
specifically, which sign to use? I know R uses a different signs from S plus.
Is it correct that the model is:
(
Is anyone aware of point.in.polygon giving different results for 32-bit vs.
64-bit R? Our OS is 64-bit Windows 7 Enterprise. I'm working with someone
else's extensive R program and the final results are close but not exactly
matching. We're thinking it might be something with the point.in.polygo
t is much harder.
> You'll notice that everyone does CDF curves for survival data (
> Kaplan-Meier = estimate of 1-CDF), but no one does histograms, which
> estimate a density. That isn't because we wouldn't like density estimates.
>
> Terry T.
>
>
>
> On 11/04/
Hi Dr. Therneau,
Thanks very much for your kind help! Does survfit() just give me the
survival curve? What I wanted is the baseline hazard estimates, i.e.,
lambda_{0} (t). How can I obtain this estimate from coxph()? Or using
basehaz()?
Thanks,
YH
On Mon, Nov 4, 2013 at 5:29 PM, Terry T
Hello,
Does anyone know how to obtain nonparametric baseline hazard estimates in
the gamma frailty model by using coxph()? I can easily get the other
parameter estimates by using coxph() but don't know how to get the baseline
hazard estimates from it.
Thanks,
YH
[[alternative HTML versio
I can easily get the other parameter estimates by using coxph() but don't
know how to get the baseline hazard estimates from it. Does anyone know how
to obtain nonparametric baseline hazard estimates in the gamma frailty
model?
Thanks,
YH
[[alternative HTML version deleted]]
Thanks very much for your help, Terry and Göran!
As pointed out by Göran, the difficult part is that it's an open set. How
to obtain a valid MLE in this case?
Thanks,
YH
On Wed, Oct 16, 2013 at 9:55 AM, Göran Broström wrote:
>
>
> On 2013-10-16 14:33, Terry Therneau wrote:
>
>>
>>
>> On
Hello,
I'm trying to use coxph() function to fit a very simple Cox proportional
hazards regression model (only one covariate) but the parameter space is
restricted to an open set (0, 1). Does anyone know how to obtain an
estimate in a particular parameter space by using coxph function? Or any
othe
Hello,
I'm trying to use coxph() function to fit a very simple Cox proportional
hazards regression model (only one covariate) but the parameter space is
restricted to an open set (0, 1). Can I still obtain a valid estimate by
using coxph function in this scenario? If yes, how? Any suggestion would
r tout renseignement.
Y. Ledroit
Charg de clientle
[8]Plateforme de Transfert de fichier
Si ce message ne s'affiche pas correctement, [9]visualisez la version en
ligne.
Si vous ne s
Hi,
Is there an easy way to add symbols to a line in a scatter plot, so that
only a few symbols are added per line (a line drawn based on a large set of
data points) , in order to distinguish several lines in one graph (not a
symbol for each data point as is the default in plot()).
Thanks a lot!
2011 at 11:56 PM, Gabor Grothendieck <
ggrothendi...@gmail.com> wrote:
> On Sun, Mar 13, 2011 at 11:32 PM, Paul Y. Peng wrote:
> > Many thanks to Jim Holtman and Gabor Grothendieck for your quick
> responses.
> > Jim's solution works beautifully for my tasks. Thanks. I
11, 2011 at 4:53 PM, jim holtman wrote:
> > You can do:
> >
> > source(readLines(yourFile)[10:20]) # lines 10-20 of the file
> >
> > On Fri, Mar 11, 2011 at 4:41 PM, Paul Y. Peng wrote:
> >> I have a text file of R commands. Some times I only want to run a
I have a text file of R commands. Some times I only want to run a few lines
of the R commands in an existing R session and wonder whether there is a
simple way to do this.
To run a few lines in a new session of R, I could use sed to pick up the
lines from the file and pipe them into R.
source() d
data.frame(group, dose, ldose, n, r, p)
SF <- cbind(success=d$r, failure=d$n - d$r)
#80 record set-up
dose2<-c(7,8,9,10,11)
doserep<-sort(rep(dose2,8))
x<-c(doserep,doserep)
log10x<-log10(x)
y_U<-c(rep(0,8), 1, rep(0, 7), 1, 1, 1, rep(0,5), rep(1, 16))
y_M<-c(rep(0,24), rep(1,4), r
data.frame(group, dose, ldose, n, r, p)
SF <- cbind(success=d$r, failure=d$n - d$r)
#80 record set-up
dose2<-c(7,8,9,10,11)
doserep<-sort(rep(dose2,8))
x<-c(doserep,doserep)
log10x<-log10(x)
y_U<-c(rep(0,8), 1, rep(0, 7), 1, 1, 1, rep(0,5), rep(1, 16))
y_M<-c(rep(0,24), rep(1,4), r
do this. Sorry for the
basic questions.
On Sat, Apr 17, 2010 at 2:49 PM, David Winsemius wrote:
>
> On Apr 17, 2010, at 1:18 PM, David Winsemius wrote:
>
>>
>> On Apr 17, 2010, at 12:45 PM, Charles C. Berry wrote:
>>
>>> On Sat, 17 Apr 2010, Akito Y. Kawahara
Hi, I am new to R, and have a quick question regarding an R script
that I received from a kind colleague.
I am trying to determine the "peak" (maximum value) of the graph that
is plotted when executing the following. There is an input file called
"rates_values.txt" which begins as:
rateValue
0.37
Hi all,
I'm using windowsXP and R 2.10.0. I downloaded "randomForest
4.5-33.tar.gz" from its archive, but how do I make it into a
installation ZIP file?
Thanks,
Ted
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
P
Hi, I would plot a Henon map:
> plot(henon(start=rnorm(2), a=1.4, b=0.3, n.sample=2000,
> n.transient=10),pch=".",cex=0.1)
but give me some errors:
Error in henon(start = rnorm(2), a = 1.4, b = 0.3, n.sample = 2000, n.transient
= 10) :
could not find function "checkVectorType"
Error in plo
Hi,
I am interested to use Kernel Density Estimation on bivariate data for
multi-class classification.
So far, I have managed to use the 'ks' package to plot the contours of the
kernel density estimates based on 8-class training dataset with only 2
variables.
However, I do not know how to m
Hi,
I am interested to use Kernel Density Estimation on bivariate data for
multi-class classification.
So far, I have managed to use the 'ks' package to plot the contours of the
kernel density estimates based on 8-class training dataset with only 2
variables.
However, I do not know how
Hi, all,
My objective is to split a dataframe named "cmbine" according to the value of
"classes". After the split, I will take the first instance from each class and
bin them into a new dataframe, "df1". In the 2nd iteration, I will take the 2nd
available instance and bin them into another new
Greetings,
When reading the random forest manual by Liaw, in the examples
"set.seed" is always used before building a forest.
Does it matter if I don't set the seed?
If I set a different seed manually each time I build a forest, will this
give better randomness?
Thank you all in advance!
--
We are starting the process for purchasing a server with Windows Server 2008
OS. We are mainly interested in storage for multiple users, but would like to
have the capability to run R simulations on the server. Are there any issues we
should be aware of for installing R on the server or with hav
Dear r-help,
How can I add a strip to show group weights using lattice package? For
example, in the following code, I'd like to using "wt" variable in a trip
to demonstrate the relative size of groups.
(Following is just the simplest form to demonstrate the question. A
conditional variable will
On Wed, Dec 17, 2008 at 3:06 AM, Martyn Plummer wrote:
> The artefacts that you see are a normal result of using bitmap graphics
> devices. I have tried to explain these below:
Thanks very much for your explanations, MP; they were quite informative!!
I recognize that others may feel differently
On Wed, Dec 17, 2008 at 12:54 AM, Prof Brian Ripley
wrote:
> Your PDF problems indicate a broken viewer. How were you viewing PDF?
You are absolutely correct about the PDF files; I've since checked the
PDF files in other viewers and have not been able to reproduce the
problem. There was definite
On my current home system, I am getting undesirable output from
graphical devices such as png() and pdf(). The graphical output is
blurry. I haven't experienced the problem on other systems. As you
will see from the attached text file (more information on this file
below), the problem does not occu
Hi everyone, when I use the two sample Kolmogorov¨CSmirnov ks2Test like this:
x=read.table("e:/x.txt")
y=rstable(1000,alpha,beta,gamma,delta)
I alway get results as follows:
Warning messages:
1: In ks.test(x = x, y = y, alternative = "two.sided") :
cannot compute correct
37 matches
Mail list logo
