Hi,
I'm having trouble with quantmod's addTA plotting functions. They seem to
work fine when run from the command line. But when run inside a function,
only the last one run is visible. Here's an example.
test.addTA - function(from = 2010-06-01) {
getSymbols(^GSPC, from = from)
library(rgdal)
my_asc=dir(~/Pulpit/dods/karol/TVDI
113_121,pattern=.asc,recursive=T,full.names=T)
for (i in 1:length(my_asc))
{
r - readGDAL(my_asc[i])
z - as.matrix(r)
vectordata[i] - mean(z)
vectordatamax[i] - max(z)
vectordatamin[i] - min(z)
vectordev[i] - sd(z,na.rm=True)
hist(z)
hello,
the following questions will without doubt reveal some fundamental
ignorance, but hopefully you can still help me out.
I'd like to bootstrap a coefficient gained on the basis of the
coefficients in a logistic regression model (the mean differences in
the predicted probabilities between
Hi R users
I was wondering on how to use lapply co when the applied function has a
conditional statement and the output is a 'growing' object.
See example below:
list1 - list('A','B','C')
list2 - c()
myfun - function(x,list2)
{
one_elem - x
cat('one_elem= ', one_elem, '\n')
random -
Hi,
I have two matrices with 4885 cols and 36 cols respectively . I would like
to extract these 36 columns from the bigger matrix, the column values are
different but the column names would be the same.Hence based on the names I
would like to perform my operation. So is there any way of
Thank you very much for your help and for your quick answers.
Lucía Cañás Ferreiro
Instituto Español de Oceanografía
Centro Oceanográfico de A coruña
Paseo Marítimo Alcalde Francisco Vázquez, 10
15001 - A Coruña, Spain
Tel: +34 981 218151 Fax: +34 981 229077
lucia.ca...@co.ieo.es
May 04, 2011; 5:50pm Cristina Silva wrote:
In lattice graphs, panels are drawn from left to right and bottom to
top. The flag as.table=TRUE changes to left to right and top to
bottom. Is there any way to change to first top to bottom and then left
to right? didn´t find anything neither in
Thanks a lot to Steve Lianoglou and Peter Savicky for their help!
Alfredo
-Messaggio originale-
Da: Steve Lianoglou [mailto:mailinglist.honey...@gmail.com]
I'd to match-merge 2 tables in such a manner that I keep all the rows in
table 1, but not the rows that are in both table 1 and
Hi Lorenzo,
On Thu, May 5, 2011 at 8:38 AM, Lorenzo Cattarino l.cattar...@uq.edu.au wrote:
Hi R users
I was wondering on how to use lapply co when the applied function has a
conditional statement and the output is a 'growing' object.
See example below:
list1 - list('A','B','C')
list2 -
The only reason the boot package will take more memory for 2000
replications than 10 is that it needs to store the results. That is
not to say that on a 32-bit OS the fragmentation will not get worse,
but that is unlikely to be a significant factor.
As for the methodology: 'boot' is support
Hi William!
So... many things to say
First, your first steps are completely unnecessary.
Just do this:
data1 - read.csv(file.choose(), header=TRUE)
data1$IND - factor(data1$IND)
M - manova(as.matrix(data1[,-1])~data1$IND, data=data1)
Though I don't understand why this doesn't work:
M -
Hello everyone,
I have a few questions about the print() fn:
1) I have the following code that does not center the character string:
print(The TITLE,quote=FALSE,justify=center)
2) How can I get R to not print the leading [1], etc. when using print()?
(Sorry, I don't know what the leading [1]
On May 04, 2011 at 8:26pm Lucia Cañas wrote:
I would like to combine lattice plot (xyplot) and standard R plot (plot
and plotCI) in an unique figure.
Hi Lucia,
Combining the two systems can be done. See:
Paul Murrell. Integrating grid graphics output with base graphics output. R
News,
Hi,
Try this,
cat(format(The TITLE, width=80, justify=centre))
HTH,
baptiste
On 5 May 2011 19:28, Dan Abner dan.abne...@gmail.com wrote:
Hello everyone,
I have a few questions about the print() fn:
1) I have the following code that does not center the character string:
print(The
Hi
I got a string that looks something like this
1 2 3 4 5 6 7 8 9 ...
and I want it to be
123456789...
So I want to remove all spaces (or whitespaces) from my string.
Anyone know a good way of doing this?
//Joel
--
View this message in context:
Stackoverflow users have attempted to help you improve your question -
please don't ignore advice this with an undeclared cross-post here:
http://stackoverflow.com/questions/5892548/problem-with-loops-in-r
For this list
PLEASE do read the posting guide
On 4 May 2011 15:32, peter dalgaard pda...@gmail.com wrote:
On May 4, 2011, at 15:11 , JP wrote:
Peter thanks for the fantastically simple and understandable explanation...
To sum it up... to find the z values of a number of pairwise wilcox
tests do the following:
# pairwise tests with
Hi Joel,
Try this:
x - 1 2 3 4 5 6 7 8 9
gsub( , , x)
Ivan
Le 5/5/2011 10:50, Joel a écrit :
Hi
I got a string that looks something like this
1 2 3 4 5 6 7 8 9 ...
and I want it to be
123456789...
So I want to remove all spaces (or whitespaces) from my string.
Anyone know a good way
Also take a look at sprintf, which offers everything the C-equivalent
has. sprintf returns a string which can be sent to the console via cat.
Andreas
baptiste auguie schrieb:
Hi,
Try this,
cat(format(The TITLE, width=80, justify=centre))
HTH,
baptiste
On 5 May 2011 19:28, Dan Abner
For those not familiar Ruby Koans is a fantastic platform for teaching many
of the basics of programming in the Ruby language. It uses unit tests
written for methods each of which describe a component of the Ruby
programming language.
http://rubykoans.com/
The koans platform has become popular
the dataset's form is changed after my post
so I repost it here
t
0
0.3403
0.4181
0.4986
0.7451
1.0069
1.5535
1.8049
2.4979
6.4903
13.5049
27.5049
41.5049
V(t)
6.053078443
5.56937391
5.45484486
5.193124598
4.31386722
3.645422269
3.587710965
3.740362689
3.699837726
2.908485019
1.888179494
ID1 ID2 t V(t)
1 1 0 6.053078443
2 1 0.3403 5.56937391
3 1 0.4181 5.45484486
4 1 0.4986 5.193124598
5 1 0.7451 4.31386722
6 1 1.0069 3.645422269
7 1 1.5535 3.587710965
8
sed?
1,$s/ //g
- Original Message -
From: Joel
To: r-help@r-project.org
Sent: Thursday, May 05, 2011 6:50 PM
Subject: [R] Remove all whitespaces
Hi
I got a string that looks something like this
1 2 3 4 5 6 7 8 9 ...
and I want it to be
123456789...
So I want
Thank you. The line of code you give certainly resolves several of the
issues.
I didn't realize that font support is such a tough matter to realize. Let me
express my gratitude to those who provide this for us in R.
On 04-05-11, Prof Brian Ripley rip...@stats.ox.ac.uk wrote:
On Wed, May 4, 2011 at 9:20 PM, Cristina Silva csi...@ipimar.pt wrote:
Hi all,
In lattice graphs, panels are drawn from left to right and bottom to top.
The flag as.table=TRUE changes to left to right and top to bottom. Is
there any way to change to first top to bottom and then left to right?
Hi:
Is this what you have in mind?
m1 - matrix(rpois(100, 10), nrow = 10, dimnames = list(NULL,
paste('V', 1:10, sep = '')))
m2 - matrix(rpois(40, 10), nrow = 10, dimnames = list(NULL,
paste('V', c(2, 5, 7, 10), sep = '')))
colnames(m1)
colnames(m2)
m1[, colnames(m2)]
HTH,
Dennis
On Wed, May
If your whitespace does not only contain regular spaces, this might also
be useful:
x.1 - \t1 2 \n3 4
write(x.1, test.txt)
x.2 - readLines(test.txt)
x.3 - gsub([[:space:]], , x.2)
x - paste(x.3, collapse=)
See ?regex for details on regular expressions in R.
HTH,
Janko
On 05.05.2011 11:02,
Yea it did...
Thanks :)
On Thu, May 5, 2011 at 3:07 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
Is this what you have in mind?
m1 - matrix(rpois(100, 10), nrow = 10, dimnames = list(NULL,
paste('V', 1:10, sep = '')))
m2 - matrix(rpois(40, 10), nrow = 10, dimnames = list(NULL,
I still need to do some repetitive statistical analysis on some outcomes
from a dataset.
Take the following as an example;
id sex hiv age famsize bmi resprate
1 M Pos 23 2 16 15
2 F Neg 24 5 18
Apologies, but I don't see a question here ... am I missing something
obvious?
Andrew
On Thu, May 05, 2011 at 01:20:33AM -0700, sterlesser wrote:
ID1 ID2 t V(t)
1 1 0 6.053078443
2 1 0.3403 5.56937391
3 1 0.4181 5.45484486
4
On 2011-05-05 0:47, Russ Abbott wrote:
Hi,
I'm having trouble with quantmod's addTA plotting functions. They seem to
work fine when run from the command line. But when run inside a function,
only the last one run is visible. Here's an example.
test.addTA- function(from = 2010-06-01) {
Hello, Derek,
see below.
On Thu, 5 May 2011, dereksloan wrote:
I still need to do some repetitive statistical analysis on some outcomes
from a dataset.
Take the following as an example;
id sex hiv age famsize bmi resprate
1 M Pos 23 2
Thank you very much.
Cristina
On 05/05/2011 10:44, Deepayan Sarkar wrote:
On Wed, May 4, 2011 at 9:20 PM, Cristina Silvacsi...@ipimar.pt wrote:
Hi all,
In lattice graphs, panels are drawn from left to right and bottom to top.
The flag as.table=TRUE changes to left to right and top to
Hi all R users
Ive got a file that contains diffrent settings in the manor of:
setting1=value1
setting2=value2
setting3=value3
setting4=value4
.
.
.
What I want to do is open the file and change the value of a specific
setting
like wanna change setting4=value4 - setting4=value5 and then save
Please read the documentation for the rms package, particularly the datadist
function.
Note that in your subject line glm should be lrm.
Frank
Komine wrote:
Hi all R users
I did a logistic regression with my binary variable Y (0/1) and 2
explanatory variables.
Now I try to draw my
try this:
a - readLines(textConnection('setting1=value1
setting2=value2
setting3=value3
setting4=value4'))
closeAllConnections()
# change values
ac - sub('setting4=value4', 'setting4=value5', a)
writeLines(ac, con='myFile.txt')
On Thu, May 5, 2011 at 8:16 AM, Joel joda2...@student.uu.se wrote:
Hi
r-help-boun...@r-project.org napsal dne 05.05.2011 14:16:04:
Joel joda2...@student.uu.se
Odeslal: r-help-boun...@r-project.org
05.05.2011 14:16
Hi all R users
Ive got a file that contains diffrent settings in the manor of:
What file, what is its structure, is it some R object
Dear Boule,
thank you for your interest in hyperSpec.
In order to look into your *problem* I need some more information.
I suggest that we solve the error off-list. Please note also that
hyperSpec has its own help mailing list:
hyperspec-h...@lists.r-forge.r-project.org
(due to the amount of
Well your question is quite general the solution would involve several steps.
Probably the easiest solution would be to read the data in as a dataframe
(using read.table()) and using the '=' as the separator of the columns. Then
change the desired values in the dataframe and save it back as a
jholtman wrote:
a - readLines(textConnection('setting1=value1
setting2=value2
setting3=value3
setting4=value4'))
closeAllConnections()
# change values
ac - sub('setting4=value4', 'setting4=value5', a)
writeLines(ac, con='myFile.txt')
Problem is that I dont know the value on all the
Hi
r-help-boun...@r-project.org napsal dne 05.05.2011 15:13:34:
Joel joda2...@student.uu.se
Odeslal: r-help-boun...@r-project.org
05.05.2011 15:13
jholtman wrote:
a - readLines(textConnection('setting1=value1
setting2=value2
setting3=value3
setting4=value4'))
On May 5, 2011, at 10:58 , JP wrote:
On 4 May 2011 15:32, peter dalgaard pda...@gmail.com wrote:
On May 4, 2011, at 15:11 , JP wrote:
Peter thanks for the fantastically simple and understandable explanation...
To sum it up... to find the z values of a number of pairwise wilcox
tests
Your code may be untested but it works - also helping me slowly to start
understanding how to write functions. Thank you.
However I still have difficulty. I also have some categorical variables to
analyse by age hiv status - i.e. my dataset expands to (for example);
id sex hiv age
Thanks once again Peter, I understand your points -- I fiddled and
googled and read some more and found an eas(ier) route:
library(coin)
t - wilcoxsign_test(a ~ b, alternative = two.sided, distribution =
exact()) # This is equivalent to paired wilcox.test
pval - pvalue(t)
sweet.zscore -
On May 5, 2011, at 1:31 AM, Maximo Polanco wrote:
I have been trying to do a nls model and gives me the error of a
nonnumeric
argument
This is my data set
TT
FFT
V
C
table(file=c:/tt2.txt,header=T)
fit.model - nls(TT~60*(1+alpha*(v/c)^beta),data=tt2,
start=list(alpha=1,
Date: Thu, 5 May 2011 01:20:33 -0700
From: sterles...@hotmail.com
To: r-help@r-project.org
Subject: Re: [R] nls problem with R
ID1 ID2 t V(t)
1 1 0 6.053078443
2 1 0.3403 5.56937391
3 1 0.4181 5.45484486
4 1 0.4986 5.193124598
5 1 0.7451
Hi all R users
I did a logistic regression with my binary variable Y (0/1) and 2
explanatory variables.
Now I try to draw my nomogram with predictive value. I visited the help of R
but I have problem to understand well the example. When I use glm fonction,
I have a problem, thus I use lrm. My
Thank you, it perfectly worked, except I had to modify some codes on the
rep(, maxLen - length(sp)) because it gave me some error, thanks
--
View this message in context:
http://r.789695.n4.nabble.com/Convert-a-presence-absence-matrix-to-a-list-of-presence-only-tp3468479p3498116.html
Sent from
Hello!
I am trying to produce a histogram of measurement data (orientation of
archaeological structures) that are a subject to measurement error.
The normal histogram just computes frequencies, but does not take into
account that a particular value is spread over a range of values (in
my case the
I'm trying to add the exact value on top of each column of an Histogram, i
have been trying with the text function but it doesn't work.
The problem is that the program it self decides the exact value to give to
each column, and ther is not like in a bar-plot that I know exactly which
values are
Dear all,
First of all apologies, I'm pretty newbie, and secondly I know this is not a
question for the forum but I'd be very grateful if you'd help me.
I had problems with normality assumptions, and I tried with weights in the
linear regression this way:
modelw
I am intesrested by your heatmap function (allowing matrix in ColSideColors
option). Can you give your complete code of your function? thanks
--
View this message in context:
http://r.789695.n4.nabble.com/heatmap-3-tp1566584p3498156.html
Sent from the R help mailing list archive at Nabble.com.
Hello all,
I am wondering if there is anyway to create distance matrix for replicated
data
for example,
I have a data like
sample pop id var1var2var3var4
1.1 1 a 1 1 0 1
1.2 1 a 0 0 1 0
1.3 1
This is probably a very simple question but I am completely stumped!I am
trying to do shapiro.wilk(x) test on a relatively small dataset(75) and each
time my variable and keeps coming out as 'NULL', and
shapiro.test(fcv)
Error in complete.cases(x) : no input has determined the number of cases
Hi,
@Steven: Since Beta distribution is a generic distribution by which i
mean that by varying the parameter of alpha and beta we can fit any
distribution.
So to check this i generated a random data from Normal distribution
like
x.norm-rnorm(n=100,mean=10,sd=10);
Now i want to estimate the
A more elegant way would be:
myString-1 2 3 4 5
myString-paste(unlist(strsplit(myString, )),collapse=)
The output will be 12345
Regards,
Som Shekhar
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hello.
I'm trying to use the functions pandit and treebase. They are in the package
apTreeshape. Once I've loaded the package, R responses:
- no function pandit/treebase.
Somebody knows why or what is the reason?
Thanks,
Arnau.
Hi Derek,
You can accomplish your loop jobs by following means:
(a) use for loop
(b) use while loop
(c) use lapply, tapply, or sapply. (i feel lapply is the elegant
way )
---For Loop-
for loops are pretty simple to use and is almost similar to any
other
Dear All,
I am trying to calculate a 95% confidence interval for the difference in two
c statistics (or equivalently D statistics). In Stata I gather that this
can be done using the lincom command. Is there anything similar in R?
As you can see below I have two datasets (that are actually two
which is the maximum large of digits that R has?, because SQL work with 50
digits I think. and I need a software that work with a lot of digits.
Thanks.
Pamela Santelices Elgueta
Estadístico
Instituto Nacional de Estadísticas
Fono: (56-2) 7962491
Paseo Bulnes 209 oficina 82
Santiago
Oops, I hope not too. Don't know why I had the brackets around B+C. My
model is actually A*B+C. And I'm not sure how to obtain the two-way
prediction of AB with C marginalized. Thanks.
Pang
-Original Message-
From: Andrew Robinson [mailto:a.robin...@ms.unimelb.edu.au]
Sent:
Beta is not as general as you think. Its support is limited to [0,1], but you
are trying to fit data that lies outside of its support. Please read about the
beta distribution from a basic stats/prob book.
Ravi.
From: r-help-boun...@r-project.org
Hello,
I am using the party module to estimate the relationship between the
probability of being a student and number of siblings (alive). However, I
need to include a number of relevant covariates. My code is below:
fm3 - mob(Student ~ age + alive + sex2 + cwa + cha + cym | Religion+Servant
+
On May 5, 2011, at 10:01 AM, dereksloan wrote:
Your code may be untested but it works - also helping me slowly to
start
understanding how to write functions. Thank you.
However I still have difficulty. I also have some categorical
variables to
analyse by age hiv status - i.e. my dataset
On May 5, 2011, at 8:20 AM, Laura Bonnett wrote:
Dear All,
I am trying to calculate a 95% confidence interval for the
difference in two
c statistics (or equivalently D statistics). In Stata I gather that
this
can be done using the lincom command. Is there anything similar in R?
Have
Em 5/5/2011 07:34, antu escreveu:
Hello all,
I am wondering if there is anyway to create distance matrix for replicated
data
for example,
I have a data like
sample pop id var1var2var3var4
1.1 1 a 1 1 0 1
1.2 1 a 0
I did read How To Ask Questions The Smart Way
http://www.catb.org/%7Eesr/faqs/smart-questions.html and I don't have
a problem with calling me stupid and/or a winer for posting this...
What's the purpose of an R Mailing list (especially r-help), after all?
IMHO it should be about R users being
Hi,
I need construct box plot graph, but I want keep Groups
order
karla = data.frame(
Groups =
factor(rep(c('CPre','SPre','C7','S7','C14','S14','C21','S21'),
11)),
Time = rep(c(0,7,14,21), 11),
Resp = valor
)
boxplot(Resp~Groups, order=T)
doesn't work.
How do this?
Hi, I wonder whether the following is possible with R, and whether anyone has
done that and can share his/her code with me. I have a correlation matrix, and
I want to create a correlation table that I can copy to Microsoft Word with a
superscript above each correlation, indicating significant
There were no problems when I repeated the test with data provided.
I simply created a vector
Hello-
Does anyone know of packages that could emulate what J. Vermunt's LEM does ?
What is the closest relative in R ?
I use both R and LEM but have trouble transforming my multiway tables in R into
a .dat file compatible with LEM.
Thanks,
David Joubert
Hi,
I am trying to compile a FORTRAN program to call from R under Windows 7 but
I am having problem in the compiling step. To demonstrate this is the
program testit.f:
--
subroutine TESTIT(x,n,m)
dimension x(n)
do 10 i=1,n
10
Dear David,
Thank you for your reply. I have come across rcorrp.cens before. However,
I'm not sure it does quite what I want it to. It seems to compare whether
one predictor is more concordant than another within the same survival
function. I want to see whether one predictor is more
Hello,
I`m tryint to plot some variable over x and y coordinate, but can`t figure
out hot to do it ( I could do a simple scatterplot3d, but there I can`t
change the viewing angle).
My dataset looks something like this (dataframe):
X Y Q95
21 2628711 1104437 0.7723994
22
On May 5, 2011, at 10:48 AM, yoav baranan wrote:
Hi, I wonder whether the following is possible with R, and whether
anyone has done that and can share his/her code with me. I have a
correlation matrix, and I want to create a correlation table that I
can copy to Microsoft Word with a
Dear R users,
May any of you tell me how to simulate data on:
y_t = a+b*y_{t-1} + u_t
where u_t~N(0,sigma^2), b1, and for some constant a.
Many thanks
Tsegaye tseg...@exeter.ac.uk
[[alternative HTML version deleted]]
__
R-help@r-project.org
You are compiling a subroutine not a program and you compile line
should read:
gfortran testit.f -c testit.o
You then reference that object code testit.o in your final loading
stage after compiling other routiens and the main program.
--
Clint BowmanINTERNET:
?arima.sim
- Phil Spector
Statistical Computing Facility
Department of Statistics
UC Berkeley
Maybe you were doing something like
fcv - read.csv('fcv.csv')
instead of
fcv - read.csv('fcv.csv')[1]
(I haven't tested this.)
Tom
On Thu, May 5, 2011 at 8:48 AM, pcc polly...@hotmail.com wrote:
This is probably a very simple question but I am completely stumped!I am
trying to do
The logspline package does density estimation in a different way than KDE, but
it does allow for interval censored data (I know this value is between a and b,
but not where in that range) using the oldlogspline function. This may be what
you are looking for.
--
Gregory (Greg) L. Snow Ph.D.
Thanks. You're right. I didn't see that. I read the ?addTA help page, which
(annoyingly) didn't mention that feature, but I didn't read the ?TA page.
(That page was mentioned as a see also, but not as a must see.)
I don't know what it means to wrap these calls in a plot call. I tried to
put the
Mikael Anderson wrote:
Hi,
I am trying to compile a FORTRAN program to call from R under Windows 7
but
I am having problem in the compiling step. To demonstrate this is the
program testit.f:
--
subroutine TESTIT(x,n,m)
dimension
with(fcv,shapiro.test(case))
--
Clint BowmanINTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.edu
Department of Ecology VOICE: (360) 407-6815
PO Box 47600FAX:(360) 407-7534
Are you really sure that you want to do that?
Read the discussion starting with this post:
http://tolstoy.newcastle.edu.au/R/e2/help/07/08/22858.html for reasons why you
probably don't (yes, the question is about bar plots not histograms, but much
of it will still apply). Near the end of the
There is a struggle in documentation that revolves around being too brief to
be useful and too verbose which then is often ignored. In general, R
proper is far less verbose than quantmod docs - so if you have trouble
with quantmod...
from ?addTA
Value:
addTA will invisibly return an S4
Hi:
Try this:
karla = data.frame(
Groups = factor(rep(c('CPre','SPre','C7','S7','C14','S14','C21','S21'), 11),
levels =
c('CPre','SPre','C7','S7','C14','S14','C21','S21')),
Time = rep(c(0,7,14,21), 11),
Resp = rnorm(88)
)
boxplot(Resp~Groups, data = karla)
Since you didn't
Thanks David,
I did notice that and I got his code to work using wilcox.test for the
continuous variables.
The problem is that when I tried to alter the code to do chisq.test on my
categorical variables there is something wrong with the syntax and I don't
know what.
Derek
--
View this message
The persp function expects z to be a matrix, so you could reshape your data so
that z is a matrix (the reshape function or package may help). Or the
wireframe function in the lattice package expects data more like what you show,
that may be the easiest solution.
--
Gregory (Greg) L. Snow
Russ,
All you have to do is replace
addTA(GSPC.EMA.3, on = 1, col = #ff)
with
plot(addTA(GSPC.EMA.3, on = 1, col = #ff))
etc.
I can sympathize with the documentation frustration, but I think
that much of the documentation in R and in many R packages is
actually very good. I
Sorry, I forgot to give you a quick example of 'wrapping in plot'
plot(addTA(...))
The objects returned by most of the charting functions in quantmod results
from the desire for the functions to be syntactically identical whether
called from inside of chartSeries (e.g. TA=addMACD() ) as they are
On May 5, 2011, at 1:08 PM, dereksloan wrote:
Thanks David,
I did notice that and I got his code to work using wilcox.test for the
continuous variables.
The problem is that when I tried to alter the code to do chisq.test
on my
categorical variables there is something wrong with the syntax
Thanks a lot,
I understand what you say but I'm having problems - maybe with the syntax or
the specific command.
You are right - I have a dataframe to store the data and want to automate
the analysis.
i.e. I want do a chisq.test with to know if alcohol intake (Y/N) differs
between sexes, then
I don't know whether I understood your question, but 1.1, 1.2 , 1.3 all are
subsample of 1 , so, rather than comparing 1000 subsample, comparison of 20
pop level makes more sense in my case.
thanks for query
-
Ananta Acharya
Graduate Student
--
View this message in context:
On May 5, 2011, at 12:40 PM, yoav baranan wrote:
Here is an example for my earlier question.
Say you have a 3x3 correlation matrix:
corrs - matrix(c(0.25,0.32,0.66,0.14,0.24,0.34,0.44,0.34,0.11),
nrow=3, ncol=3, dimnames = list(c('varA','varB', 'varC'),
c('varA','varB', 'varC')))
And
This is not mission critical, but it's bothering me. I'm getting
inconsistent results when I use the $ accessor in the gam formula
*In window #1:*
library(mgcv)
dat=data.frame(x=1:100,y=sin(1:100/50)+rnorm(100,0,.05))
str(dat)
gam(dat$y~s(dat$x))
Error in eval(expr, envir, enclos) : object
Here is an example for my earlier question. Say you have a 3x3 correlation
matrix:corrs - matrix(c(0.25,0.32,0.66,0.14,0.24,0.34,0.44,0.34,0.11), nrow=3,
ncol=3, dimnames = list(c('varA','varB', 'varC'), c('varA','varB', 'varC')))And
another matrix for the sample size of each
On May 5, 2011, at 1:45 PM, dereksloan wrote:
Thanks a lot,
I understand what you say but I'm having problems - maybe with the
syntax or
the specific command.
You are right - I have a dataframe to store the data and want to
automate
the analysis.
i.e. I want do a chisq.test with to
On May 5, 2011, at 1:08 PM, Gene Leynes wrote:
This is not mission critical, but it's bothering me. I'm getting
inconsistent results when I use the $ accessor in the gam formula
*In window #1:*
library(mgcv)
dat=data.frame(x=1:100,y=sin(1:100/50)+rnorm(100,0,.05))
str(dat)
Dear All,
I am trying to build a package for a set of functions. I am
able to build the package and its working fine. When I check it with
R CMD check
I get a following warning : no visible global function
definition for biocLite
I have
Hello R-Help
How can i exctact and store the within group mean squared difference from an
anova summary table into a varible.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
1 - 100 of 145 matches
Mail list logo