I wrote a function and I hope this helps
b<-c(1,2,3,4);lambda<-c(0.3,0.4,0.5,0.6)
m<-c(0.10,0.11,0.12,0.13);t<-c(7,8,9,10)
N=4
Y<-matrix(data=NA,nrow=4, ncol=1)
for (i in 1:N){
y=b*(1-exp(lambda*(t^m)))
Y<-matrix(data=y, nrow=4, ncol=1)
Y<-as.numeric(Y)
sim<-data.frame(t, Y);plot(s
I saw a similar question but I still don't fully understand how to implement
optim.
Can someone help me out with this?
Thanks. Keun-Hyung
> vol<-rep(c(0.03, 0.5, 2, 4, 8, 16, 32), 3)
> time<-rep(c(2,4,8),each=7)
> p.mated<-c(0.47, 0.48, 0.43, 0.43, 0.26, 0.23, "null", 0.68, 0.62, 0.64,
0.58, 0.53,
Hello..!!!
I got a problem while built R program with tcltk. I try to list the name of
all data frame in my memory by "tklistbox". when we choose the name of data
frame, then i capture the name, and catch the data frame by
tclObj(dat)<-as.data.frame(get(x,envir=.GlobalEnvir))
tclvalue(
hi i have to make my function
y=b*(1-exp(lambda*(t^m))
and my data are in excel file with names:
Year /b/lambda/m
and t is an age in another file calls pcb
could u help me please to calculate all y for all years please
i can make this step by step but i want an program for having all y and plot
the
Hi the list
I do not manage to build a package using package.skeleton.
My pacakge in called 'partition'
I have no probleme to source the file. Then I run
package.skeleton("partition")
I get :
Saving function and data...
Error in file(file,ifelse(append,"a","w")):
impossible to open the conne
Anyone know what's up w/ tolstoi.newcast.edu.au (home of the threaded
R-help archives)? I've got a "can't connect" for several days now.
Carl
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
> I have a question on how to identify the overlapped items from two vectors:
>
> >x<-c(0,1,2)
> >y<-c(1,2,2,3)
intersect(x,y)
> and plot the number of the overlapped as well as non-overlapped in a
> diagram?
Are you talking about a Venn diagram? There is a vennDiagram function in
the limma pac
On Sat, 31 May 2008, Powell, Jeff wrote:
Dear All,
Thanks to an answer which I received from a previous post, I'm now able to
create a series of
environments using the following:
nmes <- c("en1", "en2", "en3")
for(i in nmes) assign(i, new.env(parent = .GlobalEnv))
My next question is how,
> I have 2 vectors x & z
> x has 10 variables in it
> z has 75 variable in it
>
> I need to check all 75 variables in z and see if they are present in the
> vector x.
In addition to the solution using %in% which has already been suggested,
you may also want to play with something like this:
# m
Hello,
On May 31, 2008, at 4:10 PM, ss wrote:
Hi list,
I have a question on how to identify the overlapped items from two
vectors:
x<-c(0,1,2)
y<-c(1,2,2,3)
You can use the %in% operator to return a logical vector that
indicates the overlap.
> x<-c(0,1,2)
> y<-c(1,2,2,3)
> y %in%
I have two examples that I run hclust on:
a = c(0,1,1.5,1.5)
b = c(1,0,1.5,1.5)
c = c(1.5,1.5,0,0.5)
d = c(1.5,1.5,0.5,0)
ll = as.matrix(rbind(a,b,c,d))
test = as.dist(ll)
long = hclust(test)
a = c(0,0.3,1,1)
b = c(0.3,0,1,1)
c = c(1,1,0,0.5)
d = c(1,1,0.5,0)
ll = as.matrix(rbind(a,b,c,d))
test
I have 2 vectors x & z
x has 10 variables in it
z has 75 variable in it
I need to check all 75 variables in z and see if they are present in the
vector x.
So far this is what I have:
for (i in 0:75){
if (z[i]==x) (w=w+1)
}
Thanks in advance
w
--
View this message in context:
http://www.nabbl
Hello everyone!
Quite a similar question has been posed here some time ago, but there was no
explicit solution offered. So I hope that it is OK to pose it again.
I want to perform a KPSS test using the packages "urca" or "tseries". But I
neither want to use the predefined lag structures, "short"
darthgervais wrote:
>
> I have 2 vectors x & z
> x has 10 variables in it
> z has 75 variable in it
>
> I need to check all 75 variables in z and see if they are present in the
> vector x.
> So far this is what I have:
>
> for (i in 0:75){
> if (z[i]==x) (w=w+1)
> }
>
> Thanks in advance
>
Dear All,
Thanks to an answer which I received from a previous post, I'm now able to
create a series of
environments using the following:
nmes <- c("en1", "en2", "en3")
for(i in nmes) assign(i, new.env(parent = .GlobalEnv))
My next question is how, using "load", can I automatically place d
Is this ok?
set.seed(1)
# ... your code ...
legend("topright", title = "Legend", legend=c("observed density",
"normal with \nobserved mean\n& sd", "normal with \n'true' mean & sd"),
lty = 1:3, col = 1:3, box.lty = 0, text.width = 1.1, cex = 0.7)
On Sat, May 31, 2008 at 3:06 PM, Paul Johnson <[EM
Hi list,
I have a question on how to identify the overlapped items from two vectors:
>x<-c(0,1,2)
>y<-c(1,2,2,3)
and plot the number of the overlapped as well as non-overlapped in a
diagram?
Thanks much,
Allen
[[alternative HTML version deleted]]
_
Paul Johnson wrote:
Hello, everybody:
I recently encountered an example with in which the graph was placed
in a way that did not leave room for a legend. Maybe you would
describe it as "legend too big", I'm not sure. I found myself wishing
I could force in some space after the title.
Here's
Hello, everybody:
I recently encountered an example with in which the graph was placed
in a way that did not leave room for a legend. Maybe you would
describe it as "legend too big", I'm not sure. I found myself wishing
I could force in some space after the title.
Here's working example where
Paul,
On 30 May 2008 at 15:47, Paul Hewson wrote:
| Hello,
|
| We have R working with Rmpi/openmpi, but I'm a little worried.
Specifically, (a) the -np flag doesn't seem to override the hostfile (it works
fine with fortran hello world) and (b) I appear to have twice as many processes
runnin
Hello all,
I'm working on a project that uses a SOM to generate feature maps from
image data. Each image is a 100x75 pixels and RGBA, so 30,000 elements
per image.
I wanted to observe the structure of the pixel by pixel euclidean
distance (where each image is a point in 30,000 dimensional space)
[EMAIL PROTECTED] wrote:
Hi the list,
I do not manage to build a package using package.skeleton.
My pacakge in called 'partition'
I run package.skeleton, then I write the Rd file.
R CMD check partition is ok
But R CMD build --binary partition does not compile. I get :
WARNING
* some HT
Hi the list,
I do not manage to build a package using package.skeleton.
My pacakge in called 'partition'
I run package.skeleton, then I write the Rd file.
R CMD check partition is ok
But R CMD build --binary partition does not compile. I get :
WARNING
* some HTML links may not be found
in
Hello,
I am looking for advice on how to adaptively sample from a simulated population.
I have simulated a population of individuals in a matrix and need to do the
following:
1) select a random sample from that population
2) apply a function to that sample to determine if I detect any individual
Dear all:
I am searching for a regression model (e.g. y=Xb+e) in which dummy-coded events
(to time point t) on different regressors in X exhibit an effect on subsequent
responses in the vector y (to time-points t+1, t+2,… t+n). My aim is to
estimate how long the memory effect is and how strong
The zoo package has a yearqtr class that can represent quarterly
data. It is also possible to create zoo series of quarterly data. Such
series can be converted to ts class and that class can some
builtin support for quarterly data as well.
> library(zoo)
> as.yearqtr("200706", "%Y%m")
[1] "2007 Q
Robert Felty indiana.edu> writes:
>
> I am extremely puzzled by this behavior in R. I have a data frame called
> Trials in which I have results from an experiment. I am trying to do a
> subjects analysis, but getting weird results. Each row has 1 trial in it,
> which includes a column for the
I have financial data on a a set of firms, with a quarterly period
(fundamental data). The data spans 10 years, and four quarters per
year. The present file (.csv) reads the Date columns as "200706" for
the second quarter of 2007; "199809" for the third quarter of 1997.
Is there a way I can conver
28 matches
Mail list logo
