A
http://polynomial-trimonial-binomial.blogspot.com/2011/10/how-to-deal-with-polynomials.html
polynomial function is evaluated by foiling out the equation so you can
solve for 'x'.
--
View this message in context:
http://r.789695.n4.nabble.com/lm-how-are-polynomial-functions-interpreted-tp87566
Hello Petr,
The demo's don't run either, with the same errors.
Thanks for your help.
Best wishes.
--
View this message in context:
http://r.789695.n4.nabble.com/Error-in-matrix-not-ordered-vectors-or-numerical-value-and-SIAR-tp4024578p4033682.html
Sent from the R help mailing list archive
Hi all
I am trying to fit a distribution (i.e. gamma) to some data and I understand
how to use the fitdistr from the MASS package. However, what should you do
when the data has a probability associated with it that are not all equal.
e.g.
Pr(Occurrence) Size
0.00460
0.02
I'm trying out a basic plot, but something about the way I subset my data
leads to problems with the plot.
Here is the first bit of my data set
year,date,location,quadrat_juvenile,photo_location,photo_exists,genus,count,divers
2005,2005-04-30, 1 Fringing Reef,1, 1 Fringing Reef Coral Transect Po
On Nov 11, 2011, at 9:08 PM, RD235 wrote:
The code
library(Design)
f <- lrm(y~x1+x2+x1*x2, data=data)
plot(f)
produces a plot of log odds vs x2 with 0.95 confidence intervals.
How do I
get a plot of odds ratios vs x2 instead?
You would construct a dataset that had selected combinations o
The Design package is obsolete. Use its replacement rms - see
http://biostat.mc.vanderbilt.edu/Rrms. Use something like the following to
reproduce what you already have:
To plot odds ratios against a specific value of x2 and at a specific x1,
type ?contrast.rms to see examples.
Frank
RD235 wrot
The code
library(Design)
f <- lrm(y~x1+x2+x1*x2, data=data)
plot(f)
produces a plot of log odds vs x2 with 0.95 confidence intervals. How do I
get a plot of odds ratios vs x2 instead?
Thanks
--
View this message in context:
http://r.789695.n4.nabble.com/Odds-ratios-from-lrm-plot-tp4033340p403
On Nov 11, 2011, at 4:04 PM, Ben Bolker wrote:
> Jun Ji stanford.edu> writes:
>
>> Is there any encryption package of R? For instance, 3DES or AES algorithm
>> package or whatever...
>> Thanks a lot in advance.
>>
>
> library(sos)
> findFn("encryption")
>
> finds a partial AES implementati
Dear list members,
I know that there is the arules package with the implementation of the
apriori algorithm. However i want to use the "predictive apriori" instead.
These algorithm can mine as rules as i want and there is an implementation
on weka.
There is some implementation on R?
--
Att,
Fl
Hi folks,
i have some problems with my evaluation. We have collect tons of data from
23 testpersons for our new road study.
I have now a time series for each person and all the logs when he
accelerates or hits the break trying to solve five different tasks.
The dataset lools like:
#Unixtim
Hi:
This doesn't sort the data by strain level, but I think it does what
you're after. It helps if strain is either a factor or character
vector in each data frame.
h <- function(x, y) {
tbx <- table(x$strain)
tby <- table(y$strain)
# Select the strains who have more than one memb
Hi:
R-Bloggers picked up on this web site that contains several
interesting posts re usage of R, including this one on using a
Bayesian approach to multivariate mixed effects models using the
MCMCglmm package:
http://www.quantumforest.com/2011/11/coming-out-of-the-bayesian-closet-multivariate-vers
Thanks Ted!
I really appreciate your time!
Thanks for the link about the 'problem of calibration', and your suggestion to
reformulate my model. I had no idea about it before. I certainly learnt
something today.
I will try your suggestions later today and let you know how it works out.
Stefan
Use 'intersect' to get the items common in both dataframes and then use
that to extract the data in common.
On Friday, November 11, 2011, kickout wrote:
> I've scoured the archives but have found no concrete answer to my
question.
>
> Problem: Two data sets
>
> 1st data set(x) = 20,000 rows
> 2nd
On Fri, Nov 11, 2011 at 4:41 PM, Supreet kaur wrote:
> hello all R experts,
> how do I calculate the reliability between the two groups
> using the ICCs?
Possibly by using GmeanRel() from the multilevel package.
Searching for
reliability group ICC
at http://www.rseek.org
offers s
What about merge() with all=FALSE?
> x <- data.frame(a=letters[1:6], b=1:6)
> y <- data.frame(a=letters[4:9], b=11:16)
> x
a b
1 a 1
2 b 2
3 c 3
4 d 4
5 e 5
6 f 6
> y
a b
1 d 11
2 e 12
3 f 13
4 g 14
5 h 15
6 i 16
> merge(x, y, by="a", all=FALSE)
a b.x b.y
1 d 4 11
2 e 5 12
3 f 6 13
Follow-up: See at end.
On 11-Nov-11 21:16:02, Ted Harding wrote:
> On 11-Nov-11 14:51:19, Schreiber, Stefan wrote:
>> Dear list members,
>>
>> I have just a quick question:
>>
>> I fitted a non-linear model y=a/x+b to describe my data
>> (x=temperature and y=damage in %) and it works really nice
Hi
I have an organism directory that contains two folders galGal3 and hg19 and
many other files.
orgDir = '/home/mary/org'
When I try to use list.dir() function, it gives me the same answer, no
matter what is the value of full.names argument.
> list.dirs(path = indexDir, full.names = FALSE)[1]
hello all R experts,
how do I calculate the reliability between the two groups
using the ICCs?
I'll appreciate your reply,
Thanks
Sincerely,
Supreet kaur,
Biomedical research engineer,
Nationwide Childrens Hospital,
Columbus, OH
(614)355-3509
[[alternative HTML version
I've scoured the archives but have found no concrete answer to my question.
Problem: Two data sets
1st data set(x) = 20,000 rows
2nd data set(y) = 5,000 rows
Both have the same column names, the column of interest to me is a variable
called strain.
For example, a strain named "Chab1405" appear
I would suggestto Mario that, if he wishes to use colorspace and keep
AVG, he should install R 2.13.2. The colorspace library built with
that version of R does not give a false positive with AVG. To update
paackages he must remove not only colorspace but also its reverse
dependencies, the revers
Jun Ji stanford.edu> writes:
> Is there any encryption package of R? For instance, 3DES or AES algorithm
> package or whatever...
> Thanks a lot in advance.
>
library(sos)
findFn("encryption")
finds a partial AES implementation, but not much else.
Dear Hans,
[see inline below]
On 11 November 2011 22:44, Hans W Borchers wrote:
> baptiste auguie googlemail.com> writes:
>
>>
>> Dear list,
>>
>> [cross-posting from Stack Overflow where this question has remained
>> unanswered for two weeks]
>>
>> I'd like to perform a numerical integration i
Dear Ravi,
Thank you for your answer.
The integrand I proposed was a dummy example for demonstration
purposes. I experienced a similar slowdown in a real problem, where
knowing in advance the shape of the integrand would not be so easy.
Your advice is sound; I would have to study the underlying
Dear colleagues,
I'm trying to fit a multinomial logistic regression for an ordinal variable.
I see in the help pages for multinom in nnet that one should scale the
predictors from 0-1. Is that really necessary?
Also: can anyone clarify what the difference between alternative-specific and
indi
On 11-Nov-11 14:51:19, Schreiber, Stefan wrote:
> Dear list members,
>
> I have just a quick question:
>
> I fitted a non-linear model y=a/x+b to describe my data
> (x=temperature and y=damage in %) and it works really nicely
> (see example below). I have 7 different species and 8 individuals
> p
Hi Francesca,
Try something like this:
x <- c(1, 3, 7)
dati <- lapply(1:4, function(i) {datiP[datiP$city == i, x]})
dati[[1]] # datiP1
dati[[2]] # datiP2
dati[[3]] # datiP3
dati[[4]] # datiP4
if the *only* groups are 1, 2, 3, 4 (i.e., 1:4 is exhaustive), this
can be simplified:
dati <- by(dat
Hi,
It is not at all clear to me what type of model you are trying to fit.
You could consider ?manova for a multivariate analysis of variance
require(lme4)
?lmer
for longitudinal regression type models
or perhaps
OpenMx: http://openmx.psyc.virginia.edu/
for structural equation modelling in R
Not sure if this helps, but did you try Google?
http://www.jstatsoft.org/v35/i09/paper
http://cran.r-project.org/web/packages/lcmm/lcmm.pdf
http://www.warwick.ac.uk/statsdept/useR-2011/abstracts/010411-liquetbenoit.pdf
yurirouge wrote:
>
> HI,
>
> I am relatively new to R and would appreci
Dear Contributors
I would like to perform this operation using a loop, instead of repeating
the same operation many times.
The numbers from 1 to 4 related to different groups that are in the
database and for which I have the same data.
x<-c(1,3,7)
datiP1 <- datiP[datiP$city ==1,x];
datiP2 <
HI,
I am relatively new to R and would appreciate some help or directions for
this.
I am trying to model 3 longitudinal outcomes jointly and to identify some
predictors for these 3 joint outcomes (all continuous). I am trying to find
some codes that I may modify to do this but cannot seem to fin
It seems that there is a bug in the forecast::tslm function. I have forwarded
what I think is the bug (the call to get() should supply the argument
'envir=parent.frame()').
Thank you.
On Nov 11, 2011, at 12:11 PM, Joshua Wiley wrote:
> Seems like it could work---can you save your script as a
Hi Arnau,
Not aware of direct implementation. It was discussed in octave project as well
(http://octave.1599824.n4.nabble.com/Jordan-canonical-form-td2216965.html)
It is numerically ill-conditioned to compute that.
See this http://dl.acm.org/citation.cfm?id=355912
Best,
Mehmet Süzen, PhD
Man
Here's my setup:
- I'm on a Windows machine (I don't have full admin rights)
- I have a folder with an *.RData file and an .RProfile file
- I want the user to be able to start R by double clicking on the
*.RData file
Can I specify the application start up options (like --no-save --md
Dear all:
Is there any encryption package of R? For instance, 3DES or AES algorithm
package or whatever...
Thanks a lot in advance.
Best regards,
Jason
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz
Hi,
I was wondering if it is possible to harmonize the ouput of symbols()
and legend() both from the graphics package.
Let us take this example:
x<-runif(10)
y<-runif(10)
z<-runif(10)
leg<-round(seq(min(z),max(z),l=4),2) # 4 values rounded up to 2 decimals
for the legend
symbols(x,y,circl
On 11.11.2011 18:48, Gyanendra Pokharel wrote:
Following is my code, can some one help on the error at the bottom?
mh<-function(iterations,alpha,beta){
+ data<-read.table("epidemic.txt",header = TRUE)
+ attach(data, warn.conflicts = F)
+ k<-97
+ d<- (sqrt((x-x[k])^2 + (y-y[k])^2))
+
Seems like it could work---can you save your script as a .txt file and
send it as an attachment, upload it online, or make a reproducible
example?
Cheers,
Josh
On Fri, Nov 11, 2011 at 9:52 AM, wrote:
>
> I have an R script with the following applicable lines:
>
> xshort <- windo
I have an R script with the following applicable lines:
xshort <- window(s, start=st, end=ed)
. . .
xshort <- ts(xshort, frequency=1, start=1)
. . .
m1 <- m2 <- m3 <- m4 <- m5 <- m6 <- NULL
m1 <- tslm(xshort ~ trend)
I get an error:
Following is my code, can some one help on the error at the bottom?
> mh<-function(iterations,alpha,beta){
+ data<-read.table("epidemic.txt",header = TRUE)
+ attach(data, warn.conflicts = F)
+ k<-97
+ d <- (sqrt((x-x[k])^2 + (y-y[k])^2))
+ p <- 1-exp(-alpha*d^(-beta))
+ p.alpha<-1
for the googleable r-help archives, I thought I would post what I
wrote into my .Rprofile to automatically set some system information.
the most relevant aspect is the determination of mc.cores. this is
useful when users want to use the parallel package
options(uname= system("uname", intern=TRU
As the doctor says, if it hurts "don't do that".
A factor is a sequence of integers with a corresponding list of character
strings. Factors in two separate vectors can and usually do map the same
integer to different strings, and R cannot tell how you want that resolved.
Convert these columns t
I suspect that part of the hesitancy in replying to your query is that your
academic e-mail and subsequent google-ability suggest you are a student in the
economics department and that this might be homework.
In the meanwhile, you might want to look into Pfaff's book on cointegration and
time
Hello.
Is it possible to find the Jordan Form of a matrix with R?
Arnau.
Arnau Mir Torres
Edifici A. Turmeda
Campus UIB
Ctra. Valldemossa, km. 7,5
07122 Palma de Mca.
tel: (+34) 971172987
fax: (+34) 971173003
email: arnau@uib.es
UR
On Fri, Nov 11, 2011 at 2:09 AM, Prof Brian Ripley wrote:
> On Thu, 10 Nov 2011, Tyler Pirtle wrote:
>
> Hi there,
>>
>> I'm writing an R extension that has a C component that relies on two third
>> party libraries that I'm bundling
>> with the extension.
>>
>> I'd like to statically link my resu
On Fri, 11 Nov 2011, Tyler Pirtle wrote:
On Fri, Nov 11, 2011 at 2:09 AM, Prof Brian Ripley
wrote:
On Thu, 10 Nov 2011, Tyler Pirtle wrote:
Hi there,
I'm writing an R extension that has a C component
that relies on two third
party librar
Your question is a bit too general to give a useful answer. One possible
answer to your question is:
> mrv <- matrix( runif(1000), ncol=10 )
Which generates multivariate random observations, but is unlikely to be what
you are really trying to accomplish. There are many tools for generating
m
Those formulas are the standard way to convert from polar coordinates to
Euclidean coordinates. The polar coordinates are 'r' which is the radius or
distance from the center point and 'theta' which is the angle (0 is pointing in
the positive x direction).
If r is constant and theta coveres a f
'The R Inferno' Circle 8.1.66
http://www.burns-stat.com/pages/Tutor/R_inferno.pdf
On 11/11/2011 00:17, Worik R wrote:
It seems obvious to me that the empty string "" is length 0.
cheers
Worik
[[alternative HTML version deleted]]
__
R-help@r
The integrand is highly peaked. It is approximately an impulse function where
much of the mass is concentrated at a very small interval. Plot the function
and see for yourself. This is the likely cause of the problem.
Other types of integrands where you could experience problems are: integran
- Forwarded Message -
From: vrama...@neo.tamu.edu
To: "bernhard pfaff"
Sent: Friday, November 11, 2011 9:03:11 AM GMT -06:00 US/Canada Central
Subject: Use of R for VECM
Hello Fellow R'ers
I am a new user of R and I am applying it for solving Bi-Variate (Consumption
and Output) VECM wit
Some tips:
1) Excel did not, as far as I can determine, find a solution. No point seems to
satisfy
the KKT conditions (there is a function kktc in optfntools on R-forge project
optimizer.
It is called by optimx).
2) Scaling of the input vector is a good idea given the seeming wide range of
val
Thank you very much to everyone who replied!
As I mentioned - I am not a mathematician, so sorry for stupid
comments/questions.
I intuitively understand what you mean by scaling. While the solution
space for the first parameter (.alpha) is relatively compact (probably
between 0 and 2), the second o
On 10.11.2011 17:53, Kevin Burton wrote:
The problem with this documentation is two-fold. One it seems to concentrate
on building from source which I don't need. Two it doesn't address the
upgade. I have a number of packages and so I need to do what has been
suggested and install the latest ver
That line in the documentation refers to groups (i.e. columns) on each row
of data. C.1 and C.2 are your groups and if you use the beside=TRUE option
C.1 and C.2 will appear side-by-side with a single labels on the axis
(instead of stacked one atop the other. You could rearrange your data so
that c
Dear list members,
I have just a quick question:
I fitted a non-linear model y=a/x+b to describe my data (x=temperature and
y=damage in %) and it works really nicely (see example below). I have 7
different species and 8 individuals per species. I measured damage for each
individual per specie
Hi Dimitri,
Your problem has little to do with local versus global optimum. You can
convince yourself that the solution you got is not even a local optimum by
checking the gradient at the solution.
The main issue is that your objective function is not differentiable
everywhere. So, you have
Hi,
On Fri, Nov 11, 2011 at 7:26 AM, Oliver wrote:
> Christofer Bogaso gmail.com> writes:
> [...]
>> Here my question is, is there any speed reduction if I put them within a
>> function (I think there may be some speed reduction at least within
>> for-loop, because that loop needs to call that f
In http://finzi.psych.upenn.edu/R/Rhelp02a/archive/14637.html
some rudimentary R functions were given for drawing
proportional area venn diagrams with area of each intersection ~ the
count in a 2 x 2 x 2 table.
I'm interested in this, for another application: showing the
correlations among
It's the same problem as with any regression model where the predictors
are not strictly orthogonal. You can make the explained deviances per
term sum to the whole model explained deviance by dropping terms
sequentially, but then the explained deviance per term depends on the
order in which you
I won't requote all the other msgs, but the latest (and possibly a bit glitchy)
version of
optimx on R-forge
1) finds that some methods wander into domains where the user function fails
try() (new
optimx runs try() around all function calls). This includes L-BFGS-B
2) reports that the scaling i
Or even simpler (unbelievable!), see ?title:
par(mfcol=c(2,2), oma=c(0,0,1,0))
replicate(4, plot(1))
title("Hello World!", outer=TRUE)
Best,
Uwe Ligges
On 10.11.2011 17:13, Sarah Goslee wrote:
Try ?mtext for information on how to plot into the outer margin of a group of
plots.
Sarah
On T
Hans W Borchers googlemail.com> writes:
>
> Ben Bolker gmail.com> writes:
>
> >
> > Simulated annealing and other stochastic global optimization
> > methods are also possible solutions, although they may or may not
> > work better than the many-starting-points solution -- it depends
> > on
The posting guide asks for the output of sessionInfo()
And what does Sys.timezone() say (it isn't always helpful).
On Fri, 11 Nov 2011, Eduardo Mendes wrote:
Hello
Many thanks for the replies.
I am note sure whether you've got what you meant (Prof. Ripley) but here is the
output of Sys.g
Christofer Bogaso gmail.com> writes:
[...]
> Here my question is, is there any speed reduction if I put them within a
> function (I think there may be some speed reduction at least within
> for-loop, because that loop needs to call that function many times),
> relative to if I used that group of c
Hello
Many thanks for the replies.
I am note sure whether you've got what you meant (Prof. Ripley) but here is
the output of Sys.getlocale()
> Sys.getlocale()[1] "LC_COLLATE=English_United
> States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United
> States.1252;LC_NUMERIC=C;L
>
> Re: [R] Odp: Error in matrix, not ordered vectors or numerical value,
and SIAR.
>
> Thank-you Petr.
>
> I have consulted as many manuals and help pages, and search engines, and
> trying various things in SIAR, but continuous errors prevail..
>
> I have tried changing all matrices to num
Th CRAN taskview on clustering has many clustering and mixture packages
which
may be useful for your precise situation:
http://cran.r-project.org/web/views/Cluster.html
hth, Ingmar
On Thu, Nov 10, 2011 at 7:41 PM, Dan Abner wrote:
> Hello everyone,
>
> Can anyone suggest a decently documented
Ben Bolker gmail.com> writes:
>
> Simulated annealing and other stochastic global optimization
> methods are also possible solutions, although they may or may not
> work better than the many-starting-points solution -- it depends
> on the problem, and pretty much everything has to be tuned.
On Fri, 11 Nov 2011, Mario Valle wrote:
I just upgraded my Win7 32bits installation to 2.14.0 after deinstalling
2.12.x
First thing I moved the win-library from 2.12 to 2.14 and executed a
update.packages(ask='graphics',checkBuilt=TRUE)
(Swiss mirror).
This aborts with the console message:
Er
I just upgraded my Win7 32bits installation to 2.14.0 after deinstalling
2.12.x
First thing I moved the win-library from 2.12 to 2.14 and executed a
update.packages(ask='graphics',checkBuilt=TRUE)
(Swiss mirror).
This aborts with the console message:
Error in if (any(diff)) { : missing value wh
Thank you very much for your help.
I'm using the second suggestion in my program and it works very well.
Jonas
-Original Message-
From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Sent: Thu 11/10/2011 5:58 AM
To: Richter-Dumke, Jonas
Cc: r-help@r-project.org
Subject: Re: [R] mat
On Thu, 10 Nov 2011, Tyler Pirtle wrote:
Hi there,
I'm writing an R extension that has a C component that relies on two third
party libraries that I'm bundling
with the extension.
I'd like to statically link my resulting extension so that I can rely on
the bundled versions of the libraries I'm
baptiste auguie googlemail.com> writes:
>
> Dear list,
>
> [cross-posting from Stack Overflow where this question has remained
> unanswered for two weeks]
>
> I'd like to perform a numerical integration in one dimension,
>
> I = int_a^b f(x) dx
>
> where the integrand f: x in IR -> f(x) in I
Hi:
Try the log = "" argument to beanplot(). Here's an example:
library('beanplot')
v <- exp(runif(1000, 0, 10))
beanplot(v)
beanplot(v, log = '"")
HTH,
Dennis
On Thu, Nov 10, 2011 at 1:27 PM, Twila Moon wrote:
> Is it possible to force beanplot not to use a log scale? I want to be able to
>
On 11-11-10 1:28 PM, Gyanendra Pokharel wrote:
I have the following code,
gibbs<-function(m,theta = 0.25, lambda =0.55, n =1){
alpha<- 1.5
beta<- 1.5
gamma<- 1.5
x<- array(0,c(m+1, 3))
x[1,1]<- theta
x[1,2]<- lambda
x[1,3]<- n
for(t in 2:(m+1)){
x[
On 11-11-10 9:25 AM, ftonini wrote:
Hi everybody,
I started to receive a weird message in R that I have never seen
before...also I haven't found anything on google or on this forum about it.
Whenever I use the command source(...) to point to one of my scripts, I get
the following message:
Error
If you find yourself asking "is there a way..?", put that question mark in
front of the function you are wondering about and push return
?write.csv
In this case, you will see that the write function has a parameter to say
whether to write out row names ornot, which defaults to TRUE...
row.nam
The \r\n worked as the the diff can see no difference between my program's
file and the file generated from the windows computer.
For David,
I use Emacs whitespace mode to see the spaces and the \n's. The \r does not
show up in the whitespace mode. Maybe there is some way of turning it on.
but T
Hello all,
Have a little annoying issue with wreite.csv. It always writes out
observation number. Is there a way to avoid this. Thanks,
DL
--
View this message in context:
http://r.789695.n4.nabble.com/Have-write-csv-write-csv-without-observation-number-tp4030672p4030672.html
Sent from the R he
On Nov 11, 2011 at 1:06am Alex DJ wrote:
Alex,
I haven't followed this thread closely (now split), but this most recent
request for help from you is very difficult to make sense of. I think we
need access to your data set. Further, there appear to be problems with the
function you are using in pa
Thank-you Petr.
I have consulted as many manuals and help pages, and search engines, and
trying various things in SIAR, but continuous errors prevail..
I have tried changing all matrices to numeric, but the first column , of the
two latter matrices, returns as NA. Any other suggestions?
I had
82 matches
Mail list logo
