On Thu, 21 May 2009, Paolo Radaelli wrote:
Dear R-users,
I'm working with the Rpart package and trying to understand how the procedure
I presume you mean 'rpart': R package names are case sensitive.
select the best split in the case the method "class" and the
splitting index "Gini" are used
Dear all,
The "tractor.base" package has recently been added to CRAN. This
package provides functions to read, write, visualise and manipulate
magnetic resonance images. The standard Analyze, NIfTI and DICOM file
formats are supported (read-only for DICOM), and all metadata is
stored with the imag
Simon Blomberg wrote:
How about to insert a variable a2 inbetween the first and second columns
of dat:
dat2 <- cbind(dat[,1], a2=a2, dat[,2:3])
Where a2 is the new variable. This mangles the variable name for column
1, unfortunately. Surely someone else will offer a better solution.
Simon.
Hi,
I would like to know if it is possible to have a "R code" to estimate the
parameters of a mixture of bivariate (or multivariate) normals via EM
Algorithm. I tried to write it, but in the estimation of the matrix of
variance and covariance, i have some problems. I generate two bidimensional
Hello everybody.
I am new to R. Yesterday I passed the afternoon reading the
introduction and language reference, but I could'nt find a way to do a
3d plot of the density of a data table of size 2.
I am trying with:
plot(density(t(t2)))
but it mixes the two columns and calculate the density like
Hi,
Jose M. Blanco-Moreno an myself have implemented
Syrjala's test in the ecespa package. As a matter
of fact, Jose M. has implemented a Frotran
version of Syrjala's original QBasic function.
Using your data the results are very close to
your reported # statistic= 0.224 and # p-value
Hi Friends,
I am trying to import a table from SQL server to R(2.9.0), however i am
getting errors while running the below codes. Can anyone identify and let me
know where did i go wrong??? Thanks in anticipation :)
library(RODBC)
myconn <- odbcConnect("RDATABASE")
myconn
RODB Connection 6
D
>I am trying to parse XML file ( binary hex) but get an error.
> Code I am using is:
> xsd = xmlTreeParse(system.file("exampleData", "norel.xsd", package =
> "XML"), isSchema =TRUE) doc = xmlInternalTreeParse(system.
> file("exampleData", "LogCallSummary.bin", package = "XML")) Start
> tag ex
I use the following function to export some figures to .eps:
p.eps <- function(p, fname, title = NULL, width, height)
{
postscript(file=fname, onefile=FALSE, paper="special",
width=width, height=height, horizontal=FALSE)
print(p + opts(title = title))
dev.off()
}
Whenever I hav
Dear Fernando and all:
Thanks for your help. Now works. This is
a training example to learn how to estimate a
Box-Cox (right and/or left side transformations)
with R (as LIMDEP does) in order to compare these
estimations with the ones derived by applying
NLS, ones the dependent var
Hi Guys,
i think this is a relatively simple question.
I have coded the following polynomial function
y= function(x) x^3-2*x^2+1
I need to find the inverse of this, but the code i am using now isnt returning
what i want it to.
What is the general code for finding an inverse function?
Thanks
Look up packages flexmix and mclust!
Christian
On Thu, 21 May 2009, daniele riggi wrote:
Hi,
I would like to know if it is possible to have a "R code" to estimate the
parameters of a mixture of bivariate (or multivariate) normals via EM
Algorithm. I tried to write it, but in the estimation
Apparently the way to deal with this error message is to set
R_CStackLimit = (uintptr_t)-1
I tried typing this in the R console, but it says Error: object
"R_CStackLimit" not found.
So where do I type it? In one of the initialization files that R uses when
it starts up?
I can't find the answer a
Hi,
I am trying to calculate ci from lmer. However, I have the error message
below:
library(gmodels)
library(lme4)
fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy)
ci(fm2)
Error in as.vector(x, mode) :
cannot coerce type 'S4' to vector of type 'any'
In addi
> I am new to R. Yesterday I passed the afternoon reading the
> introduction and language reference, but I could'nt find a way to do a
> 3d plot of the density of a data table of size 2.
> I am trying with:
>
> plot(density(t(t2)))
>
> but it mixes the two columns and calculate the density like
On Thu, May 21, 2009 at 08:36:27PM +1000, Kon Knafelman wrote:
>
> I have coded the following polynomial function
>
> y= function(x) x^3-2*x^2+1
>
> I need to find the inverse of this, but the code i am using now isnt
> returning what i want it to.
>
This function is not injective, so the inve
On Thu, 21 May 2009, Zeljko Vrba wrote:
I use the following function to export some figures to .eps:
p.eps <- function(p, fname, title = NULL, width, height)
{
postscript(file=fname, onefile=FALSE, paper="special",
width=width, height=height, horizontal=FALSE)
print(p + opts(title
Dear R people
I ask again
1. Is there a published reference presenting the normal score
transformation? Is there a published paper (in any field) using that
transformation in the analysis of data?
And this is a new question
2. The outliers library has a function called rm.outlier
On Thu, May 21, 2009 at 12:32:28PM +0100, Prof Brian Ripley wrote:
>
> Using Acrobat Reader to view PostScript! It is a PDF viewer.
>
Ah, sorry, I explicitly convert the PS with ghostscript's ps2pdf.
>
> suspect you need to track down where conversion to PDF is happening
> and disable auto-ro
Fernando Marmolejo Ramos wrote:
Dear R people
I ask again…
1. Is there a published reference presenting the normal score
transformation? Is there a published paper (in any field) using that
transformation in the analysis of data?
Dunno, but I'd search for "Normal Scores Test", AKA van de
Is it possible to have offsets in a glm when using bigglm?
_
[[elided Hotmail spam]]
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.c
Zeljko Vrba wrote:
On Thu, May 21, 2009 at 12:32:28PM +0100, Prof Brian Ripley wrote:
Using Acrobat Reader to view PostScript! It is a PDF viewer.
Ah, sorry, I explicitly convert the PS with ghostscript's ps2pdf.
suspect you need to track down where conversion to PDF is happening
and disabl
On Thu, May 21, 2009 at 02:14:01PM +0200, Peter Dalgaard wrote:
>
> I think the trick is
>
> jade:~/> env | grep GS_
> GS_OPTIONS=-dAutoRotatePages=/None
>
Thanks, I found that myself. However, when using ps2pdf from Miktex 2.7, I
get the following error:
Unrecoverable error: typecheck in .put
this is the command i made for a normal distribution, but when i try to plot
the histograms, i dont know why the bars don't stick on the line...
nsamples<-1000
sampsize<-15
Samples<-matrix(rnorm(nsamples*sampsize,0,1),nrow=nsamples)
a<-apply(Samples,1,var)
NC14<-a*14
x<-0:40
plot(x,dchisq(x,14),
Hello,
I am a very new user to R so please have patience with me. :clap:
I am trying to evalute the "internal response" for a couple of different
cluster methods with the help of the AdjustedRandIndex, which is included in
the mclust package.
However, I do get a bit puzzled when I get a negat
You can also try:
http://www.rseek.org/
Cheers
Duncan Murdoch wrote:
On 20/05/2009 10:01 AM, cr...@binghamton.edu wrote:
For Google searches, I find that throwing in the term cran on every search
helps weed out irrelevant pages.
For example, instead of
r residuals
I type
r cran residua
> You are very picky. When I enter
>
> R residuals
>
> into Google, 8 out of the first 10 hits are for R topics. Isn't that
> good enough for you?
>
> I think this is true of most Google searches: the letter R most
often
> means the R project.
Although it does not appear to be a factor wit
Thanks! That helps a lot!
A quick follow-up question - I can't really tell what part of the
commands tell it to only look at the child nodes of . Is there any
way to also access the fields that are in the heirarchy? (ie the
S, D, C, and F)
I wouldn't necessarily want those repeated thousands
Hello!
I'm writing a script with a lot of loops and it executes really slowly over
huge amounts of data. I assume it's because I don't know how to avoid using
loops. Logical subscripts are more desirable, but I don't know how to
implement them. One example of that issue:
library(seqinr)
GCsequen
Brigid Mooney wrote:
Thanks! That helps a lot!
A quick follow-up question - I can't really tell what part of the
commands tell it to only look at the child nodes of .
xmlRoot(bri) gives us the C node.
xmlSApply(node, f) is short-hand for
sapply(xmlChildren(node), f)
so that is where we
> this is the command i made for a normal distribution, but when i try to
plot
> the histograms, i dont know why the bars don't stick on the line...
>
> nsamples<-1000
> sampsize<-15
> Samples<-matrix(rnorm(nsamples*sampsize,0,1),nrow=nsamples)
> a<-apply(Samples,1,var)
> NC14<-a*14
> x<-0:40
>
Dear R Users,
I have created a 1500 x 2 data frame - DataSeq. Each of the 1500
rows represents a data sequence.
I have another data frame iData that stores the information of these
1500 data sequences in the same order, for example, condition, gender,
etc.
If I use "subset" to select certain
I am having a slight problem with probabilities.
To calculate the final probability of an event p(F), we can take the product of
the chance that each independent event that makes p(F) will NOT occur.
So...
p(F) = 1- ( (1-p(A)) * (1-p(B)) * (1-p(C))...(1-p(x)) )
If the chance of an event within
Assuming that you get the list of indices into iData for the criteria, then
you can use that to get the appropriate rows:
indx <- which(iData >5) # or whatever your criteria is
DataSeq[indx,, drop=FALSE] # gives you a subset matrix of just the rows you
are interested in.
On Thu, May 21, 2009 at
On 2009-May-21 , at 05:40 , Marcelino de la Cruz wrote:
Jose M. Blanco-Moreno an myself have implemented Syrjala's test in
the ecespa package. As a matter of fact, Jose M. has implemented a
Frotran version of Syrjala's original QBasic function. Using your
data the results are very close to
Madan asks:
> I am trying to import a table from SQL server to R(2.9.0), however i am
> getting errors while running the below codes. Can anyone identify and
let me
> know where did i go wrong??? Thanks in anticipation :)
> ...
> NEWDATASQL1 <- sqlFetch(myconne, CampaignDataLarge)
>
> Error in od
Peter Flom wrote:
charles78 wrote
I have a stupid question on how to get the real p-values for wilcox.test and
correlation. the minmun can be reached is 2.2E-16 using the R version
2.6.2. I do not think it is the R version causing this but other issues.
Any help is highly appreciated.
There are several arbitrary precision packages available: gmp (an
interface to the GNU multi-precision library on CRAN) and
bc (an R interface to the bc arbitrary precision calculator):
http://r-bc.googlecode.com
There are also packages providing R interfaces to two computer
algebra systems and th
On 21-May-09 14:15:08, Mark Bilton wrote:
> I am having a slight problem with probabilities.
>
> To calculate the final probability of an event p(F), we can take the
> product of the chance that each independent event that makes p(F) will
> NOT occur.
> So...
> p(F) = 1- ( (1-p(A)) * (1-p(B)) * (1
Someone knows the existence of a function to sum the elements of the same
place A[i,j] B[i,j] of a matrix?thank you
--
Dr. Daniele Riggi, PhD student
University of Milano-Bicocca
Department of Statistics
Building U7, Via Bicocca degli Arcimboldi, 8
20126 Milano, Italy
cell. +39 328 3380690
mailt
Dear Daniele,
Try this:
A <- matrix(1:10,ncol=2)
B <- matrix(10:20,ncol=2)
A+B
#[,1] [,2]
#[1,] 11 21
#[2,] 13 23
#[3,] 15 25
#[4,] 17 27
#[5,] 19 29
HTH,
Jorge
On Thu, May 21, 2009 at 11:11 AM, daniele riggi wrote:
> Someone knows the existence of a function to sum the el
daniele riggi wrote:
Someone knows the existence of a function to sum the elements of the same
place A[i,j] B[i,j] of a matrix?thank you
Yes:
+
Uwe Ligges
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEAS
For this, you might want to take a look at the first example in ?Reduce.
HTH,
Jorge
On Thu, May 21, 2009 at 11:21 AM, daniele riggi wrote:
> Ok, I know this waybut if i have a lot of matrix, does a function exist to
> sum all of them elements by elements?
>
>
> 2009/5/21 Jorge Ivan Velez
>
>
>
On 5/21/2009 10:15 AM, Mark Bilton wrote:
I am having a slight problem with probabilities.
To calculate the final probability of an event p(F), we can take the product of
the chance that each independent event that makes p(F) will NOT occur.
So...
p(F) = 1- ( (1-p(A)) * (1-p(B)) * (1-p(C))...(1
Someone knows the existence of a function to sum the elements of the same
place A[i,j] B[i,j] of a matrix?thank you
--
Dr. Daniele Riggi, PhD student
University of Milano-Bicocca
Department of Statistics
Building U7, Via Bicocca degli Arcimboldi, 8
20126 Milano, Italy
cell. +39 328 3380690
mailt
Ok, I know this waybut if i have a lot of matrix, does a function exist to
sum all of them elements by elements?
2009/5/21 Jorge Ivan Velez
>
> Dear Daniele,
> Try this:
>
> A <- matrix(1:10,ncol=2)
> B <- matrix(10:20,ncol=2)
> A+B
> #[,1] [,2]
> #[1,] 11 21
> #[2,] 13 23
> #[3,] 15
I know this is easy, but I am stumped:
> gsub("0","K","8.00+00")
[1] "8.KK+KK"
> gsub("+","K","8.00+00")
Error in gsub("+", "K", "8.00+00") : invalid regular expression '+'
In addition: Warning message:
In gsub("+", "K", "8.00+00") :
regcomp error: 'Invalid preceding regular expression'
I d
If your matrices are the same size, you can just add them. If they aren't the
same size, I'm not sure what you are trying to accomplish.
> mat1 <- matrix(1:9, nrow=3)
> mat1
[,1] [,2] [,3]
[1,]147
[2,]258
[3,]369
> mat2 <- matrix(runif(9), nrow=3)
> mat2
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Mark Bilton
> Sent: Thursday, May 21, 2009 7:15 AM
> To: r-help@r-project.org
> Subject: [R] Product of 1 - probabilities
>
>
> I am having a slight problem with probabilities.
Tom La Bone wrote:
I know this is easy, but I am stumped:
gsub("0","K","8.00+00")
[1] "8.KK+KK"
gsub("+","K","8.00+00")
Error in gsub("+", "K", "8.00+00") : invalid regular expression '+'
In addition: Warning message:
In gsub("+", "K", "8.00+00") :
regcomp error: 'Invali
On 21-May-09 14:45:20, Gabor Grothendieck wrote:
> There are several arbitrary precision packages available: gmp (an
> interface to the GNU multi-precision library on CRAN) and
> bc (an R interface to the bc arbitrary precision calculator):
> http://r-bc.googlecode.com
>
> There are also packages
Dear Tom,
Use "\\" before "+":
gsub("\\+","K","8.00+00")
[1] "8.00K00"
See ?regex.
HTH,
Jorge
On Thu, May 21, 2009 at 11:45 AM, Tom La Bone wrote:
>
> I know this is easy, but I am stumped:
>
> > gsub("0","K","8.00+00")
> [1] "8.KK+KK"
>
> > gsub("+","K","8.00+00")
> Error in gsub("+", "K",
You need to escape the '+' since it is used in regular expressions:
> gsub("\\+","K","8.00+00 ")
[1] "8.00K00"
>
On Thu, May 21, 2009 at 11:45 AM, Tom La Bone wrote:
>
> I know this is easy, but I am stumped:
>
> > gsub("0","K","8.00+00")
> [1] "8.KK+KK"
>
> > gsub("+","K","8.00+00")
> Error i
Here are a few ways:
gsub("[+]", "K", "8.00+00")
gsub("\\+", "K", "8.00+00")
gsub("+", "K", "8.00+00", fixed = TRUE)
Note that with gsubfn you can replace several at once as
it is like gsubfn but can take a replacement translation list:
library(gsubfn)
gsubfn(".", list("+" = " plus", "0" = " zer
Hello everybody
How do you subset a data.frame when your boundaries are a combination of
explicit and implicit limits?
For example, I need to subset from the fourth (explicit) to the last
(implicit) column a data.frame named A.
In other languages you would do A[ , 4:]. Would anybody show me the
You can for example use ncol(A) to get the number of columns.
Sarah
On Thu, May 21, 2009 at 12:06 PM, culpritNr1 wrote:
>
> Hello everybody
>
> How do you subset a data.frame when your boundaries are a combination of
> explicit and implicit limits?
>
> For example, I need to subset from the four
For future reference, you may want to look at the symbols function instead of
points (for this example points works, but symbols gives more
options/flexibility, the my.symbols function in the TeachingDemos package gives
even more options for adding symbols to a plot).
Hope this helps,
--
Greg
> a <- matrix(c(1:4),ncol=2)
> b <- matrix(c(5:8),ncol=2)
> a
[,1] [,2]
[1,]13
[2,]24
> b
[,1] [,2]
[1,]57
[2,]68
> a+b
[,1] [,2]
[1,]6 10
[2,]8 12
Atenciosamente,
Leandro Lins Marino
Centro de Avaliação
Fundação CESGRANRIO
Rua Santa
Dear culprit
Try this:
A[ , 4:ncol(A) ]
HTH,
Jorge
On Thu, May 21, 2009 at 12:06 PM, culpritNr1 wrote:
>
> Hello everybody
>
> How do you subset a data.frame when your boundaries are a combination of
> explicit and implicit limits?
>
> For example, I need to subset from the fourth (explicit
Hi Sarah and Jorge,
ncol(). How elegant!
Thank you.
Jorge Ivan Velez wrote:
>
> Dear culprit
> Try this:
> A[ , 4:ncol(A) ]
>
> HTH,
>
> Jorge
>
>
>
> On Thu, May 21, 2009 at 12:06 PM, culpritNr1
> wrote:
>
>>
>> Hello everybody
>>
>> How do you subset a data.frame when your bounda
I would suggest looking at the documentation found in Technical Report #61 in
our department series.
http://mayoresearch.mayo.edu/mayo/research/biostat/
The formulas for priors, losses, etc in Gini splitting are such that I can't
deduce them myself from the source code (and I wrote
Patrick Burns kindly provided an article about this issue called 'The R
Inferno'. However, I will expand a little bit my question because I think it
is not clear and, if I coud improve the code it will be more understandable
to other users reading this messages when I will paste it :)
In my examp
> Several changes in print.survfit, plot.survfit and seemingly in the
> structure
> of ratetabels effect some of my syntax files.
> Is there somewhere a documentation of these changes, besides the code itself?
I agree, the Changelog.09 file is not as comprehensive as one would like.
Specific
retama wrote:
Patrick Burns kindly provided an article about this issue called 'The R
Inferno'. However, I will expand a little bit my question because I think it
is not clear and, if I coud improve the code it will be more understandable
to other users reading this messages when I will paste it
On 21-May-09 16:56:23, retama wrote:
> Patrick Burns kindly provided an article about this issue called
> 'The R Inferno'. However, I will expand a little bit my question
> because I think it is not clear and, if I coud improve the code
> it will be more understandable to other users reading this m
Have you looked at the boxcox function in the MASS package? That may do what
you want.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
> -Original Message-
> I'm looking for algorithms that assist in spreading out crowded labels, e.g.
> labels of points in a scatter plot, in order to obtain a nicer visual
> appearance and better legibility.
>
> I'm probably just stuck because I didn't find the right key words for a
> succ
Hi all,
I got a few comments offline. Here's the final of how to get graphics
remotely using ESS on a mac. It's working fine for us now.
1) change your /etc/sshd_config file. There is a line that reads:
#X11 Forwarding no
change this to:
X11 Forwarding yes
(note no # at start of line)
2) make su
Hello:
Is there a more natural way to get all elements that satisfy a condition
when there are NAs in the sample?
> x=c(1,2,NA)
> x>=2
[1] FALSE TRUENA
> x[x>=2]
[1] 2 NA ## I would expect here to get just "2"
> x[!is.na(x) & x>=2] ## seems a bit cumbersome
[1] 2
Thanks,
Szila
Hi,
Try:
which( x>=2 )
HTH,
Jorge
On Thu, May 21, 2009 at 2:28 PM, Szilard wrote:
> Hello:
>
> Is there a more natural way to get all elements that satisfy a condition
> when there are NAs in the sample?
>
> > x=c(1,2,NA)
>
> > x>=2
> [1] FALSE TRUENA
>
> > x[x>=2]
> [1] 2 NA ##
Dear list,
This might be a topic for r-devel but i may be missing something
obvious.
I don't understand the rationale in the absolute sizes of the point
symbols, and I couldn't find it documented. The example below uses
Grid to check the size of the symbols against a square of 10mm x 10mm
Hello,
I want to use seq with multiple from values and am getting unexpected
(to me) behavior. I'm wondering if this behavior is intentional or not.
> seq(2, by=3, length.out=4)
[1] 2 5 8 11
> seq(3, by=3, length.out=4)
[1] 3 6 9 12
Now if I want the combined sequence, I thought I could p
The dynIdentify and TkIdentify functions work in 2 dimensions (and I think the
thigmophope.labels does as well). The algorithm in spread.labs could be
adapted to 2 dimensions (the example shows a semi-2 dimensional approach) if
you can define what you want to happen (can labels overplot the poi
Deepayan, thank you very much for your response.
I have a general question. And please remember - I am really just a
beginner in R.
Is it truly the case that in order to build quite a basic bar chart
with value labels attached to it I have to be a true R graphics guru -
because the only way to do a
Hi,
I am running an ordered logistic regression model with an interaction, using
the polr command. I am trying to find a way to calculate the marginal
effects and their significance in R. Does anybody have any suggestion?
Thank you!
Enrico
[[alternative HTML version deleted]]
___
Hi,
How do you obtain the limits of the plotting region in a scatterplot3d
plot? `par('usr')' does not seem to give sensible values, and that
vector only has 4 elements (not the expected 6).
Alan
__
R-help@r-project.org mailing list
https://stat.ethz.
On Thu, May 21, 2009 at 12:58 PM, Dimitri Liakhovitski wrote:
> Deepayan, thank you very much for your response.
> I have a general question. And please remember - I am really just a
> beginner in R.
> Is it truly the case that in order to build quite a basic bar chart
> with value labels attache
On Wed, May 20, 2009 at 12:28 PM, Gabor Grothendieck <
ggrothendi...@gmail.com> wrote:
> There is a times class in the chron package.
Perfect! Just what I was looking for.
On Wed, May 20, 2009 at 12:19 PM, jim holtman wrote:
> If you want the hours from a POSIXct, here is one way of doing it
Dear Terry,
first of all, thank you for your immense work. At the moment, I don't
have a small reproducible example for the ratetable difficulty I
have. I will work on it. Maybe the error message I get is of some
information to you.
Error in match.ratetable(m[, rate], ratetable) :
Data has
Dear List,
I am having a problem using gsub to remove dates from a date/time string.
For example:
x<-c("5/31/2009 12:34:00","6/1/2009 1:14:00")
I would like to remove the date and have just the time.
I have tried:
gsub("[0-9+]/[0-9+]/[0-9+]","",x)
and various versions. I think my problem is
Thank you very much, Deepayan!
On Thu, May 21, 2009 at 4:49 PM, Deepayan Sarkar
wrote:
> On Thu, May 21, 2009 at 12:58 PM, Dimitri Liakhovitski
> wrote:
>
>> Deepayan, thank you very much for your response.
>> I have a general question. And please remember - I am really just a
>> beginner in R.
On May 21, 2009, at 4:13 PM, Tim Clark wrote:
Dear List,
I am having a problem using gsub to remove dates from a date/time
string.
For example:
x<-c("5/31/2009 12:34:00","6/1/2009 1:14:00")
I would like to remove the date and have just the time.
I have tried:
gsub("[0-9+]/[0-9+]/[0-9+]"
Thanks for the help and the various options! Putting the + outside the
brackets worked, but I like the strsplit option. Always nice to learn new
functions!
Aloha,
Tim
--- On Thu, 5/21/09, Marc Schwartz wrote:
> From: Marc Schwartz
> Subject: Re: [R] help with gsub and date pattern
> To
Thanks everyone
--
View this message in context:
http://www.nabble.com/Reorder-variables-in-a-dataframe-tp23647222p23660941.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailma
Dear guRus:
I am using lmer for a mixed model that includes a random intercept for a
set of effects that have the same distribution, Normal(0, sig2b). This set
of effects is of variable size, so I am using an as.formula statement to
create the formula for lmer. For example, if the set of rand
Hello,
I'm trying to use the vcd package to analyze survey data. Expert judges
ranked possible features for product packaging. Seven features were
listed, and 19 judges split between 2 cities ranked them.
The following code (1) works, but the side-by-side plots for Cities PX,
SF are shrunk too
Could anybody point me to the latest status of the most user-friendly
debugger in R?
How I wish I don't have to stare at my (long) code for ages and stuck...
Thanks!
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEA
On Thu, 21 May 2009, Zeljko Vrba wrote:
On Thu, May 21, 2009 at 02:14:01PM +0200, Peter Dalgaard wrote:
I think the trick is
jade:~/> env | grep GS_
GS_OPTIONS=-dAutoRotatePages=/None
Thanks, I found that myself. However, when using ps2pdf from Miktex 2.7, I
get the following error:
Unrec
Hi all,
Could anybody point me to a good/robust numerical optimization program
to use in R?
I am doing some MLE fitting. Thanks!
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://w
Hi all,
I am fitting a model to time series of intra-day financial data.
Could anybody point me to some hands-on books about model selection,
model specification test, goodness-of-fit test, feature selection and
statistical time series data analysis? I am looking not for
theoretical or math books
On 21-May-09 23:02:28, David Scott wrote:
> Well most people deal with that problem by not using Acrobat to
> read .pdf files. On linux you can use evince or xpdf. On windows
> just use gsview32. Those readers don't lock the .pdf.
>
> I am with Peter and generally go straight to pdf these days. Th
Dear Jim,
Thanks for your suggestion.
I have a follow-up question for fellow R Users that have followed this thread:
*I used to create two lists by some very flexible criteria to compare
from "iData" and pass the two lists to *ANOTHER FUNCTION* that further
decompose the two lists and do some cas
I use barplot to plot graph, and there is axis of y, but no x, just label,
so where to add x-axis in barplot.
--
Xiaogang Yang
Sensorweb Research Laboratory
http://sensorweb.vancouver.wsu.edu/
Washington State University Vancouver
[[alternative HTML version deleted]]
___
Here are a few more to add to the list.
x <- c("5/31/2009 12:34:00","6/1/2009 1:14:00")
# 1
read.table(textConnection(x), as.is = TRUE)[, 2]
# 2
format(as.POSIXct(x, format = "%m/%d/%Y %H:%M:%S"), "%H:%M:%S")
# 3
library(gsubfn)
strapply(x, "[0-9]*:..:..", simplify = c)
Also see R News 4/1.
Tena koe
Is this the sort of thing you had in mind?
tt <- barplot(sample(1:5))
axis(1, tt, letters[1:5])
HTH ...
Peter Alspach
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of Xiaogang Yang
> Sent: Friday, 22 May 2009 11:3
On Thu, May 21, 2009 at 5:04 PM, Stavros Macrakis wrote:
> On Wed, May 20, 2009 at 12:28 PM, Gabor Grothendieck
> wrote:
>>
>> There is a times class in the chron package.
>
> Perfect! Just what I was looking for.
>
> On Wed, May 20, 2009 at 12:19 PM, jim holtman wrote:
>>
>> If you want the ho
Dear R Users,
I have some dynamic selection rules that I want to pass around for my functions:
>rules <- paste(g$TrialList==1 & g$Session==2)
>myfunction <- function(rules) {
> index <- which(rules)
> anotherFunction(index)
> }
However, I can't find a way to pass around these selection rules
If you are willing to go down one level and work at the grid level
then you can do it without modifying the panel function.
Below gg.ls$name lists the grid object names. Within that
list look at the ones that have rect in their name. Among
those are 5 in success (the 2nd through 6th rect objects
Maybe you mean:
x[which(x>=2)]
El jue, 21-05-2009 a las 14:34 -0400, Jorge Ivan Velez escribió:
> Hi,
> Try:
>
> which( x>=2 )
>
> HTH,
>
> Jorge
>
>
> On Thu, May 21, 2009 at 2:28 PM, Szilard
> wrote:
>
> > Hello:
> >
> > Is there a more natural way to get all elements that satisfy a con
--- Begin Message ---
I use to install the R Site Search Sidebar in firefox.
http://addictedtor.free.fr/rsitesearch/
Unfortunalely you can not install it in new versions
of firefox unlees you disable "checkCompatiblity" and
"checkUpdateSecurity" on firefox (which is not recommendable).
Still I d
1 - 100 of 109 matches
Mail list logo
