Hola,
He probado con diferentes formas que me ha ido proporcionando la IA de
Claude y al final, esta es la que funciona.
Al parecer la función del paquete, tiene un bug (según Claude).
#--------------------------
calcular_aiken_ic <- function(datos, min_val = 1, max_val = 7, conf = 0.95)
{
n <- nrow(datos)
resultados <- sapply(datos, function(col) {
S <- sum(col)
V <- (S - n * min_val) / (n * (max_val - min_val))
z <- qnorm(1 - (1 - conf)/2)
se <- sqrt(V * (1 - V) / n)
lower <- max(0, V - z * se) # No puede ser menor que 0
upper <- min(1, V + z * se) # No puede ser mayor que 1
c(AikenV = V, lowerCI = lower, upperCI = upper)
})
return(t(resultados))
}
# Aplicar
resultados <- calcular_aiken_ic(df[, 2:ncol(df)])
resultados
#--------------------------
> #--------------------------> calcular_aiken_ic <- function(datos, min_val =
> 1, max_val = 7, conf = 0.95) {+ n <- nrow(datos)+ + resultados <-
> sapply(datos, function(col) {+ S <- sum(col)+ V <- (S - n *
> min_val) / (n * (max_val - min_val))+ + z <- qnorm(1 - (1 -
> conf)/2)+ se <- sqrt(V * (1 - V) / n)+ + lower <-
> max(0, V - z * se) # No puede ser menor que 0+ upper <- min(1, V + z
> * se) # No puede ser mayor que 1+ + c(AikenV = V, lowerCI =
> lower, upperCI = upper)+ })+ + return(t(resultados))+ }> > #
> Aplicar> resultados <- calcular_aiken_ic(df[, 2:ncol(df)])> resultados
> AikenV lowerCI upperCI
item1 0.9000000 0.7685216 1
item2 0.9166667 0.7955377 1
item3 0.8833333 0.7426415 1
item4 0.8750000 0.7300587 1
item5 0.8916667 0.7554547 1
item6 0.8916667 0.7554547 1
item7 0.9083333 0.7818709 1
item8 0.8750000 0.7300587 1
item9 0.9000000 0.7685216 1
item10 0.9083333 0.7818709 1
item11 0.9083333 0.7818709 1
item12 0.9000000 0.7685216 1> #--------------------------
>
Gracias,
Carlos Ortega
www.qualityexcellence.es
El dom, 25 ene 2026 a las 22:17, Jose Betancourt Bethencourt (<
[email protected]>) escribió:
> rm(list = ls())
>
> library(readxl)
> df <- read_excel("D:/Betan/Escritorio/encuestas/df.xlsx")
>
> df
> library(spheredata)
>
>
>
> attach(df)
> aikenV(df)
>
> ikenV lowerCI upperCI
> item1 1.80 NaN NaN
> item2 1.83 NaN NaN
> item3 1.77 NaN NaN
> item4 1.75 NaN NaN
> item5 1.78 NaN NaN
> item6 1.78 NaN NaN
> item7 1.82 NaN NaN
> item8 1.75 NaN NaN
> item9 1.80 NaN NaN
> item10 1.82 NaN NaN
> item11 1.82 NaN NaN
> item12 1.80 NaN NaNWarning message:
> In sqrt(4 * n * k * V * (1 - V) + (z^2)) : NaNs produced
>
> hay algún error que no veo
> saludos
>
> El dom, 25 ene 2026 a las 13:31, Carlos Ortega (<[email protected]>)
> escribió:
>
>> Hola,
>>
>> Así, funciona:
>>
>> #----------------------------------------
>> library(spheredata)
>> library(readxl)
>>
>> aikenV(df[, 2:ncol(df)])
>> #----------------------------------------
>>
>> El problema que tenías era que la primera columna "Personas" que son los
>> individuos, no tiene que entrar en el análisis, solo las columnas de las
>> mediciones.
>>
>> Gracias,
>> Carlos Ortega
>> www.qualityexcellence.es
>>
>>
>> El dom, 25 ene 2026 a las 19:21, Proyecto R-UCA (<[email protected]>)
>> escribió:
>>
>> > Buenas Jose:
>> >
>> > ¿Cómo estás cargando el fichero de datos? Sospecho que el problema es
>> ese.
>> >
>> > ¿De qué paquete estás usando la función aikenV?
>> >
>> > Un saludo.
>> >
>> > El dom, 25-01-2026 a las 09:05 -0500, Jose Betancourt Bethencourt
>> escribió:
>> > > estimados
>> > > al tratar de calcularlo
>> > > aikenV(df)
>> > >
>> > > me da el error
>> > >
>> > >
>> > >
>> > > Error in colSums(data) : 'x' must be numeric
>> > >
>> > > adjunto excel
>> > >
>> > > aprciaria ayuda
>> > > saludos
>> > > jose
>> > > --
>> > > Dr. Jose A. Betancourt Bethencourt
>> > > Universidad de Ciencias Medicas Carlos j. Finlay
>> > > _______________________________________________
>> > > R-help-es mailing list
>> > > [email protected]
>> > >
>> >
>> https://urldefense.com/v3/__https://stat.ethz.ch/mailman/listinfo/r-help-es__;!!D9dNQwwGXtA!XsqN6H5xuanagnHR3-vU-dpF24ARnSKsScgtqztxZKYeQ9ObzrpDH8THhhgQ7Gazuw9msTKXtqNGJ9Y$
>> > >
>> >
>> > _______________________________________________
>> > R-help-es mailing list
>> > [email protected]
>> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>> >
>>
>>
>> --
>> Saludos,
>> Carlos Ortega
>> www.qualityexcellence.es
>>
>> [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> R-help-es mailing list
>> [email protected]
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>
>
> --
> Dr. Jose A. Betancourt Bethencourt
> Universidad de Ciencias Medicas Carlos j. Finlay
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
