search for: stds

...<-read.csv(file="y.csv",head=FALSE,sep=";") Y=y[,1] D_1=mat[,1] D_2=mat[,2] D_3=mat[,3] D_4=mat[,4] Y_1=mat[,5] Y_2=mat[,6] Y_3=mat[,7] Y_4=mat[,8] s_F=mat[,9] s_G=mat[,10] L.minor=data.frame(Y=Y,D_1=D_1,D_2=D_2,D_3=D_3,D_4=D_4,Y_1=Y_1,Y_2=Y_2,Y_3=Y_3,Y_4=Y_4,s_F=s_F,s_G=s_G) stdS=1/(sd(s_F)) stdG=1/(sd(s_G)) L.minor.m1<-nls(Y~a_1_1*D_1+a_2_1*D_2+a_3_1*D_3+a_4_1*D_4+a_1_2*(exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)))/(1+exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)))*D_1+a_2_2*(exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F...

...turns the mean and standard deviation for each row in each sample. Sample=function(name, n){ res=replicate(n,name[,sample(colnames(name),n)],simplify=FALSE) Means=do.call(rbind,lapply(res,function(x) rowMeans(x))) rownames(Means)=paste('sample',1:n,sep="") Gmean=colMeans(Means) STDs=do.call(rbind,lapply(res,function(x) rowSds(x))) rownames(STDs)=paste('sample',1:5,sep="") Gsd=sd(STDs) return(Gmean,Gsd)} I then need to use the mean and standard deviation from each row to calculate the Power for 2 sample means. here is the function I use to do it: Power=func...

...(1:Total) result=replicate(nSample,trans[,sample(colnames(trans),nPatient)],simplify=FALSE) show(result) Means=do.call(rbind,lapply(result,function(x) rowMeans(x))) rownames(Means)=paste('sample',1:nSample,sep="") show(Means) Gmean=colMeans(Means) Gmean # load genefilter package STDs=do.call(rbind,lapply(result,function(x) rowSds(x))) rownames(STDs)=paste('sample',1:nSample,sep="") show(STDs) Gsd=sd(STDs) Gsd return(Gmean,Gsd)} -- View this message in context: http://www.nabble.com/Help-with-functions-tp20113920p20113920.html Sent from the R help mailing lis...

...to first estimate parameter values from a set of calibration measurements, and then later to use those estimates to calculate another set of values with “predict.lm”. First I have a calibration dataset of absorbance values measured from standard solutions with known concentration of Bromide: > stds abs conc 1 -0.0021 0 2 0.1003 200 3 0.2395 500 4 0.3293 800 On this small calibration series, I perform a linear regression to find the parameter estimates of the relationship between absorbance (abs) and concentration (conc): > linear1 <- lm(abs~conc, data=stds) > summary...

...am/smb/samba-3.0.0/source/smbwrapper -I. -D_LARGEFILE64_SOURCE -D_FILE_OFFSET_BITS=64 -D_GNU_SOURCE -I/usr/local/include -I/home/adam/smb/samba-3.0.0/source LIBS = -lcrypt -lresolv -lnsl -ldl LDSHFLAGS = -shared LDFLAGS = Compiling auth/pass_check.c auth/pass_check.c:38:22: afs/stds.h: No such file or directory auth/pass_check.c:39:25: afs/kautils.h: No such file or directory auth/pass_check.c: In function `afs_auth': auth/pass_check.c:53: error: `KA_USERAUTH_VERSION' undeclared (first use in this function) auth/pass_check.c:53: error: (Each undeclared identifier is r...

...number. Q2 error message I tried to display sub-model with only temp_ran using the code below but failed. Please kindly suggest the potential failure cause. code rm(list=ls()) library(MuMIn) datam <-read.csv("c:/migration/Mig_ratio_20100817.csv",header=T, row.names=1) mig.stds <-glm(datam.std$SummerM_ratio~ datam.std$temp_max++datam.std$evi_mean+datam.std$topo_var+datam.std$topo_mean+datam.std$coast+datam.std$Iso_index_0808,data=datam.std,family=gaussian) target.model <- dredge(mig.stds, subset=datam.std$temp_ran) error in eval(expr, envir, enclos), 'tem...

...code (which is not very well, obviously), each iteration assigns one element from the sample.from.vector to the current make.vector element; the element choosen from sample.from.vector is tagged so that i is not selected again. However, after the loop, make.vector contains, for example, 1758 stds and 242 devs. Each iteration results in a different proportion of stds and devs (although close to the desired, not right on). Any help would be greatly apprecitated. -Mark Orr ________________________________ Mark G. Orr Postdoctoral Research Fellow Dept. of Neuroscience RM 825 Kennedy...

...wever, upon doing this, the function stops with an error report "Matrix build from transformed variables has a constant column" I was able to fix it as follows. In the following line, taken from the same block of code producing the stop error, modify the existing line from X.to.C.stds <- sqrt(apply(X.to.C.w, 2, var)) to X.to.C.stds <- (!is.null(sweep.out)) * c(1,sqrt(apply(X.to.C.w[,-is.null(sweep.out)], 2, var))) + (is.null(sweep.out)) * sqrt(apply(X.to.C.w,2,var))) That should work.

...c/xdr.h:249: warning: this is the location of the previous definition /usr/afsws/include/rx/xdr.h:268: warning: `IXDR_PUT_U_SHORT' redefined /usr/include/rpc/xdr.h:250: warning: this is the location of the previous definition In file included from passdb/pass_check.c:140: /usr/afsws/include/afs/stds.h:94: warning: useless keyword or type name in empty declaration /usr/afsws/include/afs/stds.h:94: warning: empty declaration /usr/afsws/include/afs/stds.h:108: warning: useless keyword or type name in empty declaration /usr/afsws/include/afs/stds.h:108: warning: empty declaration In file included...

Hi all. Belove is the example from the cluster-help page wtih the output. I simply cannot figure out how to relate the estimate and robust Std. Err to the p-value. I am aware this a marginal model applying the sandwich estimator using (here I guess) an emperical (unstructered/exchangeable?) ICC. Shouldent it be, at least to some extend, comparable to the robust z-test, for rx :

Hi I am trying to print out means, STDs and histograms under two sets of factors. I can manage it for one set - see below but not for two sets. That is, I want ot print out the mean STD and Histogram for each ITEM code within each DELIVERABLE code. In addition I can only get to view the histogram for the last item. How do you get R to s...

...der Windoze on the same machine the error bars plotted in the code below appear instantaneously. has anyone else noticed this problem under Ubuntu 9.04? # make a load of means: J<- 100 mu<- .02 sigma<- .04 set.seed(1) y.bar<- rnorm(J, mean=mu, sd=sigma) # make a load of (fake) stds about the means: n.j<- round(runif(J, min=10, max=100)) sigma.alpha<- sqrt(var(y.bar) / n.j) ylim<- range(c(y.bar-2*sigma.alpha, y.bar+2*sigma.alpha)) par(las=1) plot(n.j, y.bar, cex.lab=.9, cex.axis=1, xlab="sample size", ylab="some stuff", p...

...with-ldap=/opt/openldap-2.2.4/ --with-krb5=/opt/kerberos5-1.3.1/ --with-afs --with-mysql-prefix=/opt/mysql-4.0.17/ --with-winbind --with-included-popt But i always recieve this error: Compiling auth/auth_ntlmssp.c Compiling auth/pampass.c Compiling auth/pass_check.c auth/pass_check.c:38:22: afs/stds.h: No such file or directory auth/pass_check.c:39:25: afs/kautils.h: No such file or directory auth/pass_check.c: In function `afs_auth': auth/pass_check.c:53: error: `KA_USERAUTH_VERSION' undeclared (first use in this function) auth/pass_check.c:53: error: (Each undeclared identifier is re...

...epts the test based on Cochran's Recommendations. * R says " chi squared approximation may be incorrect." Does it mean the same as what EpInfo saying " Chi square is not valid" Regards, Dr. B. Manoj Aravind. *Table for analysis* Number of STDs Identified Type of Health worker <2 2 or > ASHA 39 5 AWW 22 1 ANM 1 12 ...................................... *In R the oup...