well, Harrell,
I understand sweave or R2html could be a solution.
but please show me their applications in a large business setting. On
the contrary, I can give you many such cases using SAS.
On 9/13/06, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote:
> Wensui Liu wrote:
> > For your 1st question
Wensui Liu wrote:
> For your 1st question, you can write query against the tables in DB using
> RODBC.
>
> Being a SAS programmer, I have to say that reporting function of R is
> not as good as that of SAS.
I beg to differ. See for example
http://biostat.mc.vanderbilt.edu/StatReport
Frank Har
For your 1st question, you can write query against the tables in DB using RODBC.
Being a SAS programmer, I have to say that reporting function of R is
not as good as that of SAS.
On 9/13/06, Thorsten Muehge <[EMAIL PROTECTED]> wrote:
>
>
> Hello Colleagues,
> I programmed in SAS for 3 years and
Thorsten Muehge de.ibm.com> writes:
> 1. Is it possible to query R files by SQL internally on data frames (not on
> a database) and how is the syntax (I have the RODBC package installed).
It is possible to do similar things "conceptually" in R as in SQL---at least the
basic SQL queries (I have n
Thorsten Muehge wrote:
>
> Hello Colleagues,
> I programmed in SAS for 3 years and would like to switch to a not so costly
> software product.
>
> Hence I started to evaluate R, and my first test look promising.
>
> However I have some question:
>
> 1. Is it possible to query R files by SQL int
Hello Colleagues,
I programmed in SAS for 3 years and would like to switch to a not so costly
software product.
Hence I started to evaluate R, and my first test look promising.
However I have some question:
1. Is it possible to query R files by SQL internally on data frames (not on
a database)
Hi,
I'm writing a function f() which includes probability density function
computation, I would like to make this generic, that is , I will specify the
density when I call f(), here's what I wrote: do.call(den, args=list(d,
arg))
the question is , say, if I let den=dnorm, I have t
Hi Gunther,
Gunther Höning wrote:
> Dear list,
>
> I'm trying to find out the following in an AffyBatch.
This question is related specifically to a Bioconductor package, so
should be asked on the bioconductor listserv rather than R-help.
Best,
Jim
>
> To get the indices from a loction on a
Dear list,
I'm trying to find out the following in an AffyBatch.
To get the indices from a loction on a chip I use the function xy2i() for
the hgu133plus2 by Affymetrix.
But now I want to know the name of the probe located at this spot. How can
this be done?
And what about the locations used as Q
Dear Pref. Fox
Sorry, I didn't receive your reply. I try the new sem package. It's
great. The following is the results that I got. The fit indices are
fine.
Model Chisquare = 208 Df = 98 Pr(>Chisq) = 6.6e-10
Chisquare (null model) = 1741 Df = 120
Goodness-of-fit index = 0.9
Adjusted
rio
Canada L8S 4M4
905-525-9140x23604
http://socserv.mcmaster.ca/jfox
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Guo Wei-Wei
> Sent: Monday, September 04, 2006 12:34 AM
> To: r-help@stat.math.ethz.c
Dear all,
I met a problem while doing SEM by sem package. I got a negative
chi-square of null model. Because the theoretical value of chi-square
cannot be negative, I checked the source code of sem.R in sem package
and I found the Chi-square of null model was computed by the following
expression:
Hello
I have 2 questions:
1. How I can do the test of hartley for homogeneity of variances?
2. How I can do the test of Box for homogeneity of variances?
Thanks in advance
--
Marco Vinisio Martinez
Profesor Departamento de Matematicas y Departamento de Biologia
Pontificia Universidad Javeriana
Check out with XLSolutions Corp
[EMAIL PROTECTED]
--- jeffrey ludlow <[EMAIL PROTECTED]> wrote:
> Dear R users,
>
> The R Project website doesn't seem to have any
> links devoted to R training. Are there any R
> trainers out there?
>
> Thank you for your help!
>
> Jeffrey V.
Dear R users,
The R Project website doesn't seem to have any links devoted to R training.
Are there any R trainers out there?
Thank you for your help!
Jeffrey V. Ludlow
-
[[alternative HTML version deleted]]
__
Try
library(MASS)
y.lmr <- lm.ridge(y ~., longley)
coefs <- print(y.lmr)
On 8/21/06, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> Dear all,
>
> I am trying to calculate the fitted values using a ridge model
> (lm.ridge(), MASS library). Since the predict() does not work for lm.ridge
> object,
ng
process." - George E. P. Box
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED]
> Sent: Monday, August 21, 2006 2:40 PM
> To: r-help@stat.math.ethz.ch
> Subject: [R] question about 'coef' method
Dear all,
I am trying to calculate the fitted values using a ridge model
(lm.ridge(), MASS library). Since the predict() does not work for lm.ridge
object, I want to get the fitted_value from the coefficients information.
The following are the codes I use:
fit = lm.ridge(myY~myX,lambda=la
Hi, everyone,
I got puzzled with this "varmx.pca.fd" function which is in the fda
package. According to the description, it carries out the rotation using
the "VARIMAX criterion".
But in fact, if you use this "varmx.pca.fd" on any objects of class
"pca.fd", it will devide the variance equally for
On Sat, 19 Aug 2006, [EMAIL PROTECTED] wrote:
> Dear all,
>
> I have a question about how to get a matrix by combining a large number of
> columns from a data file. Suppose I read a file which have 1000 columns
> by:
>
> test = read.table("dat.txt", header=F)
>
> I know I could use "cbind()". I
First 4 colums of builtin data frame anscombe as a matrix:
as.matrix(anscombe[1:4])
On 8/19/06, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> Dear all,
>
> I have a question about how to get a matrix by combining a large number of
> columns from a data file. Suppose I read a file which have 1
Dear all,
I have a question about how to get a matrix by combining a large number of
columns from a data file. Suppose I read a file which have 1000 columns
by:
test = read.table("dat.txt", header=F)
I know I could use "cbind()". It's easy to do when the number of columns
is small (i.e. cbind(te
ailto:[EMAIL PROTECTED] On Behalf Of r user
> Sent: Tuesday, August 15, 2006 7:01 AM
> To: rhelp
> Subject: [R] question re: "summarry.lm" and NA values
>
> Is there a way to get the following code to include
> NA values where the coefficients are "NA"?
>
&g
On Tue, 15 Aug 2006, Petr Pikal wrote:
> Hi
>
> On 15 Aug 2006 at 7:01, r user wrote:
>
> Date sent:Tue, 15 Aug 2006 07:01:13 -0700 (PDT)
> From: r user <[EMAIL PROTECTED]>
> To: rhelp
> Subject: [R] que
user <[EMAIL PROTECTED]>, wrote:
From: Petr Pikal <[EMAIL PROTECTED]>
To: r user <[EMAIL PROTECTED]>, rhelp
Subject: Re: [R] question re: "summarry.lm" and NA values
Date sent: Tue, 15 Aug 2006 17:15:01 +0
Hi
On 15 Aug 2006 at 7:01, r user wrote:
Date sent: Tue, 15 Aug 2006 07:01:13 -0700 (PDT)
From: r user <[EMAIL PROTECTED]>
To: rhelp
Subject: [R] question re: "summarry.lm" and NA values
> Is there a way t
Is there a way to get the following code to include
NA values where the coefficients are NA?
((summary(reg))$coefficients)
explanation:
Using a loop, I am running regressions on several
subsets of data1.
reg <- ( lm(lm(data1[,1] ~., data1[,2:l])) )
My regression has 10 independent vari
Dear list members,
sorry for my incompleteness!
My problem is the following:
(R 2.3.1 on Mac OS X 10.4.7 RAM 1GByte using Mac GUI)
I have a function like this:
foo1 <- function()
{
out <- NULL
for(i in 1:10010) out <- paste(out, i, ". line\n", sep="")
return(out)
> "HansJB" == Hans-Joerg Bibiko <[EMAIL PROTECTED]>
> on Mon, 14 Aug 2006 17:56:30 +0200 writes:
HansJB> Hi, I have a tiny question concerning
HansJB> .Options$max.print
HansJB> I have to set up this value to a greater value than
HansJB> 1 because I want to concat
Hi,
I have a tiny question concerning .Options$max.print
I have to set up this value to a greater value than 1 because I
want to concatenate my output of a function to one single string (for
connivence).
I did this via .Options$max.print <- 64000 or options(max.print=64000)
Then I cal
Hello,
I am trying to predict the classes of a test data set after training an
rpart tree.
When I run:
predict(rpart_object_based_on_training_data, newdata = "testdata", type
= "class", na.action = na.pass)
I get an error message saying that a variable that is present in both
training and test
Hi Jessica,
I think you want to use names()
See:
?names
?colnames (or rownames)
If not, it's perhaps a good idea to send us a more complete example.
JeeBee.
On Fri, 04 Aug 2006 17:12:24 +0200, Jessica G. wrote:
> How to convert a rowname vector of numbers into a real column of the matrix,
>
>
Hi Mr Plate,
I have a little question
How to convert a rowname vector of numbers into a real column of the matrix,
My problem is that I applied a rowsum function on a matrix.
Then I get a matrix in which the names of the columns are the values of the
group (numbers)
Now I need to make calculati
Hi,
I am facing a problem in extrapolation of data series. It is a series of
Bond yields, I am having the yield for 1 year to 30 years. I want to find
the yield for 0.5 year and 30.5 years. I used the Langrange's Extrapolation
but the extrapolation deviates from the normal trend ( as we can see in
On Thu, 3 Aug 2006, Thomas Lumley wrote:
> On Thu, 3 Aug 2006, Benn Fine wrote:
>
> > I have ported some R code to C to make it faster.
> >
> > I can perform .Call("foobar",) once and it
> > works fine. Absolutely correct answer.
> >
> > If I put a loop inside foobar and run the main code
>
On Thu, 3 Aug 2006, Benn Fine wrote:
> I have ported some R code to C to make it faster.
>
> I can perform .Call("foobar",) once and it
> works fine. Absolutely correct answer.
>
> If I put a loop inside foobar and run the main code
> routine more than 100 times, it crashes R.
>
> Or if I cal
I have ported some R code to C to make it faster.
I can perform .Call("foobar",) once and it
works fine. Absolutely correct answer.
If I put a loop inside foobar and run the main code
routine more than 100 times, it crashes R.
Or if I call .Call("foobar") seperately more than
two tims
Dear all,
I am using the PLS package for PLSR analysis. And I have a basic question
about the standardize procedure, which I feel the PLS manual does not
explain clearly. I am hoping that I could get some help from the list.
>From the example in the "Standardization of Data Matrices" section, I c
Dear all,
I am using different multiple regression models (OLS and principal
component regression (PCR)) to make prediction of my test set. And those
models come from the same training set, except that the number of
variables or descriptors (columns of X) used in OLS is less than those
used in PCR
Dear all,
I am using different multiple regression models (OLS and principal
component regression (PCR)) to make prediction of my test set. And those
models come from the same training set, except that the number of
variables or descriptors (columns of X) used in OLS is less than those
used in PCR
Hello, I do this in such way:
DATAX <- matrix(seq(1,6,1),2,3)
DATAY <- matrix (seq(1,4,1),2,2)
rownames(DATAX) <- c("s1","s2")
rownames(DATAY) <- c("s1","s2")
colnames (DATAX) <- c("v1","v2","v3")
colnames (DATAY) <- c("respone_1","response_2")
KAL <- data.frame (N = rownames(DATAX))
KAL$Y <- DATA
Dear all,
I am trying to apply Ridge regression to my dataset, and then I would like
to predict the Y responses using the Ridge model (of certain lambda) for
new data point. The only Ridge regression functions I found is in "MASS"
library. However, there are very few functions available: lm.ridge(
On 7/29/06, Nantachai Kantanantha <[EMAIL PROTECTED]> wrote:
> Hi everyone,
>
> I would like to ask a question regarding to the data used to fit the mixed
> model.
>
> I wonder that, for the response variable data used to fit the mixed model
> (either via "spm" or "lme"), we must have several obser
Try:
?sensory
str(sensory)
dput(sensory)
lapply(sensory, class)
lapply(sensory, dim)
to see what it looks like inside. Seems that sensory is a data frame
consisting of two columns each of which is a matrix except that each
has a class of "AsIs". Thus try this (where I(...) creates objects of
cl
Dear all,
I am trying to my dataframe for the PLS analysis using the PLS package.
However I have some trouble generating the correct dataframe. The main
problem is how to use one name to represent several columns in the
dataframe.
The example dataframe in PLS package is called "sensory". I cannot
data structure. Or, you can do something like
head(star) to see the 1st 6 rows and see how the data are layed out.
I hope this helps,
Harold
-Original Message-
From: [EMAIL PROTECTED] on behalf of Nantachai Kantanantha
Sent: Sun 7/30/2006 12:33 AM
To: r-help@stat.math.ethz.ch
Subject: [R
Hi everyone,
I would like to ask a question regarding to the data used to fit the mixed
model.
I wonder that, for the response variable data used to fit the mixed model
(either via "spm" or "lme"), we must have several observations per subject
(i.e. Yij, i = 1,..,M, j = 1,.., ni) or it can b
Try to insert a space before the last bracket! R
sometimes does not accept a curly bracket as a first
character on a line of a script.
Regards,
Valentin
--- Ivan Kalafatic <[EMAIL PROTECTED]> wrote:
> I tried to make the following function:
>
> function(x, y){
>dates<-intersect(x[,1],y[,1]
Ivan Kalafatic wrote:
> I tried to make the following function:
>
> function(x, y){
>dates<-intersect(x[,1],y[,1])
>m<-matrix(NA,length(dates),3)
>m[,1]<-dates
>j<-1
>k<-1
>for(i in fdax[,1]){
> if(is.element(i,dates)){
> m[j,3]<-as.numeric(fdax[k,2])
>
"Ivan Kalafatic" <[EMAIL PROTECTED]> writes:
> I tried to make the following function:
>
> function(x, y){
>dates<-intersect(x[,1],y[,1])
>m<-matrix(NA,length(dates),3)
>m[,1]<-dates
>j<-1
>k<-1
>for(i in fdax[,1]){
> if(is.element(i,dates)){
> m[j,3]<-as.nu
I tried to make the following function:
function(x, y){
dates<-intersect(x[,1],y[,1])
m<-matrix(NA,length(dates),3)
m[,1]<-dates
j<-1
k<-1
for(i in fdax[,1]){
if(is.element(i,dates)){
m[j,3]<-as.numeric(fdax[k,2])
j<-j+1
}
k<-k+1
}
return
On Tue, 18 Jul 2006, Jim Frederick wrote:
> Hi,
>
> When I used
> x <- c(10,10,10,10,5,5,5,5)
> sum(x)
> x <- x - .Last.value/8
> the result was zero, as I expected.
> However, when I used
> x <- rnorm(101,0,10)
> sum(x)
> x <- x - .Last.value/101
> sum(x)
> I did not get zero, but -2.664535e-15.
Hi,
When I used
x <- c(10,10,10,10,5,5,5,5)
sum(x)
x <- x - .Last.value/8
the result was zero, as I expected.
However, when I used
x <- rnorm(101,0,10)
sum(x)
x <- x - .Last.value/101
sum(x)
I did not get zero, but -2.664535e-15. OK, that's fairly close to
zero, but I tried to improve the figur
On Thu, 2006-07-13 at 15:02 +0100, Pavlidis, Efthymios wrote:
> Hello,
>
> I am having the following silly problem with lm.
>
> Let X be a dataframe with X[,1] the dependent variable and X[,-1] the
> independent variables. I want to run the following
> but without including an intercept.
>
>
Hello,
I am having the following silly problem with lm.
Let X be a dataframe with X[,1] the dependent variable and X[,-1] the
independent variables. I want to run the following
but without including an intercept.
for(i in 1:100 ){
lm( X[,100-i] ) # this works fine but it returns an inter
Dear all,
I am trying to do a standard mutivariate linear regression for my training
set, and want to perform leave-one-out cross-validation to see how well
the model is. I found several CV in R, such as cv.lm in DAAG, cv.glm in
boot, and crossval in bootstrap. Which one should I use? And how can
Than you, Gavin. You helped me out a lot of problems.
Thank you very much!
Wei-Wei
2006/7/12, Gavin Simpson <[EMAIL PROTECTED]>:
> On Wed, 2006-07-12 at 00:51 +0800, Guo Wei-Wei wrote:
> > Thank you, Gavin. I think that might be what I need. But I'm a little
> > bit wandering what's the scale
On Wed, 2006-07-12 at 00:51 +0800, Guo Wei-Wei wrote:
> Thank you, Gavin. I think that might be what I need. But I'm a little
> bit wandering what's the scale of resid(mod). Is it
> scale(dist)/scale(speed), for example kilometer / (kilometer per
> hour)? or something else?
>
> Thank you very much
Thank you, Gavin. I think that might be what I need. But I'm a little
bit wandering what's the scale of resid(mod). Is it
scale(dist)/scale(speed), for example kilometer / (kilometer per
hour)? or something else?
Thank you very much!
Wei-Wei
2006/7/12, Gavin Simpson <[EMAIL PROTECTED]>:
> On Tue
On Tue, 2006-07-11 at 23:51 +0800, Guo Wei-Wei wrote:
> Dear all,
>
> I don't know what's my question is called. I have a performance
> variable A, such as sales. And I have another variable B, let's say
> establish time of firm. I want to create the third variable that is
> sales without the effe
Dear all,
I don't know what's my question is called. I have a performance
variable A, such as sales. And I have another variable B, let's say
establish time of firm. I want to create the third variable that is
sales without the effect of establish time. Maybe it can be called
partial effect proble
I need a help, I want to use p with mysql, i've installed mysql and the package
DBI but when i'm tapping :library(DBI) and there is no probleme, then i tap :
MySQL(max.con = 16, fetch.default.rec = 5000, force.reload = F) and there is
a message erreur ;
Erreur dans library(MySql) : aucun pack
Hi
[EMAIL PROTECTED] wrote:
> i created a postscipt file in R and then i downloaded a free version
> of ghostview to view it. unfortunately, i get the message
>
> fata error : dynamic memory exhausted
> when i try to view it.
>
> when i do a dir on windows xp, the file size is 149,034,475
> and
i created a postscipt file in R and then i downloaded a free version
of ghostview to view it. unfortunately, i get the message
fata error : dynamic memory exhausted
when i try to view it.
when i do a dir on windows xp, the file size is 149,034,475
and i know there about 17,000 graphs. is there
a
On Thu, 29 Jun 2006, Manoj wrote:
> Hello All,
> I am trying to work on writing the following piece of (pseudo)
> code in an optimal fashion:
>
>
> # Two data frames with some data
>
> a = data.frame(somedata)
> b = data.frame(somedata)
>
>
Hello All,
I am trying to work on writing the following piece of (pseudo)
code in an optimal fashion:
# Two data frames with some data
a = data.frame(somedata)
b = data.frame(somedata)
for(i in 1:nrow(dt) {
# Merge dates for a given da
Hello,
I have just completed my analysis of microarray data using limma. But, when I
ask for a list of topTable genes, I am getting the following error:
topTable(fit2,coef=5,adjust="fdr")
Error in array(x, c(length(x), 1), if (!is.null(names(x))) list(names(x), :
attempt to set a
Thanks a lot Spencer for your tips - I'll look into all of them.
Rainer
Spencer Graves wrote:
> It's difficult to say much at this level of generality, but I have
> four suggestions:
>
> 1. Have you tried creating a reasonable grid of starting values
> using "expand.grid" and then p
It's difficult to say much at this level of generality, but I have
four suggestions:
1. Have you tried creating a reasonable grid of starting values
using "expand.grid" and then plotting the resulting likelihood surface?
If you have more than 2 parameters, you may want to
Hi
I hope this is the right forum - if not, point me please to a better one.
I am using R 2.3.0 on Linux, SuSE 10.
I have a question concerning mle (method="BFGS").
I have a few models which I am fitting to existing data points. I
realised, that the likelihood is quite sensitive to the start v
Have you tried Bioconductor (www.bioconductor.org)? I haven't seen
any replies to your question, and I know next to nothing about
microarray data. I understand that Bioconductor specializes in
microarray and related data, has a listserve, etc.
hope this helps,
Sp
Hi,
I have a microarray data which has about 1 genes and 15 measurements. The
15 measurements contain 6 control, 3 low dose, 3 medium dose, 3 high dose. If I
treat the dose level as a continuous variable, how can I use ANCOVA to do gene
selection? (I mean to select those genes that respons
Dear Debayan,
I only use Genepix files but I will try to help you with the experience I
have concerning layout troubleshooting and limma.
First : RG$R, RG$Rb, RG$G and RG$Gb seem strange to me. As I understand
these data, it has to be a vector of the same length than the number of
spots.
It seems
Hi Everyone,
I have been trying to read some Arrayvision files( 2 channel cDNA) and
am having some
problem. My code is :
setwd('C:/work/data/limma/ndd1');
files <- c('ndd1_1.txt','ndd1_2.txt','ndd1_3.txt');
RG=read.maimages(files,"arrayvision",sep="\t");
#Normalisation
MA=normalizeWithinAr
ROTECTED] On Behalf Of singyee ling
Sent: Wednesday, May 17, 2006 8:01 AM
To: r-help@stat.math.ethz.ch
Subject:[R] question about survSplit
Dear R-users,
I use the survsplit function in the survival package to change my data into
counting-process format
and the transformed form
Dear R-users,
I use the survsplit function in the survival package to change my data into
counting-process format
and the transformed format is as follow:
(a)
start stop event DP age
0 5 01 20
5 1001 20
1025 11
om: MARK LEEDS
> > To: [EMAIL PROTECTED]
> > Sent: Sunday, May 14, 2006 9:49 PM
> > Subject: simple .bat question but i can't find out where to post. my
> > apologiesfor posting a non R question
> >
> >
> > I have what I think is a simple .bat questi
4/06, MARK LEEDS <[EMAIL PROTECTED]> wrote:
>
> - Original Message -
> From: MARK LEEDS
> To: [EMAIL PROTECTED]
> Sent: Sunday, May 14, 2006 9:49 PM
> Subject: simple .bat question but i can't find out where to post. my
> apologiesfor posting a non R question
&
- Original Message -
From: MARK LEEDS
To: [EMAIL PROTECTED]
Sent: Sunday, May 14, 2006 9:49 PM
Subject: simple .bat question but i can't find out where to post. my
apologiesfor posting a non R question
I have what I think is a simple .bat question but i can't figure ou
Hello
I am looking for informations about the R/Java package provided by the
Omegahat project.
I've already tried to contact the Omega project,but I'm steel without any
answer.
So, I wonder if somebody has already use this R/Java interface on
Windows.
In this case, I would appreciatte any help. H
Hi all.
I have a doubt with weighted linear regression. I've noted that supplying
integer weights gives different residuals df than giving 'double' weights.
Here an example:
x <- 1:5
y <- -1 + 2*x + rnorm(length(x))*0.1
y <- c(y, x + rnorm(length(x))*0.1)
dat <- data.frame(x=rep(x,2), y = y)
#int
zhongmiao wang gmail.com> writes:
> I am studying the effect of schools on student achievement growth over
> time. School effect is random. The effects of schools in prior years
> are assumed to be persistent till the current year. Thus, the total
> school effect in the second year is like J=J2+J
Dear Adrian and Ian (and r-helpers),
To answer my own question: The difference is simply due to rounding; that
is, there are actually small differences in the R^2 values that account for
the differences in the BICs.
Sorry for the false alarm,
John
> -Original Message-
> From: John Fox
Hello:
I am studying the effect of schools on student achievement growth over
time. School effect is random. The effects of schools in prior years
are assumed to be persistent till the current year. Thus, the total
school effect in the second year is like J=J2+J1. Does anyone know how
to model this
Dear Adrian and Ian (and r-helpers),
I encountered a curious result in developing an example using the bicreg
function in the BMA package: I noticed that pairs of models with equal R^2
and equal numbers of predictors had nevertheless different BIC values.
Looking at the bicreg function, the defini
See argument ``fixed'' in help(arima).
cheers,
Rolf Turner
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
How to force some AR coefficients to be ZERO and only allow the function
"ar" in R to change the other coefficients?
Thanks a lot!
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listin
Hello everyone,
I would highly appreciate your assistance with the following:
I would like to add a line to the coherency plot of a multivariate time
series, denoting the critical value above which the null hypothesis of
having no coherency between the two series may be rejected.
Is there a mecha
Andris Jankevics wrote:
> But I have a sligthy different results with my real data. Where can the
> problem be?
I think you have to supply some details for anyone to be able to
answer. At least what you did (the code), what you got (the results)
and what you expected to get.
--
Bjørn-Helge Me
I am sorry. It was my fault. My example is wrong. I need also rearrange a
validation data set too.
But I have a sligthy different results with my real data. Where can the
problem be?
Andris Jankevics
On Otrdiena, 18. Aprīlis 2006 17:55, Andris Jankevics wrote:
> Hello useRs,
>
> I am new user t
Hello useRs,
I am new user to R and also statistics. Why predicted results in this example
are different? Is the order of variables in X matrix important?
library (pls)
set.seed (1)
Y1 <- c(1,2,3,4,5,6,7,8,9,10)
Y2 <- c(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0)
X1 <- rnorm(10,sd=0.2)
X2 <- rnorm(
Hi useRs,
I have been running a regression of the following kind:
> summary(lm(dx[2:2747] ~ 0 + (dx[1:2746]>15)))
Call:
lm(formula = dx[2:2747] ~ 0 + (dx[1:2746] > 15))
Residuals:
Min1QMedian3Q Max
-46.35871 -3.15871 0.04129 3.04129 30.04129
Coefficient
On Wed, 12 Apr 2006, James Anderson wrote:
> Hi,
> Just a silly question about R source package:
> I downloaded some source package of R and found there are some source
> codes written in C, is this true for all the R packages? but it seems
> that not all the functions described in the package
The minimum you need is to have a single file, DESCRIPTION, plus
R and man directories/folders which contain your .R and .Rd
files respectively. A package need not use C and I suspect most
don't.
On 4/12/06, James Anderson <[EMAIL PROTECTED]> wrote:
> Hi,
> Just a silly question about R source p
Hi,
Just a silly question about R source package:
I downloaded some source package of R and found there are some source codes
written in C, is this true for all the R packages? but it seems that not all
the functions described in the package are contained in the source package
written in C.
Thank you both very much.
On 3/30/06, Spencer Graves <[EMAIL PROTECTED]> wrote:
>
> Dear Prof. Ripley:
>
> Thanks for the clarification.
>
> spencer graves
>
> Prof Brian Ripley wrote:
>
> > I don't think any of those cover such models, which are not GLMs.
> >
> > One can of co
Dear Prof. Ripley:
Thanks for the clarification.
spencer graves
Prof Brian Ripley wrote:
> I don't think any of those cover such models, which are not GLMs.
>
> One can of course run anything in R:
>
> Evelyn Hall: I would like to know how (if) I can extract some of the
>
I don't think any of those cover such models, which are not GLMs.
One can of course run anything in R:
Evelyn Hall: I would like to know how (if) I can extract some of the
information from the summary of my nlme.
Simon Blomberg: This is R. There is no if. Only how.
-- Evelyn Hall and Simon `Y
Have you considered "lmer" with the mle4 and Matrix packages and
"glmmPQL" in library(MASS)? If yes, PLEASE do read the posting guide!
"www.R-project.org/posting-guide.html", then gives us the simplest
example you can devise of what you want, explaining what you tried and
the obstacl
Andris Jankevics wrote:
> Hello useRs,
>
> I can't figure out how can I store a data frame of values of X and D from
> this loop:
>
> Z <- c(1:10)
>
> for (i in 1:(length(Z)-2))
> {D <- x[1:(2+i)];
> X <- paste("x",sep="",i)
> print (X)
> print (D)
> }
Assign
301 - 400 of 977 matches
Mail list logo