Note what are supposed to be quotation marks around "ID" in your post are "a
circumflex" instead. Maybe the problem is that the quotation marks that are
supposed to be around ID are not recognized by your R version. May you do
not use the proper font encoding. Use straight quotation marks or
apost
It likely means that your x and y are differently long. That is, affect1 and
adh1scr do not contain the same number of values in that instance. That
precludes them from being plotted against each other. abline and lowess
would fail for the same reason.
x<-c(1,2,3)
y<-c(2,4)
plot(y~x)
complete.case
Hi RUsers,
I am currently tying to use RExcel via R installed in my remote server
(CentOS 5.5).
Could you please help me set up a R remote server?
All I know about the server is its IPAddress.
Apologies for my fundamental question.
Thank you so much!
Taka
_
Dear R People: (particularly those who have build Rcmdr Plugin packages):
I'm building a new Plugin and keep getting the following error:
Error in if (is.null(where) || where >= n) rbind(object1, object2)
else if (where < :
missing value where TRUE/FALSE needed
Error in library(RcmdrPlugin
H
Dube, Jean-Pierre wrote:
>
> To whom it may concern,
>
> I am trying to maximize a log-likelihood function using optim. This is a
> simple problem with only 18 parameters. To conserve memory, I am using
> sparse matrices (SLAM) for some of the data matrices used in the
> computation of the lik
Something like this could work for you (note that a few other packages would be
faster, such as data.table), but plyr is intuitive.
require(plyr)
# example data set
dat <- data.frame(col1 = rep(c("a","b"), each=5),
col2 = rep(c("prog1","prog2","prog3","prog4","prog5"), 2),
col3 = rnorm(10), co
To whom it may concern,
I am trying to maximize a log-likelihood function using optim. This is a
simple problem with only 18 parameters. To conserve memory, I am using sparse
matrices (SLAM) for some of the data matrices used in the computation of the
likelihood. However, optim appears to co
Hi,
I am working on a heatmap and clustering all the rows with the method of
'ward'.
But I notice that my rows are in a weird order which is not what I want /
expect.
Could anyone help on this? My questions are with the figure below. Thanks
in advance
http://r.789695.n4.nabble.com/file/n360381
Hello,
I am working on a project to create some scatter plots. I have syntax
for 26 plots, and 22 of them display as they should. But here, for
example, is a sample of the command syntax I am using:
good <- complete.cases(affect1,adh1scr)
plot
(jitter(affect1,2.0),jitter(adh1scr,1.0),xli
Hi R Help list
I'm looking to visualise US foreign aid 1946-2009 and I have the dataset for
this. The trouble is it's a bit too complex and I need to simply it
I want to merge all of the rows with the same country together and add up the
individual totals to make one total figure per country
Dear R users,
I hope this email finds you well,
My name is Mariam and I am currently using R in my thesis
project. I is about modeling investors' sentiment.
My R skills are very modest and I am trying to solve a Garch in mean
equation using the "rgarch" package.
The main issue I am facing is with
the examples work great but it gives me this error on this file
http://r.789695.n4.nabble.com/file/n3603737/il1ra_L.csv il1ra_L.csv
> Age.co<-c(dat2b$Dr_Age)
> Sex.co<-as.factor(dat2b$Sex)
> casecont.co<-as.factor(dat2b$Self_T1D)
> stat.co<-as.factor(dat2b$status)
>
> m <- matchControls(Sex.co~A
Dear R-users,
I computed a simple coxph model and plotted survival over time with
associated confidence intervals for 2 covariate levels (males and
females).
M1 <- coxph(survobject~sex, data=surv)
M1
survsex <- survfit(survobject~sex,data=surv)
summary(survsex)
plot(survsex, conf.int=T, col=c("b
Sorry, Patrick and other victims. My bad~~
The data file can be accessed at
http://zhanglab.ccmb.med.umich.edu/TM-fold/featureSelection_rec.vectors
The code is attached below. Error happens at "data.cforest.varimp <-
varimp(data.cforest, conditional = TRUE)": "Error in
model.matrix.default
I am still thinking about this problem. The solution could look something
like this (it's net yet working):
k<-lapply(h, function (x) x*0) # I keep the same format as h, but set all
values to 0
years<-c(1997:1999) # I define the years
for (t in 1:length(years))
{
year = as.characte
I am new to R and am having difficulty merging 2 data sets, both of which have
the same 30 variables and about 2,000 records. I've looked high and lowâPaul
Johnson's help page, through some of the 5,420 (yikes!) results for âmergeâ
on the R list archives, and google searchesâand I'm ge
I made different assumptions than Josh.
xf<-c("W",NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA,NA)
xg<-c(NA,"k","h",NA,"g","r","j",NA,"v","d",NA,"v",NA,"z","r","r","i")
xf
xg
unlist(apply(cbind(xf,xg), 1, function(x) x[!is.na(x)]))
as.vector(unlist(apply(cbind(xf,xg), 1, function(x) x[!is.na(x)]))
Hi,
You did not specify what assumptions you can make about xf and xg
(such as will they have identical lengths and is it possible both
could contain nonmissing values in the same element?), but this seems
a straightforward approach in your little example:
index <- is.na(xf)
xf[index] <- xg[index
Dear R People:
Suppose I have the following two character vectors:
xf
[1] "W" NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
> xg
[1] NA "k" "h" NA "g" "r" "j" NA "v" "d" NA "v" NA "z" "r" "r" "i"
>
I want to end up with
"W" "k" "h" ...
What is the best way to achieve t
Without seeing the model I'm not clear on the cause. In general it is not
very useful to compute 1 d.f. tests for parameters that represent only part
of a meaningful hypothesis.
Note that anova.lrm computes Wald statistics, which are subject to the
Hauck-Donner effect (which may not be the cause
caret has many useful preprocessing facilities for numeric predictors,
such as removing near-zero-variance predictors, removing correlated
predictors, removing linear combinations, centering/scaling, and
Box-Cox transformations. However, these don't work for factor data,
which is common in many le
Folks:
This is sufficiently "tricky" that a word of explanation might be
helpful to reveal the magic.
>
> Rates <- c(0.006280048, 0.330934659)
> Index.Matrix <- matrix(c(1, 2, 1, 2), 2, 2)
>
>> Rates
> [1] 0.006280048 0.330934659
>
>> Index.Matrix
> [,1] [,2]
> [1,] 1 1
> [2,] 2 2
On Jun 16, 2011, at 4:52 PM, Jeremy Beaulieu wrote:
> Hi all~
>
> I have a bit of stupid question. And I promise, I have struggled through this
> and it sort of pains me to post this question...
>
> I have the following vector which is are the estimated parameter values
> output from "optim".
Hi all~
I have a bit of stupid question. And I promise, I have struggled through this
and it sort of pains me to post this question...
I have the following vector which is are the estimated parameter values output
from "optim".
$Rates
[1] 0.006280048 0.330934659
I want to match the order in
Hi all,
I have the following script which fills the values which are less than
the mean of a given timeseries.
If you look closely, the colored regions are "out of line".
Any suggestions how I can rectify this?
Thanks
Muhammad
# -
#rm(list=ls())
x <- abs(rnorm(100))
tt <- 1:100
m <- m
Hi,
As far as I know secondary y-axis and multiple pages are not possible
in ggplot2 (there are workarounds for the latter in the ggplot2 list
archives). For the subtitle, you could implement it with grid.text and
grid viewports,
library(gridExtra)
library(ggplot2)
grid.arrange( qplot(1,1), sub
Greetings to the help mailing list.
I am in the process of translating a large graphic from xyplot to ggplot2 (13
columns by about rows).
I have been unsuccessful trying to understand how to place the following text
strings after three days of:
* Perusing Wickham's "gg2plot" book;
* Searching hi
Dear List,
I wanted to calculate the asymmetrical confidence interval based on the sample
statistic and standard error that available from the published report (complex
survey-based).
The calculation details can be seen from pages 17-18 of the document at the
following link: http://www.oas.sa
You could do this using permutation tests. If everything is balanced and
orthogonal then you can permute the predictors, otherwise you can follow these
main steps:
1. Formulate the test of interest as a full and reduced model test and find the
F statistic (or other, but I will assume F) for th
I don't think that this approach is appropriate here. Each iteration after the
1st the lm/predict combination will assume that the new data is exact when in
fact it is an estimate with some error involved. To properly do this you need
to take into account that variability. There is a time ser
On Thu, 16-Jun-2011 at 09:19AM -0400, Jinrui Xu wrote:
> Thanks for your feedback.
> I think the problem is not because of many levels. There is only 1
> column with two levels as class labels in my input data.
>
> Below is my code. The commandline "data.cforest.varimp <-
> varimp(data.cforest
On Thu, Jun 16, 2011 at 3:43 PM, Prof Brian Ripley
wrote:
> On Thu, 16 Jun 2011, Sarah Goslee wrote:
>
>> Hi Frank,
>>
>> It's not inaccuracy, exactly, and it's intentional.
>>
>> The default behavior is to use a round endcap on a line, and with such
>> a large linewidth that endcap is arcing belo
On Thu, 16 Jun 2011, Sarah Goslee wrote:
Hi Frank,
It's not inaccuracy, exactly, and it's intentional.
The default behavior is to use a round endcap on a line, and with such
a large linewidth that endcap is arcing below and above the starting
points.
This being R, there's a par for that: putt
Hi,
If 100 frames can work but 730 cannot, I think this might be due to
the limit of ImageMagick. If that is the case, I do not have any
solutions from R's aspect, instead, you may want to try
GraphicsMagick, which is similar to ImageMagick. I'm not sure if it is
better than ImageMagick in terms o
Hi Frank,
It's not inaccuracy, exactly, and it's intentional.
The default behavior is to use a round endcap on a line, and with such
a large linewidth that endcap is arcing below and above the starting
points.
This being R, there's a par for that: putting lend=1 in your lines()
statements will g
Using a line width > 1 results in not only a thicker line but also some fuzz
in the other direction, as shown in this example program. You will see that
the thick vertical black lines extend below the gray scale horizontal lines.
Does anyone know whether this is intended or is it a bug? The appl
Kristian,
> Another update: I ditched the pre-installed version of openmpi and
> installed openmpi from source (./configure --without-openib).
A very reasonable approach. Remembering my own struggles with SuSE years
ago, I did not dare to propose that procedure.
> Compiling and running your tes
Hi,
my problem sounds simple but I've been trying and looking for a solution,
but just can not find out how to solve it.
I have a dataset and want to do a GAM with binomial link function, random
effects (for individuals as there are severeal measurements for everyone)
and a time-spline.
For thi
Bill,
Thank you for the correction; yes, Nevilles' is an algorithm to produce the
same Lagrange interpolating polynomial. And thank you for the polynomF
package!
Don Slowik
--
View this message in context:
http://r.789695.n4.nabble.com/functions-for-polynomial-and-rational-interplation-tp359533
Hi All - I am using the daisy function from the cluster library to create a
dissimilarity matrix. I'm going to use that matrix to run a cluster
analysis. My participants are identified with the variable, hhid. However,
when I try to keep hhid in the dataset that I use to create the
dissimilarity
Dear members,
I'm fitting linear model using "lm" which has numerous auto-regressive terms as
well as other explanatory variables. In order to calculate prediction
intervals, i've used a for-loop as the auto-regressive parameters need to be
updated each time so that a new forecast and corresp
Hi,
I am working under Linux and have an executable R script like this:
#! /usr/bin/env Rscript
Non- Parametric equivalent to parametric anova is Kruskal-Wallis test - I am
not sure R has any function for this - I am new to R (I think R should have it
already) - SAS has this option.
Sharif
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org
Petr Savicky-2, thanks you. It seems the code works very well. I use rate=1,
so "common" is between 0 to 1. Thanks again.
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-generate-bivariate-exponential-distribution-tp3596944p3602539.html
Sent from the R help mailing list arch
Hi,
If you truly have an array, this is option that should be much faster
than a loop:
index <- which(is.na(dat))
dat[index] <- dat[index - 1]
the only catch is that when there previous value is NA, you may have
to go through the process a few times to get them all. One way to
automate this wou
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
On 06/16/2011 01:47 PM, Rob James wrote:
> Ben,
>
> Thanks for this. Very helpful and clearly others have tripped over the
> same problem
> I would have supposed that the solution was to ask lrm (or glm) to use
> LR rather than Wald, but I don't see s
Dear list,
I have read several posts on this topic. I would use the same
methodology as proposed
by Simon Wood in this post:
http://r.789695.n4.nabble.com/variance-explained-by-each-term-in-a-GAM-td836513.html
My first question is:
Does anyone know a scientific source (paper, book,...) that exp
Do you mean prompts (the > that indicates R is waiting for input)?
Your issue is not clear. Please post a copy/paste of your session, as
well as your sessionInfo().
On Thu, Jun 16, 2011 at 7:16 AM, Rosario Garcia Gil
wrote:
> Hello
>
> I have a annoying problem with R (which I am running in MAC).
Dear R.oo package users,
while testing some functionality of the R.oo, I found that during the first
construction of a object from a class, the constructor is twice called, but
only one object is finalized.
In all subsequent creation processes, the constructor is (expectedly) called
once.
Here
Hello
I have a annoying problem with R (which I am running in MAC).
Every time I enter a function and I press enter to run it, I get many
additional "<´s", it is really bothering me. I have tried to google to find out
how to fix the problem but I did not succeed. Anyone with similar experience
Jannis,
Beautiful! Thanks for the useful responses from others as well.
Regards,
Phil.
On Thu, 16 Jun 2011 12:49:58 +0100 (BST), Jannis wrote:
Try something like:
par(mar=c(4,2,1,12),tcl=0.2,mgp=c(1,0,0))
plot(1:10,1:10)
for (i in 1:7)
{
par(new=TRUE)
plot(1:10,(i+(1:10)+rnorm(10,0
I got an array similar to the one below, and want to replace all NAs with the
previous value.
99 8.2 b
NA 8.3 x
NA 7.9 x
98 8.1 b
NA 7.7 x
99 9.3 b
...
i.e. the first two NAs should be replaced to 99, whereas the last one should
be 98.
I would like to apply a function to reach row, checking if th
Hi everybody,
I just tried to run R on one of my projects
but it did not want to run:
R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-pc-linux-gnu (64-bit)
Loading required package: utils
Error in loadNamespace(i
Hi all,
What is the best package to use to evaluate correlations between events
by time (species co-occurrence).
Chron?
I would like to determine fi there is any solid relationship between the
occurrence of critter A given critter B is present at a given time frame.
Bruce
Greg Snow imail.org> writes:
> The individual tests on coefficients in logistic regression
> are generally based on a Wald test statistic.
> Unfortunately there is a bit of a paradox possible in this
> case where the coefficient is highly
> significant, but due to a flattening of the likelihoo
Those are in fact the coefficients for p and q they are estimating, though
their M is different. Who knows what they did with that.
The data source was:
Title: Diffusion models of mobile telephony in Greece
Source: Telecommunications policy [0308-5961] Michalakelis yr:2008 vol:32
iss:3-4 pg:234
I am not sure what you actually want to accomplish. If you try to correlate
numeric/rank data with text data, I have no clue how that could be achieved
(other than with text length, or presence/absence of a comment). If you try
to figure out how you prevent the cor() function to fail when it encoun
Hi
>
> Dear R users,
> I am using package alphahull version 0.2-0 with R-2.13.0. Here is my
problem.
> I have produced a 2-D plot of points and need to draw convex hulls
around
> some subsets of these points. So I just need to add the convex hulls to
> the already existing plot. That is why I
On 06/16/2011 03:04 AM, siddharth arun wrote:
>
> I will break my query into 2 parts:
>
> 1. My main idea is to give some data as an input and if any seasonality
> is there in the input, it should be captured automatically and the
> corresponding frequency automatically generated. please comment
Dear R users,
I am using package alphahull version 0.2-0 with R-2.13.0. Here is my problem.
I have produced a 2-D plot of points and need to draw convex hulls around some
subsets of these points. So I just need to add the convex hulls to the already
existing plot. That is why I cannot use the pl
On 16 June 2011 at 07:17, Frans Marcelissen wrote:
|
| Hello,
| For one reason or another, I cannot install Rcmdr on my R 2.12.2 ubuntu 11.4
system.
| I run R as root When I try to install Rcmdr, then whatever mirror I select, I
get the message:
Rcmd is an Ubuntu package too, so just do
$ s
On Thu, 16 Jun 2011, Mr Rhelp wrote:
Hi all,
I am trying to re-compile some "unmaintained" (it seems) packages, namely rSoNIA
and dynamicnetwork from:
http://csde.washington.edu/~skyebend/rsonia/rsoniaDemo/
These packages predates R 2.10.0 so they need to be recompile.
Well, re-installed.
A
On 15.06.2011 15:58, Franc Lucas wrote:
Hello,
I am wondering how to get the quotation marks into a variable expression. I
can't escape it with the backslash \ ...
Example:
I can access my data frame via
TABLE$"2011-01-02"$columnD
Now I want to do this automatically
bbolker wrote:
>
> biologie.uni-marburg.de> writes:
>
>> you answerd to Nancy Shackelford about Clarks 2Dt function.
>> Since the thread ended just after your reply,
>> I would like to ask, if you have an idea how to use this function in R
>
> I got started on your problem, but I didn't finish
bbolker wrote:
>
> biologie.uni-marburg.de> writes:
>
>> you answerd to Nancy Shackelford about Clarks 2Dt function.
>> Since the thread ended just after your reply,
>> I would like to ask, if you have an idea how to use this function in R
>
> I got started on your problem, but I didn't finish
On 16.06.2011 14:46, Ben Rhelp wrote:
Hi all,
I am trying to re-compile some "unmaintained" (it seems) packages, namely rSoNIA
and dynamicnetwork from:
http://csde.washington.edu/~skyebend/rsonia/rsoniaDemo/
These packages predates R 2.10.0 so they need to be recompile.
After split the single
On Thu, 16 Jun 2011, Joel wrote:
Found what was wrong with my code
Please follow the posting guide, and give the context asked for there.
changed:
inFile <- file('gwas_data_chr10.gen')
to:
inFile <- file('gwas_data_chr10.gen','r')
Dunno why its important but it is :P
You didn't open the f
The script command "coda" is not necessary, unless you want to read in in
coda-format output. Necessary are the other commands: "history & save
(odc)". "save (txt)" lets you read the on-screen output of winbugs later as
a text-file.
Another way is also to take the saved last values of the chains a
Found what was wrong with my code
changed:
inFile <- file('gwas_data_chr10.gen')
to:
inFile <- file('gwas_data_chr10.gen','r')
Dunno why its important but it is :P
--
View this message in context:
http://r.789695.n4.nabble.com/Read-file-line-by-line-tp3602349p3602399.html
Sent from the R help m
Hi all,
I am trying to re-compile some "unmaintained" (it seems) packages, namely
rSoNIA
and dynamicnetwork from:
http://csde.washington.edu/~skyebend/rsonia/rsoniaDemo/
These packages predates R 2.10.0 so they need to be recompile.
After split the single big file in /man in each packages into
Hi
I got a file that looks like this: (i have shorten it alot the real file is
over 200k rows long)
chr10 rs7909677 101955 A G 0 1 0 1 0 0...
chr10 rs2436025 238506 C G 1 0 0 1 0 0...
chr10 rs11253562 148946 C T 0 1 0 0 1 ...
chr10 rs1105116 230788 G T 0 0 1 0 0 1...
chr10 rs4881551 149076 A G 0
On 11-06-16 7:50 AM, christiaan pauw wrote:
Hi everybody.
Was the problem below ever solved? I have the same problem on Windows 7 with
R 2.13.0 and LyX 2 - everything freshly installed. My Rweave.bat and
MakeSweave.R files are similar to below (I followed
http://wiki.lyx.org/uploads/LaTeX/Sweave
Hi everybody.
Was the problem below ever solved? I have the same problem on Windows 7 with
R 2.13.0 and LyX 2 - everything freshly installed. My Rweave.bat and
MakeSweave.R files are similar to below (I followed
http://wiki.lyx.org/uploads/LaTeX/SweaveR/lyx_with_r_and_sweave_instructions.pdf
)
I a
Try something like:
par(mar=c(4,2,1,12),tcl=0.2,mgp=c(1,0,0))
plot(1:10,1:10)
for (i in 1:7)
{
par(new=TRUE)
plot(1:10,(i+(1:10)+rnorm(10,0,0.4)),yaxt='n',col=i,type='l')
axis(4,line=(i-1)*1.5,col=i-1,col.ticks=i-1)
}
You may need to tweek around the options a bit and check that the t
On 16.06.2011 10:25, Animesh wrote:
I am trying to optimize a simple function, but I want to optimize it for
values of k between 0 and 1.
I read the help file for constrOptim function but it is always giving me
the error
constrOptim((0.4),simulation,NULL,ui=rbind(c(0,0),c(1,1)))
Your param
Another update: I ditched the pre-installed version of openmpi and
installed openmpi from source (./configure --without-openib).
Compiling and running your test programme gives me the following error:
unger@zytosrv01dmi:~/R_projects/OS> ./mtest
[zytosrv01dmi:6094] *** An error occurred in MPI_Sen
get(a)[[b]][[c]]
with
b <- "2011-01-02"
Sent from my iPad
On Jun 15, 2011, at 9:58, "Franc Lucas" wrote:
>
> Hello,
> I am wondering how to get the quotation marks into a variable expression. I
> can't escape it with the backslash \ ...
> Example:
> I can access my data frame via
On 06/16/2011 09:23 AM, Jim Lemon wrote:
> On 06/16/2011 07:12 PM, amrita gs wrote:
>> Hi everyone,
>>
>> I want to do data analysis using the values retrieved from a MySQL
>> database. Is there a way to save the results of data analysis in R.
>>
> Hi amrita,
> If you mean "save the output of
On 06/16/2011 07:12 PM, Philip Rhoades wrote:
People,
I want to create seven plots on the one graph - each with it's own Y
axis and scale - I want the plots overlaid because it is the shape of
the graphs that I want to compare, not so much the values. I have found
examples where people show how
On 06/16/2011 07:12 PM, amrita gs wrote:
Hi everyone,
I want to do data analysis using the values retrieved from a MySQL
database. Is there a way to save the results of data analysis in R.
Hi amrita,
If you mean "save the output of an analysis", there are several ways.
1) Enter "sink() b
Dear Philip,
If the values are not important, then you don't need different Y-axes. Why not
standardise the seven datasets so they have a common scale?
Best regards,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur-
People,
I want to create seven plots on the one graph - each with it's own Y
axis and scale - I want the plots overlaid because it is the shape of
the graphs that I want to compare, not so much the values. I have found
examples where people show how to put two Y axes on the left of the plot
Hi everyone,
I want to do data analysis using the values retrieved from a MySQL
database. Is there a way to save the results of data analysis in R.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.
On Thu, 16 Jun 2011, Frans Marcelissen wrote:
Hello,
For one reason or another, I cannot install Rcmdr on my R 2.12.2 ubuntu 11.4
system.
I run R as root When I try to install Rcmdr, then whatever mirror I select, I
get the message:
install.packages('Rcmdr',dependencies=TRUE)
Installing pack
Hi everyone,
I retrieved the values of form elements such as radiobuttons, text
area,textbox etc into MySQL database. I retieved the values from MySQL db to
R using RODBC package. If i want to do a corelational analysis using the
values in textarea and radiobutton how is it possible as textarea m
On 06/15/2011 09:51 PM, bjmjarrett wrote:
Hi all,
I have scoured the archives of this forum but nothing quite seems to fit the
bill...
I would like to plot a graph displaying two variables (y axes) that share
date as the x axis. However, the date values for each variable are not the
same - for
Dear Sir,
I would like to ask whether there is any tool in R, developed for the function
of non parametric Anova, where the non parametric analysis is able to compute
the p-value for interaction as well (similar to Anova)? If not, would sir like
to suggest any of the other statistical software
Hello,
For one reason or another, I cannot install Rcmdr on my R 2.12.2 ubuntu 11.4
system.
I run R as root When I try to install Rcmdr, then whatever mirror I select, I
get the message:
install.packages('Rcmdr',dependencies=TRUE)
Installing package(s) into /usr/local/lib/R/site-library
(as
I am trying to optimize a simple function, but I want to optimize it for
values of k between 0 and 1.
I read the help file for constrOptim function but it is always giving me
the error
> constrOptim((0.4),simulation,NULL,ui=rbind(c(0,0),c(1,1)))
"Error in ui %*% theta : non-conformable arguments
Hello everyone,
I use the following command lines to get important variable from training
dataset.
data.controls <- cforest_unbiased(ntree=500, mtry=3)
data.cforest <- cforest(V1~.,data=rawinput,controls=data.controls)
data.cforest.varimp <- varimp(data.cforest, conditional = TRUE)
I got er
Dear all,
I am estimating a bivariate probit model using the package VGAM:
fit1 = vglm(cbind(y1, y2) ~ x1+x2, binom2.rho, data=dat, trace=TRUE)
I would like to estimate this via the method vglm.fit (in an analogous way
as lm.fit for lm or glm.fit for glm). Unfortunately my trials did no
Use the `xlim` parameter of the `plot` function.
IFor example, if X1 and X2 are the variables that contain the values of your
x-axis,
do xlim1 <- range(X1,X2),
than do plot(, xlim = xlim1,...)
Hope this helps.
2011/6/15 bjmjarrett
> Hi all,
>
> I have scoured the archives of this foru
I will break my query into 2 parts:
1. My main idea is to give some data as an input and if any seasonality is
there in the input, it should be captured automatically and the
corresponding frequency automatically generated. please comment whether this
is possible with auto.arima()
2. If frequenc
93 matches
Mail list logo
