ahimsa campos arceiz wrote:
> Sorry I have a very simple question:
>
> I used somers2 function from Design package:
>
>
>>z<- somers2(x,y, weights=w)
>
>
> results are:
>
>
>>z
>
> CDxynMissing
> 0.88 0.76 5000.00
>
> Now I want to call only the value of C to be used
Sorry I have a very simple question:
I used somers2 function from Design package:
> z<- somers2(x,y, weights=w)
results are:
>z
CDxynMissing
0.88 0.76 5000.00
Now I want to call only the value of C to be used in further analyses, but I
fail to do it. I have tried:
> z$C
N
On Wed, 2006-01-04 at 13:28 +0800, Vincent Deng wrote:
> Dear R-helpers
>
> I have a stupid question about cbind function. Suppose I have a
> dataframe like this
> Frame:
> A 10
> C 20
> B 40
>
> and a numeric matrix like this
> Matrix:
> A 1
> B 2
> C 3
>
> cbind(Frame[,2],Matrix[,1]) simply bi
Dear R-helpers
I have a stupid question about cbind function. Suppose I have a
dataframe like this
Frame:
A 10
C 20
B 40
and a numeric matrix like this
Matrix:
A 1
B 2
C 3
cbind(Frame[,2],Matrix[,1]) simply binds these two columns without
checking the order, I mean, the result will be
A 10 1
B 2
On Tue, 2006-01-03 at 22:56 -0500, Mark Leeds wrote:
> I don't want to bother anyone
> with specific questions
> because I am a R newbie and
> I see that there is TON ( emphasis on
> Ton ) of documentation
> out there but could
> someone just tell me the
> best placed to look/read
> for learning a
I replicated your 'false convergence' using R 2.2.0:
> sessionInfo()
R version 2.2.0, 2005-10-06, i386-pc-mingw32
attached base packages:
[1] "methods" "stats" "graphics" "grDevices" "utils" "datasets"
[7] "base"
other attached packages:
nlme MASS
"3.1-66" "7.2-23"
On 3 January 2006 at 22:56, Mark Leeds wrote:
| out there but could
| someone just tell me the
| best placed to look/read
| for learning about ( for R-2-2.1 in Windows )
|
| .Rprofile
| .REnviron.
| .Rdata
?Startup
| .First function ( analogous to the one in Splus ).
?.First
| Analog of Sp
I don't want to bother anyone
with specific questions
because I am a R newbie and
I see that there is TON ( emphasis on
Ton ) of documentation
out there but could
someone just tell me the
best placed to look/read
for learning about ( for R-2-2.1 in Windows )
.Rprofile
.REnviron.
.Rdata
.First f
Hi All,
Sorry if this is a repost (a quick browse didn't give me the answer).
I wonder if there are packages that can do the feature selection and
classification at the same time. For instance, I am using SVM to classify my
samples, but it's easy to get overfitted if using all of the features. Th
Gabor Grothendieck wrote on 1/3/2006 2:37 PM:
>Looking at the first few queries,
>see how easy it is to take the top few in Stata whereas in R one would
>have a complex use of order. Its not hard in R to write a function
>that would make it just as easy but its not available off the top
>of one'
I'm not certain what you are asking. I just got 10 hits for
'RSiteSearch("Glimmix")'. Seven of them mentioned SAS PROC GLIMMIX:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/65945.html
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/65954.html
http://finzi.psych.upenn.edu/R/Rhelp02a
Jason
the lmer() function in the Matrix package is what you will need.
-Original Message-
From: [EMAIL PROTECTED] on behalf of Jason Marshal
Sent: Tue 1/3/2006 8:10 AM
To: r-help@stat.math.ethz.ch
Cc:
Subject:[R] Including random effects in logistic regression.
I'm
>Hello,
>Unlike most posts on the R mailing list I feel qualified to comment on
>this one. For about 3 months I have been trying to learn use R, after
>having used various versions of SPSS for about 10 years.
I think it is far too simplistic to ascribe non-use of R to laziness. This
ma
On Tue, 2006-01-03 at 18:57 -0600, Eberhard F Morgenroth wrote:
> I have a list as follows
>
> P <- list(A = c("CS", "CX"),
> B = 1:4,
> Y = c(4, 9))
>
> I now would like to prepare a new list where the rows of the new list
> provide all possible combinations of the element
I have a list as follows
P <- list(A = c("CS", "CX"),
B = 1:4,
Y = c(4, 9))
I now would like to prepare a new list where the rows of the new list
provide all possible combinations of the elements in the orginal list.
Thus, the result should be the following
CS 1
Hi,
I'm looking for ideas or packages with relevant algorithms for
calculating the connectivity across a grid, where connectivity is
defined as the minimum amount of cross-sectional area along a continuous
path. The upper boundary of the cross-sectional area is a fixed
elevation, and the lower bo
I'm working with a scatterplot of data, using plot() with log="xy" to get
log-log axes. But I can't get the regression line to plot correctly. I use
abline(lm(log(Y)~log(X))) and ge a line that looks like the correct slope,
but the Y-intercept is messed up. I haven't changed the y-axis other tha
I'm trying to analyse some data using logistic regression in R, but I
want to include random effects in the model. The glm function
appears not to have options for including random effects, and the lme
and nlme documentation indicates that these functions are for
continuous, not dichotomous, r
On Tue, 3 Jan 2006, Matt Williamson wrote:
> Hello List,
> I was wondering if a package or piece of code exists that will allow all
> possible subsets regression model selection within program R. I have
> already looked at step(AIC) which does not test differing combinations
> of variables within
Hello List,
I was wondering if a package or piece of code exists that will allow all
possible subsets regression model selection within program R. I have
already looked at step(AIC) which does not test differing combinations
of variables within a model as far as I can tell. In addition I tried
to
Have you tried nlme?
I tried to something similar.
Here is the code that I used for a negative binomial random effects model
library(nlme)
mydata<-read.table("C:\\Plx\\plx.all\\plxall.txt",header=TRUE)
loglike =
function(PLX_NRX,PD4_42D,GAT_34D,VIS_42D,MSL_42D,SPE_ROL,XM2
On Tue, 3 Jan 2006, Verena Hoffmann wrote:
> Hello!
>
> I want to compare two Kaplan-Meier-Curves by using the Logrank-Test:
>
> logrank(Surv(time[b], status[b]) ~ group[b])
>
> This way I only get the value of the test-statistic, but not the p-value.
>
> Does anybody know how I can get the p-valu
Thanks to all. I didn't
Realize that you
Got so many packages
Automatically.
I've used S+ for roughly
10 years on and off and
I am starting to switch over
( was finally forced to because my new company preferred
me to use R for cost. I am the only user )
and it's unbelievable what has been done
On Tue, 3 Jan 2006, Mark Leeds wrote:
> I'm sorry to bother everyone with a stupid
> question but, when I am at an R prompt in Windows,
> is there a way to see what packages
> you already have installed from the R site so
> that you can just do library(name_of_package)
> and it will work.
>
> I've
On Tue, 3 Jan 2006, Mark Leeds wrote:
> I'm sorry to bother everyone with a stupid
> question but, when I am at an R prompt in Windows,
> is there a way to see what packages
> you already have installed from the R site so
> that you can just do library(name_of_package)
> and it will work.
>
> I'
On Tue, 2006-01-03 at 16:07 -0500, Mark Leeds wrote:
> I'm sorry to bother everyone with a stupid
> question but, when I am at an R prompt in Windows,
> is there a way to see what packages
> you already have installed from the R site so
> that you can just do library(name_of_package)
> and it will
"Mark Leeds" <[EMAIL PROTECTED]> writes:
> I'm sorry to bother everyone with a stupid
> question but, when I am at an R prompt in Windows,
> is there a way to see what packages
> you already have installed from the R site so
> that you can just do library(name_of_package)
> and it will work.
>
>
Hello!
I want to compare two Kaplan-Meier-Curves by using the Logrank-Test:
logrank(Surv(time[b], status[b]) ~ group[b])
This way I only get the value of the test-statistic, but not the p-value.
Does anybody know how I can get the p-value?
Thanks in advance!
Verena Hoffmann
_
I'm sorry to bother everyone with a stupid
question but, when I am at an R prompt in Windows,
is there a way to see what packages
you already have installed from the R site so
that you can just do library(name_of_package)
and it will work.
I've looked at help etc but I can't find
a command like t
From: Stephen Choularton
>
> Hi
>
> I am trying to use randomForest for classification. I am using this
> code:
>
> > set.seed(71)
> > rf.model <- randomForest(similarity ~ ., data=set1[1:100,],
> importance=TRUE, proximity=TRUE)
> Warning message:
> The response has five or fewer unique valu
Hi
I am trying to use randomForest for classification. I am using this
code:
> set.seed(71)
> rf.model <- randomForest(similarity ~ ., data=set1[1:100,],
importance=TRUE, proximity=TRUE)
Warning message:
The response has five or fewer unique values. Are you sure you want to
do regression? in:
On 1/3/06, Thomas Lumley <[EMAIL PROTECTED]> wrote:
> On Tue, 3 Jan 2006, Peter Dalgaard wrote:
> > One thing that is often overlooked, and hasn't yet been mentioned in
> > the thread, is how much *simpler* R can be for certain completely
> > basic tasks of practical or pedagogical relevance: Calcu
As others have pointed out, since R is more of a programming language than a
statistical package, yes, it is _harder_ to learn. I would say its easier
to learn than C++, harder to learn than VBA, and on par with learning Java,
but that's all debatable.
One thing that makes R slightly more intimid
Accepting this stacked representation for the
moment try this. When reordering the dates do it
in reverse order. Then loop over all codes
applying the zoo function na.locf to the the
prices for that code. locf stands for last
observation carried forward. Since our dates
are reversed it will bri
A candidate for the fortunes package?
(Perhaps the highest honor one can receive: being "verbified" :-) )
> And the fear of
> getting Ripleyed on the mailing list also makes me think, read, and
> improve before submitting half baked questions to the list.
>
> -- Eric Kort
Cheers,
Be
On Tue, 3 Jan 2006, Peter Dalgaard wrote:
> One thing that is often overlooked, and hasn't yet been mentioned in
> the thread, is how much *simpler* R can be for certain completely
> basic tasks of practical or pedagogical relevance: Calculate a simple
> derived statistic, confidence intervals from
Your 2-million loop is overkill, because apparently in the (vast)
majority of cases you don't need to loop at all. You could try
something like this:
1. Split the price by id, e.g.
price.list <- split(price,id)
For each id,
2a. When price is not NA, assign it to next price _without_ using a
for loo
Berton Gunter writes
> U
>
> I cannot say how easy or hard R is to learn, but in response to the
UCLA
> commentary:
>
> > However, I
> > feel like R
> > is not so much of a statistical package as much as it is a
statistical
> > programming environment that has many new and cutting e
Wensui Liu wrote:
>Another big difference between R and other computing language such as
>SPSS/SAS/STATA.
>You can easily get a job using SPSS/SAS/STATA. But it is extremely difficult
>to find a job using R. ^_^.
>
>
Actually in finance it is getting easier all the time for
knowledge of R to be
U
I cannot say how easy or hard R is to learn, but in response to the UCLA
commentary:
> However, I
> feel like R
> is not so much of a statistical package as much as it is a statistical
> programming environment that has many new and cutting edge
> features.
Please note: the first se
>>> "Ben Fairbank" <[EMAIL PROTECTED]> 1/3/2006 12:42 pm >>> wrote
<<<
One implicit point in Kjetil's message is the difficulty of learning
enough of R to make its use a natural and desired "first choice
alternative," which I see as the point at which real progress and
learning commence with any ne
One implicit point in Kjetil's message is the difficulty of learning
enough of R to make its use a natural and desired "first choice
alternative," which I see as the point at which real progress and
learning commence with any new language. I agree that the long learning
curve is a serious problem,
Another big difference between R and other computing language such as
SPSS/SAS/STATA.
You can easily get a job using SPSS/SAS/STATA. But it is extremely difficult
to find a job using R. ^_^.
On 03 Jan 2006 17:53:40 +0100, Peter Dalgaard <[EMAIL PROTECTED]>
wrote:
>
> Patrick Burns <[EMAIL PROTECTE
Try this:
Pstats <- function(x) c(Max = max(x),
Min = min(x),
AMean = mean(x),
AStdev = sd(x),
Samples = length(x),
quantile(x, 1:9/10, na.rm = TRUE))
res <- with(areas, by(AdRes, N_Type, Pstats))
do.call("rbind", res)
Also, check out summaryBy in the doBy package at
http://geneti
Looks like I may have found a function that addresses my needs. Bootcov in
Design handles bootstrapping from clustered data and will save the
coefficients. I'm not entirely sure it handles clusters the way I'd like,
but I'm going through the code. If it doesn't, it looks easily
re-writeable. As
This most often indicates a problem with dispersion estimate. See the
cautionary tale in MASS4 chapter 7. If you have a reliable dispersion
estimate that low for genuine counts, they are either not independent or
not Poisson (for example, limited), and one would want to find out what is
going
I am running R 2.1.1 in a Microsoft Windows XP environment.
I have a matrix with three vectors (columns) and ~2 million rows. The
three vectors are date_, id, and price. The data is ordered (sorted) by code
and date_.
(The matrix contains daily prices for several thousand stocks
Patrick Burns <[EMAIL PROTECTED]> writes:
> I have had an email conversation with the author of the
> technical report from which the quote was taken. I am
> formulating a comment to the report that will be posted
> with the technical report.
>
> I would be pleased if this thread continued, so I
On 3 Jan 2006 at 7:35, Thomas Lumley wrote:
> On Mon, 2 Jan 2006, Philippe Grosjean wrote:
> >
> > That said, I think one should interpret Mitchell's paper in a different
> > way. Obviously, he is an unconditional and happy Stata user (he even
> > wrote a book about graphs programming in Stata). H
I have written a few different summary functions. I want to calculate
the statistics by groups and I am having trouble getting the output as a
dataframe. I have attached one example with a small dataset that
calculates summary stats and percentiles, I have others that calculate
upper confidence lim
On 1/3/2006 10:30 AM, begert wrote:
> Dear R- Users,
>
>is there a way to determine the size of
>an rgl window (rgl.open()) either in advance or
>afterwards, (without using the mouse, of course) ?
>
>Intuitively, one would assume to set the size by:
>
>library("rgl");
>pa
I'm a french medicine student and i work on oncology. I work about
treatment oh breast cancer. I have 3 sub group of patient. I made some
kaplan meyer survival curve, and i made a cox model.
On the survival curve, on the last observations there is a crossmatch of
the different survival curve at 15
>>> "John Fox" <[EMAIL PROTECTED]> 1/3/2006 9:35 am >>> as always,
raises some excellent points. I have some responses, interspersed
It's not reasonable to argue with someone's experience -- that is, if
people
tell me that they found R harder to learn than SAS, say, then I believe
them
-- bu
Dear Peter et al.,
It's not reasonable to argue with someone's experience -- that is, if people
tell me that they found R harder to learn than SAS, say, then I believe them
-- but that's not my experience in teaching relatively inexperienced
students to use statistical software. A few points:
(1)
"Brian Perron" <[EMAIL PROTECTED]> writes:
> Hello all:
>
> I am interested in computing what the multilevel modeling literature
> calls a multiple membership model. More specifically, I am working
> with a data set involving clients and providers. The clients are the
> lower-level units who a
I am trying to use Poisson regression to model count data. My results are
suggestive of under dispersion (0.79). How close to one does one want the
measure of dispersion to be before one accepts the results of the analysis?
I know that there is no definitive answer to my question, but I would li
On Tue, 3 Jan 2006, Brian Perron wrote:
> Hello all:
>
> I am interested in computing what the multilevel modeling literature
> calls a multiple membership model. More specifically, I am working with
> a data set involving clients and providers. The clients are the
> lower-level units who ar
Souds like a model with cross-classified random effects. Lme4 can handle
this easily.
Shige
On 1/3/06, Brian Perron <[EMAIL PROTECTED]> wrote:
>
> Hello all:
>
> I am interested in computing what the multilevel modeling literature calls
> a multiple membership model. More specifically, I am work
In fact it's just as easy in Internet Explorer: right-click + Open in
New Window, or Shift-Click, followed by Ctrl+D. Or, right-click + Add
to Favorites.
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of
> Charles Annis, P.E.
> Sent: Monday, January 0
On Mon, 2 Jan 2006, Philippe Grosjean wrote:
>
> That said, I think one should interpret Mitchell's paper in a different
> way. Obviously, he is an unconditional and happy Stata user (he even
> wrote a book about graphs programming in Stata). His claim in favor of
> Stata (versus SAS and SPSS, and
I have had an email conversation with the author of the
technical report from which the quote was taken. I am
formulating a comment to the report that will be posted
with the technical report.
I would be pleased if this thread continued, so I will know
better what I want to say. Plus I should be
Dear R- Users,
is there a way to determine the size of
an rgl window (rgl.open()) either in advance or
afterwards, (without using the mouse, of course) ?
Intuitively, one would assume to set the size by:
library("rgl");
par3d(viewport=c(0,0,500,500));
#rgl.open();
for ex
On 1/3/06 10:24 AM, "Constantine Tsardounis" <[EMAIL PROTECTED]>
wrote:
> Hello, Happy New Year!...
>
> I am encountering a problem trying to work on the data that I load in R.
>
> I have loaded to R a series of stock data using
> (csv files are named e.g. IBM.R)
>
> length.R <- length(list.
Hello, Happy New Year!...
I am encountering a problem trying to work on the data that I load in R.
I have loaded to R a series of stock data using
(csv files are named e.g. IBM.R)
length.R <- length(list.files(".", pattern=".R")) # the number of
files with one #column in the director
Is this happening with the example as you show it, or are you trying to
print mydata?
There is a bug in the print method for TSdata objects, which I have
fixed and was intending to put on CRAN in a few days. This bug does give
the infinite recursion error, but would only happen when you print t
Hello all:
I am interested in computing what the multilevel modeling literature calls a
multiple membership model. More specifically, I am working with a data set
involving clients and providers. The clients are the lower-level units who are
nested within providers (higher-level). However
If I use cex.lab = 2 and cex.axis = 2 the yaxis label in a plot exceed
the plot area. How do I get the plot itself smaller to get space for the
label so I still use cex.lab = 2 and cex.axis = 2?
Kare
--
###
Kare Edvardsen <[EMAIL PROTECTED]>
Norwegian In
Clark,
I agree with Vitor that working in R might be easier, but it seems that you
are working in the excel VBA environment and there may be good reason why
you want to do so that we don't know about. If so, why use Rexcel function
to read the file into excel when you can use VBA code to open the
> On Tue, 3 Jan 2006 07:29:27 +,
> hadley wickham (hw) wrote:
>> A solution would be a content-management system that produced the HTML
>> of the site from some other form of input. Only the output HTML would
>> need to be mirrored. Care to put together such a thing, and import
Check the way you imported the data / the SQLite documentation. The
\r\n that you see (you're on Windows, right?) is used to indicate the
end of the data lines in the source file - \r is a carriage return,
and \n is a new line character.
> -Original Message-
> From: [EMAIL PROTECTED]
> [m
Hello Liu,
this might be caused by NA entries in your SQLite table. Have a look at the
following code:
(test <- data.frame(matrix(c(1:10, NA, NA), ncol=2, nrow=6)))
con <- dbConnect(SQLite(), dbname = "test.db")
dbWriteTable(con, "test", test, type="BLOB", overwrite=TRUE)
d1 <- dbReadTable(con,
On 1/3/06 6:46 AM, "Mohammed Asifulla - CTD , Chennai" <[EMAIL PROTECTED]>
wrote:
> Hi,
>
> I am new-comer to statistics and R-Project. I would like to know if these
> features can be attained in R-Project.Please help.
>
> 1) beta 1 and Beta 2, or gamma one and gamma two for skewness and kur
On 1/1/06 2:07 PM, "Erin Hodgess" <[EMAIL PROTECTED]> wrote:
> Dear R People:
>
> Could someone direct me to some documentation on the
> difference between S3 and S4 classes, please?
>
> For example, why would a person use one as opposed to another?
> Maybe pros and cons of each?
The Biocon
Hi,
I am new-comer to statistics and R-Project. I would like to know if these
features can be attained in R-Project.Please help.
1) beta 1 and Beta 2, or gamma one and gamma two for skewness and kurtosis,
respectively, including standard errors and tests for significance (relative
to values for
On 1/3/06, Peter Flom <[EMAIL PROTECTED]> wrote:
> >>> "Rau, Roland" <[EMAIL PROTECTED]> >>> wrote
> <<<
> IMO this is a very good proposal but I think that the main problem is
> not the "translation" of one function in SPSS/Stata/SAS to the
> equivalent in R.
> Remembering my first contact with R
the help page on "setOldClass" might help you. In particular the section
"Register or Convert?".
Matthias
hadley wickham schrieb:
>I'm trying to create an extension to data.frame with more complex row
>and column names, and have run into some problems:
>
>
>
>>setClass("new-data.frame", repre
>>> "Rau, Roland" <[EMAIL PROTECTED]> >>> wrote
<<<
IMO this is a very good proposal but I think that the main problem is
not the "translation" of one function in SPSS/Stata/SAS to the
equivalent in R.
Remembering my first contact with R after using SPSS for some years (and
having some experience
Dear All,
I have the following error message when I fitted lmer to a binary data with
the "AGQ" option:
Error in family$mu.eta(eta) : NAs are not allowed in subscripted assignments
In addition: Warning message:
IRLS iterations for PQL did not converge
Any help?
Thanks in advance,
Abderrahi
Roland,
Yes, indeed, you are perfectly right. The problem is that R richness
means R complexity: many different data types, "sub-languages" like
regexp or the formula interface, S3/S4 objects, classical versus lattice
(versus RGL versus iplots) graphs, etc. During translation of R in
French, I
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Gabor
> Grothendieck
> Sent: Monday, January 02, 2006 4:59 PM
> To: Philippe Grosjean
> Cc: Kort, Eric; Kjetil Halvorsen; R-help@stat.math.ethz.ch
> Subject: Re: [R] A comment about R:
>
>
> Probably
I'm trying to create an extension to data.frame with more complex row
and column names, and have run into some problems:
> setClass("new-data.frame", representation("data.frame"))
[1] "new-data.frame"
Warning message:
old-style ('S3') class "data.frame" supplied as a superclass of
"new-data.frame"
On 12/30/05, Âíµ¤ <[EMAIL PROTECTED]> wrote:
>
> Hello,
>I am a new R user and I need R code for GLARMA.
> I will be really thankful if you help me.
You should really spell out an acronym like GLARMA.
RSiteSearch("GLARMA") does't give anything. ou could have a look at
package sspir.
Kjetil
This has come up before: it needs a bug fix which Paul Gilbert has already
implemented (but not yet released).
Please use an informative subject line, and don't SHOUT at us. (All caps
is regarded as shouting, and BTW the package bundle is dse not DSE.)
On Tue, 3 Jan 2006, Sumanta Basak wrote:
Hello,
I am a new R user and I need R code for GLARMA.
I will be really thankful if you help me.
Yours sincerely,
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-p
84 matches
Mail list logo