Re: [R] problem in saving the history file '~/.Rhistory'

2019-06-26 Thread William Dunlap via R-help
Did you (or a startup script) set the environment variable R_HISTORY to be '~/somefile'? R does not expand the tilde in this case so you need to let the shell do it. Perhaps ~ doesn't expand for the root user. (IMO, it is crazy to run R, or almost any other program, as root.) Bill Dunlap TIBCO

[R] problem in saving the history file '~/.Rhistory'

2019-06-26 Thread lejeczek via R-help
hi guys, I'm on Centos 7.6 with R 3.6 and I see this: >  Save workspace image? [y/n/c]: y Warning message: problem in saving the history file '~/.Rhistory'  from strace: rt_sigaction(SIGWINCH, {SIG_DFL, [], SA_RESTORER, 0x7fb9dd34b5d0}, {0x7fb9dc383820, [], SA_RESTORER|SA_RESTART, 0x7fb9dd34b5d

Re: [R] Help with maximum likelihood estimation

2019-06-26 Thread peter dalgaard
Bricks fly fine with sufficient thrust, but you have lff with a mu argument that never gets used, so the negative log-likelihood is constant and mle() cannot minimize it. You need to read up on the definition of (log-) likelihood and write a proper one for your problem. -pd > On 26 Jun 2019,

[R] IRT discrimination value (ltm and psych package)

2019-06-26 Thread shreepad khandve
Hello Sir, I am learning R and its syntax and I have successfully converted irt.item.diff.rasch into python code and pass the inputs as per the function made in R. In R :- function (items) { ncases <- nrow(items) item.mean <- colMeans(items, na.rm = TRUE) item.mean[item.mean < (1/n

Re: [R] How to select max data according to week?

2019-06-26 Thread SITI AISYAH BINTI ZAKARIA
Dear All, Thank you very much for helping me.. i will try my best by using all suggestion.. thank you again.. hopefully is work. - Original Message - From: "Eric Berger" To: "Marc Schwartz" Cc: "Bert Gunter" , "R-help" , "SITI AISYAH BINTI ZAKARIA" Sent: Thursday, June 20, 2019 3:2

[R] Help with maximum likelihood estimation

2019-06-26 Thread avadhoot velankar
I am analyzing animal movement pattern using levy flight pattern and want to fit power function to observed data and estimate exponent using Maximum Likelihood Estimation. I am using lff<-function(mu){1-1/mean(log(x))} library(stats4) mle(lff, start = list(mu = 1)) where x is the observed data.