On Tue, 13 Jan 2009, Stefano Leonardi wrote:
Thanks for the answers.
Still I am not totally convinced about the interpretation of intercept as a
mean of fitted values for group belonging to first level of each factor (those
having 0 in all columuns in matrix.models, except the first column) be
There are real examples; they are all fairly obscure. It can't be a big
problem because the standard formal argument name for a data frame in modelling
and graphics functions is 'data'. That's actually a more serious problem than
the function called data() -- having local and global variable
The original Olsson's paper
(http://www.citeulike.org/user/ctacmo/article/553309) did mention that
the greatest biases and numeric problems were encountered when the two
variables had opposite skewness. Your example is even more extreme:
tetrachoric and polychoric correlations do not like zero coun
I also have this problem! I would be great to hear a solution.
Mike Williamson-9 wrote:
>
> I am creating a CDF plot function more user-friendly than any default
> r
> function. Depending upon the bimodality of the data (it is often
> bimodal),
> or any other strange data trends, the poi
Since you have a time series of times this may not be the best application
for zoo. Read R News 4/1 and try this:
Lines <- "Spec_code DateTime
Eptfus 6/4/03 22:44
Eptfus 6/4/03 22:44
Eptfus 6/9/03 21:59
Eptfus 6/10/03 3:24
Eptfus 6/13/03 0
Thanks Jim,
What happens if I want one solid line and four custom lines
("92","11","12","62).
I tried c("10","92","11","12","62) and it said that you can't use 0's. I
also tried:
ownlines<-c("92","11","12","62)
lty=c(1,ownlines)
which resulted in the same thing.
jholtman wrote:
>
> You c
Hi all,
As I understand the zoo package will aggregate dates and times, but
unclear how to tackle this problem.
I need to accomplish the following:
1. Generate a scatter plot of bat activity with dates on the X-axis
and time on the Y-axis.
include sunset and sunrise curves as an added part
On Tue, 13 Jan 2009, joseph.g.bo...@gsk.com wrote:
The first program generates an error message and does not execute the
regression of y on x.
x<-1:10;
y<-rnorm(10) + x;
prac <- function( model, wghts ){ lm(model, weights = wghts) }
prac(model = y~x, wghts = rep(1
I would like to be able to extract the significance levels (or standard
errors) of the variance-covariance parameters in nlme objects (i.e. from
the varClasses and/or corClasses). Is there a way to do it, or do I have
to reverse the intervals function calculations ?
Thanks.
Rob
See ?emptyspace in the plotrix package.
On Tue, Jan 13, 2009 at 5:02 PM, John Fort wrote:
> I have written code which puts text to the upper right of points on a plot.
> Some of my point are very close together and the new text is writing over
> the previous text.
>
> Has anyone written a clev
David Winsemius wrote:
> help.search did not offer any clues. Using Baron's search page, I got
> reference to r-help entries from 2004, but looking at the current
> documentation , I see no such function. I am wondering if the name of the
> function was changed in later versions of the package
I have written code which puts text to the upper right of points on
a plot. Some of my point are very close together and the new text is
writing over the previous text.
Has anyone written a clever "smart-text" function that writes the
text in a blank spot near the point (perhaps on the
Dear Dorothee,
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Dorothee
> Sent: January-13-09 8:49 PM
> To: r-help@r-project.org
> Subject: Re: [R] polychoric correlation: issue with coefficient sign
>
>
> Thank you so much f
> Hello all,
> I try the follow test with an odf file.
>
> *INPUT: teste.odt*
> >=
> x <- matrix(rnorm(6), nc=2)
> xst <- tableStyles(x, useRowNames = TRUE, header = NULL)
> odfTable(x, styles=xst, colnames = NULL, useRowNames = TRUE)
> @
Here is one issue: you are asking for the row names to be
Dear List:
After reading some of the e-mails related to this topic I found one that
converts a ftable to latex nicely.
library(Hmisc)
n=500
sex <- factor(sample(c("m","f"), n, rep=TRUE))
treatment <- factor(sample(c("Drug","Placebo"), n, rep=TRUE))
symptom <- factor(sample( c('H','S','G'), n
Thank you so much for all your answers! And sorry for being scarce on the
details.
My dataset has 12 variables (6 ordinal coded from 1 to 5, and 6 binary) and
384 cases without missing value. High values mean 'positive' attitude toward
the object of study.
I probably went too fast in my earlier i
You can create your own dashed lines:
lty
The line type. Line types can either be specified as an integer
(0=blank, 1=solid (default), 2=dashed, 3=dotted, 4=dotdash,
5=longdash, 6=twodash) or as one of the character strings "blank",
"solid", "dashed", "dotted", "dotdash", "longdash", or "twodash",
I am creating a CDF plot function more user-friendly than any default r
function. Depending upon the bimodality of the data (it is often bimodal),
or any other strange data trends, the points can end up gathering in just
about any corner of the plot. So, when I add a legend, whether I choose
Dear R-Users
I have 2 questions:
Firstly, If I create a matplot and legend for multiple vectors and then tag
another vector on using matlines (e.g. a 'total' of all vectors), is there
anyway to add the new line to the legend without recreating it? I have
created the plot this way because I w
The first program generates an error message and does not execute the
regression of y on x.
x<-1:10;
y<-rnorm(10) + x;
prac <- function( model, wghts ){ lm(model, weights = wghts) }
prac(model = y~x, wghts = rep(1, 10))
But the next program works:
x<
David Winsemius wrote:
On Jan 13, 2009, at 5:54 PM, Stefano Leonardi wrote:
David Winsemius wrote:
It would be so *if* you had estimated a saturated model, but you did
not. Your fit
of Y|A values is adjusted for the B values.
Thank you.
What do you exactly mean by "Y|A values is adjusted fo
On Jan 13, 2009, at 5:54 PM, Stefano Leonardi wrote:
David Winsemius wrote:
It would be so *if* you had estimated a saturated model, but you
did not. Your fit
of Y|A values is adjusted for the B values.
Thank you.
What do you exactly mean by "Y|A values is adjusted for the B values"?
Anot
David Winsemius wrote:
It would be so *if* you had estimated a saturated model, but you did
not. Your fit
of Y|A values is adjusted for the B values.
Thank you.
What do you exactly mean by "Y|A values is adjusted for the B values"?
And Yes. It works! Using the interaction I got what I expe
On Jan 13, 2009, at 5:15 PM, Stefano Leonardi wrote:
Peter Dalgaard wrote:
Actually, notice that you are averaging identical values, so the
"mean"
in the tapply is slightly misleading.
Notice also that the intercept may be defined even when _no_
observations have zero entries in the design m
On 13-Jan-09 22:18:00, glenn wrote:
> Very simple one sorry:
>
> How do I join 2 lists please
>
> Thanks
> glenn
c(list1,list2)
Example:
A<-list(a=1,b=2,c=3)
B<-list(d=4,e=5,f=6)
c(A,B)
# $a
# [1] 1
# $b
# [1] 2
# $c
# [1] 3
# $d
# [1] 4
# $e
# [1] 5
# $f
# [1] 6
Ted.
---
Dear R-Users,
I put a small sample data set and script.
Aim: combine two partly overlapping series to one by prediction.
Problem: only overlapping data points are predicted.
Question: How do I predict data for rows 1-9 and 14-16?
Thanks in advance for your advince,
Tim
### CODE ###
x <- read
Very simple one sorry:
How do I join 2 lists please
Thanks
glenn
[[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 http://www.R-project.org/
Peter Dalgaard wrote:
Actually, notice that you are averaging identical values, so the "mean"
in the tapply is slightly misleading.
Notice also that the intercept may be defined even when _no_
observations have zero entries in the design matrix. This is the usual
case in linear regression, for
Dear list,
I would like to run Lasso (least absolute shrinkage and selection
operator) on a bimodal dataset. Can I do that using packages like lasso
or lars ?
Thanks in advance.
Stella Sim
DISCLAIMER:\ This email contains confidential informatio...{{dropped:11}}
Hello all,
I try the follow test with an odf file.
*INPUT: teste.odt*
>=
x <- matrix(rnorm(6), nc=2)
xst <- tableStyles(x, useRowNames = TRUE, header = NULL)
odfTable(x, styles=xst, colnames = NULL, useRowNames = TRUE)
@
and I get this error:
*OUTPUT error:*
Error: chunk 1 (label=teste)
Err
xtable package and depending on what word processor you're using export as
latex or html. Fancier approach is xtable combined with Sweave.
--- On Tue, 1/13/09, Jason Rupert wrote:
> From: Jason Rupert
> Subject: [R] Tables for journal/conference publications from within R
> To: r-help@r-pr
Hi Harold: Below works on your data set but check it a lot because I am
a little worried that
I could have missed something. Hopefully someone can send a a little
clearer way.
dat <- data.frame(id = c(1,1,2,2,2), var1 = c(10,10,20,20,25), var2 =
c('foo', 'foo', 'foo', 'foobar', 'foo'))
print(
Gregor Gorjanc wrote:
hsl.gov.uk> writes:
I'm running OpenBUGS model via the R2WinBUGS package interface, under
Windows. Is it possible to terminate running models, short of using the
Windows Task Manager to forcibly exit the program?
If you use OpenBUGS, then I guess you can not since R
> look at zoo and ts, and it all depends on what you want to do.
I researched a bit.
I am looking at the application of different prediction methods
(predict.lm etc.) for time series.
I have two series of data.
One contains measurements of 30 years.
The other just 5 years. Both overla
on 01/13/2009 01:17 PM Doran, Harold wrote:
> Suppose I have a dataframe as follows:
>
> dat <- data.frame(id = c(1,1,2,2,2), var1 = c(10,10,20,20,25), var2 =
> c('foo', 'foo', 'foo', 'foobar', 'foo'))
>
> Now, if I were to subset by id, such as:
>
>> subset(dat, id==1)
> id var1 var2
> 1 1
Nice. I was thinking maybe length(table(x) == 1), but this works great
> -Original Message-
> From: Carlos J. Gil Bellosta [mailto:c...@datanalytics.com]
> Sent: Tuesday, January 13, 2009 2:55 PM
> To: Doran, Harold
> Cc: r-help@r-project.org
> Subject: Re: [R] Comparing elements for equa
Hello,
You could build your output dataframe along the following lines:
foo <- function(x) length( unique(x) ) == 1
results <- data.frame(
freq = tapply( dat$id, dat$id, length ),
var1 = tapply( dat$var1, dat$id, foo ),
var2 = tapply( dat$var2, dat$id, foo )
)
Best reg
On Jan 13, 2009, at 10:28 AM, Gabor Grothendieck wrote:
If efficiency is not of concern then this is easy to understand:
for(i in 1:5) for(j in 1:5) if (i < j) { ... }
or?:
> xx <- expand.grid(i=1:5, j=1:5)
> subset(xx, i < j)
i j
6 1 2
11 1 3
12 2 3
16 1 4
17 2 4
18 3 4
21 1 5
22 2 5
Hi Brigid,
You will probably get some more informed answers in a bit, but to give
you some quick things to try...
There's no size limit to an object like you are referring to. When
you're having problems with small objects or simple operations, it
sounds like you've used up the memory for ru
Suppose I have a dataframe as follows:
dat <- data.frame(id = c(1,1,2,2,2), var1 = c(10,10,20,20,25), var2 =
c('foo', 'foo', 'foo', 'foobar', 'foo'))
Now, if I were to subset by id, such as:
> subset(dat, id==1)
id var1 var2
1 1 10 foo
2 1 10 foo
I can see that the elements in var1 ar
Try writing out the dataframe in smaller pieces. If you have 800K row
and 18 columns, then this is about 14M items and if you have 8 bytes
per items (assuming all numeric; more for character), this is 115MB of
space that would probably be required to construct the output data (or
more). Your erro
On Tue, 13 Jan 2009, Matthieu Stigler wrote:
Hello
I'm trying to run a fortran code which use LAPACK subroutines. I think I
should use some points shown in the manual 5.5 Creating shared objects but it
is too technical for me :-(...
Be sure to look over section 6, esp. 6.17
The headers f
gregor rolshausen biologie.uni-freiburg.de> writes:
>
> ok. sorry for being blurry.
>
> I have x,y data, that probably fits a asymptotic curve (asymptote at N).
> now I want to fit a curve onto the data, that gives me the N. therefore
> I thought to fit an e-function, namely N(1-e^(y/x)) onto
hsl.gov.uk> writes:
> I'm running OpenBUGS model via the R2WinBUGS package interface, under
> Windows. Is it possible to terminate running models, short of using the
> Windows Task Manager to forcibly exit the program?
If you use OpenBUGS, then I guess you can not since R2WinBUGS just passes
on 01/13/2009 08:43 AM joe1985 wrote:
> Hello
>
> I have a problem that i ant make a Bar plot like the one i have tried to
> illustrate below (made in paint);
> http://www.nabble.com/file/p21437080/LG5%2Bgraf%2Bredigeret.jpg
>
> http://www.nabble.com/file/p21437080/LG5%2Bgraf%2Bredigeret.JPG
>
It's fairly clear from the documentation that approxfun() will not
extrapolate.
help.search("extrapolate")
library(Hmisc)
?approxExtrap
Some sort of minimization approach:
> approxExtrap(x=c(0,5,10,15,20), y=c(16,45,77,101,125),xout=c(-4,0,4))
$x
[1] -4 0 4
$y
[1] -7.2 16.0 39.2
> approxE
on 01/13/2009 11:25 AM Peter Dalgaard wrote:
> Marc Schwartz wrote:
>
>>> DF.fitted
>> Y A B F.lm
>> 1 21.86773 0 a 23.52957
>> 2 25.91822 0 a 23.52957
>> 3 20.82186 0 a 23.52957
>> 4 42.97640 1 a 36.18023
>> 5 36.64754 1 a 36.18023
>> 6 30.89766 1 a 36.18023
>> 7 47.43715 0 b
Marc Schwartz wrote:
>
>> DF.fitted
> Y A B F.lm
> 1 21.86773 0 a 23.52957
> 2 25.91822 0 a 23.52957
> 3 20.82186 0 a 23.52957
> 4 42.97640 1 a 36.18023
> 5 36.64754 1 a 36.18023
> 6 30.89766 1 a 36.18023
> 7 47.43715 0 b 46.50615
> 8 48.69162 0 b 46.50615
> 9 47.87891 0 b
On Tue, 13 Jan 2009, jim holtman wrote:
That is certainly true because I have seen differences due to the
sharing of values. I also look at what 'gc()' shows at the memory
being used. Does this provide a reasonable estimate of the total
space being used?
Yes, but not just by your workspace,
My apologies if this is a bit of a 'newbie' question.
I am using R v 2.8.0 in Windows and am a bit confused about the memory
size/allocation.
A script I wrote faulted out with the error: "Error: cannot allocate vector
of size 5.6 Mb"
After this error, I still have:
> memory.size()
[1] 669.3517
... which is the same as
setdiff(union(x, y), intersect(x, y))
-Christos
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of
> dav...@rhotrading.com
> Sent: Tuesday, January 13, 2009 10:56 AM
> To: Juliet Hannah; r-help@r-proj
on 01/13/2009 09:32 AM Stefano Leonardi wrote:
> Hallo,
> yesterday I was puzzled when I discovered that I
> probabliy miss something in the interepretation of intercept
> in two-way lm models.
>
> I thought that the intercept, using the default contr.treatment
> contrasts, represents the mean of
Hello everyone,
I have a data set that looks like the following:
Year Days to the beginning of YearValue
1 30 100
1 60200
1..
Dear all,
is there in R a function to run 2D kernel density estimation (similar to e.g.
'kde2d()') in a space that is split internally?
I think about a store floor space - basically a rectangle but divided by
shelves etc. into "alleys". The "points" are people in the store and I would
like to d
On Tue, 13 Jan 2009, David Winsemius wrote:
I remember have the same consternation using GLIM with binomial models on
grouped and ungrouped data, but I was counseled by my betters only to
consider differences in models. The differences in deviance are the same up
to rounding error.
859.8018
That is certainly true because I have seen differences due to the
sharing of values. I also look at what 'gc()' shows at the memory
being used. Does this provide a reasonable estimate of the total
space being used?
On Tue, Jan 13, 2009 at 11:14 AM, Prof Brian Ripley
wrote:
> On Tue, 13 Jan 2009
Hello,
The symmetric set difference of A and B is the set of elements in A or B
but not in A intersection B, i.e., ( (A U B) \ (A intersection B) ).
The asymmetric set difference of A and B is the set of elements of A
except those in B, i.e., (A \ B).
Best regards,
Carlos J. Gil Bellosta
http:/
On 1/13/2009 10:41 AM, Juliet Hannah wrote:
Sorry if this is a silly question. What does asymmetric refer to from
?intersect, and are there any
aspects of the result that may be different from expected (reason for
exclamation point) ?
There are two definitions for a set difference: the symmetr
Hi, i'm an italian student that use for the first time the mailing list. I need
an help to use a function for elaborate a multinomial logit analisys. I'm
making a paper on the USA's commons and their way of organize their service.
The various tipologies of service are 4: public, contracting in
>
> What is the problem that you are trying to solve?
>
>From the data I provided: x=c(0,5,10,15,20) y=c(16,45,77,101,125); I
want to obtain the value of x when y=0.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
There is also a symmetric set difference = union(setdiff(x, y),
setdiff(y, x))
denoted x \Delta y
-- David
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Juliet Hannah
Sent: Tuesday, January 13, 2009 9:41 AM
To: r-help@r-project.o
hi, I have two questions:
#first (SPSS vs. R):
I just compared the output of different PCA routines in R (pca, prcomp,
princomp) with results from SPSS. the loadings of the variables differ
vastly! in SPSS the variables load constantly higher than in R.
I made sure that both progr. use the cor
Because there is also a symmetric difference: see
http://en.wikipedia.org/wiki/Symmetric_difference
Sorry if this is a silly question. What does asymmetric refer to from
?intersect, and are there any
aspects of the result that may be different from expected (reason for
exclamation point) ?
As
On Tue, Jan 13, 2009 at 10:23 AM, jim holtman wrote:
> How about this:
>
> > data(ToothGrowth)
> > ls()
> [1] "ToothGrowth"
> > data <- function(x){invisible(NULL)}
> > data(ToothGrowth)
> > ls()
> [1] "data"
> >
>
Yep, that sure does cause a problem alright. Is it the case that that
problems ari
On Tue, 13 Jan 2009, jim holtman wrote:
Here is a function I use to see how big the objects in my workspace are:
Be careful with the caveats spelled out in ?object.size. Especially
for character data such summations can be way off.
my.ls <-
+ function (pos = 1, sorted = F)
+ {
+ .r
If efficiency is not of concern then this is easy to understand:
for(i in 1:5) for(j in 1:5) if (i < j) { ... }
On Tue, Jan 13, 2009 at 8:20 AM, Niccolò Bassani
wrote:
> Dear R users,I come to you with a quite silly question I think, but I hope
> you can answer me...
> That is, I've got some pr
I remember have the same consternation using GLIM with binomial models
on grouped and ungrouped data, but I was counseled by my betters only
to consider differences in models. The differences in deviance are the
same up to rounding error.
> 859.8018 - 711.3479
[1] 148.4539
> 168.8-20.3
[1]
Hello
I have a problem that i ant make a Bar plot like the one i have tried to
illustrate below (made in paint);
http://www.nabble.com/file/p21437080/LG5%2Bgraf%2Bredigeret.jpg
http://www.nabble.com/file/p21437080/LG5%2Bgraf%2Bredigeret.JPG
LG5+graf+redigeret.JPG
Where each line represents a
On 1/13/2009 9:01 AM, Powers, Randall - BLS wrote:
Hello All,
I'm actually the system administrator of a UNIX system where several
users use R version 2.6.0. I have a user who is trying to use the SURVEY
package, and when he does, he gets the message:
survey' is not a valid package -- installe
This doesn't address your original question, but you might want to look at
using the Sweash.sh script available from
http://lib.stat.cmu.edu/R/CRAN/contrib/extra/scripts/. Download the file,
move it to /usr/local/bin (e.g., mv Sweave.sh /usr/local/bin) and make sure
it's executable (chmod +x /usr/l
Sorry if this is a silly question. What does asymmetric refer to from
?intersect, and are there any
aspects of the result that may be different from expected (reason for
exclamation point) ?
As toy data, if needed, here are the examples from the help page.
(x <- c(sort(sample(1:20, 9)),NA))
(y <-
R does understand the class of the object so that it knows it should
not try to call the dataframe 'data' if the function 'data' is being
called. It is probably best to avoid what would be termed "reserved"
words just to prevent mistakes in the future. How about calling them
"Data"?
On Tue, Jan
Hello All,
I'm actually the system administrator of a UNIX system where several
users use R version 2.6.0. I have a user who is trying to use the SURVEY
package, and when he does, he gets the message:
survey' is not a valid package -- installed < 2.0.0?
When I run the library() command, I get
Hallo,
yesterday I was puzzled when I discovered that I
probabliy miss something in the interepretation of intercept
in two-way lm models.
I thought that the intercept, using the default contr.treatment
contrasts, represents the mean of the group of observations
having zero in all column of the
approxfun returns a function; that is not an error message:
> x=c(0,5,10,15,20)
> y=c(16,45,77,101,125)
>
> approx(x,y,method="linear")
$x
[1] 0.000 0.4081633 0.8163265 1.2244898 1.6326531 2.0408163
2.4489796 2.8571429 3.2653061
[10] 3.6734694 4.0816327 4.4897959 4.8979592 5.306
The data.table package may be more in line with what you are after, but xts
and zoo can also do what you need in this particular example:
> a <- xts(c('a1','a2','a3'), timeBasedSeq(20090101/20090103))
> colnames(a) <- 'foo'
> b <- xts(c('b1'), as.Date('2009-01-04'))
> colnames(b) <- 'foo'
> a
Here is a function I use to see how big the objects in my workspace are:
> my.ls <-
+ function (pos = 1, sorted = F)
+ {
+ .result <- sapply(ls(pos = pos, all.names = TRUE), function(..x)
object.size(eval(as.symbol(..x
+ if (sorted) {
+ .result <- rev(sort(.result))
+ }
+
How about this:
> data(ToothGrowth)
> ls()
[1] "ToothGrowth"
> data <- function(x){invisible(NULL)}
> data(ToothGrowth)
> ls()
[1] "data"
>
On Tue, Jan 13, 2009 at 9:53 AM, Ista Zahn wrote:
> From: baptiste auguie
> To: Dimitris Rizopoulos
> Date: Tue, 13 Jan 2009 09:38:09 +
> Subject: Re
On 13-Jan-09 13:20:54, Niccolò Bassani wrote:
> Dear R users,I come to you with a quite silly question I think,
> but I hope you can answer me...
> That is, I've got some problems in using the if and while conditions
> in a loop.
> Substantially, I want to select rows in a dataset depending on an
>
Withut looking up your reference, are you not comparing grouped and
ungrouped deviances? And polr() does not say anything about accepting
a model (or not), only about the comparison between two models.
'Deviances' are in comparison with some 'saturated' model, and I would
say that M&N are com
Dear all,
Is there a way we can find the total object.size of
all the objects in our R script?
The reason we want to do this because we want to know
how much memory does our R script require overall.
Rprofmem(), doesn't seem to do it.
and Unix 'top' command is dynamic and
it doesn't give the ex
Generally extraction of subsets is better done by subset().
myDat <- subset(dataset, subset= (variable==i) | (variable==j) ) #, or
myDat <- dataset[which( (dataset$variable==i) | (dataset
$variable==j) ), ]
#note the need to add the name of the dataframe in the second version.
Both of these
From: baptiste auguie
To: Dimitris Rizopoulos
Date: Tue, 13 Jan 2009 09:38:09 +
Subject: Re: [R] indexing question
> you can also look at subset,
>
>
>my.data.frame <- data.frame(a=rnorm(10),
>> b=factor(sample(letters[1:4], 10, replace=T)))
>>str(my.data.frame)
>>my.
Thanks for your most reasonable reply, Henrique.
- Gundala Viswanath
Jakarta - Indonesia
On Tue, Jan 13, 2009 at 8:01 PM, Henrique Dallazuanna wrote:
> Try this:
>
> unlist(ifelse(q %in% names(x), x[q], NA))
>
> On Tue, Jan 13, 2009 at 8:49 AM, Gundala Viswanath
> wrote:
>>
>> Dear all,
>>
>>
Which R version do you have? I'm asking this because my 2.7.0 gives a
different error message:
> x[[q]]
Error in x[[q]] : recursive indexing failed at level 2
Anyway, as Wacek said, x[[q]] is equivalent to
x[["some"]][["more"]][["not_there"]] -- and you don't have an element called
"more" in x[[
On 08/01/2009, Greg Snow wrote:
> If you want to just linearly interpolate, then use the functions approx or
> approxfun from the stats package (one of those that is loaded by default).
I have read the guide for approx and approxfun functions. Below is my data.
x=c(0,5,10,15,20)
y=c(16,45,77,101,
Dear Jason,
I would suggest xtable (from the xtable package) and LaTeX / Sweave.
HTH,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
Cel biometrie, met
help.search did not offer any clues. Using Baron's search page, I got
reference to r-help entries from 2004, but looking at the current
documentation , I see no such function. I am wondering if the name of
the function was changed in later versions of the package and you are
mixing old doc
Dear all,
I've replicated the cheese tasting example on p175 of GLM's by McCullagh
and Nelder. This is a 4 treatment (rows) by 9 ordinal response (cols)
table.
Here's my simple code:
cheese
library(MASS)
options(contrasts = c("contr.treatment", "contr.poly"))
y = c
Yeah.
addtable2plot in the plotrix is really close, but, working from within R, I
would really like to be able format tables in such a way that they look like
they are ready for journal or conference paper publication.
Is addtable2plot in the plotrix the right tool for that or are there o
On 1/13/2009 8:20 AM, Niccolò Bassani wrote:
Dear R users,I come to you with a quite silly question I think, but I hope
you can answer me...
That is, I've got some problems in using the if and while conditions in a
loop.
Substantially, I want to select rows in a dataset depending on an index
vari
Shouldn't your loop be:
for (i in 1:4){
for (j in (i+1):5{
...
}
}
On Tue, Jan 13, 2009 at 8:20 AM, Niccolò Bassani
wrote:
> Dear R users,I come to you with a quite silly question I think, but I hope
> you can answer me...
> That is, I've got some problems in using the if and while cond
Dear R users,I come to you with a quite silly question I think, but I hope
you can answer me...
That is, I've got some problems in using the if and while conditions in a
loop.
Substantially, I want to select rows in a dataset depending on an index
variable (suppose it ranges from 1 to 5), so to mak
Vincenzo,
there is no function 'mcmcFmodel' in Geneland. Perhaps you intend to
use 'MCMC'?
See ?MCMC for a proper use.
Additionally you might profit from reading the Geneland Manual on
http://folk.uio.no/gillesg/Geneland/Geneland.html
HTH
Christoph
___
Dear Jim and Dorothee,
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On
> Behalf Of Jim Lemon
> Sent: January-13-09 5:17 AM
> To: Dorothee
> Cc: r-help@r-project.org
> Subject: Re: [R] polychoric correlation: issue with coefficient sign
>
ok. sorry for being blurry.
I have x,y data, that probably fits a asymptotic curve (asymptote at N).
now I want to fit a curve onto the data, that gives me the N. therefore
I thought to fit an e-function, namely N(1-e^(y/x)) onto the data and
get the N from the fitted curves' equation.
in the
Hello
I'm trying to run a fortran code which use LAPACK subroutines. I think I
should use some points shown in the manual 5.5 Creating shared objects
but it is too technical for me :-(... Could anyone help me for the
procedure to do:
-which part of the manual is relevant for this type of que
gregor rolshausen wrote:
hello,
I want to fit a curve to a simple x,y dataset - my problem is, that I
want to fit it for the following term:
n(1-e^x/y) - so I get the n constant for my data...
Not an R problem in the first place, but the question arises what
"n(1-e^x/y)" means, its is jus
This is an ugly hack but add this line to your sourced script:
this.dir <- dirname(parent.frame(2)$ofile)
It must be at top level, i.e. not in a function in the script.
On Tue, Jan 13, 2009 at 6:53 AM, ppaarrkk wrote:
>
> Is it possible for an R script to read a text file, say, from the directo
Wacek Kusnierczyk wrote:
>
> x[[q]] is equivalent to x[['some']][['more']][['not_there']]. since
> x[['some']] is an atomic integer vector, it won't collaborate with
> [['more']]
>
æsj, one more lousy statement. x[['some']] is an atomic *double* vector
*with no element named 'more'*, hence x[
1 - 100 of 125 matches
Mail list logo
