Salaam Batur gmail.com> writes:
> This maybe unrelated R-question, but I am sure some of you might have a same
> problem like I do. I'd like to run some tests on files which need to be
> opened by a software tool called Business Objects. Right now, I am having
> problem with open the data files w
Hi,
I was wondering if anyone could provide me with help in entering the attached
dataset into R? I've been having a hard time in trying to figure out how to
assemble it into both a frequency table and a bar graph within R. I've been
trying to present the way I have the data arranged in Excel
Dear Millo
Thank you for the prompt and honest answer.
Please accept my appreciation for developing the 'plm' package and for the
excellent documentation associated with it. It was a great place for me to
start, and it made my initial forays into panel data analysis a lot easier.
You have giv
Thomas, thank you for pointing this out!
On another note, do you think coxph() will also take offset(logweight) as part
of the formula? Or just use the argument weight=logweight in coxph()?
Thanks!
John
From: Thomas Lumley
Sent: Sun, February 27, 2011 12:1
Hi Greg,
When you are running R from a script, you want to divert the output
somewhere. Depending exactly how you run it, R may nicely push the
text output to a file, but you should specify where you want the
graphics to go yourself. For example:
## use the PDF graphics device
pdf(file = "boxpl
Dear all,
This maybe unrelated R-question, but I am sure some of you might have a same
problem like I do. I'd like to run some tests on files which need to be
opened by a software tool called Business Objects. Right now, I am having
problem with open the data files with my lower version of Busines
(English Below)Estimados,
Acabamos de confirmar el lugar y la fecha del primer curso de R en Santiago,
Chile. El curso se realizará el día Miércoles 9 y Jueves 10 de Marzo en las
dependencias de la Facultad de Recursos Naturales de la Universidad Andres
Bello, ubicada en República 440 (Metro R
Hi Brant,
My version of Orthodont doesn't seem to have as many levesl of age as
yours but the general idea is
> library(MEMSS)
> data(Orthodont)
> library(reshape) # could use reshape function in stats package instead of cast
> c.orth <- as.data.frame(cast(Orthodont, Subject + Sex ~ age,
> value=
Thanks to you all!
Now I got it!
--
View this message in context:
http://r.789695.n4.nabble.com/Random-Forest-Cross-Validation-tp3314777p3327384.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https:
I will confess to being very wet behind the ears when it comes to scripting, so
with that as a backdrop, here is my question. I can create a box and whisker
plot without any problem from the command line, but incorporating that process
into a bash script has proved problematic. Below is the port
R-helpers:
I would like to measure the correlation coefficient between the repeated
measures of a single variable that is measured over time and is unbalanced. As
an example, consider the Orthodont dataset from package nlme, where the model
is:
fit <- lmer(distance ~ age + (1 | Subject), dat
You can compute the logarithm of it easily enough
> lchoose(54323456, 2345)
[1] 25908.4
Now, what did you want to do with it?
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Jim Silverton
Sent: Monday, 28 February 2011 10:38 AM
T
One way to do it is to use the 'abind' package
> NCurvas <- 10
> NumSim <- 15
> dW <- replicate(NumSim, matrix(rnorm(NCurvas * 3), NCurvas, 3),
+ simplify = FALSE)
> library(abind)
> DW <- do.call(abind, c(dW, rev.along = 0))
> dim(DW)
[1] 10 3 15
-Original Message-
From: r-h
Warning: This is not a helpful answer. Actually, it's a question: Why
do you want to do this? Replacing missing values with row or column
averages and then analyzing the data as if the missing values were not
there is a dangerous thing to do it can produce biased estimates and
understate the true e
Hello. I have the following object which is a list of length NumSim with
each entry being a matrix of dimensions Ncurvas x 3:
dW =
replicate(NumSim,cbind(rnorm(Ncurvas),rnorm(Ncurvas),rnorm(Ncurvas)),simplify=F)
I would like to transform it into an array of dimension Ncurvas x 3 x
NumSim. Does a
any idea how to get R to compute a combination like (54323456, 2345) ?
--
Thanks,
Jim.
[[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:
Take a look at this exchange:
http://r.789695.n4.nabble.com/Partial-correlations-and-p-values-td908641.html
Dimitri
On Sun, Feb 27, 2011 at 2:01 PM, Neeti Sinha wrote:
> Hello,
> I wanted to ask if anybody knows how to calculate part correlation
> coefficient in R. I have to calculate partial c
Here is one solution; mine differs since there should be at least one
item in the range which would be itself:
tm gr
1 12345 1
2 42352 3
3 12435 1
4 67546 2
5 24234 2
6 76543 4
7 31243 2
8 13334 3
9 64562 3
10 64123 3
> d$ct <- ave(d$tm, d$gr, FUN = function(x){
+ # de
Hi Daniel,
If your data is stored in a matrix, the following should work (and be
fairly efficient):
#
dat <- matrix(rnorm(100), nrow = 10)
dat[sample(1:10, 3), sample(1:10, 3)] <- NA
## create an index of missing values
index <- which(is.na(dat), arr.ind = TRUE)
## calculate the row m
have nobody any idea?
i have already try with tapply(d,gr, ... ) but i have problems with the
choose of the function ... also i am not really sure if that is the right
direction with tapply ...
it'll be really great when somebody comes with new suggestion..
10x
--
View this message in context
Hello,
I have some dataset, which i read it from external file using the (data <-
read.csv("my file location")) and read as a dataframe
> is(data)
[1] "data.frame" "list" "oldClass" "vector"
but i have also converted this into a matrix and tried to apply my code but
didnt work.
Anywa
Dear Peter,
perfect *as usual*, many thanks!
Cheers,
Marius
On 2011-02-26, at 17:55 , Peter Ehlers wrote:
> On 2011-02-25 15:57, Marius Hofert wrote:
>> Dear expeRts,
>>
>> I would like to use LaTeX-like symbols in keys via plotmath. Below is a
>> minimal
>> example. I get:
>> Error in fun(k
I am partial to Gary King's book:
Unifying Political Methodology: The Likelihood Theory of Statistical Inference
(University of Michigan Press, 1998)
Cheers
David Cross
d.cr...@tcu.edu
www.davidcross.us
On Feb 27, 2011, at 2:19 PM, Ning Cheng wrote:
> Dear List,
> This problem is more a st
On Sun, Feb 27, 2011 at 9:16 AM, Mike Marchywka wrote:
> I can probably find more "negative" ones if you are interested LOL :)
> I was hoping to find another interesting technical conversation
> to which I could contribute a few thoughts but subjective matters
> do come up with statistical analysi
dear R experts---I have been experimenting with the foreach package
(with doMC) for a while.
my first impression is that it is a very easy way to acquire parallel
processing capabilities. (thanks, revolution R.) the only two
gotchas were about installation (it required an exit and restart), and
Hi Eric,
On Sun, Feb 27, 2011 at 10:06 AM, eric wrote:
> Basic question but still learning
>
> How do I plot two lines (f$equity and f$bh.equity) on one of the three
> graphs under mfrow ? I tried putting brackets around the first plot and
> lines command but that didn't work.
>
> par(mfrow=
Check out Casella and Berger's Statistical Inference. ISBN 978-81-315-0394-2
or http://en.wikipedia.org/wiki/Maximum_likelihood as an online reference.
--Mark J. Lamias
From: Ning Cheng
To: r-help@r-project.org
Sent: Sunday, February 27, 2011 3:19 PM
Subject: [R] MLE estimation
Dear List,
Th
Dear List,
This problem is more a statistic one than a R one.
Any one can recommend me some references website or online paper on
maximum likelihood estimation?I'm now working on that,while still
doubt how to prove that the estimated parameters are normal
distributed.
Thanks for your time and hel
ok, solved it.
In case anyone else is interested:
1st step: import data with convert.factors = FALSE
2nd step: replace only for variable v2 numeric values to labels:
wvs1981$v2 <- factor(wvs1981$v2, labels =
labels(attributes(wvs1981)$stata.info$label.table$v2))
Original-Nachricht
> Date: Sun, 27 Feb 2011 08:27:00 -0800
> From: markkne...@gmail.com
> To: gunter.ber...@gene.com
> CC: r-help@r-project.org
> Subject: Re: [R] Hello!
>
> On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter wrote:
> > Are you a fan of James Joyce? Is the Caps key on your keyboard broken?
> >
> >
Thanks for the idea, Ben. It works, but unfortunately the numbers assigned by
as.numeric() do not correspond to the original coding.
Any other suggestions are very much appreciated.
Best regards, Tobias
Original-Nachricht
> Datum: Sun, 27 Feb 2011 15:13:47 +
> Von: Ben B
Hello,
I wanted to ask if anybody knows how to calculate part correlation
coefficient in R. I have to calculate partial correlation as well as
part-correlation. I found different libraries for partial correlation
but could not find anything related to part correlation.
Thanks for the help in advanc
Greetings,
I am trying to model a time series using arima(). For getting the
model order components(p, d, q and P,D,Q) I am using procedure
discussed in [1] in section 3.2 . It is most likely hit and trial
method based on lower AIC value.
I want to know what is the correct way to find model order
Basic question but still learning
How do I plot two lines (f$equity and f$bh.equity) on one of the three
graphs under mfrow ? I tried putting brackets around the first plot and
lines command but that didn't work.
par(mfrow=c(3,1))
{plot(f$Date,f$equity, col="blue", type="l", main="equity")
l
Hi Jannis,
thank you very much for your help. After conversion to factor everything
is working fine.
Michael
Am 26.02.2011 20:32, schrieb Jannis:
Michael,
i have no experience with the Bayesian models you use, but I could
imagine that the problem is that this classifiedPositive variable i
On 27/02/2011 05:46, Robert A LaBudde wrote:
I think you are over-concerned with the term "pseudo-replication". All
this means is that the error source is nested, and not a full replicate.
What you haven't done, and need to do, is to describe your experiment
in terms of the real variables and
Dear all, I am trying to do a n-fold cross-validation for a regularized
discrimant function analysis using rda from the package klaR. However, I have
problems to predict the groups from the test/validation sample. The exmaples of
the R documantation and some online webpage also do not work. Does
Dear Sacha,
That is the most piteous soliloquy I have read in some time. Here
then a response to soothe your dizzied mind and wrest your rest from
the ornery sprites that so threaten to plague it now.
## Some sample data
water <- as.data.frame(matrix(rnorm(120), ncol = 12))
## Adding your names
That worked! Gracias Jorge
Nico
On 2/27/2011 12:58 PM, Jorge Ivan Velez wrote:
Hi Nicolas,
Try
popn[!rownames(popn) %in% rownames(fish), ]
HTH,
Jorge
On Sun, Feb 27, 2011 at 3:29 PM, Nicolas Gutierrez <> wrote:
Hi!
I have 2 data.frames: "fish" and "popn":
>fish
xloc
A quick look says you are passing N_Male in which is a vector of
values, but you are trying to use it as if it were a character value
in:
y <- DATA[[y]]
try calling it with
comp.plot(i,"N_male",DATA=water)
Also get out of the habit of using 'attach'; it can lead to other issues.
On Sun, Feb 27
Hi Nicolas,
Try
popn[!rownames(popn) %in% rownames(fish), ]
HTH,
Jorge
On Sun, Feb 27, 2011 at 3:29 PM, Nicolas Gutierrez <> wrote:
> Hi!
>
> I have 2 data.frames: "fish" and "popn":
>
> >fish
>
> xloc yloc id birth size weight energy gonad
> 20 15 15 54 -60 107.9 63.0 15952.9 8
Hi Grant,
I don't have a solution, but just to be clearer on your situation:
One row from sheet 2 looks like this:
BC1 BC2 1 01/02/87 33 3 1 03/03/87 44 3
?
Are you using only the first 6 columns for the data to be replicated, and
using the other columns as some sort of indicators on when a seq
dear list!
I KNOW this has been answered a million times before. But, as some might
remember from their "freelance" times as a statistic consultant,
sometimes you're blinded by the facts. I KNOW I have seen this problem
solved, but due to the mental blockade i have been suffering from the
last
I think paste all column together, then using a %in% b.
Maybe others have better solution.
On Sun, Feb 27, 2011 at 9:29 PM, Nicolas Gutierrez wrote:
> Hi!
>
> I have 2 data.frames: "fish" and "popn":
>
>>fish
>
> xloc yloc id birth size weight energy gonad
> 20 15 15 54 -60 107.9 63.0
Hi!
I have 2 data.frames: "fish" and "popn":
>fish
xloc yloc id birth size weight energy gonad
20 15 15 54 -60 107.9 63.0 15952.9 8.0
21 15 15 32 -60 105.1 61.4 15538.8 7.8
914 43 96 -60 118.9 69.4 17573.2 8.8
71 324 64 -60 121.6 71.0 17976.0 9.0
I actually thought that Bert's comment was tastefully funny. I, too, do not
enjoy reading all lowercase emails or emails that use ridiculous abbreviations.
In any case, I am not sure from where Mark gets the idea that the list is way
too negative. In my earlier days of learning S-Plus and R,
Hi:
Using package sos:
# install.packages('sos') # if necessary
library(sos)
findFn('conditional logistic regression')
the following appear to be reasonable candidates to start investigating:
* clogit() in package survival
* clogistic() in package Epi
* clogistCalc()in package saws
Hi Mark,
I agree that comment was completely unnecessary, but Berton Gunter does
contribute (not defending him) especially, in the "caution make sure you
know what you are doing and consult your local statistician comments." I
would look for the nuggets of wisdom in his comments and ignore the
"ex
Dear List,
I am having trouble with a tricky merging task. I have one data sheet that
has dates (continuous) that radio collared individuals were monitored via
telemetry. I have a different sheet containing data from instances where
individuals were recaptured and associated body condition data
On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter wrote:
> Are you a fan of James Joyce? Is the Caps key on your keyboard broken?
>
> -- Bert
Are your snide comments adding anything to the conversation?
Do you allow for the possibility of an ESL speaker writing an email to
this list?
Do you know th
Hi,
> I think you'll need to post a complete example to illustrate the problem.
> I'd suggest posting it to the R-sig-mac list, since it seems to be a problem
> specific to that platform.
>
> Duncan Murdoch
ok, thank you. I will post it there.
Stefan Richter
_
Hi,
> On 11-02-27 7:31 AM, Stefan Richter wrote:
>> Hi,
>>
>> I have a tktoplevel window and in it a tkbutton:
>>
>> AUS<- function()
>> {
>> foo(parameters);
>> }
>> AUSButton.but<- tkbutton(tt,text="OK",command=AUS)
>>
>> The function foo(...) does a time-consuming calculation, and durin
Hello everyone,
I've got a pretty simple exercise to solve with genalg. The goal is to
maximize the revenue.
Even tough genalg is not respecting rather simple constraints.
What's wrong with my code?
library(genalg)
evaluate <- function(X=c()) {
returnVal = NA;
a<-c(X*c(2,0,0,4,8,1,0),X*c(0,2,1,0
Hi:
There has to be a better way to do this, but one option is to use pmin() and
pmax().
# Function to apply to pairs of columns:
ordPair <- function(x, y) data.frame(pmin(x, y), pmax(x,y))
u <- with(odataframe, cbind(ordPair(X1, X2), ordPair(X3, X4), ordPair(X5,
X6))
u
pmin.x..y. pmax.x..y. pm
On 11-02-27 10:09 AM, Stefan Richter wrote:
Hi,
On 11-02-27 7:31 AM, Stefan Richter wrote:
Hi,
I have a tktoplevel window and in it a tkbutton:
AUS<- function()
{
foo(parameters);
}
AUSButton.but<- tkbutton(tt,text="OK",command=AUS)
The function foo(...) does a time-consuming calcul
gmx.at> writes:
>
> Dear R community,
>
> I would like to import data saved with Stata and then
> run a Probit model using R.
>
> My data comes from the World Values Surveys and
> in the Probit model I want to control for countries.
>
[snip]
[snip]
Why don't you set convert.factors
Hello,
I am training to use the changeLOS package. Using data provided in this package
(los.data), I want to compare transition probability P01 and P03 like the
Kaplan-Meier Method.Can someone help me ?
Thank you.
Jan
data(los.data)
my.observ <- prepare.los.data(x=los.data)
my.model <- msmod
I would like to invite those of you in the St. Louis, Missouri area to
join the newly-created St. Louis R user group.
The aim is to create an interdisciplinary venue where practitioners
from various fields in industry, academia, and government can exchange
knowledge and experience using R for data
On 11-02-27 7:31 AM, Stefan Richter wrote:
Hi,
I have a tktoplevel window and in it a tkbutton:
AUS<- function()
{
foo(parameters);
}
AUSButton.but<- tkbutton(tt,text="OK",command=AUS)
The function foo(...) does a time-consuming calculation, and during this
calculation there is printe
Hi,
I have a tktoplevel window and in it a tkbutton:
AUS <- function()
{
foo(parameters);
}
AUSButton.but <- tkbutton(tt,text="OK",command=AUS)
The function foo(...) does a time-consuming calculation, and during this
calculation there is printed some progress information to the console.
Thank you Uwe for the clarification.
Alan Kelly
On 24.02.2011 17:48, Alan Kelly wrote:
> Deal all, as MCMClogit does not allow for the specification of several
> chains, I have run my model 3 times with different random number seeds and
> differently dispersed multivariate normal priors.
> For
Hi all -
I've been working on a lite ORM and database abstraction package for R.
Formatting complex queries by hand has always been an error-prone hassle, so
I've tried to do away with that as much as possible, instead, using R
objects to represent elements of a database system (statements, clause
Dear R community,
I would like to import data saved with Stata and then run a Probit model using
R.
My data comes from the World Values Surveys and in the Probit model I want to
control for countries.
So far I figured out that I should put "convert.factors = FALSE" when using
stata.get() in o
To provide more details:
1) the package I use is the RSNNS (as stated in the topic)
2) for input data to be split I fed in ts() object.. maybe this is a wrong
move.
Does anybody knows what is the type of object that can be fed into the
train() function from the RSNNS package?
The input dat
On 11-02-26 8:26 AM, mathijsdevaan wrote:
Mean doesn't work either... I understand that the message "replacement has 0
items, need 37597770" implies that the function is not returning any values,
but I don't understand why then this is not the case in the example.
DF = data.frame(read.table(text
Look at ?readline for input from the user. For reading Excel files, a good
place to start would be page 26 of the manual:
http://cran.r-project.org/doc/manuals/R-data.pdf
--
Robert Tirrell | r...@stanford.edu | (607) 437-6532
Program in Biomedical Informatics | Butte Lab | Stanford University
O
tsvi sabo wrote:
Hi,i would like to know what is the best way to write a procedure in R,it seems
that when i run a script it doesn't wait for previouse code lines to be
excuted,i would like to know how to force R to hold on a code line until it is
processed,for example:i want to choose from a
On Feb 27, 2011, at 01:12 , Pasicolan wrote:
> Hello Everyone,
>
> *Question: *How do you calculate the sum of squares for a two-way
> _unbalanced_ ANOVA?
There is no unique way of doing it. You can either do a decomposition as
f1 : SS(Y~f1)
f2 : SS(Y~f1+f2) - SS(Y~f1)
res: SStot - SS(Y~
Hi all,
While useing the R package "Spatstat" to detect the spatial point pattern of
my data, I met a problem. When I computes simulation envelopes with a
fitted point process model: Poisson cluster (Thomas) process
(kappa=2.751010e-05; sigma=5.634274e+01; mu=4.943639e+02), I cannot get the
high
69 matches
Mail list logo
