similar to: lmer with only random effect

Hello List and thanks in advance for all of your help, I am trying implement a permutation test of a multinomial logistic regression ('multinom' within the nnet package). In the end I want to compare the parameter estimate from my data to the distribution of randomized parameter estimates. I have figured out how to permute my dependent variable (MNNUM) x number of times, apply

Dear R gurus, I have the following code, but I still not know how to estimate and extract confidence intervals (95%CI) from resampling. Thanks! ~Adriana #data penta<-c(770,729,640,486,450,410,400,340,306,283,278,260,253,242,240,229,201,198,190,186,180,170,168,151,150,148,147,125,117,110,107,104,85,83,80,74,70,66,54,46,45,43,40,38,10) x<-log(penta+1) plot(ecdf(x),

I looked at the documentation of source() and summary(), and I could not find the reason why calling something like: > summary(resamps) from the command line, works (it prints the summary) whereas calling summary(resampls) from a file that I source with source("my_file.r") does not print anything. How can I get summary(resamps) to print when I source a file with this command?

Using the following MATLAB snippet: fnames={'chirp','perfect','block','filter'} for k=1:length(fnames) fn=fnames{k}; myfile=fopen([fn, '.fl'], 'r', 'ieee-le'); x=fread(myfile, Inf, 'float32', 0, 'ieee-le'); fclose(myfile); X=specgram(x,2048,1,kaiser(2048,14)); spectrogram_floor=-96;

Hi, I have 35 data files for reading. I would like get a program for performing reading of 35 files at once. All are of the type: Dados1.raw, Dados2.raw and so on. If the files have the same number of columns, I can read with the following commands: rm(list=ls()) filenames = list.files(path="~/Silvano/Arq", pattern="Dados+.*raw") names = substr(filenames, 1, 7) for(i in

Say I train a model in caret, e.g.: RFmodel <- train(X,Y,method='rf',trControl=myCtrl,tuneLength=1) How can I save this to disk and load it later in R? How about an object of the class "resamples"? resamps <- resamples( list( RF = RFmodel, SVM = SVMmodel, KNN = KNNmodel, NN = NNmodel )) Thanks,

Hello to the R world... I have some problems regarding a GLM - repeated measures analysis. I want to test overall differences between AgeClass and Treatment (between subject) with OpenR1+OpenR2+OpenR3 (repeated measures, within subject). The table looks kind like this: AgeClass Treatment OpenR1 OpenR2 OpenR3 1 1 0 0 12.63 1 1 12.67 3.83 45.67 1 1 38.46 65.38 75.21 1 1 14.46 0 17.96 1 2 27.83

Hello, One more question.. I have the data.frame "pop": xloc yloc gonad ind Ene W Area 1 23 20 516.74 1 0.02 20.21 1 2 23 20 1143.20 1 0.02 20.21 1 3 23 20 250.00 1 0.02 20.21 1 4 22 15 251.98 1 0.02 18.69 2 5 22 15 598.08 1 0.02 18.69 2 6 21 19 250.00 1 0.02 20.21 3 7 22 20 251.98 1

Forgot to send to R-help From: Nordlund, Dan (DSHS/RDA) Sent: Thursday, February 07, 2013 2:09 PM To: 'James Jong' Subject: RE: [R] Sourcing my file does not print command outputs James, Your code seems to have ‘…’ sitting on a line all by itself (maybe should be at the end of the preceding comment? Anyway, when I eliminated that problem and sourced the script using the following call

Hello, I'm trying to create a matrix (95x55) with data from a data.frame pop: xloc yloc go ind Ene totW 1 23 20 516 1 0.02 20.21 2 23 20 1143 1 0.02 20.21 3 23 20 250 1 0.02 20.21 4 22 15 251 1 0.02 18.69 5 22 15 598 1 0.02 18.69 6 21 19 250 1 0.02 20.21 7 22 20 251 1

Hi All, I have a data.frame: >pop consumed ind recruits gonad cell 1 516.74 1 0.02 20.21 0.25 2 1143.20 1 0.02 20.21 0.50 3 250.00 1 0.02 20.21 0.25 4 251.98 1 0.02 18.69 0.25 5 598.08 1 0.02 18.69 0.25 6 437.38 1 0.02 18.69 0.25 7 250.00 1 0.02 17.39 0.25 8 250.00 1 0.02 17.39 0.25 9

** Disclaimer: I'm looking for general suggestions ** I'm sorry, but can't send out the file I'm using, so there is no reproducible example. I'm using read.table and it's taking over 30 seconds to read a tiny file. The strange thing is that it takes roughly the same amount of time if the file is 100 times larger. After re-reviewing the data Import / Export manual I think

Hello List! We're having trouble making call deflection on ISDN PRI. We would like to transfer a call to an external extension but keeping the callerid of the caller so it can be presented to the receiver of the transferred call. At the time we're using Zaptel 1.4.5.1, Asterisk 1.4.11 and Digium hardware TE420B. We've ordered the service (CD) from the phone company. The

Hi Simon, Thanks for reposting the example, and looking at it more carefully, I think it is very similar to my first proposal. This was met with some resistance here because it dumps loop information in the vector preheader. Doing it this early, we want to emit this in the vectoriser, puts a restriction on (future) optimisations that transform vector loops to honour/update/support this intrinsic

Invitation accepted, I am happy to help out with reviews, like I did with the previous VP patches. And of course agreed that things should be well defined, and that we shouldn't paint ourselves in a corner, but I don't think that this is the case. And it's not that I am in a rush, but I don't think this change needs to be predicated on a big change landing first like the LV

> You have similar problems with https://reviews.llvm.org/D79100 The new revision D79100<https://reviews.llvm.org/D79100> solves your comment 1), and I don't think your comments2) and 3) apply as there are no vendor specific intrinsics involved at all here. Just to quickly discuss the optimisation pipeline, D79100<https://reviews.llvm.org/D79100> is a small extension for the

Hi everyone, This is a graphics question, concerning presentation rather than actually producing figures. After much fiddling I've managed to create bar plots showing error bars. I've put in an x-axis, which goes along all the bars, but doesn't appear between the y-axis and the first bar. I've searched all the resources I have access to but can't find anything addressing this

Hi R Users, When I use package lme4 for mixed model analysis, I can't distinguish the significant and insignificant variables from all random independent variables. Here is my data and result: Data: Rice<-data.frame(Yield=c(8,7,4,9,7,6,9,8,8,8,7,5,9,9,5,7,7,8,8,8,4,8,6,4,8,8,9), Variety=rep(rep(c("A1","A2","A3"),each=3),3),

Hi all! I'm going mad about getting my shares mapped to my win2k and winXP machines from my Samba PDC Version 3.0.2 on Debian. Win-machines doesn't seem to execute logon.bat at all, except when i log on as root. I've set the ownership and chmod 0775 all logon.bat files to eatch user on my network. Am i just blind or are there something obvious/very wrong with my config? my smb.conf

Hi, I abandoned that approach and followed Eli's suggestion, see somewhere earlier in this thread, and emit an intrinsic that represents/calculates the active mask. I've just uploaded a new revision for D79100 that implements this. Cheers. ________________________________ From: Simon Moll <Simon.Moll at EMEA.NEC.COM> Sent: 18 May 2020 13:32 To: Sjoerd Meijer <Sjoerd.Meijer at