Displaying 20 results from an estimated 1200 matches similar to: "Help to make a for for index"
2004 Aug 25
3
integrate function
Is possible to integrate this diferential equation:
dN/dt = Nr(1-(N/K))
in R using the integrate() function?
Or any other diferential equation?
If yes, how?
If no, anybody know any software on linux that make this?
Inte
Ronaldo
--
O problema de ter os dois p??s bem firmes no ch??o ?? que voc?? n??o vai conseguir
tirar as cal??as.
--An??nimo
--
|> // | \\
2005 Nov 16
3
Difficulties with for() {while(){}}
Hi,
I have the follow function:
function() {
## Init of function
...
for(i in test) {
...
while(j <= test2) {
...
}
}
}
The problem is that sometimes, naturally, the while is not possible to be
resolved, and so the program abort.
In this case I need that program return to the init of function and run again.
How I can make this? Abort the while, abort the for
2005 Sep 09
2
best way to fit a model
Hi,
I have some data that have this behaviour:
|
|*******
| *
| *
| *
| *
|----------------
What is the best and simpler way to fit this in R?
Thanks
Ronaldo
--
Ela pilotava um Continenal 2001 com igni????o autom??tica Magiclic...
--
|> // | \\ [***********************************]
| ( ?? ?? ) [Ronaldo Reis J??nior ]
|>
2006 Mar 15
1
manipulating weeks dates
Hi,
I have these vectors:
> WEEK <- rep(c(1:52),2)
> YEAR <- rep(c(2000,1999),c(52,52))
How to make a vector of Date with weeks in years? I try as.date from survival
package, but it dont work with weeks, just only with days, months etcs.
Thanks
Ronaldo
--
"Realmente minha cidade e muito facultativa"
--Elivelton, ao rep?rter da Jovem Pan que falava das muitas
2003 Oct 14
1
[OFF] Dataset for extra Crawley Chapter
Hi,
anybody have the dataset used in Gamma Errors chapter of the Crawley's books
(An Introduction to Data Analysis using S-Plus).
specifically the functionalresponse and the Density datasets.
Thanks
Ronaldo
--
For every problem there is one solution which is simple, neat, and wrong.
-- H. L. Mencken
--
|> // | \\ [***********************************]
| ( ? ? ) [Ronaldo Reis
2005 Oct 24
3
lme and lmer syntax
Hi,
I have this:
lme(y~x1+x2,random=~1|x1/x2)
How to make this random effect using lmer?
I try this:
lmer(y~x1+x2+(1|x1/x2)
But it dont work.
Any idea?
Thanks
Ronaldo
--
System halted!
--
|> // | \\ [***********************************]
| ( ?? ?? ) [Ronaldo Reis J??nior ]
|> V [UFV/DBA-Entomologia ]
| / \ [36570-000 Vi??osa -
2002 Jun 21
3
generating points over a function, its possible???
Hi,
I have a new (for me) situation.
I have a function for one linear model:
y = 9.7909-0.035*x
I have the variance of system
var = 13.01403
Is possible to create a random vector points of Y (simulating the real data)
over this curve using the function and the variance? I need this to a
discipline.
I have this:
| \
| \ -> y = 9.7909-0.035*x and var = 13.014
| \
| \
| \
| \
2003 Oct 14
2
different results depending of variable position.
Hi,
I make an analysis and depending of the order of the variables, the
significance change, look.
m1 <- glm((infec/ntot)~idade+sexo+peso,family=binomial,weights=ntot)
> anova(m1,test="F")
Analysis of Deviance Table
Model: binomial, link: logit
Response: (infec/ntot)
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev F Pr(>F)
2003 Feb 04
2
testing slope
Hi all,
I try to test a linear slope using offset.
I have:
> m2 <- glm(Y~X*V)
> summary(m2)
Call:
glm(formula = Y ~ X * V)
Deviance Residuals:
Min 1Q Median 3Q Max
-2.01688 -0.56028 0.05224 0.53213 3.60216
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.3673 0.8476 1.613 0.119788
X
2005 May 09
4
Error in F test on version 2.1.0
Hi,
I make a upgrade to R 2.1.0 and in some analysis I give an error:
anova(model,test="F")
Analysis of Deviance Table
Model: binomial, link: logit
Response: landing/total
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev F Pr(>F)
NULL 16 105.079
trat 1 93.149 15 11.930 93.15
Warning message:
2004 Jun 03
1
[OFF] program to estimate the best fit
Hi,
exist in R or any other program for linux that estimate the best fit for data
using severals functions? Somethink like tablecurve
Thanks
Ronaldo
--
Um menino-prodigio e uma crianca cujos pais tem muita
imaginacao.
-- Jean Costeau
--
|> // | \\ [***********************************]
| ( ?? ?? ) [Ronaldo Reis J??nior ]
|> V [UFV/DBA-Entomologia
2006 Jan 23
1
simple problem
Hi,
look this:
> summary(fam??lia)
Anacardiaceae Annonaceae Bombacaceae Cecropiaceae
2 4 1 3
Chrysobalanaceae Clusiaceae Euphorbiacaea Fabacea
1 1 4 3
Fabaceae Flacourtiaceae Humiriaceae indeterminada
2003 Mar 21
2
Trying to make a nested lme analysis
Hi,
I''m trying to understand the lme output and procedure.
I''m using the Crawley''s book.
I''m try to analyse the rats example take from Sokal and Rohlf (1995).
I make a nested analysis using aov following the book.
> summary(rats)
Glycogen Treatment Rat Liver
Min. :125.0 Min. :1 Min. :1.0 Min. :1
1st Qu.:135.8
2003 Mar 06
1
Problems with variable types.
Hi all,
I have problems in a dataframe variables types.
Look:
from a loop function:
for(...){
...
dados.fin <- rbind(dados.fin, c(L=j, A=j^2,
Nsp=nsps,
N=length(amosfin$SP),
AmT="am",NAm=nam,
AMST=amst))
dados.fin <- rbind(dados.fin, c(L=j, A=j^2,
2004 Mar 05
1
Problems with SJava instalation
Hi,
I'm try to use the SJava package. The install is OK. In R I have this error
message:
--------------------------
> library(SJava)
Error in dyn.load(x, as.logical(local), as.logical(now)) :
unable to load shared library
"/opt/lib/R/site-library/SJava/libs/SJava.so":
libRSNativeJava.so: cannot open shared object file: No such file or
directory
Error in
2002 Aug 02
1
doubt about contrib...
Hi,
I have a little doubt about contrib packages.
I make some little functions in a R package format, this functions are for my
personal use (envelopes for glm models and a plot bar routine).
Some friends ask me about the liberation of this functions in a web site, R
contrib etc.
what is the correct procedure for send a package to R contrib?
The functions is not create for me, is adaptaded to R
2003 Jun 20
1
glmm and overall goodness of fit
Hi,
exist in R any glmm function that have any tools for test for overall goodness
of fit?
Thanks
Ronaldo
--
O papel da impressora ? sempre mais forte na parte picotada.
--
| // | \\ [***********************************]
|> ( ? ? ) [Ronaldo Reis J?nior ]
| V [UFV/DBA-Entomologia ]
|> / \ [36571-000 Vi?osa - MG ]
|
2004 Nov 09
1
glm.nb stop on Error.
Hi,
I make an analysis sequence on R. In some cases the function glm.nb fail to
ajust the model. Its Ok. The problem is that this error stop the program. I
need treat this error and not stop the program. Something like this:
...
model <- glm.nb(y~x,maxit=1000)
if(glm.nb fail) {
teste[i] <- 0
}
else {
teste[i] <- anova(modelo)$"P(>|Chi|)"[2]
}
...
I try this:
...
2006 Apr 04
1
Problem with Crawley book example
Hi,
I try to run the example of Crawley's Book on the page 661, but it fail, look
> repmeasures <-
read.table("../Packages/Crawley/data/repmeasures.txt",header=T)
> attach(repmeasures)
> rep <- as.factor(rep)
> library(nlme)
> model <- lme(height~seed,random=~time|rep/seed)
Erro em lme.formula(height ~ seed, random = ~time | rep/seed) :
iteration limit
2005 Sep 20
1
Change the mirror
Hi,
Please, change the brazilian mirror http://www.termix.ufv.br/CRAN/ to
http://www.insecta.ufv.br/CRAN/ in R homepage.
Thanks
ROnaldo
--
Se dois homens no mesmo trabalho concordam o tempo todo, um deles ?? demais. Se
discordam sem parar, ent??o os dois s??o dispens??veis
--Darryl F. Zanuck
--
|> // | \\ [***********************************]
| ( ?? ?? ) [Ronaldo Reis J??nior