 | | | rhelp | | Dear all
When fitting an "ols.model " the confidence interval at 95% doesn 't cover
the plotted data points because it is very narrow.
Does this mean that the model is 'overfitted ' or is there a speHi
I have a (numeric) matrix X of the type
> X <- matrix(c(1 9) ncol 3 nrow 3)
> colnames(X) <- c( "A " "B " "C ")
> rownames(X) <- c( "D " "E " "F ")
Also a have a data frame Y like
> Y <- as.data.Hi
I am trying to get the basis matrix and penalty matrix
for natural cubic splines. I found R function "ns " can
generate the B-spline basis matrix for a natural cubic
spline and "bsplinepen "
canHey
I am estiamting a partially linear model y X\beta+f(\theta) where the f(\theta) is estiamted using wavelets.
Has anyone heard of methods to test if the betas are significant or to ad
I 'm trying to compile R on AIX 5.3. I 've gotten "configure " to work but
"make " generates the following errors
Target "R " is up to date.
Target "R " is up to date.
Target "R " is up to date.
Target "Does anybody know what is relation between 'T ' value calculated by
'wilcox_test ' function (coin package) and more common 'W ' value?
Jarek Tuszynski
>library(coin)
> ### Tritiated Water Diffusion Ac
Dear r-list
I am trying to visually seperate the two main clusters of a dendrogram.
The idea is to use
'edgePar list(lty 3) ' for 'dend1[[1]] ' and
'edgePar list(lty 1) ' for 'dend1[[2]] '
I have noHas anybody implemented code to extract coefficients for a Butterworth
low-pass filter? I know Matlab has it implemented in the signal toolbox. I
want to make use of a 10 point Butterworth low-pass fiHI
Suppose I have the following data structure.
LRT tp
1 1.50654010 522
2 0.51793929 522
3 0.90340299 522
4 1.20293325 522
5 1.05578774 523
6 0.01617942 523
7 0.68183543 523In response to a thread where a statistical question
unrelated to R was asked Uwe Ligges wrote
> PLEASE do read the posting guide!
> http //www.R-project.org/posting-guide.html
>
> Please ask a lhy all
I wish to switch in a part in my code that use "read.table " to "scan " actually i use this
tab <-scan( "data.dat " what integer() skip 1)
dim(tab) <-c(75 length(tab)/75)
tab <-t(tab)
It gives mR-help
I would like to place nine (3X3) plots per page. I am not properly
implement mfrow(3 3) in the script below
jpeg( "xyplot.jpg ") #names output file
my_args <- commandArgs() #sets up to take I don 't care what you guys say -- that email was for me!! That money is
MINE baby!!!!!!
Hi all
not to invent the wheel a second time I wonder if someone out
there has programmend Snee 's CADEX algorithm (Validation
of regression models Methods and examples Technometrics 19
415-428).Dear statisticians
I would like to save results for a "for loop " in a vector previously created.
My result would be of class "matrix ".
I tried the following script
Script
#Creation of a previousHi all
I would like to test the relevance of a vector field (i.e. if the
vectors are organized or not).
To do so I would like to use an autocorrelation test so that I have
two questions
- is th-----BEGIN PGP SIGNED MESSAGE-----
Hash SHA1
Hello
I 'm having problems with this line of code
X.lm <- lapply(names(d) function(x) lm(d[ "cls "] ~ d[x] data d))
d[x] is what is giving trouble hOn Thu 15 Dec 2005 14 21 17 -0500 Strickland Matthew wrote
SM > Hi
SM >
SM > My data consists of a set of point locations (x y).
SM >
SM > I would like to know if there is a procedure for bivariate kHi All
Using glmmPQL you can extract the full table of estimates SE p-values
etc using as an example
mymodel <-glmmPQL(mymodel here)
summary(mymodel)[[18]]
How can I pull this table out of a lHi
I have a question on write.table
I have a data.frame called t7 as below
> dim(t7)
[1] 14015184 6
> t7[1 5 ]
uci uce par line graphical.forms stems
1 0 0 0 0 active You can try `locfit ' though it does local likelihood rather than
garden-variety kernel density estimation. Here 's an example
library(locfit)
data(cldem)
den.fit <- locfit(~ x1 + x2 data cltrainHi
My data consists of a set of point locations (x y).
I would like to know if there is a procedure for bivariate kernel
density estimation in R that returns the density estimates at the
observed Hi
?qt says
ncp non-centrality parameter delta currently for 'pt() ' and
'dt() ' only for 'ncp < 37.62 '.
'pt ' doesn 't complain when I give ncp > 37.62 (pnt.c seems to have a
work
Dear all
I have a 10 columns matrix which has 2^10 1024 unique rows with values 0 and
1.
What I would like to do is a little complicated in a simple statement for a
subset (say 1000 rows) I waDear Friends
I 'm trying to learn to use parallel computation using
snow & rpvm
Now I have two boxes
1) Pentium 4 pvm 3.4.5 R 2.2.0
with FreeBSD 5.4 (box uffbsd)
2) PowerPC pvm 3.4.5 R 2.1The name conflicts in Epi and ROC packages (2 'ROC ' functions are the
problem) cause the following code
to work once but not twice
library(MASS) data(cats)
x cats[ 2]
y ifelse(cats[ 1]Hello there
I would like to simulate X --Normal (20 5)
Y-- Normal (40 10)
and the correlation between X and Y is 0.6. How do I do it in R?
Thank you very much
Lisa WangHi all
I would like to test the uniformity (or no uniformity) of a vector field.
To do so I would like to use a ks.test (Kolmogorov-Smirnov) but don 't
know what value to give for the second parameHello
I need a little initial help with the MNP package.
I am trying to analyze a discrete choice experiment with the following
settings
- choice sets consist of 3 alternatives
- choice alternati > library(foreign)
> da <-read.dta(file.choose())
> da
startdat starttim enddate endtime days hoursmin secused
1 2005-01-11 2 2005-12-15 20.19 NA NA 9
> attributes(da)
$d |
|
|
