HI all,
I am a beginner to R and trying to do some bootstrapping exercise using
meboot package.
would like to test
1. The joint hypothesis of regression coefficients say beta1 =0 and beta2 =
1 and
2. Individual hypothesis of beta2=1.
However, the code engineered by meboot is applicable for
Hi Pablo,
Prof John Fox is a frequent contributor to this forum, and no-one is better
placed than him to respond on your particular point, but I would use the sem
and the polycor packages. There's an example in the documentation of sem which
uses polychoric correlations.
Incidentally you
Hello guys
I have the following questions.
I am currently running various regressions using the lm() command.
As a result I get the following:
http://r.789695.n4.nabble.com/file/n4673551/Bildschirmfoto_2013-08-12_um_11.40.56.png
My question is why R always adds the es although for example they
Dear All,
please let me write about a slightly off-topic, but highly R and
community-related event, as I am more than happy to announce the birth of
the first Hungarian RUG in Budapest started right now:
http://www.meetup.com/Budapest-Users-of-R-Network
We will meet on the 30th of August (around
sir i have not received vol 126 issue 11 and 13. plz send me again
thanks
[[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
Hi,
I would like to understand why the residual standard error, and the degrees
of freedom are changing when I define custom contrasts, which are not
orthogonal.
For example:
y - rnorm(40)
x - factor(rep(1:10,4))
summary(lm(y~x)) #standard model: Residual standard error: 1.103 on 30
degrees of
Hi Listers,
I have generated a random normal distribution using rnorm
y - rnorm(180,116.7,5.4)
How do I get a list of the x values that correspond to each y value that
was generated. If I do a call to y, then I get the y-values but I also need
the x values so I can generate a plot in another
Dear R Users
This is a summary of the things I tried with read.table.ffdf and fixed-width
files. I would like to thank Jan Wijffels and Jan van der Laan for their
suggestions and the time they spent on my problem!
My objective was to import a file with 6'079'455 lines and 32 variables using
What do you expect the 'x' values to be? Are you just looking for the
sequence 1:180 for plotting? 'y' is just a vector of length 180 with
values and has no 'x' associated with it.
On Mon, Aug 12, 2013 at 9:15 AM, Andrew Halford andrew.half...@gmail.comwrote:
Hi Listers,
I have generated a
TMiller thomas.mueller at student.unisg.ch writes:
Hello guys
I have the following questions.
I am currently running various regressions using the lm() command.
As a result I get the following:
[snip]
(it's generally better to just cut-and-paste text into your
message ...)
My
Thanks, that's exactly what I wanted.
Am 12.08.2013 um 15:03 schrieb arun smartpink...@yahoo.com
:
Hi,
You could change ?options()
options(scipen)
#$scipen
#[1] 0
print(1e5)
#[1] 1e+05
op- options()
options(scipen=20)
print(1e10)
#[1] 100
Using the example from
Thanks to Bert and Göran for your responses.
To answer Göran's comment, yes I did plot the Schoenfeld residuals using
plot.cox.zph and the lines look horizontal (slope = 0) to me, which makes
me think that it contradicts the results of cox.zph.
What alternatives do I have if I assume
I'm using lavaan indeed
On Mon, Aug 12, 2013 at 6:21 AM, Jose Iparraguirre
jose.iparragui...@ageuk.org.uk wrote:
Hi Pablo,
Prof John Fox is a frequent contributor to this forum, and no-one is
better placed than him to respond on your particular point, but I would use
the sem and the
Hello,
I really need help, I am completely new in using R and many things were
possible to figure out but not my last problem.
I created a dotchart with the dotchart2 command. On the y axis are my 16
groups and plotted is the mean of the data from each group. And now I want
to add the SD for
Dear All,
I am trying to fit a regression (growth curve) through my data using nls().
I want to learn how this works, so I made up my own simple data:
PGC = function(K, P0, r, t) { (K*P0*exp(r*t)) / (K+P0*(exp(r*t)-1)) }
x - seq(1:20)
y-PGC(6, 0.1, 0.4, x)
Plotting this data yields a nice
Greetings,
I have a user who is running an R program on two different Linux systems. For
the most part, they are very similar in terms of hardware and 64bit OS.
However, they perform significantly different. Under one box the program uses
upwards of 20GB of ram but fluctuates around 15GB of ram
I have this
mcfa-
impulsivity=~imp1r+imp2+imp3+imp4
physical=~phys1+phys2+phys3+phys4
risky=~risk1+risk2+risk3+risk4
selfish=~self1+self2+self3+self4
simple=~simp1+simp2+simp3+simp4r
temper=~temp1+temp2+temp3+temp4
control=~impulsivity+physical+risky+selfish+simple+temper
where: the model is a
On 08/12/2013 10:18 AM, Stackpole, Chris wrote:
Greetings, I have a user who is running an R program on two different
Linux systems. For the most part, they are very similar in terms of
hardware and 64bit OS. However, they perform significantly different.
Under one box the program uses upwards
Alan,
You might want to have a look to the R2STATS package on CRAN. It is a
GUI for GLM and GLMM written in GTK (with the use of the nice RGtk2 and
gWidgets packages by Michael Lawrence and John Verzani).
Don't expect any gain in performance on large datasets though. But at
least the use
Perhaps this clear warning on the manual page for nls (?nls) has
something to do with it:
Warning
Do not use nls on artificial zero-residual data.
The nls function uses a relative-offset convergence criterion
that compares the numerical imprecision at the current parameter
estimates to the
Hi,
Just for fun, I wanted to try function closures.
Here is simple function definition.
The following generateFiles function return function fileGen which will
return a a regular file in directory , every time a call it (fileGen).
generateFiles - function(path, pattern=*.xml)
{
Change your 'fileIndex - fileIndex + 1' to 'fileIndex - fileIndex + 1'
so the fileIndex in the environment of the call to generateFiles() is altered
(instead of creating and forgetting about a variable called fileIndex
in the function created by generateFiles.
generateFiles - function(path,
Thanks Bill.
Can you please tell me why did you change exp pattern=*.xml to
pattern=\\.xml$.
I know that second one also makes sense.
On Mon, Aug 12, 2013 at 9:59 AM, William Dunlap wdun...@tibco.com wrote:
Change your 'fileIndex - fileIndex + 1' to 'fileIndex - fileIndex + 1'
so the
Try changing the assignment to use -
fileindex - fileindex + 1
On Mon, Aug 12, 2013 at 12:44 PM, Ashwani Rao ashwamegh1...@gmail.comwrote:
Hi,
Just for fun, I wanted to try function closures.
Here is simple function definition.
The following generateFiles function return
Because dir()'s pattern argument should be a 'regular expression'
and *.xml as a regular expression matches things you probably
do not want:
file.create(c(one.pxml, two.xml, wxml.three))
[1] TRUE TRUE TRUE
dir(pattern=*.xml)
[1] one.pxml two.xmlwxml.three
dir(pattern=\\.xml$)
Thanks for your help.
Ashwin
On Mon, Aug 12, 2013 at 10:38 AM, William Dunlap wdun...@tibco.com wrote:
Because dir()'s pattern argument should be a 'regular expression'
and *.xml as a regular expression matches things you probably
do not want:
** **
Hi,
I'm studying SVMs and found that if I run SVM in R, Weka, Python their results
are differ. So, to eliminate possible pitfalls, I decided to use standard iris
dataset and wrote implementation in R, Weka, Python for the same SVM/kernel. I
think the choice of kernel does not matter and only
Hi,
You could change ?options()
options(scipen)
#$scipen
#[1] 0
print(1e5)
#[1] 1e+05
op- options()
options(scipen=20)
print(1e10)
#[1] 100
Using the example from ?lm()
ctl - c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt -
Interesting...interesting...
Thank you very much for this. I'll play with this tomorrow see what
happens.
Regards,
Alan.
On 12 August 2013 17:15, Yvonnick Noel yvonnick.n...@uhb.fr wrote:
Alan,
You might want to have a look to the R2STATS package on CRAN. It is a GUI
for GLM and GLMM
Dear All,
Could anyone suggest a quick way to combine all columns in a data frame into a
vector?
For example, I have a data frame of 205 columns with character data types, many
data values are repeated in all the columns. Actually, I would like to
retrieve all the unique values from this data
Sohail:
1. Are they character or factor?
2. ?unlist
unique(unlist(yourframe))
-- Bert
On Mon, Aug 12, 2013 at 1:23 PM, Khan, Sohail skha...@nshs.edu wrote:
Dear All,
Could anyone suggest a quick way to combine all columns in a data frame into
a vector?
For example, I have a data frame
Thanks Bert,
All are character values.
-Sohail
-Original Message-
From: Bert Gunter [mailto:gunter.ber...@gene.com]
Sent: Monday, August 12, 2013 4:35 PM
To: Khan, Sohail
Cc: greatest.possible.newbie; r-help@r-project.org
Subject: Re: [R] combine all data frame columns into a vector.
Hi,
May be this help:
dat1- structure(list(V1 = c(h, f, s, n, r, x, h, t,
u, g), V2 = c(p, j, r, r, i, x, f, b, n,
d), V3 = c(c, o, s, d, f, r, b, p, q, b
), V4 = c(i, g, j, d, y, f, s, q, s, z),
V5 = c(m, j, h, f, b, b, k, j, g, i),
V6 = c(m, w, m, s, o, z, l, h, e, d),
V7 =
Thanks Arun and Bert.
Both options work.
-Sohail
-Original Message-
From: arun [mailto:smartpink...@yahoo.com]
Sent: Monday, August 12, 2013 4:51 PM
To: Khan, Sohail
Cc: R help
Subject: Re: [R] combine all data frame columns into a vector.
Hi,
May be this help:
dat1-
Hi,
I am trying to do a chi sqaure on a set of values, and my different groups
are not even. Is there away to add arbetrary symbols or #s to make the
matrices even? Or do I need to do a different type of pvalue analysis?
S-1:86
B-1:15
V-1:45
table(S)
S
1 2 3 4 5 6 7 8 9 10 11 12 13
Your post is incomprehensible. Please do read the posting guide:
http://www.R-project.org/posting-guide.html
That said, this might help:
http://stats.stackexchange.com/questions/9099/how-to-perform-diallel-analysis-in-r
Kevin
On Sun, Aug 11, 2013 at 1:02 AM, waqas shafqat waqas1...@gmail.com
Hello everyone,
I'm having trouble pulling out paired samples from a data set... I have the
following:
reps-c(4,15) #the variable reps is a list of all paired samples
data
SameName
Individual
Age
Gender
1
4
80
M
Hi,
The question is not clear so not sure this is what you wanted.
dat1- read.table(text=
SameName Individual Age Gender
1 4 80 M
2 15 56 F
3 1 75 F
4 15 56 F
5 2 58 F
6 4 80 M
,sep=,header=TRUE,stringsAsFactors=FALSE)
reps-c(4,15)
dat1$Newcol-as.numeric(dat1$Individual%in%
38 matches
Mail list logo