Thank you, Steven. It helps!
Best,
Dana
Steven McKinney wrote:
>
> Hi Dana
>
>> -Original Message-
>> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED]
>> On Behalf Of Dana77
>> Sent: Wednesday, December 03, 2008 3:24 PM
>> To: r-help@r-project.org
>> Subject: [R] Help with reading
Hi,
I obtained custom code for the infamous HL test, which I know is frowned
upon here.
In any event, I'm feeding it a matrix and a data frame. How can I set the
function to ignore missing values?
I tried:
HL <- function(p, y, NA.rm=T)
{
options(na.action=na.exclude)
y <- y[order(p)];
p <- p
On Fri, Dec 5, 2008 at 12:11 AM, Sean Davis <[EMAIL PROTECTED]> wrote:
> On Thu, Dec 4, 2008 at 5:21 AM, Rainer M Krug <[EMAIL PROTECTED]> wrote:
>> Did you make any progress in using R in the SGE / OpenMPI environment?
>> I am asking, because I am struggling with the same problem. Using qsub
>> do
I have problems converting my dataset from long to wide format. Previous
attempts using reshape package and aggregate function were unsuccessful as they
took too long. Apparently, my simplified solution also lasted as long.
My complete codes is given below. When sample.size = 1, the execu
Hello,
I'm doing bayesian network analyses with the deal package. I am at a loss
for how to interpret output from the analysis (i.e. what is a good score, what
is a bad score, which stats tell me what about the network edges/nodes).
Here is an example node with its posterior scores for all
Hi William,
Thanks for your suggestion. I learn alot from it, some concepts are new which
need some digestion. > Subject: Re: [R] How to optimize this codes ?> Date:
Thu, 4 Dec 2008 11:49:30 -0800> From: [EMAIL PROTECTED]> To: [EMAIL PROTECTED]>
CC: [EMAIL PROTECTED]> > > -Original Messa
Making data, especially data that have been the subject of published
papers, widely available, can be a useful spinoff from the R project,
another gift to the scientific community beyond the provision of
computing and analytic tools. Nowadays, in a complete publication of
a scientific resu
... Of course what I should have said is:
lmer( y ~ ... + (-1+A|Grp) + (-1+B|Grp) + (-1+C|Grp))
Cheers,
Bert Gunter
Genentech
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://ww
On Thu, 4 Dec 2008, Johannes Huesing wrote:
Hi all,
I can't come over a problem with formula. Suppose I have a coxmod model
with the following formula:
somemod$formula
Surv(lebzeit, tot == 1) ~ sex + (alter >= 65) + diff3k + zelltyp_k_c +
q_nuc_3k + kar_k80_g80 + stadium
You need to give
Folks:
Suppose I have 3 random effects, A,B, and C. Using the older lme() function
(in nlme) it was possible (using the pdMat classes) to specify that they are
uncorrelated with identical variances. Is it possible to do this with lmer?
My understanding is that if I specify them as
lmer( y ~ ...
Hi:
I'm trying to estimate a latent variable model in mnl discrete choice
framework using R. I need to do first a uni dimensional integral
within each observation (row) in the database and then sum over
observations. I'm stacked in the point shown below. Apparently I have
a dimensionality problem
Skotara wrote:
Dear all,
I apologize for my basic question.
I try to calculate an anova for repeated measurements with 3 factors
(A,B,C) having 2, 2, and 7 levels.
or with an additional fourth between subjects factor D.
Everything works fine using
aov(val ~ A*B*C + Error(subject/ (A*B*C) ) )
On Thu, Dec 4, 2008 at 5:21 AM, Rainer M Krug <[EMAIL PROTECTED]> wrote:
> Did you make any progress in using R in the SGE / OpenMPI environment?
> I am asking, because I am struggling with the same problem. Using qsub
> does not give me an interactive session which I can use (wherefore I
> don't s
Hi list,
Well, this time Ive a doubt with mapping generation.
I was already able to read and plot shapefiles, plot point on this map. All
this without any problems, but now I want to do something that I think, only
Golden Software Surfer is capable of.
I would like to plot a kriging result on t
Dear R-help,
I was comparing SAS (I do not know what version it is) and R (version
2.6.0 (2007-10-03) on Linux) survival analyses with time-dependent
covariates. The results differed significantly so I tried to understand
on a short example where I went wrong. The following example shows that
Hi all,
I can't come over a problem with formula. Suppose I have a coxmod model
with the following formula:
> somemod$formula
Surv(lebzeit, tot == 1) ~ sex + (alter >= 65) + diff3k + zelltyp_k_c +
q_nuc_3k + kar_k80_g80 + stadium
and I want to drop the stadium explanatory variable from the m
> -Original Message-
> From: William Dunlap
> Sent: Thursday, December 04, 2008 9:59 AM
> To: '[EMAIL PROTECTED]'
> Cc: 'R help'
> Subject: Re: [R] How to optimize this codes ?
>
>[R] How to optimize this codes ?
>Daren Tan daren76 at hotmail.com
>Thu Dec 4 17:02:49 CET 2008
>
agrep{base} is very useful, but it would be very useful to also have the
offset of the first match returned, e.g.
agrepx("foo","sdxforzoo",max.distance=1) => 4 (there is room for quibbling
over the definition of first match...). As far as I can tell, the only way
to do that now is to do multiple ag
Dear all,
I apologize for my basic question.
I try to calculate an anova for repeated measurements with 3 factors
(A,B,C) having 2, 2, and 7 levels.
or with an additional fourth between subjects factor D.
Everything works fine using
aov(val ~ A*B*C + Error(subject/ (A*B*C) ) ) or
aov(val ~ (D
Just an update to this - it is possible to connect to an 11G database
via using the oracle 10G (libsqora.so.10.1) driver.
Quoting "Prof Brian Ripley" <[EMAIL PROTECTED]>:
On Fri, 28 Nov 2008, Simon Collins wrote:
Hi
I'm presently trying to connect to Oracle through RODBC / UnixODBC
on
You can get all the rownames by calling write.table once on something
like
# if all list elements have equal length
t(data.frame(MyList))
# otherwise, let NA extend where needed
t(do.call(cbind, lapply(MyList, ts)))
You can take care of NAs as needed on your end.
HTH
-Original Message---
csiro.au> writes:
>
>
> Hi List,
>
> I am unable to read in a 7.8Gb ascii grid using readAsciiGrid {maptools} - R
runs out of memory. I have 4Gb ram
> and 4Gb swap, so things are getting tight.
This function is written for small rasters only, and is not suited to data of
this size.
>
> I
Yes, that is basically truncation, but the original poster said they wanted
something fast and did not give detail. If they just wanted some data to put
into a glm model to demonstrate, then this would work.
Another idea that I had is instead of dropping all the extreme values, is take
the ext
Checkout the portfolio package.
On Thu, Dec 4, 2008 at 9:44 AM, <[EMAIL PROTECTED]> wrote:
> Hi all
>
> Is there an R package that enables me to generate "squarified treemap" plots?
>
> Thanks
> Rory
>
> Rory Winston
> RBS Global Banking & Markets
> 280 Bishopsgate, London, EC2M 4RB
> Office: +4
[R] How to optimize this codes ?
Daren Tan daren76 at hotmail.com
Thu Dec 4 17:02:49 CET 2008
How to optimize the for-loop to be reasonably fast for
sample.size=1 ?
You may want to change sample.size=1000 to have an idea what I am
achieving.
set.seed(143)
A <- mat
There are the multcomp and multtest packages as well as various tools in other
packages (estimable in gmodels, TukeyHSD in stats, etc.). Do any of these do
what you want? If not, then some more detail would help us help you.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain
V interesting point Greg.
But are you not just suggesting left and right truncation?
It strikes me that if the data are Poission then a mixture is likely to be
better - something akin to zero-deflated.
Neg binomail works for greater variance == mix a gamma and poisson, but I'm
unsure what to mix t
I'm trying to establish the occurrence of a specific day of the week in a
given month.
Namely I want to know if a specific Tuesday was the 2nd,3rd Tuesday in a
month.
Is there a good way to this in R?
thank you,
Alina Sheyman
[[alternative HTML version deleted]]
___
What are you planning to do with the file after it is created?
The best way depends on what you want to do. The write.table function may not
be the best choice, there are several functions for writing data to files
including dump, write, cat, save, and probably others, some of these will write
2008/12/4 João Cruz <[EMAIL PROTECTED]>:
> I need help with these R simulations
Try asking 'jrrbc', who has asked the exact same questions on this forum:
http://www.artofproblemsolving.com/Forum/viewtopic.php?p=1338520#1338520
where he does at least admit this is Master's course homework. Oh
ha
thank you all, Jagat, Bert and Prof. Ripley for your kind help.
The PROBLEM was:
give descripitve names to the columns in
JJ<-data.frame( c( as.character(rep( gender,3))) , c( F76,6- F83, F90) )
Possible SOLUTIONs are:
1.JJ <- data.frame( gender=c( as.character(rep( gender,3))) ,
J.value=c(
Accurate probably depends on what exactly you want to accomplish with the
simulation, but one simple approach would be to generate more data than you
need from a regular Poisson distribution, then throw away the extreme values on
both sides.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Cent
Your teacher did not assign this homework to you because he/she does not know
the answer, rather he/she has some vague hope that by You doing Your homework
You may actually learn something. The chances of You learning something
(supposedly the purpose behind you taking the class) is greatly inc
I'm trying to establish the occurrence of a specific day of the week in a
given month.
Namely I want to know if a specific Tuesday was the 2nd,3rd Tuesday in a
month.
Is there a good way to this in R?
thank you,
Alina Sheyman
[[alternative HTML version deleted]]
___
on 12/04/2008 10:38 AM Lars Beckmann wrote:
> Hey,
>
> two questions regarding setwd and list.files:
>
> I try to use the function list.files:
>
>> datapath<-'/nfs/home/m/cremer/gaw/data/traws_chrs'
>
>> fl<-list.files(datapath, pattern="\\.tped")
>
> but I get the following error message:
>
hi there,
I have a list called MyList.
MyList[[1]]
[1] "SRY"
[2] "RPS4Y1"
[3] "ZFY"
Hey,
two questions regarding setwd and list.files:
I try to use the function list.files:
>datapath<-'/nfs/home/m/cremer/gaw/data/traws_chrs'
>fl<-list.files(datapath, pattern="\\.tped")
but I get the following error message:
Warning message:
In list.files(datapath, pattern = "\\.tped") :
li
Does R has some option on calculating the Guttman's scale?
thank you.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.or
How to optimize the for-loop to be reasonably fast for sample.size=1 ?
You may want to change sample.size=1000 to have an idea what I am achieving.
set.seed(143)
A <- matrix(sample(0:1, sample.size, TRUE), ncol=10, dimnames=list(NULL,
LETTERS[1:10]))
B <- list()
for(i in 1:10) {
Hi,
Thank you Stefan,
I think, this package will do what I want.
Edwin
> Edwin Sendjaja schrieb:
> > Hi,
> >
> > Can someone help me how to create a table for Min, Max, Mean, Deviance.
>
> Try the summaryBy function of the doBy package. The package has an
> excellent vignette (=documentation).
On Thu, 4 Dec 2008, Spiderschwein wrote:
Say F_n is my ecdf, generated from my sample x.
With sample(x, 2000, replace = TRUE) I get a uniformly distributed sample
from x.
Which is not the same as uniformly from unique(x).
But I need a (re)sample, which is F_n distributed.
Now I'm confused,
Hi,
If you have Python running on your system you could create this script
(graphics.py):
import win32api
import win32con
import win32gui
def windowEnumerationHandler(hwnd, windowlist):
windowlist.append((hwnd, win32gui.GetWindowText(hwnd)))
mywindows = []
win32gui.EnumWindows(windowEnumerat
Hi,
I am wondering if there is a way to change the value of the "record" option
in a graphics device that is already open (and accepts this option). I
don't want to open a new device with, for example "dev.new(record=T)", but
just want to change the settings of the current device. This can be
Dear R help,
I use the package plm e the function plm() to analyse a panel data and
estimate a dynamic model.
How can I estimate a model without intercept?
Thanks,
Andrea Ferroni
[[alternative HTML version deleted]]
__
R-help@r-project.org
Hello,
Anyone who knows a fast and accurate algorithm for generating draws from an
underdispersed Poisson distribution. Or even better, if there is a package
containing such an implementation.
Thanks
Rene
__
R-help@r-project.org mailing list
https:/
I'm forwarding to this list as I've got no
answer in [R-sig-Geo] and perhaps this
issue is not exclusive of readOGR()
I'm using ubuntu 8.04 andR version 2.8.0 (2008-10-20)
Original Message
Subject: readOGR and non-ascii characters
Date: Wed, 03 Dec 2008 14:22:34 +0100
From: Ag
I need help with these R simulations
___
The generation of pseudo-random numbers with uniform distribution [0.1]
can be done through a method called Congruencial Committee, which consists
the use of a recursive relationship of the type:
xi +1 = (a.xi + c) mod
Hello,
I want to plot replicates of an experiment which are already grouped by a
different factor than the replicate number into a single panel using
the xyplot function from the lattice lib.
Example:
My data.frame consists of the following collumns:
x, y, factor_for_y, replicate_no, se
Say F_n is my ecdf, generated from my sample x.
With sample(x, 2000, replace = TRUE) I get a uniformly distributed sample
from x.
But I need a (re)sample, which is F_n distributed.
Now I'm confused, if that makes a difference or not...
Regards,
Martin
Dimitris Rizopoulos-4 wrote:
>
> say 'x'
Edwin Sendjaja schrieb:
> Hi,
>
> Can someone help me how to create a table for Min, Max, Mean, Deviance.
>
>
Try the summaryBy function of the doBy package. The package has an
excellent vignette (=documentation).
hth
Stefan
__
R-help@r-project.org
Hello,
I am having two difficulties with xspline:
1) Passing lty=2 (via ...) appears to have no effect
2) Apparently depending on the current plot axes, I may or may not
get the failure: add_point - reached MAXNUMPTS (25200)
This happens both on Linux/R 2.8.0 (details below) and Wi
Hi all
Is there an R package that enables me to generate "squarified treemap" plots?
Thanks
Rory
Rory Winston
RBS Global Banking & Markets
280 Bishopsgate, London, EC2M 4RB
Office: +44 20 7085 4476
***
The Royal
If I understand your situation correctly, you may be able to make use of
the "strata" and "sampsize" arguments in randomForest() to get bootstrap
samples that resemble the original data distribution. They allow you to
specify stratified samples using the "strata" variable.
Best,
Andy
From: Ragh
Hi Benjamin,
thanks for the hint.
In barplot.default I found "grey.colors(nrow(height))" :-)
(I had the same problem with the help file for grey.colors)
Ciao,
Antje
Nutter, Benjamin schrieb:
I believe the defaults in barplot are found using
grey.colors{GrDevices}.
?grey.colors
(For some re
Hi
use the logarithmic representation for your problem.
The Brobdingnag package uses such a form in a (more-or-less)
user-transparent manner.
HTH
rksh
Marc Jekel wrote:
Dear R Fans,
I have a simple probem but cannot find any reference to the soultion.
I want to do calculations with small n
pufftissue pufftissue wrote:
Hi,
When I use logistic regression, each variable has a p value associated with
it. Do I only include the variables that have a statistically significant p
value (<0.05), or are there situations when I should include variables when
their p values are high? I had he
jestinah mahachie wrote:
Dear all,
I have been using rpart to get trees. May someone help me with how I can do
some logistic regression modelling to derive odds ratios on the branches.
Also help on how to derive estimates on prediction error.
--
Kind Regards,
Jestinah
The odds ratios from
ExpeRts,
I'm trying to compare three survival curves using the function "survdiff" in
the survival package. Following is my code and corresponding error message.
> survdiff(Surv(st_months, status) ~ strata(BOR), data=mydata)
Error in survdiff(Surv(st_months, status) ~ strata(BOR), data = mydata)
Puff -
There are many strategies, ideas, and literature on this topic. A great
introduction that leads to many of the references that are interesting
is Frank Harrell's book, "Regression Modeling Strategies". I would
highly recommend it.
pufftissue pufftissue wrote:
Hi,
When I use logi
Dear R Fans,
I have a simple probem but cannot find any reference to the soultion.
I want to do calculations with small numbers (for max likelihood estimations).
The lowest value R is storing by default is 1*10^-323, a smaller numer like
1*10^-324 is stored as a 0. How can I circumvent this pro
Gabor Grothendieck wrote:
> R does not support tail recursion so not using it would
> not matter.
>
>
this is a feature, i guess.
another issue is that the recursive calls receive substrings as
arguments, which means copying the whole remaining content, and the
copies would not be deallocated
Well, if it is transparent then you will not see it: I presume you mean
translucent (aka semi-transparent). With 50% opacity you could use
red <- rgb(0x9E, 0x47, 0x70, 127, max=255)
polygon(xvalue1,yvalue1,col=red,border=red)
See ?rgb.
You can also specify a translucent colour as e.g. "#9E477
R does not support tail recursion so not using it would
not matter.
On Thu, Dec 4, 2008 at 5:04 AM, Wacek Kusnierczyk
<[EMAIL PROTECTED]> wrote:
> John Fox wrote:
>> By coincidence, I have a version of strsplit() that I've used to
>> illustrate recursion:
>>
>> Strsplit <- function(x, split){
>>
Hi list,
Does any one know how I can make the red colour transparant?
#
par()
#Set up blank plot to remove plot border type lty="n"
xvalue1<- c(1,2,2,1)
xvalue2<- c(1,3,3,1)
yvalue1<- c(0,0,22,22)
yvalue2<- c(0,0,44,44)
plot(c(min(xvalue2),max(xvalue2)),c(min(yvalue2),max(y
I believe the defaults in barplot are found using
grey.colors{GrDevices}.
?grey.colors
(For some reason my machine won't pull up the help files for grey.colors
from the command line, but I can still access it through the html help).
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMA
John Fox wrote:
> Dear Wacek,
>
> "Wrong" is a bit strong, I think -- limited to single-pattern characters is
> more accurate.
nothing is ever wrong if seen from an appropriate perspective. for
example, there is nothing wrong in that many core functions in r deparse
some, but not all, of the arg
Dear Wacek,
"Wrong" is a bit strong, I think -- limited to single-pattern characters is
more accurate. Moreover, it isn't hard to make the function work with
multiple-character matches as well:
Strsplit <- function(x, split){
if (length(x) > 1) {
return(lapply(x, Strsplit, split)) #
Hello,
Can anybody help me to find out which colors are used automatically when
calling barplot (e.g. 3 series beside each other will get different gray values).
I want to apply a legend but I don't know the colors used...
Antje
__
R-help@r-project.
say 'x' is your sample, then try
sample(x, 2000, replace = TRUE)
I hope it helps.
Best,
Dimitris
Spiderschwein wrote:
Good day everyone,
I have about 500 values to generate my ecdf.
Now I want to sample 2000 values i.i.d. distributed according to my ecdf.
Is it possible to sample from an
Good day everyone,
I have about 500 values to generate my ecdf.
Now I want to sample 2000 values i.i.d. distributed according to my ecdf.
Is it possible to sample from an ecdf?
If yes, I would appreciate your help.
Martin
--
View this message in context:
http://www.nabble.com/Resampling-fro
Did you make any progress in using R in the SGE / OpenMPI environment?
I am asking, because I am struggling with the same problem. Using qsub
does not give me an interactive session which I can use (wherefore I
don't see a way of using RMPISNOW and getMPIcluster()) and the
spawning (makeMPIcluster(
John Fox wrote:
> By coincidence, I have a version of strsplit() that I've used to
> illustrate recursion:
>
> Strsplit <- function(x, split){
> if (length(x) > 1) {
> return(lapply(x, Strsplit, split)) # vectorization
> }
> result <- character(0)
> if (nchar(x) == 0) r
Dear all,
> I have been using rpart to get trees. May someone help me with how I can do
> some logistic regression modelling to derive odds ratios on the branches.
> Also help on how to derive estimates on prediction error.
>
> --
> Kind Regards,
>
> Jestinah
>
>
>
> Home Address:
> Brusselsepoo
On 12/4/2008 12:59 AM, Gundala Viswanath wrote:
> Dear all,
>
> How can I plot the data below with scatter plot
> so that the next to the marker it contain the label for "species"
>
> __DATA__
> class true_hits false_hits species
> big 2350180 lion
> big 375 20 dog
We are happy to inform you that the online abstract submission and
registration for `useR! 2009' is now available online from
http://www.R-project.org/useR-2009/
This meeting of the R user community will take place at Agrocampus,
Rennes, France, July 8 to 10, 2000. The conference schedule comp
Edwin Sendjaja wrote:
Hi,
Can someone help me how to create a table for Min, Max, Mean, Deviance.
Hi Edwin,
You might find that the describe function in the prettyR package does
what you want.
Jim
__
R-help@r-project.org mailing list
https://st
76 matches
Mail list logo
