Ben Bolker ufl.edu> writes:
> Perhaps read.fwf is what you want? (fwf stands for
> "fixed width format"). You would have to work out the
> field widths, but it would seem to be pretty straightforward).
> [A little bit of extra verbiage will make gmane happier ...]
>
> Ben Bolker
Tru
I am using R 2.6.2. under windows
On Wed, Feb 27, 2008 at 12:43 PM, <[EMAIL PROTECTED]> wrote:
> >From: Arun Kumar Saha <[EMAIL PROTECTED]>
> >Date: 2008/02/27 Wed AM 01:03:26 CST
> >To: "[EMAIL PROTECTED]" <[EMAIL PROTECTED]>
> >Subject: [R] Loading user defined functions autometically each time
Hi, new to R and ggplot2. I've been trying to get a facet plot in which the
order of the facets is as I require, rather than ordered numerically,
alphabetically, by Roman numerals, mean (answers to these were posted here
after much searching). Here's some test code to demonstrate what I get.
seri
On Tue, 26 Feb 2008, Andre Nathan wrote:
> I know about stem, but the data set has 1 million points, so it's not
> very useful here. I want to avoid binning just to have an idea about the
> shape of the distribution, before deciding how I'll bin it.
Ideas:
1) use a much smaller sample of the dat
Hi all,
I wrote some user defined function for my own. Now I want to get a mechanism
so that every time I start R, those function will automatically be loaded in
R without manually copying pasting. Can gurus here pls tell me how to do
that? Or I have to build my own packages bundled with those fun
On Tue, 26 Feb 2008, Judith Flores wrote:
> Hi there,
>
> I found one reference to add a reactangle behind a
> plot using plot(...,add=T), I tried this but didn't
> obtain the desired result.
>
> If a I have the following code:
>
> plot(x,y)
> rect(xleft, ybottom, xright,ytop,col='green)
>
> T
While we're at it, and since you're new...
Your example will be much easier for r-help folks to read if you do
it like this:
j <- 0.4
for(i in 1:20) {
j=j+0.1
cp[i] <- pnorm(-j*3)*10^6
ratio[i] <- j
}
table <- data.frame(ratio,cp)
table
But the loop is unnecessary. Try
j <- seq(0
Ultimately my aim is to get a plot of density faceted by 2 factors with a
horizontal boxplot overlaid on each density plot in the grid to indicate
summary stats. So I've been experimenting with creating boxplots and density
plots. Here's some representative data.
series = c('C2','C4','C8','C10','
Dear Sir/Madam:
I am newbie of R. I a currently using multdrc object to generate fitting curve
and IC50. My 384 well format raw data contains multi dose response curves. My
script goes through sets of data then produce curve and ic50.
Here is my sudo code:
For (plateid in platelist)
{
Input d
If the goal is to get a sense of the 'shape' of the overall distribution
of 'x', then why not use:
plot(density(x))
?
HTH,
Marc Schwartz
Peter Alspach wrote:
> Andre
>
> If I understand you correctly, you could try a barplot() on the result
> of table().
>
> HTH ..
>
> Peter Alspach
>
Has anyone computed a Backwards lasso with positive coefficients? I
would like the R commands or guidance on fitting such a model.
The backwards lasso and positively constrained lasso are both described
in Least Angle Regression, Efron etal 2002.
__
R
Peter Dunn wrote:
> Hi all
>
> I'm using Sweave and LaTeX, and love how they interact.
>
> But here's a different interaction I'm not sure how to achieve; I
> hope someone can help. I use a simple example, of course, to
> demonstrate.
>
> Suppose in my LaTeX document I have this:
>
> Here
Have you looked at the 'dlm' package? It has a vignette to help
you learn to use it. Also, I've heard that a book about that package is
scheduled to appear in the next few months.
I have looked at the Kalman functions in the 'stats' package but
have not found documentation that s
Just replot the points over it again:
points(x, y)
Also there is an example using lattice graphics in the Example
section of:
library(zoo)
?xyplot.zoo
On Tue, Feb 26, 2008 at 7:16 PM, Judith Flores <[EMAIL PROTECTED]> wrote:
> Hi there,
>
> I found one reference to add a reactangle behind a
>
Thibaut Jombart wrote:
> Duncan Murdoch wrote:
>
>>> Maybe a clue: we can use
>>> cat(readLines(as.character(?plot)),sep="\n")
>>>
>>> to display the help (here, of plot) directly to the screen. So we could
>>> use something like:
>>>
>>> <>=
>>> ?plot
>>> @
>>>
>>>
>>> <>=
Hi all
I'm using Sweave and LaTeX, and love how they interact.
But here's a different interaction I'm not sure how to achieve; I
hope someone can help. I use a simple example, of course, to
demonstrate.
Suppose in my LaTeX document I have this:
Here is a linear equation:
\beg
Andre
If I understand you correctly, you could try a barplot() on the result
of table().
HTH ..
Peter Alspach
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Andre Nathan
> Sent: Wednesday, 27 February 2008 1:34 p.m.
> To: roger koenker
> C
Hi,
I am trying to replace (lattice) standard xyplot with xYplot variant from
Hmisc package to be able to add error bars to my plots. However, this does
not work, e.g:
library(lattice)
d <- data.frame(
SKU=gl(3, 1, 21, labels=c("a", "b", "c")),
Weekday=gl(7, 3, 21),
QC
Hi,
When I install new packages from CRAN, I frequently
find that some packages were missing from the download
queue. For example, on one of my computer with R2.6.2,
I can not find package glmpath from the download
queue. On my other computer with R2.5.1, I could still
find that particular package
I know about stem, but the data set has 1 million points, so it's not
very useful here. I want to avoid binning just to have an idea about the
shape of the distribution, before deciding how I'll bin it.
Andre
On Tue, 2008-02-26 at 16:20 -0600, roger koenker wrote:
> take a look at
>
> ?ste
> "Valentin" == Valentin Bellassen <[EMAIL PROTECTED]> writes:
> Hello, I have a data frame with 3 vectors $x, $y, and
> $type. I would like to plot $x~$y and having different
> colors for the corresponding points, one for each level of
> $type. Would someone know how to do tha
Hi there,
I found one reference to add a reactangle behind a
plot using plot(...,add=T), I tried this but didn't
obtain the desired result.
If a I have the following code:
plot(x,y)
rect(xleft, ybottom, xright,ytop,col='green)
The rectangle appear on top of the plot.
Any help will be gre
Hi,
looking for a function or syntax to estimate the score test in logistic
regression for the null hypothesis b1=0 in the model logit(p)=b0+ b1*x
+b2*z. Data comes from the binomial distribution (n,p).
Thanks,
ben
__
R-help@r-project.org mailing lis
Hello,
I am attempting to fit a non-linear model (Von Bertalanffy growth model)
to fish length-at-age data with the purpose of determining if any of the
three parameters differ between male and female fish. I believe that I
can successfully accomplish this goal assuming an additive error
struc
Hi everyone!
I have an array containing the following fields for over hundred compounds:
cpd, activity, fixterm, energy1, energy2, energy3, ...
I want to run a multiple linear regression on all entries of an array.
Therefore I tried to do this with a for loop. (Maybe there is a direct way of
c
Hi,
I am newbie of R. I a currently using multdrc object to generate fitting
curve and IC50. My 384 well format raw data contains multi dose response
curves. My script goes through set of data then produce curve and ic50.
Here is my sudo code:
For (plateid in platelist)
{
Input data (
Thanks, Erik,
This is a partial copy of my code. I want to get rid of the left column of
integers and I want to control the number of decimal places after the
decimal point.
> j<-0.4
> for(i in 1:20){
+ j=j+0.1;cp[i]<-pnorm(-j*3)*10^6;ratio[i]<-j}
> table<-data.frame(ratio,cp)
> table
rati
Just a couple of notes.
* What you are attempting to fit here is a linear model, or, more precisely
13112 simple linear regressions. Why not just use 'lm'? 'glm' is for fitting
generalized linear models and using it for this special case is going to be
computationally costly. [NB the concept
On Tue, 26 Feb 2008, John Fox wrote:
> Dear John and Stephen,
>
> As I recall, importing from Excel, Access, or dBase files works properly
> only under Windows, and consequently the menu item appears only in Windows
> systems. I frankly don't recall why this is the case, since the RODBC
> package
It's hard to provide a useful answer without know what your data looks
like, but maybe something like:
library(ggplot2)
qplot(timep, mn, data=summ, geom="line", colour = tx)
Hadley
On Tue, Feb 26, 2008 at 4:24 PM, Judith Flores <[EMAIL PROTECTED]> wrote:
> Sorry, I just realized I didn't type in
On Tue, Feb 26, 2008 at 4:10 PM, Andre Nathan <[EMAIL PROTECTED]> wrote:
> Hello
>
> I need to plot a histogram, but insted of using bars, I'd like to plot
> the data points. I've been doing it like this so far:
>
> h <- hist(x, plot = F)
> plot(y = x$counts / sum(x$counts),
>x = x$br
Rob Tibshirani propose to use lasso with Cox
regression for variable selection in his 1997 paper
"The lasso method for variable selection in the Cox
model" published in Statistics In Medicine 16:385. I
understand the lars() function in lars package
implemented lasso, but it does not do lasso with C
Thanks very much.
--- John Fox <[EMAIL PROTECTED]> wrote:
> Dear John and Stephen,
>
> As I recall, importing from Excel, Access, or dBase
> files works properly
> only under Windows, and consequently the menu item
> appears only in Windows
> systems. I frankly don't recall why this is the
> cas
Dear John and Stephen,
As I recall, importing from Excel, Access, or dBase files works properly
only under Windows, and consequently the menu item appears only in Windows
systems. I frankly don't recall why this is the case, since the RODBC
package is used. The code for reading these kinds of file
Hello,
I tried the R-Sig-Mac list with this query, but had no takers. I hope
that the following isn't too far off the mark for this list. Many
thanks in advance if someone can help me out!
Quick summary: I can't get RSPerl working on a PPC G5 with the
pre-compiled binary for Mac OS X, but
I'm on Windows XP R 2.6.2
The command
Data > Import > From Excel, Access etc
Supply data.frame name
Select xls file and then select the sheet seem to work
for me
What is not showing?
--- stephen sefick <[EMAIL PROTECTED]> wrote:
> I am using R 2.6.2 on Mac OS X 10.4.9- I would like
> to impo
Sorry, I just realized I didn't type in the correct
names of the variables I am working with, this is how
it should be:
plot(1,1,type="n")
for (i in summ$tx) {
points(summ$timep[summ$tx==i],summ$mn[summ$tx==i])
lines(summ$timep[summ$tx==i],summ$mn[summ$tx==i])
}
Thank you,
Judith
On 27/02/2008, at 11:01 AM, John Kane wrote:
> Can you give a working example of what is happening
> and explain what is x?
>
> With a simple x vector of x <- rnorm(20, 5, 2)
> I don't get anything like what you seem to be getting.
>
> My code
> ===
take a look at
?stem
There is still a place for handtools in the age of integrated
circuits. Of course, avoiding binning isn't really desirable.
url:www.econ.uiuc.edu/~rogerRoger Koenker
email[EMAIL PROTECTED]Department of Economics
vox: 217-333-4558
Hello
I need to plot a histogram, but insted of using bars, I'd like to plot
the data points. I've been doing it like this so far:
h <- hist(x, plot = F)
plot(y = x$counts / sum(x$counts),
x = x$breaks[2:length(x$breaks)],
type = "p", log = "xy")
Sometimes I want to have a look
Can you give a working example of what is happening
and explain what is x?
With a simple x vector of x <- rnorm(20, 5, 2)
I don't get anything like what you seem to be getting.
My code
===
x <- rnorm(20, 5, 2)
table<-data.frame(x, scientific=F, digi
Sorry, it should have been:
col = c("red", "blue")[c(rep(1, 100), rep(2, 15))]
-Christos
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Christos Hatzis
> Sent: Tuesday, February 26, 2008 4:46 PM
> To: 'SNN'; r-help@r-project.org
> Subject: Re: [
I am not sure if I completely understand but would
something like this work? (not tested)
points( col=i)
lines ( col=i=1)
--- Judith Flores <[EMAIL PROTECTED]> wrote:
> Dear R-experts,
>
>I want to assign different colors to groups of
> points, lines and segments in a plot.
>
> I
Those are parameter to 'print'; what you want is something like:
> x <- data.frame(a=runif(10))
> print(x)
a
1 0.713705394
2 0.715496609
3 0.629578524
4 0.184360667
5 0.456639418
6 0.008667156
7 0.260985437
8 0.270915631
9 0.689128652
10 0.302484280
> print(x,scientific=F, d
If your samples are in the specified order (i.e. first 100 from group A and
remaining from group B) you can try the following in your plot call:
plot(..., col=c("red", "blue")[c(rep(100, 1), rep(15, 2))])
-Christos
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECT
Without knowing what your 'x' is, it's hard to see what is happening
that you don't like.
Your data.frame function creates a data.frame containing columns
scientific and digits, equal to FALSE and 4 for all rows, respectively.
Is that what you want?
cvandy wrote:
> Hi!
> I'm an R newbie and
Hi Martin
What do you think of the following code ?
A <- 1
imputeMyObj <- function(x){
xImp <- 3 #Calcul for Imputation of x
nam<-deparse(substitute(x))
assign(nam,xImp,envir=parent.frame())
}
imputeMyObj(A)
A
I guess it is not usual, but do you think it is 'bad' programming or
will have
?rbind
On 2/26/08, stephen sefick <[EMAIL PROTECTED]> wrote:
> I have not been able to find anything to do what I want, so I am going
> to tact to the left. I have twp continuous time series for two years
> with the same fourteen variables. I would like to simply append the
> second year to the
stephen sefick schreef:
> I have not been able to find anything to do what I want, so I am going
> to tact to the left. I have twp continuous time series for two years
> with the same fourteen variables. I would like to simply append the
> second year to the first. They both have the same column
I'm trying to learn how to use the lrm() function by simulating data using
an old dataset but it's giving me an error I don't understand:
(nst$regular already exists)
nst$regular<-as.ordered(nst$regular)
nst$age<-rnorm(n=942,mean=43.20488,sd=17.03)
nst$age<-round(age,digits=0)
regform<-regular~
Hi!
I'm an R newbie and this should be a trivial problem, but I can't make it
work and cannot find what I'm doing wrong in the literature.
I entered the the command:
table<-data.frame(x, scientific=F, digits=4)
table
This prints a column of x with 16 useless decimal places after the decimal
point.
Hi,
I have matrix of 300,000*115 (snps*individual). I ran the PCA on the
covariance matrix which has a dimention oof 115*115. I have the first 100
individuals from group A and the rest of 15 individuals from group B. I need
to plot the data in two and 3 dimentions with respect to PC1 and PC2 and
Dear R-experts,
I want to assign different colors to groups of
points, lines and segments in a plot.
I have a vector(dat$tx) that is going to have a
different number of levels every time the user wants
to plot a set of data. The number of lines is going to
be equal to the number of levels of
Camarda, Carlo Giovanni wrote:
> Dear R-users,
>
> is there any way to sum up the elements of the "diagonals" of a matrix
> without using a for-loop? While there is a simple way over rows and
> columns, I don't see a straightforward multiplication for the diagonals,
> am I too demanding? Or, more l
Thanks for the tip, Erik. I did check, but didn't see any $... elements to
choose from. For comparison, I checked another object class (lm) with the
str function, and it clearly showed $coefficients, $residuals, etc..
So it would seem that the usual approach doesn't work with adfstat
objects...
-
Try str(yourobject) to see the structure of the object. If there were a
value called, say, p.value, you could then do:
myvar <- yourobject$p.value
This is just general advice, I don't know the specifics of the package
you mention, but try this first.
Best,
Erik Iverson
Matt33 wrote:
> Hi,
>
Hi,
I'm using the ADF.test function in the uroot package to obtain an
adfstat-class object. I'm wondering how I can extract the values (test
statistic, p value, etc.) from this class, since it doesn't seem to have
usual values. I get the following summary, but I'm not sure how to do
anything with
Dear all,
thanks for your suggestions. I like the idea of including directly
the LaTeX file corresponding to the targeted topic, however, my
understanding from the reading of ?help is that these LaTeX files
are not always available, depending on the build of R.
I found a solution that works well
Does this do what you want?
> d <- col(mat) - row(mat)
> tapply( mat, d, sum )
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECT
Hello LIST,
In using Rulefit I;ve bee nabel to fit a model using rulefit without
incident. However, when trying to use intnull and interact things go
astray.
> rf=rulefit(x,"N", cat.vars="H",
> not.used=c("G","T"),huber=0.9,path.speed=1);
[snip]
RuleFit model 2/26/2008 2:17p
ave(abs(error))
Dear R-users,
is there any way to sum up the elements of the "diagonals" of a matrix
without using a for-loop? While there is a simple way over rows and
columns, I don't see a straightforward multiplication for the diagonals,
am I too demanding? Or, more likely, I'm lack some algebra trick? Is
the
See the description of the subset argument in ?xyplot
subset: logical or integer indexing vector (can be specified in terms
of variables in 'data'). Only these rows of 'data' will be
used for the plot. If 'subscripts' is 'TRUE', the subscripts
will provide ind
I have not been able to find anything to do what I want, so I am going
to tact to the left. I have twp continuous time series for two years
with the same fourteen variables. I would like to simply append the
second year to the first. They both have the same column headings
etc. Just like tappin
Try this:
xyplot(y ~ x , data = junk.frm, type = c("g", "p","smooth"), pch=20,
subset=z=="D" & y < 2)
On 26/02/2008, David Afshartous <[EMAIL PROTECTED]> wrote:
>
> All,
>
> I'm having problems w/ a simple attempt to subset an xyplot.
>
> The first plot below is a plot of y versus x for certa
All,
I'm having problems w/ a simple attempt to subset an xyplot.
The first plot below is a plot of y versus x for certain values of a third
categorical variable z. Now I'd like to further restrict this to certain
values of variable y. Neither of the two attempts below work. Any
suggestions m
On Mon, Feb 25, 2008 at 1:46 PM, Peter Dalgaard
<[EMAIL PROTECTED]> wrote:
> Ben Bolker wrote:
> > stian mail.rockefeller.edu mail.rockefeller.edu> writes:
> >
> >> Hi,I am wondering how to conduct Kenward-Roger correction in
> >> the linear mixed model using R. Any idea?
> >>
> >> Thanks a
Dennis Fisher plessthan.com> writes:
>
> Colleagues,
>
> I am trying to read a file written by Fortran. Several lines of the
> file are pasted below:
Perhaps read.fwf is what you want? (fwf stands for
"fixed width format"). You would have to work out the
field widths, but it would seem
Colleagues,
I am trying to read a file written by Fortran. Several lines of the
file are pasted below:
101 06e+050280.9777.484208.18147.20.246350.4020.7308717.882600 1 1
101 0.500280.9777.484208.18147.20.246350.4020.7308717.882600 0 0
101 100280.9777.484208.18147.20.246350.4020.730
On Tue, Feb 26, 2008 at 10:58 AM, Werner Wernersen
<[EMAIL PROTECTED]> wrote:
> Hello,
>
> I am very happy that I have Sweave and R to write my
> papers. But I still have to do some tables by hand
> since I have not found out how I can customize the
> latex tables produced by R further (I mainly us
Using eval-parse for this looks like overkill. You should just be able to
do something straightforward like:
for (i in 1:120)
assign(paste("book", i, sep=""), read.xls(paste("Book", i, ".xls",
sep=""), sheet=1, from=4, colClasses="numeric"))
which would put your spreadsheets in variables bo
Hi
My name is Vladimir Samaj. I am a student of Univerzity of Zilina. I am
trying to implement Kalman Filter into my school work. I have some problems
with understanding of R version of Kalman Filter in package stats( functions
KalmanLike, KalmanRun, KalmanSmooth,KalmanForecast).
1) Can you tell
As Thomas Lumley has frequently noted, eval(parse...)) should generally be
avoided.
You probably want do.call(), as in
## warning: UNTESTED
lapply(numbs,
function(n)do.call(read.xls,list(file=paste("Book",n,".xls",sep=""),sheet=1,
...))) ## ... contains the other explicitly names arguments
N
Michael -
>
> 3. Equally, advice about a better approach would also be v. welcome.
>
As an alternative approach, you may want to save the sheet names in a
list and then use the lapply function to read your data into a list.
This would greatly reduce the 'clutter' of objects in your global
> ### stripping off eval, looking for clues
>
>
>
> parse(paste('read.xls("Book',nam,'.xls",sheet=1,from=4,colClasses="numer
> ic")',sep=''))
>
> Error in file(file, "r") : unable to open connection
>
> In addition: Warning message:
>
> In file(file, "r") :
>
> cannot open fi
And also scree.plot psy package.
Bruno Falissard
INSERM U669, PSIGIAM
"Paris Sud Innovation Group in Adolescent Mental Health"
Maison de Solenn
97 Boulevard de Port Royal
75679 Paris cedex 14, France
tel : (+33) 6 81 82
R-helpers
I have 120 small Excel sheets to read and I am using
library(xlsReadWrite): one example below.
I had hoped to read sheets by looping over a list of numbers in their
name (eg Book1.xls, Book2.xls, etc).
I thought I had seen examples which used eval-parse-paste in this way.
Of course the second this showed up on R-help, I got rid of the formal
argument 'mydata' and just used ..., which seems to have solved the problem.
Erik
Erik Iverson wrote:
> Hello -
>
> I would like to create a wrapper to the 'save' function in the base
> package, but have a small problem wit
Hello,
I am trying to perform a glm analysis on a 68x13113 matrix (named
data.spect). The first column corresponds to the predictor
(data.spect[,1]) and the rest to the response variables
(data.spect[,2:13113]). When I try this code
glmObject <- glm(data.spect[,2:13113]~data.spect[,1])
I get t
Hello -
I would like to create a wrapper to the 'save' function in the base
package, but have a small problem with the name of the object that is
getting saved in the file. Below is a simple example illustrating my
problem.
## BEGIN SAMPLE R CODE
#
Hello,
I am trying to use the projection.matrix( ) function and am following the
example given. I have my data formatted very similar to the test.census
example.
> str(AsMi05mat)
`data.frame': 1854 obs. of 6 variables:
$ Tag : num 501 502 503 504 505 506 507 508 509 510 ...
Werner Wernersen yahoo.de> writes:
> I am very happy that I have Sweave and R to write my
> papers. But I still have to do some tables by hand
> since I have not found out how I can customize the
> latex tables produced by R further (I mainly use
> xtable()). Like for instance, I have a table whi
Hi,
i am trying to do a small job but it's bugging and i don't know what's
happening.
let's say i have four matrix, M01_a, M02_a, M01_b, M02_b that looks
like this
col1 col2
1 20
2 25
3 26
4 30
5 35
6 64
7 75
8 80
9 85
Hello,
I am very happy that I have Sweave and R to write my
papers. But I still have to do some tables by hand
since I have not found out how I can customize the
latex tables produced by R further (I mainly use
xtable()). Like for instance, I have a table which
needs an extra row every few rows as
Is there a way to get Stangle to print out Sexpr expressions as well as
code chunks?
I am using weaver and when trying to debug Sweave files, it is sometimes
difficult to figure out where the error is. The errors are hardest to
find when they are in Sexpr expressions because stepping through th
see fa.parallel in the psych package
At 7:27 AM -0800 2/26/08, Karen Douglas wrote:
>I am seeking information on whether anyone has written code to perform
>Horn's parallel analysis (a procedure that informs the selection of the
>proper number of components in PCA) in R.
>
>
>
>Thank you in advan
ade4 has the "testdim" function which implements the method proposed in:
*S. Dray*. On the number of principal components: A test of
dimensionality based on measurements of similarity between matrices.
/Computational Statistics and Data Analysis/, 52:2228-2237, 2008.
Cheers,
Karen Douglas wro
Dear Karen,
You should find the code you want in the package paran:
http://cran.r-project.org/web/packages/paran/index.html
And you may be interested in looking at the Multivariate Cran Task view which
attempts to catalogue these kinds of things:
http://cran.r-project.org/web/views/Multivariate.
I am seeking information on whether anyone has written code to perform
Horn's parallel analysis (a procedure that informs the selection of the
proper number of components in PCA) in R.
Thank you in advance for any help you can provide. Please respond
off-list at the email address below.
Karen
*hi
I aam facing following problem
my action executing twice when i try to load the tileslayout page as
actionforward success
*
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-h
Try this:
merge(tsus_same_mean, as.data.frame(tsuPIDCount), by.x="PID", by.y=0)
On 26/02/2008, Karin Lagesen <[EMAIL PROTECTED]> wrote:
>
> I have managed to create a data frame like this:
>
>
> > tsus_same_mean[1:10,]
> PIDGrpDist PercAlnPercId
> 1 12638 Acid
Dear All,
I have a question regarding predict.rpart. I use
rpart to build classification and regression trees and I deal with data with
relatively large number of input variables (predictors). For example, I build an
rpart model like this
rpartModel <- rpart(Y ~ X, method="class",
minsplit =1, m
I am using R 2.6.2 on Mac OS X 10.4.9- I would like to import excel
files into R with Rcmdr, but This option does not present itself in
the file menu. Any help would be greatly appreciated.
Stephen
--
Let's not spend our time and resources thinking about things that are
so little or so large t
I have managed to create a data frame like this:
> tsus_same_mean[1:10,]
PIDGrpDist PercAlnPercId
1 12638 Acidobacteria 0.0 1.000 1.000
2 87 Actinobacteria 0.0 0.970 0.970
3 92 Actinobacteria 0.008902000 1.000 0.991
Try this:
foo <- function(data, ...)
{
patt <- grep("var[0-9]$", names(data), value=T)
res <- paste(sapply(data[, patt], as.character), sapply(data[,
paste(patt, "lab", sep="_")], as.character), sep=": ")
cbind(a, matrix(res, ncol=length(patt), dimnames=list(NULL,
paste(patt, "new", sep="_"
}
Hi,
if the columns always follow the same order (indx,var#, var#_lab ...),
then You could use the column numbers and do the following:
1) CN<-colnames(df)[#] (#=colnum of 'var#')
2) df$NEW<-(here the expression to create the new variable)
3) colnames(df)[ncol(df)]<-c(paste(CN,"new",sep="_"))
Th
Dear All,
I have a question regarding predict.rpart. I use
rpart to build classification and regression trees and I deal with data with
relatively large number of input variables (predictors). For example, I build an
rpart model like this
rpartModel <- rpart(Y ~ X, method="class",
minsplit =
Patrick Giraudoux univ-fcomte.fr> writes:
>
> Dear listers,
>
> Here we have a strange result we can hardly cope with. We want to
> compare a null mixed model with a mixed model with one independent
> variable.
>
> > lmmedt1<-lme(mediane~1, random=~1|site, na.action=na.omit, data=bdd2)
> >
ian white a écrit :
Patrick,
The likelihoods of two models fitted using REML cannot be compared
unless the fixed effects are the same in the two models.
Many thanks for this reminder. Shame on me: it recalls me that this
subject may have been already largely discussed on this list. Now, I
You can use '\include' also:
\documentclass{article}
\begin{document}
\include{/usr/local/lib/R/library/graphics/latex/plot}
%file.path(system.file(package="graphics"), "latex", "plot")
\end{document}
I don't know if there is a way of input file.path in include by Sweave.
On 25/02/2008, Jean
Patrick,
The likelihoods of two models fitted using REML cannot be compared
unless the fixed effects are the same in the two models.
On Tue, 2008-02-26 at 14:38 +0100, Patrick Giraudoux wrote:
> Dear listers,
>
> Here we have a strange result we can hardly cope with. We want to
> compare a n
1 - 100 of 115 matches
Mail list logo
