Displaying 20 results from an estimated 639 matches for "m3".
Did you mean:
mp3
2007 Feb 11
3
merge words=data name
I would like to merge two parts of words to get a name of the data. First
M3$N (invariable) and second is a number from 0001 to 3003 -
M3$N0001,M3$N0002,...,M3$N3003. For example if I do it like this:
my.data <- paste("M3$N",2456,sep="")
I get
> my.data
[1] "M3$N2456"
But I want to get something equivalent to
my.data<- M3$N2456
Is the...
2012 Sep 25
1
REML - quasipoisson
...quasipoisson with scale = 1
#fitting is identical to the poisson case
#F1, F3, and F4 unchanged but F2 is now undefined
m2<-gam(y~s(x1)+s(x2),family="quasipoisson",method="REML",scale=1)
F1+F3-F4
m2$gcv
#reml score = formula with F2 omitted
##quasipoisson with unknown scale
m3<-gam(y~s(x1)+s(x2),family="quasipoisson",method="REML",scale=-1)
phiq<-m3$scale
#1st term
S1q<-m3$smooth[[1]]$S[[1]]*m3$sp[1]
S2q<-m3$smooth[[2]]$S[[1]]*m3$sp[2]
Sq<-matrix(0,19,19)
for (i in 2:10)
{
for (j in 2:10)
{
Sq[i,j]=S1q[i-1,j-1]
Sq[i+9,j+9]=S2q[i-1,j-1]...
2013 Oct 11
3
[LLVMdev] Generate code for ARM Cortex m0, m3, and m4.
Hi,
I am trying to cross compile code for ARM Cortex m0, m3, and m4.
For m0, I use:
-target armv6--eabi -mcpu=cortex-m0
That seems to work. For m3 and m4, I use the following which does not work
(fatal error: error in backend: CPU: 'cortex-m3' does not support ARM mode):
-target armv7m--eabi -mcpu=cortex-m3
and
-target armv7em--eabi -mc...
2007 Jul 26
3
substituting dots in the names of the columns (sub, gsub, regexpr)
Dear R users,
I have the following two problems, related to the function sub, grep,
regexpr and similia.
The header of the file(s) I have to import is like this.
c("y (m)", "BD (g/cm3)", "PR (Mpa)", "Ks (m/s)", "SP g./g.", "P
(m3/m3)", "theta1 (g/g)", "theta2 (g/g)", "AWC (g/g)")
To get rid of spaces and symbols in the names of the columns,
I use read.table(... check.names=TRUE) and I get:
str <- c(&...
2012 Nov 09
6
(sin asunto)
Saludos, tengo un problema que no puedo resolver dentro del R
Estoy creando una funcion que a partir de un objeto compuesto por diferentes matrices como el que esta a continuacion:
[[1]]
M1 M2 M3 M4
sp1 2 0 1 8
sp2 4 5 2 4
sp3 0 0 4 0
sp4 5 7 5 0
sp5 0 4 0 0
[[2]]
M3 M2 M4 M1
sp1 1 0 8 2
sp2 2 5 4 4
sp3 4 0 0 0
sp4 5 7 0 5
sp5 0 4 0 0
[[3]]
M4 M3 M1 M2
sp1 8 1 2 0
sp2 4 2 4 5
sp3 0 4 0 0
sp4 0 5 5 7
sp5 0 0 0 4
[[4]]
M3...
2009 Aug 17
1
[Fwd: Re: R code to reproduce (while studying) Bates & Watts 1988]]]
...= df.Chloride,
> start = list(
> t1 = 35,
> t2 = 0.91,
> k = 0.22))
So my error was to use nls instead that gnls. Thanks a lot, Kevin.
> summary(m2)
> plot(m2)
> lag.plot(resid(m2), do.lines=FALSE)
> acf(resid(m2))
>
> m3 <- update(m2, corr=corAR1(.67))
> summary(m3)
> plot(m3)
> lag.plot(resid(m3), do.lines=FALSE)
> acf(resid(m3))
>
> The residual plots for model m3 still show structure, unlike in Bates
> & Watts, so maybe this is not the correct model?
>
> Kevin
Actually is not e...
2011 Dec 03
1
partial mantel tests in ecodist with intential NA values.
...not square error.
is it possible to perform this test in ecodist?
many thanks
Nevil Amos
> library(ecodist)
> library(ncf)
> x<-sample(1:1000,20)
> y<-sample(1:1000,20)
> z<-sample(1:1000,20)
> M1<-as.matrix( distance(x))
> M2 <-as.matrix( distance(y ))
> M3<-as.matrix( distance(z ))
> D1<-(lower(M1))
> D2<-(lower(M2))
> D3<-(lower(M3))
> mantel(D1 ~ D2+D3, nperm=1000)
mantelr pval1 pval2 pval3 llim.2.5% ulim.97.5%
0.09014696 0.10300000 0.89800000 0.18400000 0.01857311 0.18468621
> partial.mantel.test(M1,...
2012 Oct 01
0
[Fwd: REML - quasipoisson]
...ntical to the poisson case
> #F1, F3, and F4 unchanged but F2 is now undefined
>
> m2<-gam(y~s(x1)+s(x2),family="quasipoisson",method="REML",scale=1)
> F1+F3-F4
> m2$gcv
> #reml score = formula with F2 omitted
>
> ##quasipoisson with unknown scale
> m3<-gam(y~s(x1)+s(x2),family="quasipoisson",method="REML",scale=-1)
> phiq<-m3$scale
>
> #1st term
> S1q<-m3$smooth[[1]]$S[[1]]*m3$sp[1]
> S2q<-m3$smooth[[2]]$S[[1]]*m3$sp[2]
> Sq<-matrix(0,19,19)
> for (i in 2:10)
> {
> for (j in 2:10)
>...
2012 Jul 18
2
[LLVMdev] Setting up a cross-compiler for cortex-m3
On 18 July 2012 14:33, salvatore benedetto
<salvatore.benedetto at gmail.com> wrote:
> but I still haven't figure out how to build for cortex-m3
>
> clang -march=armv7-m -mfloat-abi=soft <something missing?> testReference.cpp -c
-march should have done the trick.
You can also try -mcpu=cortex-m3,
or try -ccc-host-triple armv7m-none-gnueabi (or -eabi),
and possibly -ccc-gcc-name code sourcery's gcc,
or possibly --sysroot...
2011 Mar 08
1
SEM error
Dear All,
I am new for R and SEM. I try to fit the model with Y (ordinal outcome), X
(4 categorical data), M1-M3 (continuous), and 2 covariates (Age&sex) as a
diagram.
library(polycor)
model.ly <-specify.model()
1: x -> m1, gam11, NA
2: x -> m2, gam12, NA
3: x -> m3, gam13, NA
4: age -> m1, gam14, NA
5: age -> m2, gam15, NA
6: age -> m3, gam16, NA
7: sex -> m1, gam17, NA
8:...
2010 Oct 24
1
Optimize parameters of ODE Problem which is solved numeric
...651463 1.436980 7
4845 6 1.692710 1.267359 4
4552 5 1.686448 1.220539 6
4282 6 1.579868 1.086296 6
75441 6 2.978580 1.338100 7
I want to solve the following system of ode (ord. differential equations) numerically (f.e. with euler)
dm1/dt <- a*m1+b*m2+d*m3+e*m4
dm2/dt <- a*m1+b*m2+d*m3+e*m4
dm3/dt <- a*m1+b*m2+d*m3+e*m4
dm4/dt <- a*m1+b*m2+d*m3+e*m4
with following initial values:
m1<- PID$sato[1]
m2<- PID$hrs[1]
m3<- PID$fim[1]
m4<- PID$health[1]
a,b,d,e are free coeffient.
The parameters a,b,d,e are not fix. My goal is t...
2012 Mar 21
2
Type II and III sum of squares (R and SPSS)
...elow I report the list of commands used in R ("fattA" is the beteween
factor and "fB" is the within factor):
> a1b1<-c(10,9,8,7)
> a1b2<-c(7,6,4,5)
> a1b3<-c(3,2,3,4)
> a2b1<-c(9,9,8,7)
> a2b2<-c(8,7,9,7)
> a2b3<-c(7,8,8,6)
>
> M3<-matrix(0,8,4)
> M3[,1]<-cbind(a1b1,a2b1)
> M3[,2]<-cbind(a1b2,a2b2)
> M3[,3]<-cbind(a1b3,a2b3)
> M3[,4]<-rep(c(1,2),each=4)
>
> colnames(M3)<-c("b1","b2","b3","fattA")
>
> M3<-as.data.frame(M3)
>
>...
2012 Nov 06
1
Filling matrix elements with a function
Hi all,
I have a matrix simulating migration in a spatial model. I want to be
able to define movement (the values of m1, m2 and m3) as only != 0
between adjacent patches and contingent on certain conditions as
defined in the function.
Here is the code:
WET<-function(t) {everglades$precipitation[t]} #simply reads
precipitation data from a csv, value is numeric
AB<-function(WET,t) {ifelse(WET(t)>0, 0.0001, 0)}
BC&...
2012 Jul 18
0
[LLVMdev] Setting up a cross-compiler for cortex-m3
On Wed, Jul 18, 2012 at 3:52 PM, Renato Golin <rengolin at systemcall.org> wrote:
> On 18 July 2012 14:33, salvatore benedetto
> <salvatore.benedetto at gmail.com> wrote:
>> but I still haven't figure out how to build for cortex-m3
>>
>> clang -march=armv7-m -mfloat-abi=soft <something missing?> testReference.cpp -c
>
> -march should have done the trick.
>
> You can also try -mcpu=cortex-m3,
>
> or try -ccc-host-triple armv7m-none-gnueabi (or -eabi),
I've tried with all of the followin...
2012 Jul 22
2
[LLVMdev] Setting up a cross-compiler for cortex-m3
On Sun, Jul 22, 2012 at 9:42 PM, Anton Korobeynikov
<anton at korobeynikov.info> wrote:
>> Any suggestions?
> Try to specify CPU explicitly.
Already did.
clang++ -ccc-gcc-name arm-none-linux-gnueabi-g++ -ccc-host-triple
thumbv7m-none-gnueabi testReference.cpp -c -mcpu=cortex-m3
fatal error: error in backend: CPU: 'cortex-m3' does not support ARM
mode execution!
S.
2017 Jun 06
2
integrating 2 lists and a data frame in R
> On Jun 6, 2017, at 4:01 AM, Jim Lemon <drjimlemon at gmail.com> wrote:
>
> Hi Bogdan,
> Kinda messy, but:
>
> N <- data.frame(N=c("n1","n2","n3","n4"))
> M <- data.frame(M=c("m1","m2","m3","m4","m5"))
> C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"), I=c(100,300,400))
> MN<-as.data.frame(matrix(NA,nrow=length(N[,1]),ncol=length(M[,1])))
> names(MN)<-M[,1]
> rownames(MN...
2017 Jun 06
4
integrating 2 lists and a data frame in R
Dear Bert,
thank you for your response. here it is the piece of R code : given 3 data
frames below ---
N <- data.frame(N=c("n1","n2","n3","n4"))
M <- data.frame(M=c("m1","m2","m3","m4","m5"))
C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"), I=c(100,300,400))
how shall I integrate N, and M, and C in such a way that at the end we have
a data frame with :
- list N as the col...
2017 Jun 06
0
integrating 2 lists and a data frame in R
Hi Bogdan,
Kinda messy, but:
N <- data.frame(N=c("n1","n2","n3","n4"))
M <- data.frame(M=c("m1","m2","m3","m4","m5"))
C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"), I=c(100,300,400))
MN<-as.data.frame(matrix(NA,nrow=length(N[,1]),ncol=length(M[,1])))
names(MN)<-M[,1]
rownames(MN)<-N[,1]
C[,1]<...
2009 Dec 30
1
Fwd: Negbin Error Warnings
...sage ----------
From: Clara Brück <clara_brueck@web.de>
Date: 2009/12/30
Subject: Re: [R] Negbin Error Warnings
To: milton ruser <milton.ruser@gmail.com>
Dear Milton,
Thanks for your email. I ran the NBR-model and then tried to compute the
first differences of some
expected values:
m3.x1 <-setx(m3, advisory=1)
m3.x2 <-setx(m3, advisory=0)
m3.s1 <-sim (m3,x=m3.x1, x1=m3.x2)
However, I then get 50 error warnings, and there all like this:
In rpois(k, (mu * rgamma(k, theta))/theta) : NAs produced
Bests
Clara
> -----Ursprüngliche Nachricht-----
> Von: "milton...
2013 Oct 12
0
[LLVMdev] Generate code for ARM Cortex m0, m3, and m4.
Hi Jan,
For Cortex-M0, you should probably use the armv6m string in the target
triple. For M3 and M4 you need to use the thumbv7m arch string, -mthumb
won't be necessary.
Amara
On 11 October 2013 19:23, Jan Hoogerbrugge <
jan.hoogerbrugge at biface-tools.com> wrote:
> Hi,
>
> I am trying to cross compile code for ARM Cortex m0, m3, and m4.
>
> For m0, I use:
>...